Showing 1–43 of 43 results for author: Hayasaki, K

Revealing Physical Properties of a Tidal Disruption Event: iPTF16fnl

Authors: T. Mageshwaran, Gargi Shaw, Sudip Bhattacharyya, Kimitake Hayasaki

Abstract: Tidal disruption event (TDE) iPTF16fnl shows a relatively low optical flare with observationally very weak X-ray emission and the spectroscopic property that the helium emission line from the source dominates over the hydrogen emission line at early times. We explore these observed signatures by calculating spectral emission lines with the publicly available code, CLOUDY. We estimate five physical… ▽ More Tidal disruption event (TDE) iPTF16fnl shows a relatively low optical flare with observationally very weak X-ray emission and the spectroscopic property that the helium emission line from the source dominates over the hydrogen emission line at early times. We explore these observed signatures by calculating spectral emission lines with the publicly available code, CLOUDY. We estimate five physical parameters by fitting the observed optical UV spectra on multiple days to a theoretical model of a steady-state, slim disk with a spherical outflow. The resultant key parameters among them are black hole mass $M_{\bullet} = (6.73 \pm 0.44) \times 10^5 M_{\odot}$, stellar mass $M_{\star} = (2.59 \pm 0.17) M_{\odot}$, and wind velocity $v_{\rm w} = 7447.43 \pm 183.9~{\rm km~s^{-1}}$. The disk-wind model also estimates the radiative efficiency to be $0.01\lesssimη\lesssim0.02$ over the observational time, resulting in the disk being radiatively inefficient, and the disk X-ray luminosity is consistent with the observed low luminosity. In our CLOUDY model, the filling factor of the wind is also estimated to be 0.8, suggesting that the wind is moderately clumpy. We reveal that the helium-to-hydrogen number density ratio of the wind lies between 0.1 and 0.15, which is nearly the same as the solar case, suggesting the tidally disrupted star is originally a main sequence star. Because the optical depth of the helium line is lower than the hydrogen line by two orders of magnitude, the helium line is significantly optically thinner than the hydrogen line. Consequently, our results indicate that the helium line luminosity dominates the hydrogen line luminosity due to the optical depth effect despite a small helium-to-hydrogen number density ratio value. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 15 pages, 8 figures, submitted to Physical Review D

Evolution of Tidal Disruption Event Disks with Magnetically Driven Winds

Authors: T. Mageshwaran, Kimitake Hayasaki, Takeru K. Suzuki

Abstract: We present a time-dependent, one-dimensional magnetically driven disk-wind model based on magnetohydrodynamical (MHD) equations in the context of tidal disruption events (TDEs). We assume that the disk is geometrically thin and gas-pressure dominant and explicitly take into account magnetic braking as well as turbulent viscosity by an extended alpha viscosity prescription. We find a particular win… ▽ More We present a time-dependent, one-dimensional magnetically driven disk-wind model based on magnetohydrodynamical (MHD) equations in the context of tidal disruption events (TDEs). We assume that the disk is geometrically thin and gas-pressure dominant and explicitly take into account magnetic braking as well as turbulent viscosity by an extended alpha viscosity prescription. We find a particular wind solution for a set of basic equations that satisfies the necessary and sufficient conditions for vertically unbound MHD flows. The solution demonstrates that the disk evolves with losing mass due to wind and accretion from the initial Gaussian density distribution. We confirm that the mass accretion rate follows the power law of time $t^{-19/16}$ at late times if the wind is absent, which corresponds to the classical solution of Cannizzo et al. (1990). We find that the mass accretion rate is steeper than the $t^{-19/16}$ curve if the disk wind is present. This is because the wind removes a significant fraction of the mass and angular momentum. Mass accretion is further induced by magnetic braking, known as a wind-driven accretion mechanism, resulting in more rapid decay with time in both the mass accretion and loss rates. The ultraviolet (UV) luminosity is the highest among the optical, UV, and X-ray luminosities from early to late evolutionary phases, suggesting optical and X-ray emissions from the disk are observationally insignificant due to magnetically driven winds in TDEs. Our model predicts that late-time bolometric light curves steeper than $t^{-19/16}$ in UV-bright TDEs potentially serve as compelling indicators for magnetically driven winds. △ Less

