BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Kimitake Hayasaki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
Tidal disruption events are an excellent probe for supermassive black holes in distant inactive galaxies because they show bright multi-wavelength flares lasting several months to years. AT2019dsg presents the first potential association with neutrino emission from such an explosive event.
We show results from the radiation hydrodynamics (RHD) simulations of tidal disruption of a star on a parabolic orbit by a supermassive black hole (SMBH) based on a three-dimensional smoothed particle hydrodynamics code with radiative transfer. We find that such a tidally disrupted star fragment and form clumps soon after its tidal disruption. The fragmentation results from the endothermic processes of ionization and dissociation that reduce the gas pressure, leading to local gravitational collapse. Radiative cooling is less effective because the stellar debris is still highly optically thick in such an early time. Our simulations reveal that a solar-type star with a stellar density profile of n=3 disrupted by a 10^6 solar mass black hole produces $\sim20$ clumps of masses in the range of 0.1 to 12 Jupiter masses. The mass fallback rate decays with time, with pronounced spikes from early to late time. The spikes provide evidence for the clumps of the returning debris, while the clumps on the unbound debris can be potentially freely-floating planets and brown dwarfs. This ionization and dissociation induced fragmentation on a tidally disrupted star are a promising candidate mechanism to form low-mass stars to planets around an SMBH.
We study the mass fallback rate of tidally disrupted stars on marginally bound and unbound orbits around a supermassive black hole (SMBH) by performing three-dimensional smoothed particle hydrodynamic (SPH) simulations with three key parameters: the star is modeled by a polytrope with two different indexes ($n=1.5$ and $3$). The stellar orbital properties are characterized by five orbital eccentricities ranging from $e=0.98$ to $1.02$ and five different penetration factors ranging from $β=1$ to $3$, where $β$ represents the ratio of the tidal disruption to pericenter distance radii. We derive analytic formulae for the mass fallback rate as a function of the stellar density profile, orbital eccentricity, and penetration factor. Moreover, two critical eccentricities to classify tidal disruption events (TDEs) into five different types: eccentric ($e<e_{\rm crit,1}$), marginally eccentric ($e_{\rm crit,1}\lesssim{e}<1$), purely parabolic ($e=1$), marginally hyperbolic ($1<e<e_{\rm crit,2}$), and hyperbolic ($e\gtrsim{e_{\rm crit,2}}$) TDEs, are reevaluated as $e_{\rm crit,1}=1-2q^{-1/3}β^{k-1}$ and $e_{\rm crit,2}=1+2q^{-1/3}β^{k-1}$, where $q$ is the ratio of the SMBH to stellar masses and $0<k\lesssim2$. We find the asymptotic slope of the mass fallback rate varies with the TDE type. The asymptotic slope approaches $-5/3$ for parabolic TDEs, is steeper for marginally eccentric TDEs, and is flatter for marginally hyperbolic TDEs. The mass fallback rates of the marginally eccentric TDEs are much larger than the parabolic TDE case, while those of marginally hyperbolic TDEs are much smaller. Marginally unbound TDEs could be an origin of a very low density gas disk around a dormant SMBH.
We study high-energy neutrino emissions from tidal disruption remnants (TDRs) around supermassive black holes. The neutrinos are produced by the decay of charged pions originating in ultrarelativistic protons that are accelerated there. In the standard theory of tidal disruption events (TDEs), there are four distinct phases from debris circularization of stellar debris to super- and sub-Eddington to radiatively inefficient accretion flows (RIAFs). In addition, we consider the magnetically arrested disk (MAD) state in both the super-Eddington accretion and RIAF phases. We find that there are three promising cases to produce neutrino emissions: the super-Eddington accretion phase of the MAD state and the RIAF phases of both the non-MAD and MAD states. In the super-Eddington MAD state, the enhanced magnetic field makes it possible to accelerate the protons to $E_{p,max}~0.35 PeV (M_bh/10^{7.7}M_\odot)^{41/48}$ with the other given appropriate parameters. The neutrino energy is then $E_{ν,pk}~67 TeV (M_bh/10^{7.7}M_\odot)^{41/48}$ at the peak of the energy spectrum. For $M_bh\gtrsim10^{7.7} M_\odot$, the neutrino light curve is proportional to $t^{-65/24}$, while it follows the standard $t^{-5/3}$ decay rate for $M_bh<10^{7.7} M_\odot$. In both cases, the large luminosity and characteristic light curves diagnose the super-Eddington MAD state in TDEs. In the RIAF phase of the non-MAD state, we find $E_{p, max}~0.45 PeV (M_bh/10^7M_\odot)^{5/3}$ and $E_{ν,pk}~0.35 PeV (M_bh/10^7M_\odot)^{5/3}$, and its light curve is proportional to $t^{-10/3}$. This indicates that one can identify whether the existing RIAFs are the TDE origin or not. TDRs are potentially a population of hidden neutrino sources invisible in gamma rays.
We study the circularization of tidally disrupted stars on bound orbits around spinning supermassive black holes by performing three-dimensional smoothed particle hydrodynamic simulations with Post-Newtonian corrections. Our simulations reveal that debris circularization depends sensitively on the efficiency of radiative cooling. There are two stages in debris circularization if radiative cooling is inefficient: first, the stellar debris streams self-intersect due to relativistic apsidal precession; shocks at the intersection points thermalize orbital energy and the debris forms a geometrically thick, ring-like structure around the black hole. The ring rapidly spreads via viscous diffusion, leading to the formation of a geometrically thick accretion disk. In contrast, if radiative cooling is efficient, the stellar debris circularizes due to self-intersection shocks and forms a geometrically thin ring-like structure. In this case, the dissipated energy can be emitted during debris circularization as a precursor to the subsequent tidal disruption flare. The possible radiated energy is up to ~2*10^{52} erg for a 1 Msun star orbiting a 10^6 Msun black hole. We also find that a retrograde (prograde) black hole spin causes the shock-induced circularization timescale to be shorter (longer) than that of a non-spinning black hole in both cooling cases. The circularization timescale is remarkably long in the radiatively efficient cooling case, and is also sensitive to black hole spin. Specifically, Lense-Thirring torques cause dynamically important nodal precession, which significantly delays debris circularization. On the other hand, nodal precession is too slow to produce observable signatures in the radiatively inefficient case. We also discuss the relationship between our simulations and the parabolic TDEs that are characteristic of most stellar tidal disruptions.
Galaxy mergers produce supermassive black hole binaries, which emit gravitational waves prior to their coalescence. We perform three-dimensional hydrodynamic simulations to study the tidal disruption of stars by such a binary in the final centuries of its life. We find that the gas stream of the stellar debris moves chaotically in the binary potential and forms accretion disks around both black holes. The accretion light curve is modulated over the binary orbital period owing to relativistic beaming. This periodic signal allows to detect the decay of the binary orbit due to gravitational wave emission by observing two tidal disruption events that are separated by more than a decade.
We propose a new scenario for the evolution of a binary of primordial black holes (PBHs). We consider a dynamical friction by ambient dark matter, scattering of dark matter particles with a highly eccentric orbit besides the standard two-body relaxation process to refill the loss cone, and interaction between the binary and a circumbinary disk, assuming that PBHs do not constitute the bulk of dark matter. Binary PBHs lose the energy and angular momentum by these processes, which could be sufficiently efficient for a typical configuration. Such a binary coalesces due to the gravitational wave emission in a time scale much shorter than the age of the universe. We estimate the density parameter of the resultant gravitational wave background. Astrophysical implication concerning the formation of intermediate-mass to supermassive black holes is also discussed.
We study the warping and tearing of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on an eccentric orbit. The circumbinary disk is significantly misaligned with the binary orbital plane, and is subject to the time-dependent tidal torques. In principle, such a disk is warped and precesses, and is torn into mutually misaligned rings in the region, where the tidal precession torques are stronger than the local viscous torques. We derive the tidal-warp and tearing radii of the misaligned circumbinary disks around eccentric SMBH binaries. We find that in disks with the viscosity parameter, alpha, larger than a critical value depending on the disk aspect ratio, the disk warping appears outside the tearing radius. This condition is expressed as alpha > sqrt{H/3r} for H/r ~<0.1, where H is the disk scale height. If alpha < sqrt{H/3r}, only the disk tearing occurs because the tidal warp radius is inside the tearing radius, where most of disk material is likely to rapidly accrete onto SMBHs. In warped and torn disks, both the tidal-warp and the tearing radii most strongly depend on the binary semi-major axis, although they also mildly depend on the other orbital and disk parameters. This strong dependence enables us to estimate the semi-major axis, once the tidal warp or tearing radius is determined observationally: For the tidal warp radius of 0.1 pc, the semi-major axis is estimated to be ~10^{-2} pc for 10^7 Msun black hole with typical orbital and disk parameters. We also briefly discuss the possibility that central objects of observed warped maser disks in active galactic nuclei are supermassive black hole binaries.
We study a warping instability of a geometrically thin, non-self-gravitating, circumbinary disk around young binary stars on an eccentric orbit. Such a disk is subject to both the tidal torques due to a time-dependent binary potential and the radiative torques due to radiation emitted from each star. The tilt angle between the circumbinary disk plane and the binary orbital plane is assumed to be very small. We find that there is a radius within/beyond which the circumbinary disk is unstable to radiation-driven warping, depending on the disk density and temperature gradient indices. This marginally stable warping radius is very sensitive to viscosity parameters, a fiducial disk radius and the temperature measured there, the stellar luminosity, and the disk surface density at a radius where the disk changes from the optically thick to thin for the irradiation from the central stars. On the other hand, it is insensitive to the orbital eccentricity and binary irradiation parameter, which is a function of the binary mass ratio and luminosity of each star. Since the tidal torques can suppress the warping in the inner part of the circumbinary disk, the disk starts to be warped in the outer part. While the circumbinary disks are most likely to be subject to the radiation-driven warping on a AU to kilo-AU scale for binaries with young massive stars more luminous than 10^4Lsun, the radiation driven warping does not work for those around young binaries with the luminosity comparable to the solar luminosity.
We study a warping instability of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on a circular orbit. Such a circumbinary disk is subject to not only tidal torques due to the binary gravitational potential but also radiative torques due to radiation emitted from an accretion disk around each black hole. We find that a circumbinary disk initially aligned with the binary orbital plane is unstable to radiation-driven warping beyond the marginally stable warping radius, which is sensitive to both the ratio of vertical to horizontal shear viscosities and the mass-to-energy conversion efficiency. As expected, the tidal torques give no contribution to the growth of warping modes but tend to align the circumbinary disk with the orbital plane. Since the tidal torques can suppress the warping modes in the inner part of circumbinary disk, the circumbinary disk starts to be warped at radii larger than the marginally stable warping radius. If the warping radius is of the order of 0.1 pc, a resultant semi-major axis is estimated to be of the order of 10^-2 pc to 10^-4 pc for 10^7 Msun black hole. We also discuss the possibility that the central objects of observed warped maser disks in active galactic nuclei are binary supermassive black holes with a triple disk: two accretion disks around the individual black holes and one circumbinary disk surrounding them.
We study the basic properties of accretion flows onto binary supermassive black holes, including the cases in which a circumbinary disk is misaligned with the binary orbital plane, by means of three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a circular binary system with a misaligned circumbinary disk normally produces a double peaked mass-accretion-rate variation per binary orbit. This is because each black hole passes across the circumbinary disk plane and captures gas twice in one orbital period. Even in misaligned systems, however, a single peaked mass-accretion-rate variation per binary orbit is produced, if the orbital eccentricity is moderately large (e\lesssim0.3). The number of peaks in mass accretion rates can be understood simply in terms of the orbital phase dependence of the distance between each binary black hole and its closest inner edge of the circumbinary disk. In the cases of eccentric binary black holes having different masses, the less massive black hole can get closer to the circumbinary disk than the massive one, thus tidally splitting gas from its inner edge, but the created gas flows are comparably captured by both black holes with a short time delay. As a consequence, the combined light curve shows periodic occurrence of double-peaked flares with a short interval. This may account for the observed light variations of OJ287.
We study accretion processes for tidally disrupted stars approaching supermassive black holes on bound orbits, by performing three dimensional Smoothed Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that there is a critical value of the orbital eccentricity below which all the stellar debris remains bound to the black hole. For high but sub-critical eccentricities, all the stellar mass is accreted onto the black hole in a finite time, causing a significant deviation from the canonical $t^{-5/3}$ mass fallback rate. When a star is on a moderately eccentric orbit and its pericenter distance is deeply inside the tidal disruption radius, there can be several orbit crossings of the debris streams due to relativistic precession. This dissipates orbital energy in shocks, allowing for rapid circularization of the debris streams and formation of an accretion disk. The resultant accretion rate greatly exceeds the Eddington rate and differs strongly from the canonical rate of $t^{-5/3}$. By contrast, there is little dissipation due to orbital crossings for the equivalent simulation with a purely Newtonian potential. This shows that general relativistic precession is crucial for accretion disk formation via circularization of stellar debris from stars on moderately eccentric orbits.
When binary black holes are embedded in a gaseous environment, a rotating disk surrounding them, the so-called circumbinary disk, will be formed. The binary exerts a gravitational torque on the circumbinary disk and thereby the orbital angular momentum is transferred to it, while the angular momentum of the circumbinary disk is transferred to the binary through the mass accretion. The binary undergoes an orbital decay due to both the gravitational wave emission and the binary-disk interaction. This causes the phase evolution of the gravitational wave signal. The precise measurement of the gravitational wave phase thus may provide information regarding the circumbinary disk. In this paper, we assess the detectability of the signature of the binary-disk interaction using the future space-borne gravitational wave detectors such as DECIGO and BBO by the standard matched filtering analysis. We find that the effect of the circumbinary disk around binary black holes in the mass range $6M_sun\le{M}\lesssim3\times10^3M_sun$ is detectable at a statistically significant level in five year observation, provided that gas accretes onto the binary at a rate greater than $\dot{M}\sim1.4\times10^{17} [gs^{-1}] j^{-1}(M/10M_sun)^{33/23}$ with 10% mass-to-energy conversion efficiency, where j represents the efficiency of the angular momentum transfer from the binary to the circumbinary disk. We show that $O(0.1)$ coalescence events are expected to occur in sufficiently dense molecular clouds in five year observation. We also point out that the circumbinary disk is detectable, even if its mass at around the inner edge is by over 10 orders of magnitude less than the binary mass.
We propose the new scenario for X-ray outbursts in Be/X-ray binaries that normal and giant outbursts are respectively caused by radiatively inefficient accretion flows (RIAFs) and Bondi-Hoyle-Lyttleton (BHL) accretion of the material transferred from the outermost part of a Be disk misaligned with the binary orbital plane. Based on simulated mass-transfer rates from misaligned Be disks, together with simplified accretion flow models, we show that mass-accretion rates estimated from the luminosity of the normal X-ray outbursts are consistent with those obtained with advection-dominated accretion flows, not with the standard, radiative-cooling dominated, accretion. Our RIAF scenario for normal X-ray outbursts resolves problems that have challenged the standard disk picture for these outbursts. When a misaligned Be disk crosses the orbit of the neutron star, e.g., by warping, the neutron star can capture a large amount of mass via BHL-type accretion during the disk transit event. We numerically show that such a process can reproduce the X-ray luminosity of giant X-ray outbursts. In the case of very high Be disk density, the accretion flow associated with the disk transit becomes supercritical, giving rise to the luminosity higher than the Eddington luminosity.
We study tidal disruption and subsequent mass fallback for stars approaching supermassive black holes on bound orbits, by performing three dimensional Smoothed Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that the mass fallback rate decays with the expected -5/3 power of time for parabolic orbits, albeit with a slight deviation due to the self-gravity of the stellar debris. For eccentric orbits, however, there is a critical value of the orbital eccentricity, significantly below which all of the stellar debris is bound to the supermassive black hole. All the mass therefore falls back to the supermassive black hole in a much shorter time than in the standard, parabolic case. The resultant mass fallback rate considerably exceeds the Eddington accretion rate and substantially differs from the -5/3 power of time.
We study the interaction between the Be star and the pulsar in the TeV binary PSR B1259-63/LS 2883, using 3-D SPH simulations of the tidal and wind interactions in this Be-pulsar system. We first run a simulation without pulsar wind nor Be wind, taking into account only the gravitational effect of the pulsar on the Be disk. In this simulation, the gas particles are ejected at a constant rate from the equatorial surface of the Be star, which is tilted in a direction consistent with multi-waveband observations. We run the simulation until the Be disk is fully developed and starts to repeat a regular tidal interaction with the pulsar. Then, we turn on the pulsar wind and the Be wind. We run two simulations with different wind mass-loss rates for the Be star, one for a B2V type and the other for a significantly earlier spectral type. Although the global shape of the interaction surface between the pulsar wind and the Be wind agrees with the analytical solution, the effect of the pulsar wind on the Be disk is profound. The pulsar wind strips off an outer part of the Be disk, truncating the disk at a radius significantly smaller than the pulsar orbit. Our results, therefore, rule out the idea that the pulsar passes through the Be disk around periastron, which has been assumed in the previous studies. It also turns out that the location of the contact discontinuity can be significantly different between phases when the pulsar wind directly hits the Be disk and those when the pulsar wind collides with the Be wind. It is thus important to adequately take into account the circumstellar environment of the Be star, in order to construct a satisfactory model for this prototypical TeV binary.
If the activity of active galactic nuclei (AGNs) is predominantly induced by major galaxy mergers, then a significant fraction of AGNs should harbor binary massive black holes in their centers. We study the mass function of binary massive black holes in nearby AGNs based on the observed AGN black-hole mass function and theory of evolution of binary massive black holes interacting with a massive circumbinary disk in the framework of coevolution of massive black holes and their host galaxies. The circumbinary disk is assumed to be steady, axisymmetric, geometrically thin, self-regulated, self-gravitating but non-fragmenting with a fraction of Eddington accretion rate, which is typically one tenth of Eddington value. The timescale of orbital decay is {then} estimated as ~10^8yr for equal mass black-hole, being independent of the black hole mass, semi-major axis, and viscosity parameter but dependent on the black-hole mass ratio, Eddington ratio, and mass-to-energy conversion efficiency. This makes it possible for any binary massive black holes to merge within a Hubble time by the binary-disk interaction. We find that (1.8+-0.6%) for the equal mass ratio and (1.6+-0.4%) for the one-tenth mass ratio of the total number of nearby AGNs have close binary massive black holes with orbital period less than ten years in their centers, detectable with on-going highly sensitive X-ray monitors such as Monitor of All-sky X-ray Image and/or Swift/Burst Alert Telescope. Assuming that all binary massive black holes have the equal mass ratio, about 20% of AGNs with black hole masses of 10^{6.5-7}M_sun has the close binaries and thus provides the best chance to detect them.
We study an accretion flow during the gravitational-wave driven evolution of binary massive black holes. After the binary orbit decays due to an interaction with a massive circumbinary disk, the binary is decoupled from the circumbinary disk because the orbital-decay timescale due to emission of gravitational wave becomes shorter than the viscous timescale evaluated at the inner edge of circumbinary disk. During the subsequent evolution, the accretion disk, which is truncated at the tidal radius because of the tidal torque, also shrinks as the orbital decay. Assuming that the disk mass changed by this process is all accreted, the disk becomes radiatively inefficient when the semi-major axis is several hundred Schwarzschild radii. The high-energy radiations, in spite of a low bolometric luminosity, are emitted from an accretion disk around each black hole long before the black hole coalescence as well as the gravitational wave signals. The synchrotron process can notably produce potentially observable radio emissions at large distances if there is a strong, dipole magnetic field around each black hole. In unequal mass-ratio binaries, step-like light variations are seen in the observed light curve because the luminosity is higher and its duration time are shorter in the radio emission by the disk around the secondary black hole than those of the primary black hole. Such a precursor would be unique to not a single black hole system but a binary black hole system, and implies that binary black holes finally merge without accretion disks.
It is still unknown how the BBH evolves after its semi-major axis reached to the sub-parsec/parsec scale where the dynamical friction with the neighboring stars is no longer effective (the so-called the final parsec problem). In this paper, we propose a new mechanism by which the massive BBH can naturally coalesce within a Hubble time. We study the evolution of the BBH with triple disks which are composed of an accretion disk around each black hole and one circumbinary disk surrounding them. While the circumbinary disk removes the orbital angular momentum of the BBH via the binary-disk resonant interaction, the mass transfer from the circumbinary disk to each black hole adds some fraction of its angular momentum to the orbital angular momentum of the BBH. We find that there is a critical value of the mass-transfer rate where the extraction of the orbital angular momentum from the BBH is balanced with the addition of the orbital angular momentum to the BBH. The semi-major axis of the BBH decays with time whereas the orbital eccentricity of the BBH grows with time, if the mass transfer rate is smaller than the critical one, and vice versa. Its evolutionary timescale is characterized by the product of the viscous timescale of the circumbinary disk and the ratio of the total black hole mass to the mass of the circumbinary disk. Since a minimum value of the critical mass-transfer rate is larger than the Eddington accretion rate of massive black holes with masses in the 10^{6}M_{sun} to $10^{9}M_{sun} range, it is promising that the critical mass-transfer rate is larger than the mass transfer rate. Most of massive BBHs, therefore, enable to merge within a Hubble time by the proposed mechanism, which helps to solve the final parsec problem.
We study the accretion flows from the circumbinary disks onto the supermassive binary black holes in a subparsec scale of the galactic center, using a smoothed particles hydrodynamics (SPH) code. Simulation models are presented in four cases of a circular binary with equal and unequal masses, and of an eccentric binary with equal and unequal masses. We find that the circumblack-hole disks are formed around each black holes regardless of simulation parameters. There are two-step mechanisms to cause an accretion flow from the circumbinary disk onto supermassive binary black holes: First, the tidally induced elongation of the circumbinary disk triggers mass inflow towards two closest points on the circumbinary disk from the black holes. Then, the gas is increasingly accumulated on these two points owing to the gravitational attraction of black holes. Second, when the gas can pass across the maximum loci of the effective binary potential, it starts to overflow via their two points and freely infalls to each black hole. In circular binaries, the gas continues to be supplied from the circumbinary disk (i.e. the gap between the circumbinary disk and the binary black hole is always closed.) In eccentric binaries, the mass supply undergoes the periodic on/off transitions during one orbital period because of the variation of periodic potential. The gap starts to close after the apastron and to open again after the next periastron passage. Due to this gap closing/opening cycles, the mass-capture rates are eventually strongly phase dependent. This could provide observable diagnosis for the presence of supermassive binary black holes in merged galactic nuclei.