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 $\beta=1$ to $3$, where $\beta$ 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}\beta^{k-1}$ and $e_{\rm crit,2}=1+2q^{-1/3}\beta^{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 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 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_{\nu,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_{\nu,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 rates of tidal disruption of stars by intermediate-mass to supermassive black holes on bound to unbound orbits by using high-accuracy direct N-body experiments. The approaching stars from the star cluster to the black hole can take three types of orbit: eccentric, parabolic, and hyperbolic orbits. Since the mass fallback rate shows a different variability depending on these orbital types, we can classify tidal disruption events (TDEs) into three main categories: eccentric, parabolic, and hyperbolic TDEs. Respective TDEs are characterized by two critical values of the orbital eccentricity: the lower critical eccentricity is the one below which the stars on eccentric orbits cause the finite, intense accretion, and the higher critical eccentricity above which the stars on hyperbolic orbits cause no accretion. Moreover, we find that the parabolic TDEs are divided into three subclasses: precisely parabolic, marginally eccentric, and marginally hyperbolic TDEs. We analytically derive that the mass fallback rate of the marginally eccentric TDEs can be flatter and slightly higher than the standard fallback rate proportional to $t^{-5/3}$, whereas it can be flatter and lower for the marginally hyperbolic TDEs. We confirm by N-body experiments that only few eccentric, precisely parabolic, and hyperbolic TDEs can occur in a spherical stellar system with a single intermediate-mass to supermassive black hole. A substantial fraction of the stars approaching to the black hole would cause the marginally eccentric or marginally hyperbolic TDEs.

When a planet inspirals into its host star, it releases gravitational energy which is converted into an expanding bubble of hot plasma. We study the radiation from the bubble and show that it includes prompt optical-infrared emission and a subsequent radio afterglow. The prompt emission from M31 and Large Magellanic Cloud is detectable by optical-near infrared transient surveys with a large field of view. The subsequent radio afterglows are detectable for $10^{3-4}$~years. The event rate depends on uncertain parameters in the formation and dynamics of giant planets. Future observation of the rate will constrain related theoretical models. If the event rate is high (> a few events per year), the circumstellar disk must typically be massive as suggested by recent numerical simulations.

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 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.

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 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 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.

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.

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

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 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.

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.

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.

Circumbinary disks are considered to exist in a wide variety of astrophysical objects, e.g., young binary stars, protoplanetary systems, and massive binary black hole systems in active galactic nuclei (AGNs). However, there is no definite evidence for the circumbinary disk except for some in a few young binary star systems. In this Letter, we study possible oscillation modes in circumbinary disks around eccentric and circular binaries. We find that progarde, nonaxisymmetric waves are induced in the inner part of the circumbinary disk by the tidal potential of the binary. Such waves would cause variabilities in emission line profiles from circumbinary disks. Because of prograde precession of the waves, the distance between each component of the binary and the inner edge of the circumbinary disk varies with the beat period between the precession period of the wave and the binary orbital period. As a result, light curves from the circumbinary disks are also expected to vary with the same period. The current study thus provides a new method to detect circumbinary disks in various astrophysical systems.

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.

Hierarchical structure formation inevitably leads to the formation of supermassive binary black holes (BBHs) with a sub-parsec separation in galactic nuclei. However, to date there has been no unambiguous detection of such systems. In an effort to search for potential observational signatures of supermassive BBHs, we performed high-resolution smoothed particle hydrodynamics (SPH) simulations of two black holes in a binary of moderate eccentricity surrounded by a circumbinary disk. Building on our previous work, which has shown that gas can periodically transfer from the circumbinary disk to the black holes when the binary is on an eccentric orbit, the current set of simulations focuses on the formation of the individual accretion disks, their evolution and mutual interaction, and the predicted radiative signature. The variation in mass transfer with orbital phase from the circumbinary disk induces periodic variations in the light curve of the two accretion disks at ultraviolet wavelengths, but not in the optical or near-infrared. Searches for this signal offer a promising method to detect supermassive BBHs.

Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of how the evolution changes when the star is in close proximity to a companion. Within this work we complete a full population synthesis study of Be binaries for the Galaxy. The results for the first time match aspects of the observational data, most notably the observed upper limit to the period distribution. We conclude that greater detailed studies on the evolution of Be stars within X-ray binaries needs to be completed, so that rapid binary evolution population synthesis packages may best evolve these systems.

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.

We investigate the accretion flows onto the supermassive binary black holes (SMBBHs) from the circumbinary disk with the equal mass, eccentric binary on the subparsec scale, using Smoothed Particle Hydrodynamics (SPH) code. We find that the material can be supplied from circumbinary disk, which leads to the formation of two accretion disks around the SMBBHs. The mass accretion rates significantly modulate with the binary orbital motion. These could provide the observable diagnosis of the existence of the supermassive binary black holes (e.g. OJ287) on the subparsec scale in merged galactic nuclei.

We numerically study the long-term evolution of the accretion disc around the neutron star in a coplanar Be/X-ray binary with a short period and a moderate eccentricity. From three dimensional Smoothed Particle Hydrodynamics simulations, we find that the disc evolves through three distinct phases, each characterized by different mass accretion patterns. In the first "developing phase", the disc is formed and develops towards a nearly Keplerian disc. It has a relatively large, double-peaked mass-accretion rate with the higher peak by the direct accretion at periastron, which is followed by the lower peak by the accretion induced by a one-armed spiral wave. In the second "transition phase", the disc is approximately Keplerian and grows with time. The mass-accretion rate increases as the disc grows. In the second phase, there is a transition in the mass accretion rate from a double peaked to a single peaked pattern. In the final quasi-steady state, the mass-accretion rate is on average balanced with the mass-transfer rate from the Be disc and exhibits a regular orbital modulation. In the quasi-steady state, the mass-accretion rate has a single peak by the wave-induced accretion as in a later stage of the transition phase. The orbital modulation of X-ray maxima could provide not only a circumstantial evidence for the persistent disc but also an observational diagnosis of the disc evolutionary state.

we study the long-term evolution of the accretion disk around the neutron star in Be/X-ray binaries. We confirm the earlier result by Hayasaki & Okazaki (2004) that the disk evolves via a two-stage process, which consists of the initial developing stage and the later developed stage. The peak mass-accretion rate is distributed around apastron after the disk is fully developed. This indicates that the modulation of the mass accretion rate is essentially caused by an inward propagation of the one-armed spiral wave. The X-ray luminosity peak around the apastron could provide circumstatial evidence for an persistent disk around the neutron star in Be/X-ray binaries.

We compare the accretion flow onto the neutron star induced by Roche lobe overflow with that by the overflow from the Be disk, in a zero eccentricity, short period binary with the same mass transfer rate, performing three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a persistent accretion disk is formed around the neutron star in both cases. The circularization radius of the material transferred via Roche lobe overflow is larger than that of the material transfered from the Be disk. Thus, the growth of the accretion disk in the former case becomes significantly slower than in the latter case. In both cases, the mass accretion rate is very small and varies little with orbital phase, which is consistent with the observed X-ray behaviour of Be/X-ray binaries with circular orbits (e.g. XTE J1543-568).

We study non-axisymmetric structure of accretion discs in Be/X-ray binaries, performing three dimensional Smoothed Particle Hydrodynamics simulations for a coplanar system with a short period and a moderate eccentricity. We find that ram pressure due to the phase-dependent mass transfer from the Be-star disc excites a one-armed, trailing spiral structure in the accretion disc around the neutron star. The spiral wave is transient; it is excited around the periastron passage, when the material is transferred from the Be disc, and is gradually damped afterwards. The disc changes its morphology from circular to eccentric with the development of the spiral wave, and then from eccentric to circular with the decay of the wave during one orbital period. It turns out that the inward propagation of the spiral wave significantly enhances the mass-accretion rate onto the neutron star. Thus, the detection of an X-ray luminosity peak corresponding to the peak in enhanced mass-accretion rate provides circumstantial evidence that an accretion disc is present in Be/X-ray binaries.

The non-axisymmetric structure of accretion disks around the neutron star in Be/X-ray binaries is studied by analyzing the results from three dimensional (3D) Smoothed Particle Hydrodynamics (SPH) simulations. It is found that ram pressure due to the phase-dependent mass transfer from the Be-star disk excites a one-armed, trailing spiral structure in the accretion disk around the neutron star. The spiral wave has a transient nature; it is excited around the periastron, when the material is transferred from the Be disk, and is gradually damped afterwards. It is also found that the orbital phase-dependence of the mass-accretion rate is mainly caused by the inward propagation of the spiral wave excited in the accretion disk.

We study the accretion on to the neutron star in Be/X-ray binaries, using a 3D SPH code and the data imported from a simulation by \citet{oka2} for a coplanar system with a short period ($P_{\rm orb}=24.3 \rm{d}$) and a moderate eccentricity $(e=0.34)$, which targeted the Be/X-ray binary 4U 0115+63. For simplicity, we adopt the polytropic equation of state. We find that a time-dependent accretion disc is formed around the neutron star regardless of the simulation parameters. In the long term, the disc evolves via a two-stage process, which consists of the initial developing stage and the later developed stage. The developed disc is nearly Keplerian. In the short term, the disc structure modulates with the orbital phase. The disc shrinks at the periastron passage of the Be star and restores its radius afterwards. The accretion rate on to the neutron star is also phase dependent, but its peak is broader and much lower than that of the mass-transfer rate from the Be disc, unless the polytropic exponent is as large as 5/3. Our simulations show that the truncated Be disk model for Be/X-ray binaries is consistent with the observed X-ray behaviour of 4U 0115+63.

We study accretion onto the neutron star in Be/X-ray binaries, using a 3D SPH code and the data imported from a high resolution simulation by Okazaki et al.\ (2002) for a coplanar system with a short period ($P_{\rm orb}=24.3 \rm d$) and moderate eccentricity $(\rm{e=0.34})$. We find that a time-dependent accretion disk is formed around the neutron star in Be/X-ray binaries. The disk shrinks after the periastron passage of the Be star and restores its radius afterwards. Our simulations show that the truncated Be disk model for Be/X-ray binaries is consistent with the observed X-ray behavior.