Submitted 24 December, 2023; originally announced December 2023.

Comments: 20 pages, 13 figures, submitted to Physical Review D

X-ray spectral variations of Circinus X-1 observed with NICER throughout an entire orbital cycle

Authors: Mayu Tominaga, Masahiro Tsujimoto, Ken Ebisawa, Teruaki Enoto, Kimitake Hayasaki

Abstract: Circinus X-1 (Cir X-1) is a neutron star binary with an elliptical orbit of 16.6~days. The source is unique for its extreme youth, providing a key to understanding early binary evolution. However, its X-ray variability is too complex to reach a clear interpretation. We conducted the first high cadence (every 4 hours on average) observations covering one entire orbit using the NICER X-ray telescope… ▽ More Circinus X-1 (Cir X-1) is a neutron star binary with an elliptical orbit of 16.6~days. The source is unique for its extreme youth, providing a key to understanding early binary evolution. However, its X-ray variability is too complex to reach a clear interpretation. We conducted the first high cadence (every 4 hours on average) observations covering one entire orbit using the NICER X-ray telescope. The X-ray flux behavior can be divided into stable, dip, and flaring phases. The X-ray spectra in all phases can be described by a common model consisting of a partially covered disk black body emission and the line features from a highly-ionized photo-ionized plasma. The spectral change over the orbit is attributable to rapid changes of the partial covering medium in the line-of-sight and gradual changes of the disk black body emission. Emission lines of the H-like and He-like Mg, Si, S, and Fe are detected, most prominently in the dip phase. The Fe emission lines change to absorption in the course of the transition from the dip phase to the flaring phase. The estimated ionization degree indicates no significant changes, suggesting that the photo-ionized plasma is stable over the orbit. We propose a simple model in which the disk black body emission is partially blocked by such a local medium in the line-of-sight that has spatial structures depending on the azimuth of the accretion disk. Emission lines upon the continuum emission are from the photo-ionized plasma located outside of the blocking material. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ, 16 pages, 7 figures

Impact of scale-height derivative on general relativistic slim disks in tidal disruption events

Authors: T. Mageshwaran, Kimitake Hayasaki

Abstract: We construct a numerical model of steady-state, general relativistic (GR) super-Eddington accretion flows in an optically thick, advection-dominated regime, motivated by tidal disruption events wherein super-Eddington accretion assumes a pivotal role. Our model takes into account the loss of angular momentum due to radiation and the scale-height derivative in the basic equations of the GR slim dis… ▽ More We construct a numerical model of steady-state, general relativistic (GR) super-Eddington accretion flows in an optically thick, advection-dominated regime, motivated by tidal disruption events wherein super-Eddington accretion assumes a pivotal role. Our model takes into account the loss of angular momentum due to radiation and the scale-height derivative in the basic equations of the GR slim disk. For comparison purposes, we also provide a new analytical solution for a radiation-pressure-dominant GR slim disk, which neglects the angular momentum loss due to radiation and the scale-height derivative. We find that the radiation pressure enhances by incorporating the scale height derivative into the basic equations. As a result, the surface density near the disk's inner edge decreases, whereas the disk temperature and scale height increase, brightening the disk spectrum in the soft X-ray waveband. Notably, an extremely high mass accretion rate significantly enhances the effect of the scale-height derivative, affecting the entire disk. In contrast, the inclusion of the radiation-driven angular momentum loss only slightly influences the disk surface density and temperature compared with the case of the scale-height derivative inclusion. The X-ray luminosity increases significantly due to scale height derivative for $\dot{M}/\dot{M}_{\rm Edd} \gtrsim 2$. In addition, the increment is higher for the non-spinning black hole than the spinning black hole case, resulting in a one-order of magnitude difference for $\dot{M}/\dot{M}_{\rm Edd}\gtrsim100$. We conclude that incorporating the scale-height derivative into a GR slim disk model is crucial as it impacts the disk structure and its resultant spectrum, particularly on a soft-X-ray waveband. △ Less

Submitted 2 August, 2023; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 21 pages, 15 figures, 1 table, accepted for publication in Physical Review D

Disk-wind-driven Expanding Radio-emitting Shell in Tidal Disruption Events

Authors: Kimitake Hayasaki, Ryo Yamazaki

Abstract: We study the evolution of a non-relativistically expanding thin shell in radio-emitting tidal disruption events (TDEs) based on a one-dimensional spherically symmetric model considering the effect of both a time-dependent mass loss rate of the disk wind and the ambient mass distribution. The analytical solutions are derived in two extreme limits: one is the approximate solution near the origin in… ▽ More We study the evolution of a non-relativistically expanding thin shell in radio-emitting tidal disruption events (TDEs) based on a one-dimensional spherically symmetric model considering the effect of both a time-dependent mass loss rate of the disk wind and the ambient mass distribution. The analytical solutions are derived in two extreme limits: one is the approximate solution near the origin in the form of the Taylor series, and the other is the asymptotic solution in which the ambient matter is dominant far away from the origin. Our numerical solutions are confirmed to agree with the respective analytical solutions. We find that no simple power-law of time solution exists in early to middle times because the mass loss rate varies over time, affecting the shell dynamics. We also discuss the application of our model to the observed radio-emitting TDE, AT2019dsg. △ Less

Submitted 6 July, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 28 pages, 6 figures, Accepted for publication in ApJ

Establishing accretion flares from massive black holes as a source of high-energy neutrinos

Authors: Sjoert van Velzen, Robert Stein, Marat Gilfanov, Marek Kowalski, Kimitake Hayasaki, Simeon Reusch, Yuhan Yao, Simone Garrappa, Anna Franckowiak, Suvi Gezari, Jakob Nordin, Christoffer Fremling, Yashvi Sharma, Lin Yan, Erik C. Kool, Daniel Stern, Patrik M. Veres, Jesper Sollerman, Pavel Medvedev, Rashid Sunyaev, Eric C. Bellm, Richard G. Dekany, Dimitri A. Duev, Matthew J. Graham, Mansi M. Kasliwal , et al. (4 additional authors not shown)

Abstract: The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an… ▽ More The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6$σ$ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ~10% of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares. △ Less

Submitted 3 April, 2024; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society (2024), Volume 529, Issue 3, 2559-2576

Broadband X-Ray Observations of the 2018 Outburst of the Changing-Look Active Galactic Nucleus NGC 1566

Authors: Arghajit Jana, Neeraj Kumari, Prantik Nandi, Sachindra Naik, Arka Chatterjee, Gaurava K. Jaisawal, Kimitake Hayasaki, Claudio Ricci

Abstract: We study the nature of the changing-look Active Galactic Nucleus NGC 1566 during its June 2018 outburst. During the outburst, the X-ray intensity of the source rises up to ~25-30 times compared to its quiescent state intensity. We perform timing and spectral analysis of the source during pre-outburst, outburst, and post-outburst epochs using semi-simultaneous observations with the XMM-Newton, Nucl… ▽ More We study the nature of the changing-look Active Galactic Nucleus NGC 1566 during its June 2018 outburst. During the outburst, the X-ray intensity of the source rises up to ~25-30 times compared to its quiescent state intensity. We perform timing and spectral analysis of the source during pre-outburst, outburst, and post-outburst epochs using semi-simultaneous observations with the XMM-Newton, Nuclear Spectroscopic Telescope Array (NuSTAR), and Neil Gehrels Swift Observatories. We calculate variance, normalized variance, and fractional rms amplitude in different energy bands to study the variability. The broad-band 0.5-70 keV spectra are fitted with phenomenological models, as well as physical models. A strong soft X-ray excess is detected in the spectra during the outburst. The soft excess emission is found to be complex and could originate in the warm Comptonizing region in the inner accretion disc. We find that the increase in the accretion rate is responsible for the sudden rise in luminosity. This is supported by the q-shape of the hardness-intensity diagram that is generally found in outbursting black hole X-ray binaries. From our analysis, we find that NGC 1566 most likely harbours a low-spinning black hole with the spin parameter a* ~ 0.2. We also discuss a scenario where the central core of NGC 1566 could be a merging supermassive black hole. △ Less

Submitted 2 August, 2021; v1 submitted 23 July, 2021; originally announced July 2021.

Comments: 18 pages, 10 figures, Accepted for publication in MNRAS

On the origin of late-time X-ray flares in UV/optically-selected tidal disruption events

Authors: Kimitake Hayasaki, Peter G. Jonker

Abstract: We propose a model to explain the time delay between the peak of the optical and X-ray luminosity, \dt hereafter, in UV/optically-selected tidal disruption events (TDEs). The following picture explains the observed \dt in several TDEs as a consequence of the circularization and disk accretion processes as long as the sub-Eddington accretion. At the beginning of the circularization, the fallback de… ▽ More We propose a model to explain the time delay between the peak of the optical and X-ray luminosity, \dt hereafter, in UV/optically-selected tidal disruption events (TDEs). The following picture explains the observed \dt in several TDEs as a consequence of the circularization and disk accretion processes as long as the sub-Eddington accretion. At the beginning of the circularization, the fallback debris is thermalized by the self-crossing shock caused by relativistic precession, providing the peak optical emission. During the circularization process, the mass fallback rate decreases with time to form a ring around the supermassive black hole (SMBH). The formation timescale corresponds to the circularization timescale of the most tightly bound debris, which is less than a year to several decades, depending mostly on the penetration factor, the circularization efficiency, and the black hole mass. The ring will subsequently evolve viscously over the viscous diffusion time. We find that it accretes onto the SMBH on a fraction of the viscous timescale, which is $2$ years for given typical parameters, leading to X-ray emission at late times. The resultant \dt\,is given by the sum of the circularization timescale and the accretion timescale and significantly decreases with increasing penetration factor to several to $\sim10$ years typically. Since the X-ray luminosity substantially decreases as the viewing angle between the normal to the disk plane and line-of-sight increases from $0^\circ$ to $90^\circ$, a low late-time X-ray luminosity can be explained by an edge-on view. We also discuss the super-Eddington accretion scenario, where \dt\,is dominated by the circularization timescale. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: 21 pages, 4 figures, accepted for publication in ApJ

Neutrinos from tidal disruption events

Authors: Kimitake Hayasaki

Abstract: 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. 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. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: Published in Nature Astronomy, News & Views

Journal ref: Nat Astron (2021)

Exploring the origin of stars on bound and unbound orbits causing tidal disruption events

Authors: Shiyan Zhong, Kimitake Hayasaki, Shuo Li, Peter Berczik, Rainer Spurzem

Abstract: Tidal disruption events (TDEs) provide a clue to the properties of a central supermassive black hole (SMBH) and an accretion disk around it, and to the stellar density and velocity distributions in the nuclear star cluster surrounding the SMBH. Deviations of TDE light curves from the standard occurring at a parabolic encounter with the SMBH depends on whether the stellar orbit is hyperbolic or ecc… ▽ More Tidal disruption events (TDEs) provide a clue to the properties of a central supermassive black hole (SMBH) and an accretion disk around it, and to the stellar density and velocity distributions in the nuclear star cluster surrounding the SMBH. Deviations of TDE light curves from the standard occurring at a parabolic encounter with the SMBH depends on whether the stellar orbit is hyperbolic or eccentric (Hayasaki et al. 2018) and the penetration factor ($β$, tidal disruption radius to orbital pericenter ratio). We study the orbital parameters of bound and unbound stars being tidally disrupted by comparison of direct $N$-body simulation data with an analytical model. Starting from the classical steady-state Fokker-Planck model of Cohn & Kulsrud (1978), we develop an analytical model of the number density distribution of those stars as a function of orbital eccentricity ($e$) and $β$. To do so fittings of the density and velocity distribution of the nuclear star cluster and of the energy distribution of tidally disrupted stars are required and obtained from $N$-body data. We confirm that most of the stars causing TDEs in a spherical nuclear star cluster originate from the full loss-cone region of phase space, derive analytical boundaries in eccentricity-$β$ space, and find them confirmed by $N$-body data. Since our limiting eccentricities are much smaller than critical eccentricities for full accretion or full escape of stellar debris, we conclude that those stars are only very marginally eccentric or hyperbolic, close to parabolic. △ Less

Submitted 1 December, 2023; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 28 pages, 10 figures, 1 table. Accepted for publication in ApJ

Tidal Disruption Flares from Stars on Marginally Bound and Unbound Orbits

Authors: Gwanwoo Park, Kimitake Hayasaki

Abstract: 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 eccentri… ▽ More 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. △ Less

Submitted 20 June, 2020; v1 submitted 13 January, 2020; originally announced January 2020.

Comments: 42 pages 11 figures, accepted for publication in ApJ

Ionization and dissociation induced fragmentation of a tidally disrupted star into planets around a supermassive black hole

Authors: Kimitake Hayasaki, Matthew R. Bate, Abraham Loeb

Abstract: 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 process… ▽ More 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. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 13 pages, 4 figures, submitted for publication, comments welcome

Neutrino Emissions from Tidal Disruption Remnants

Authors: Kimitake Hayasaki, Ryo Yamazaki

Abstract: 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… ▽ More 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. △ Less

Submitted 26 November, 2019; v1 submitted 28 August, 2019; originally announced August 2019.

Comments: 41 pages, 2 figures, published in ApJ, 886 114 (2019)

Classification of Tidal Disruption Events Based on Stellar Orbital Properties

Authors: Kimitake Hayasaki, Shiyan Zhong, Shuo Li, Peter Berczik, Rainer Spurzem

Abstract: 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 ty… ▽ More 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. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 37 pages, 2 tables, 7 figures, accepted for publication in ApJ

Optical-infrared flares and radio afterglows by Jovian planets inspiraling into their host stars

Authors: Ryo Yamazaki, Kimitake Hayasaki, Abraham Loeb

Abstract: 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… ▽ More 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. △ Less

Submitted 7 December, 2016; v1 submitted 6 December, 2015; originally announced December 2015.

Comments: 8 pages, 2 tables, MNRAS in press

Journal ref: MNRAS, vol.466, p.1421-1427 (2017)

Detection of Gravitational Wave Emission by Supermassive Black Hole Binaries Through Tidal Disruption Flares

Authors: Kimitake Hayasaki, Abraham Loeb

Abstract: 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 ho… ▽ More 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. △ Less

Submitted 17 October, 2016; v1 submitted 20 October, 2015; originally announced October 2015.

Comments: 22pages, 7 figures, accepted for publication on Oct. 21, 2016 in Scientific Reports

Warping and tearing of misaligned circumbinary disks around eccentric SMBH binaries

Authors: Kimitake Hayasaki, Bong Won Sohn, Atsuo T. Okazaki, Taehyun Jung, Guangyao Zhao, Tsuguya Naito

Abstract: 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 t… ▽ More 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. △ Less

Submitted 31 January, 2015; originally announced February 2015.

Comments: 17 pages, 6 figures, submitted to JCAP. arXiv admin note: substantial text overlap with arXiv:1406.2317

Circularization of Tidally Disrupted Stars around Spinning Supermassive Black Holes

Authors: Kimitake Hayasaki, Nicholas C. Stone, Abraham Loeb

Abstract: 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… ▽ More 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. △ Less

Submitted 7 September, 2016; v1 submitted 21 January, 2015; originally announced January 2015.

Comments: 23 pages, 18 figures, 1 appendix, accepted for publication in MNRAS (with significant improvement)

Radiation-Driven Warping of Circumbinary Disks Around Eccentric Young Star Binaries

Authors: Kimitake Hayasaki, Bong Won Sohn, Atsuo T. Okazaki, Taehyun Jung, Guangyao Zhao, Tsuguya Naito

Abstract: 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 v… ▽ More 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. △ Less

Submitted 29 October, 2014; originally announced October 2014.

Comments: 28 pages, 5 figures, 1 table, accepted for publication in ApJ. arXiv admin note: substantial text overlap with arXiv:1406.2317

Warped Circumbinary Disks in Active Galactic Nuclei

Authors: Kimitake Hayasaki, Bong Won Sohn, Atsuo T. Okazaki, Taehyun Jung, Guangyao Zhao, Tsuguya Naito

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2014; originally announced June 2014.

Comments: 23 pages, 4 figures, 2 tables, accepted for publication in ApJ

Origin of Two Types of X-Ray Outbursts in Be/X-Ray Binaries. I. Accretion Scenarios

Authors: Atsuo T. Okazaki, Kimitake Hayasaki, Yuki Moritani

Abstract: 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 w… ▽ More 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. △ Less

Submitted 22 November, 2012; originally announced November 2012.

Comments: 16 pages, 10 figures, accepted for publication in PASJ

Binary Black Hole Accretion Flows From a Misaligned Circumbinary Disk

Authors: Kimitake Hayasaki, Hideki Saito, Shin Mineshige

Abstract: 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 bi… ▽ More 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. △ Less

Submitted 21 November, 2012; originally announced November 2012.

Comments: 12pages, 5figures, submitted to PASJ

Finite, Intense Accretion Bursts from Tidal Disruption of Stars on Bound Orbits

Authors: Kimitake Hayasaki, Nicholas Stone, Abraham Loeb

Abstract: 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, a… ▽ More 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. △ Less

Submitted 4 October, 2012; originally announced October 2012.

Comments: 15 pages, 15 figures, submitted to MNRAS

Tidal disruption flares from stars on eccentric orbits

Authors: Kimitake Hayasaki, Nicholas Stone, Abraham Loeb

Abstract: 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 stella… ▽ More 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. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Comments: 4 pages, 2 figures, Conference proceedings to appear in "Tidal Disruption and AGN outburst workshop" (EPJ Web of Conferences)

Modeling high-energy light curves of the PSR B1259-63/LS 2883 binary based on 3-D SPH simulations

Authors: J. Takata, A. T. Okazaki, S. Nagataki, T. Naito, A. Kawachi, S. -H. Lee, M. Mori, K. Hayasaki, M. S. Yamaguchi, S. P. Owocki

Abstract: 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 t… ▽ More 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. △ Less

Submitted 13 March, 2012; v1 submitted 8 March, 2012; originally announced March 2012.

Comments: 32 pages, 6 figues. Accepted for publication in ApJ

Gravitational wave diagnosis of a circumbinary disk

Authors: Kimitake Hayasaki, Kent Yagi, Takahiro Tanaka, Shin Mineshige

Abstract: 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 under… ▽ More 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. △ Less

Submitted 13 January, 2012; originally announced January 2012.

Comments: 10 pages, 2 figures, submitted to ApJ

Hydrodynamic Interaction between the Be Star and the Pulsar in the TeV Binary PSR B1259-63/LS 2883

Authors: Atsuo T. Okazaki, Shigehiro Nagataki, Tsuguya Naito, Akiko Kawachi, Kimitake Hayasaki, Stanley P. Owocki, Jumpei Takata

Abstract: 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… ▽ More 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. △ Less

Submitted 15 May, 2011; v1 submitted 7 May, 2011; originally announced May 2011.

Comments: 10 pages, 5 figures, accepted for publication in PASJ

Radiatively inefficient accretion flows induced by gravitational-wave emission before massive black hole coalescence

Authors: Kimitake Hayasaki

Abstract: 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 circumbina… ▽ More 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. △ Less

Submitted 7 December, 2010; v1 submitted 1 September, 2010; originally announced September 2010.

Comments: 13 pages, 3 figures, accepted for publication in ApJL

Mass Function of Binary Massive Black Holes in Active Galactic Nuclei

Authors: Kimitake Hayasaki, Yoshihiro Ueda, Naoki Isobe

Abstract: 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 ci… ▽ More 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. △ Less

Submitted 17 August, 2010; v1 submitted 20 January, 2010; originally announced January 2010.

Comments: 22 pages, 11 figures, accepted for publication in PASJ. The draft was significantly revised. The major differences from the previous version are as follows: (1)The circumbinary disk is assumed to be a steady, axisymmetric, geometrically thin, self-gravitating, self-regulated but non-fragmenting. (2)The stellar scattering process is taken account of in the merging process of binary black holes

Rapid Merger of Binary Primordial Black Holes: an Implication for GW150914

Authors: Kimitake Hayasaki, Keitaro Takahashi, Yuuiti Sendouda, Shigehiro Nagataki

Abstract: 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… ▽ More 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. △ Less

Submitted 3 June, 2016; v1 submitted 9 September, 2009; originally announced September 2009.

Comments: 7 pages, no figure. Accepted for publication in PASJ

A new approach for probing circumbinary disks

Authors: Kimitake Hayasaki, Atsuo T. Okazaki

Abstract: 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 disk… ▽ More 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. △ Less

Submitted 1 January, 2009; v1 submitted 4 September, 2008; originally announced September 2008.

Comments: Accepted for publication in ApJ Letters (3/Dec/2008), 4pages, 1figure

A new mechanism for massive binary black hole evolution

Authors: Kimitake Hayasaki

Abstract: 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 ar… ▽ More 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. △ Less

Submitted 10 September, 2008; v1 submitted 22 May, 2008; originally announced May 2008.

Comments: Accepted for publication in PASJ(Publications of the Astronomical Society of Japan), 10pages, 8figures. High Resolution Version is Available at http://www2.yukawa.kyoto-u.ac.jp/~kimitake/bbhs.html

A supermassive binary black hole with triple disks

Authors: Kimitake Hayasaki, Shin Mineshige, Luis C. Ho

Abstract: 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… ▽ More 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. △ Less

Submitted 7 April, 2008; v1 submitted 19 August, 2007; originally announced August 2007.

Comments: Accepted for publication in the Astrophysical Journal, 16 pages, 11 figures. High Resolution Version is Available at http://www2.yukawa.kyoto-u.ac.jp/~kimitake/bbhs.html

Populating the Galaxy with close Be binaries

Authors: P. D. Kiel, J. R. Hurley, J. R. Murray, K. Hayasaki

Abstract: 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 h… ▽ More 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. △ Less

Submitted 7 March, 2007; originally announced March 2007.

Comments: 3 pages, no figs 1 table. Published in proceedings of Active OB stars: Laboratories for Stellar and Circumstellar Physics (Aug 29th to Sep 2, 2005; Hokkaido Uni, Sapporo, Japan)

Journal ref: Kiel P.D., Hurley J.R., Murray J.R. and Hayasaki K., 2007, in the Proceedings of Active OB stars: Laboratories for Stellar and Circumstellar Physics, ASP Conference Series Vol. 361, eds. S. Stefl, S.P. Owocki S.P. and A.T. Okazaki, p448

SPH Simulations of Accretion Flows onto the Supermassive Binary Black Holes from the Circumbinary Disks

Authors: Kimitake Hayasaki, Shin Mineshige, Hiroshi Sudou

Abstract: 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 significantl… ▽ More 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. △ Less

Submitted 19 December, 2006; originally announced December 2006.

Comments: Conference proceedings to appear in "The Central Engine of Active Galactic Nuclei", ed. L. C. Ho and J.-M. Wang (San Francisco: ASP)

Binary Black Hole Accretion Flows in Merged Galactic Nuclei

Authors: Kimitake Hayasaki, Shin Mineshige, Hiroshi Sudou

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 19 December, 2006; v1 submitted 6 September, 2006; originally announced September 2006.

Comments: 16 pages, 27 figures, 2 tables, accepted for publication in PASJ. "High Resolution Version is Available at "http://www2.yukawa.kyoto-u.ac.jp/~kimitake/bbhs.html" Three observational references are added. Grammatical errors and typos are corrected

Long-term evolution of accretion discs in Be/X-ray binaries

Authors: K. Hayasaki, Atsuo T. Okazaki

Abstract: 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… ▽ More 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. △ Less

Submitted 8 August, 2006; originally announced August 2006.

Comments: 10 pages, 7 figures, Accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.372:1140-1148,2006

SPH simulations of accretion flow via Roche lobe overflow and via mass transfer from Be disk

Authors: Kimitake Hayasaki, Atsuo T. Okazaki, James. R. Murray

Abstract: 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… ▽ More 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). △ Less

Submitted 17 November, 2005; originally announced November 2005.

Comments: 3 pages, 5 figures, to appear in "Active OB Stars: Laboratories for Stellar & Circumstellar Physics", ASP Conf. Ser. 2005, S. Stefl, S. P. Owocki & A. Okazaki, eds

Long-term evolution of accretion disks around the neutron star in Be/X-ray binaries

Authors: Kimitake Hayasaki, Atsuo T. Okazaki

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 17 November, 2005; originally announced November 2005.

Comments: 3 pages, 3 figure, to appear in "Active OB Stars: Laboratories for Stellar & Circumstellar Physics", ASP Conf. Ser. 2005, S. Stefl, S. P. Owocki & A. Okazaki, eds

One-armed spiral structure of accretion discs induced by a phase-dependent mass transfer in Be/X-ray binaries

Authors: Kimitake Hayasaki, Atsuo T. Okazaki

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2005; originally announced March 2005.

Comments: 5 pages, 7 figures. Accepted for publication in MNRAS Letters

Non-axisymmetric Structure of Accretion Disks in Be/X-ray Binaries

Authors: Kimitake Hayasaki, Atsuo T. Okazaki

Abstract: 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. T… ▽ More 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. △ Less

Submitted 9 December, 2004; originally announced December 2004.

Comments: 4pages, 6figures, Conference Proceeding for the "Interacting Binaries: Accretion, Evolution and Outcome" worshop held in Cefalu' (Italy) in July 2004

Accretion disc formation around the neutron star in Be/X-ray binaries

Authors: Kimitake Hayasaki, Atsuo T. Okazaki

Abstract: 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 accretio… ▽ More 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. △ Less

Submitted 11 February, 2004; v1 submitted 10 November, 2003; originally announced November 2003.

Comments: 13pages, 62figures, accepted to MNRAS. Revised version adds in 4 new figures, in which we have improved the English, kindly pointed out by Manson Katherine. For associated movies, see http://astro3.sci.hokudai.ac.jp/~kimi/movie.html

Journal ref: Mon.Not.Roy.Astron.Soc. 350 (2004) 971

Accretion onto the neutron star in Be/X-ray binaries

Authors: Kimitake Hayasaki, Atsuo. T. Okazaki

Abstract: 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 aft… ▽ More 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. △ Less

Submitted 15 October, 2003; originally announced October 2003.

Comments: 2 pages, 8 figures, To appear in "Young Neutron Stars and Their Environments" (IAU Symposium 218, ASP Conference Proceedings), eds F. Camilo and B. M. Gaensler