Source author record

Ken Ohsuga

Ken Ohsuga appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2022arXiv

Component of energy flow from supercritical accretion disks around rotating stellar mass black holes

By performing two-dimensional axisymmetric general relativistic radiation magnetohydrodynamics simulations with spin parameter $a^*$ varying from -0.9 to 0.9, we investigate the dependence on the black hole spin of the energy flow from supercritical accretion disk around stellar mass black hole. It is found that optically and geometrically thick disks form near the equatorial plane, and a part of the disk matter is launched from the disk surface in all models. The gas ejection is mainly driven by the radiative force, but magnetic force cannot be neglected, when $|a^*|$ is large. The energy outflow efficiency (total luminosity normalized by $\dot{M}_{\rm in} c^2 $; $\dot{M}_{\rm in}$ and $c$ are the mass accretion rate at the event horizon and the light speed) is larger for rotating black holes than for non-rotating black holes. This is $0.7\%$ for $a^*=-0.7$, $0.3\%$ for $a^*=0$, and $5\%$ for $a^*=0.7$ for $\dot{M}_{\rm in} \sim 100L_{\rm Edd}/c^2$ ($L_{\rm Edd}$ is Eddington luminosity). Also, although the energy is mainly released by radiation when $a^* \sim 0$, the Poynting power increases with $|a^*|$ and exceeds the radiative luminosity for models with $a^* \geq 0.5$ and $a^* \leq -0.7$. The more the black hole rotates, the larger the power ratio of the kinetic luminosity to the isotropic luminosity tends to be. This implies that objects with large (small) power ratio may have rapidly (slowly) rotating black holes. Among ultraluminous X-ray sources, IC342 X-1, is a candidate with a rapidly rotating black hole.

preprint2022arXiv

Investigating the Disk-Jet Structure in M87 through Flux Separation in the Linear and Circular Polarization Images

For testing different electron temperature ($T_{\rm e}$) prescriptions in general relativistic magnetohydrodynamics (GRMHD) simulations through observations, we propose to utilize linear polarization (LP) and circular polarization (CP) images. We calculate the polarization images based on a semi-Magnetically Arrested Disk (MAD) GRMHD model for various $T_{\rm e}$ parameters, bearing M87 in mind. We find an LP-CP separation in the images of the low-$T_{\rm e}$ disk cases at 230 GHz; namely, the LP flux mainly originates from downstream of the jet and the CP flux comes from the counter-side jet, while the total intensity is maximum at the jet base. This can be understood as follows: although the LP flux is generated through synchrotron emission widely around the black hole, most of the LP flux from the jet base does not reach the observer, since it undergoes Faraday rotation ($\propto T_{\rm e}^{-2}$) when passing through the outer cold disk and is thus depolarized. Hence, only the LP flux from the downstream (not passing the cold dense plasmas) can survive. Meanwhile, the CP flux is generated from the LP flux by Faraday conversion ($\propto T_{\rm e}$) in the inner hot region. Stronger CP flux is thus observed from the counter-side jet. Moreover, the LP-CP separation is more enhanced at a lower frequency such as 86 GHz but is rather weak at 43 GHz, since the media in the latter case is optically thick for synchrotron self-absorption so that all the fluxes should come from the photosphere. The same is true for cases with higher mass accretion rates and/or larger inclination angles.

preprint2022arXiv

Structure of the super-Eddington outflow and itsimpact on the cosmological scale

It is one of the biggest issues in black hole (BH) astrophysics how to precisely evaluate BH feedback to its environments. Aiming at studying the unique gas dynamics of super-Eddington flow around supermassive black hole (SMBH) seeds at high redshift, we carried out axisymmetric two dimensional radiation hydrodynamic simulations by a nested simulation-box method. Here we divide the simulation box into the inner zone at $(2 - 3 \times 10^3) r_{\rm{Sch}}$ (with $r_{\rm Sch}$ being the Schwarzschild radius) and the outer zone at $(2\times 10^{3} - 3\times 10^6) r_{\rm{Sch}}$, with smooth connection of the physical quantities, such as gas density, velocity, and radiation energy. We start the calculation by injecting mass through the outer boundary of the inner zone at a constant rate of $\dot{M}_{\rm{inj}}=10^3L_{\rm{Edd}}/c^2$, where $L_{\rm{Edd}}$ is the Eddington luminosity and $c$ is the speed of light. A powerful outflow is generated in the innermost region and it propagates from the inner zone to the outer zone. The outflows are characterized by a velocity of 0.02$c$ (0.7$c$) and density of $10^{-17}$ ($10^{-19}$) g cm$^{-3}$ for near the edge-on (face-on) direction. The outflow is gradually accelerated as it travels by accepting radiation-pressure force. The final mass outflow rate at the outermost boundary is $\dot{M}_{\rm{out}}\sim 0.3 \times \dot{M}_{\rm{inj}}$. By extrapolating the outflow structure to a further larger scale, we find that the momentum and energy fluxes at $r \sim 0.1$ pc are $\sim 10-100 L_{\rm{Edd}}/c $ and $\sim 0.1-10 L_{\rm{Edd}}$, respectively. Moreover, we find that the impacts are highly anisotropic in the sense that larger impacts are given towards the face-on direction than in the edge-on direction. These results indicate that the BH feedback will more efficiently work on the interstellar medium than that assumed in the cosmological simulations.

preprint2021arXiv

Outflow from super-Eddington flow: where it originates from and how much impact it gives?

It is widely believed that super-Eddington accretion flow can produce powerful outflow, but where it originates from and how much mass and energy are carried away to which directions? To answer to these questions, we newly perform a large-box, two-dimensional radiation hydrodynamic simulation, paying special attention lest the results should depend on adopted initial and boundary conditions. We could achieve a quasi-steady state in an unprecedentedly large range, $r=2~r_{\rm S}$-$600~r_{\rm S}$ (with $r_{\rm S}$ being the Schwarzschild radius) from the black hole. The accretion rate onto the central $10 ~M_{\odot}$ black hole is $\dot{M}_{\rm BH} \sim 180 ~L_{\rm Edd}/c^{2}$, whereas the mass outflow rate is ${\dot M}_{\rm outflow} \sim 24 ~L_{\rm Edd}/c^2$ (where $L_{\rm Edd}$ and $c$ are the Eddington luminosity and the speed of light, respectively). The ratio (${\dot M}_{\rm outflow}/{\dot M}_{\rm BH} \sim 0.14$) is much less than those reported previously. By careful inspection we find that most of outflowing gas which reach the outer boundary originates from the region at $R\lesssim140~r_{\rm S}$, while gas at $140~r_{\rm S}$-$230 ~r_{\rm S}$ forms failed outflow. Therefore, significant outflow occurs inside the trapping radius $\sim 450 ~r_{\rm S}$. The mechanical energy flux (or mass flux) reaches its maximum in the direction of $\sim 15^\circ$ ($\sim 80^\circ$) from the rotation axis. The total mechanical luminosity is $L_{\rm mec}\sim 0.16~L_{\rm Edd}$, while the isotropic X-ray luminosity varies from $L_{\rm X}^{\rm ISO}\sim 2.9~L_{\rm Edd}$, (for a face-on observer) to $\sim 2.1~L_{\rm Edd}$ (for a nearly edge-on observer). The power ratio is $L_{\rm mec}/L_{\rm X}^{\rm ISO}\sim 0.05$-$0.08$, in good agreement with the observations of Ultra-Luminous X-ray sources surrounded by optical nebulae.

preprint2020arXiv

Polarization Imaging of M87 Jets by General Relativistic Radiative Transfer Calculation based on GRMHD Simulations

The spectacular images of the M87 black hole taken by the Event Horizon Telescope (EHT) have opened a new era of black hole research. One of the next issues is to take polarization images around the central black hole (BH). Since radio emission is produced by synchrotron process, polarization properties should vividly reflect the magnetic field structure at the jet base and thus provide good information regarding the magnetic mechanism of jet formation. With this kept in mind we perform general relativistic (GR) radiative transfer calculations of polarized light based on the GR magnetohydrodynamic (MHD) simulation data of accretion flow and outflow in M87, to obtain their linear and circular polarization images in the horizon-scale. We found that the linear polarization components from the jet base and inner accretion flow should experience Faraday rotation when passing through magnetized plasmas around the BH, thus sensitively depending on the BH spin. Through the comparison with intensity image at 1.3mm by EHT and the polarized degree and the rotation measure (RM) at 1.3mm with the Submillimeter Array, the model with the spin parameter of a=0.9M_BH (with M_BH being the BH mass) is favored over other models with a=0.5M_BH or 0.99M_BH, though we need further systematic studies for confirmation. We also find in low-temperature models clear ring-like image in the circular polarization map, which arises because of Faraday conversion of the linearly polarized synchrotron emission and is thus indicative of magnetic field. This occurs only when the emission region is threaded with well-ordered magnetic fields and hence no clear images are expected in high-temperature disk models, in which disk emission is appreciable. We will be able to elucidate the field configuration through the comparison between the simulated polarization images and future polarimetry with EHT and other VLBI observations.

preprint2020arXiv

Pulsed fraction of super-critical column accretion flows onto neutron stars: modeling of ultraluminous X-ray pulsars

We calculate the pulsed fraction (PF) of the super-critical column accretion flows onto magnetized neutron stars (NSs), of which the magnetic axis is misaligned with the rotation axis, based on the simulation results by Kawashima et al.(2016, PASJ, 68, 83). Here, we solve the geodesic equation for light in the Schwarzschild spacetime in order to take into account the light bending effect. The gravitational redshift and the relativistic doppler effect from gas motions of the accretion columns are also incorporated. The pulsed emission appears since the observed luminosity, which exceeds the Eddington luminosity for the stellar-mass black holes, periodically changes via precession of the column caused by the rotation of the NS. The PF tends to increase as $θ_{\rm obs}$ approaching to $θ_{\rm B}$, where $θ_{\rm obs}$ and $θ_{\rm B}$ are the observer's viewing angle and the polar angle of the magnetic axis measured from the rotation axis. The maximum PF is around 50 %. Also, we find that the PF becomes less than 5 % for $θ_{\rm obs} \lesssim 5^\circ$ or for $θ_{\rm B} \lesssim 5^\circ$. Our results are consistent with observations of ultraluminous X-ray pulsars (ULXPs) with few exceptions, since the ULXPs mostly exhibit the PF of $\lesssim$ 50 %. Our present study supports the hypothesis that the ULXPs are powered by the super-critical column accretion onto NSs.

preprint2020arXiv

The thermal-radiative wind in low mass X-ray binary H 1743-322: II. iron line predictions from Monte Carlo radiation transfer

We show the best current simulations of the absorption and emission features predicted from thermal-radiative winds produced from X-ray illumination of the outer accretion disc in binary systems. We use the density and velocity structure derived from a radiation hydrodynamic code as input to a Monte-Carlo radiation transport calculation. The initial conditions are matched to those of the black hole binary system H1743-322 in its soft, disc dominated state, where wind features are seen in Chandra grating data. Our simulation fits well to the observed line profile, showing that these physical wind models can be the origin of the absorption features seen, rather than requiring a magnetically driven wind. We show how the velocity structure is the key observable discriminator between magnetic and thermal winds. Magnetic winds are faster at smaller radii, whereas thermal winds transition to a static atmosphere at smaller radii. New data from XRISM (due for launch Jan 2022) will give an unprecedented view of the physics of the wind launch and acceleration processes, but the existence of static atmospheres in small disc systems already rules out magnetic winds which assume self-similar magnetic fields from the entire disc as the origin of the absorption features seen.

preprint2020arXiv

The thermal-radiative wind in the neutron star low mass X-ray binary GX 13+1

We fit the observed high ionisation X-ray absorption lines in the neutron star binary GX13+1 with a full simulation of a thermal-radiative wind. This uses a radiation hydrodynamic code coupled to Monte Carlo radiation transfer to compute the observed line profiles from Hydrogen and Helium-like iron and Nickel, including all strong Kα and Kβ transitions. The wind is very strong as this object has a very large disc and is very luminous. The absorption lines from Fe Kα are strongly saturated as the ion columns are large, so the line equivalent widths (EWs) depend sensitively on the velocity structure. We additionally simulate the lines including isotropic turbulence at the level of the azimuthal and radial velocities. We fit these models to the Fe xxv and xxvi absorption lines seen in the highest resolution Chandra third order HETGS data. These data already rule out the addition of turbulence at the level of the radial velocity of ~500 km/s. The velocity structure predicted by the thermal-radiative wind alone is a fairly good match to the observed profile, with an upper limit to additional turbulence at the level of the azimuthal velocity of ~100 km/s. This gives stringent constraints on any remaining contribution from magnetic acceleration.

preprint2016arXiv

A radiation-hydrodynamic model of accretion columns for ultra-luminous X-ray pulsars

Prompted by the recent discovery of pulsed emission from an ultra-luminous X-ray source, M82 X-2 ("ULX-pulsar"), we perform a two-dimensional radiation-hydrodynamic simulation of a super-critical accretion flow onto a neutron star through a narrow accretion column. We set an accretion column with a cone shape filled with tenuous gas with density of $10^{-4} {\rm g}~ {\rm cm}^{-3}$ above a neutron star and solve the two dimensional gas motion and radiative transfer within the column. The side boundaries are set such that radiation can freely escape, while gas cannot. Since the initial gas layer is not in a hydrostatic balance, the column gas falls onto the neutron-star surface, thereby a shock being generated. As a result, the accretion column is composed of two regions: an upper, nearly free-fall region and a lower settling region, as was noted by Basko \& Sunyaev (1976). The average accretion rate is very high; ${\dot M}\sim 10^{2-3} L_{\rm E}/c^2$ (with $L_{\rm E}$ being the Eddington luminosity), and so radiation energy dominates over gas internal energy entirely within the column. Despite the high accretion rate, the radiation flux in the laboratory frame is kept barely below $L_{\rm E}/(4πr^2)$ at a distance $r$ in the settling region so that matter can slowly accrete. This adjustment is made possible, since large amount of photons produced via dissipation of kinetic energy of matter can escape through the side boundaries. The total luminosity can greatly exceed $L_{\rm E}$ by several orders of magnitude, whereas the apparent luminosity observed from the top of the column is much less. Due to such highly anisotropic radiation fields, observed flux should exhibit periodic variations with the rotation period, provided that the rotation and magnetic axes are misaligned.

preprint2016arXiv

Formation of Overheated Regions and Truncated Disks around Black Holes; Three-dimensional General Relativistic Radiation-magnetohydrodynamics Simulations

Using three-dimensional general relativistic radiation magnetohydrodynamics simulations of accretion flows around stellar mass black holes, we report that the relatively cold disk ($\gtrsim 10^{7}$K) is truncated near the black hole. Hot and less-dense regions, of which the gas temperature is $ \gtrsim 10^9$K and more than ten times higher than the radiation temperature (overheated regions), appear within the truncation radius. The overheated regions also appear above as well as below the disk, and sandwich the cold disk, leading to the effective Compton upscattering. The truncation radius is $\sim 30 r_{\rm g}$ for $\dot{M} \sim L_{\rm Edd}/c^2$, where $r_{\rm g}, \dot M, L_\mathrm{Edd}, c$ are the gravitational radius, mass accretion rate, Eddington luminosity, and light speed. Our results are consistent with observations of very high state, whereby the truncated disk is thought to be embedded in the hot rarefied regions. The truncation radius shifts inward to $\sim 10 r_{\rm g}$ with increasing mass accretion rate $\dot{M} \sim 100 L_{\rm Edd}/c^2$, which is very close to an innermost stable circular orbit. This model corresponds to the slim disk state observed in ultra luminous X-ray sources. Although the overheated regions shrink if the Compton cooling effectively reduces the gas temperature, the sandwich-structure does not disappear at the range of $\dot{M} \lesssim 100L_{\rm Edd}/c^2$. Our simulations also reveal that the gas temperature in the overheated regions depends on black hole spin, which would be due to efficient energy transport from black hole to disks through the Poynting flux, resulting gas heating.

preprint2015arXiv

Comparison between RHD simulation of supercritical accretion flows and steady model with outflows

We apply our two-dimensional (2D), radially self-similar steady-state accretion flow model to the analysis of hydrodynamic simulation results of supercritical accretion flows. Self-similarity is checked and the input parameters for the model calculation, such as advective factor and heat capacity ratio, are obtained from time-averaged simulation data. Solutions of the model are then calculated and compared with the simulation results. We find that in the converged region of the simulation, excluding the part too close to the black hole, the radial distribution of azimuthal velocity $v_ϕ$, density $ρ$ and pressure $p$ basically follows the self-similar assumptions, i.e. they are roughly proportional to $r^{-0.5}$, $r^{-n}$, and $r^{-(n+1)}$, respectively, where $n\sim0.85$ for the mass injection rate of $1000L_\mathrm{E}/c^2$, and $n\sim0.74$ for $3000L_\mathrm{E}/c^2$. The distribution of $v_r$ and $v_θ$ agrees less with self-similarity, possibly due to convective motions in the $rθ$ plane. The distribution of velocity, density and pressure in $θ$ direction obtained by the steady model agrees well with the simulation results {within the calculation boundary of the steady model}. Outward mass flux in the simulations is overall directed toward polar angle of 0.8382 rad ($\sim 48.0^\circ$) for $1000L_\mathrm{E}/c^2$, and 0.7852 rad ($\sim 43.4^\circ$) for $3000L_\mathrm{E}/c^2$, and $\sim$94\% of the mass inflow are driven away as outflow, while outward momentum and energy fluxes are focused around the polar axis. Part of these fluxes lie in the region that are not calculated by the steady model, and special attention should be paid when the model is applied.

preprint2015arXiv

Radiation Hydrodynamic Simulations of Line-Driven Disk Winds for Ultra Fast Outflows

Using two-dimensional radiation hydrodynamic simulations, we investigate origin of the ultra fast outflows (UFOs) that are often observed in luminous active galactic nuclei (AGNs). We found that the radiation force due to the spectral lines generates strong winds (line-driven disk winds) that are launched from the inner region of accretion disks (~30 Schwarzschild radii). A wide range of black hole masses ($M_{\rm BH}$) and Eddington ratios ($\varepsilon$) was investigated to study conditions for causing the line-driven winds. For $M_{\rm BH} = 10^6-10^9 M_\odot$ and $\varepsilon = 0.1-0.7$, funnel-shaped disk winds appear, in which dense matter is accelerated outward with an opening angle of 70-80 deg and with 10% of the light speed. If we observe the wind along its direction, the velocity, the column density, and the ionization state are consistent with those of the observed UFOs. As long as the obscuration by the torus does not affect the observations of X-ray bands, the UFOs could be statistically observed in about 13-28% of the luminous AGNs, which is not inconsistent with the observed ratio (~40%). We also found that the results are insensitive to the X-ray luminosity and the density of the disk surface. Thus, we can conclude that the UFOs could exist in any luminous AGNs, such as narrow-line Seyfert 1s (NLS1s) and quasars with $\varepsilon > 0.1$, in which fast line-driven winds are associated.

preprint2014arXiv

Discovery of Dramatic Optical Variability in SDSS J1100+4421: A Peculiar Radio-Loud Narrow-Line Seyfert 1 Galaxy?

We present our discovery of dramatic variability in SDSS J1100+4421 by the high-cadence transient survey Kiso Supernova Survey (KISS). The source brightened in the optical by at least a factor of three within about half a day. Spectroscopic observations suggest that this object is likely a narrow-line Seyfert 1 galaxy (NLS1) at z=0.840, however with unusually strong narrow emission lines. The estimated black hole mass of ~ 10^7 Msun implies bolometric nuclear luminosity close to the Eddington limit. SDSS J1100+4421 is also extremely radio-loud, with a radio loudness parameter of R ~ 4 x 10^2 - 3 x 10^3, which implies the presence of relativistic jets. Rapid and large-amplitude optical variability of the target, reminiscent of that found in a few radio- and gamma-ray loud NLS1s, is therefore produced most likely in a blazar-like core. The 1.4 GHz radio image of the source shows an extended structure with a linear size of about 100 kpc. If SDSS J1100+4421 is a genuine NLS1, as suggested here, this radio structure would then be the largest ever discovered in this type of active galaxies.

preprint2014arXiv

Radiation Drag Effects in Black Hole Outflows from Super-critical Accretion Disks via Special Relativistic Radiation Magnetohydrodynamics Simulations

By performing 2.5-dimensional special relativistic radiation magnetohydrodynamics simulations, we study the super-critical accretion disks and the outflows launched via the radiation force. We find that the outflow is accelerated by the radiation flux force, but the radiation drag force prevents the outflow velocity from increasing. The outflow velocity saturates around 30-40% of the light speed around the rotation axis, since then the flux force balances with the drag force. Our simulations show that the outflow velocity is kept nearly constant in the regime of \dot{M}_{BH} ~ 100-1000 L_{Edd}/c^2, where \dot{M}_{BH} is the mass accretion rate, L_{Edd} is the Eddington luminosity, and c is the light speed. Such a faster outflow is surrounded by a slower outflow of ~ 0.1c. This velocity is also determined by force balance between the radiation flux force and the radiation drag. The radiation drag works to collimate the slower outflow in cooperation with the Lorentz force, although the faster outflow is mainly collimated by the Lorentz force. The kinetic energy is carried by the slower outflow rather than by the faster outflow. The total kinetic luminosity of the outflow as well as the photon luminosity is ~ L_{Edd}, almost independent of the mass accretion rate.

preprint2014arXiv

Radiation hydrodynamics simulations of wide-angle outflows from super-critical accretion disks around black holes

By performing two-dimensional radiation hydrodynamics simulations with large computational domain of 5000 Schwarzschild radius, we revealed that wide-angle outflow is launched via the radiation force from the super-critical accretion flows around black holes. The angular size of the outflow, of which the radial velocity (v_r) is over the escape velocity (v_esc), increases with an increase of the distance from the black hole. As a result, the mass is blown away with speed of v_r > v_esc in all direction except for the very vicinity of the equatorial plane, theta=0-85^circ, where theta is the polar angle. The mass ejected from the outer boundary per unit time by the outflow is larger than the mass accretion rate onto the black hole, ~150L_Edd/c^2, where L_Edd and c are the Eddington luminosity and the speed of light. Kinetic power of such wide-angle high-velocity outflow is comparable to the photon luminosity and is a few times larger than the Eddington luminosity. This corresponds to ~10^39-10^40 erg/s for the stellar mass black holes. Our model consistent with the observations of shock excited bubbles observed in some ultra-luminous X-ray sources (ULXs), supporting a hypothesis that ULXs are powered by the super-critical accretion onto stellar mass black holes.

preprint2014arXiv

Two dimensional numerical simulations of Supercritical Accretion Flows revisited

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{θϕ}$ component of the viscous stress and consider various values of the viscous parameter $α$. We find that when $T_{θϕ}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $α$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $ρ(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $α$ ranging from $α=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $α$, the flow is marginally stable (when $α$ is small) or unstable (when $α$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

preprint2013arXiv

Recurrent Outbursts and Jet Ejections Expected in Swift J1644+57: Limit-Cycle Activities in a Supermassive Black Hole

The tidal disruption event by a supermassive black hole in Swift J1644+57 can trigger limit-cycle oscillations between a supercritically accreting X-ray bright state and a subcritically accreting X-ray dim state. Time evolution of the debris gas around a black hole with mass $M=10^{6} {\MO}$ is studied by performing axisymmetric, two-dimensional radiation hydrodynamic simulations. We assumed the $α$-prescription of viscosity, in which the viscous stress is proportional to the total pressure. The mass supply rate from the outer boundary is assumed to be ${\dot M}_{\rm supply}=100L_{\rm Edd}/c^2$, where $L_{\rm Edd}$ is the Eddington luminosity, and $c$ is the light speed. Since the mass accretion rate decreases inward by outflows driven by radiation pressure, the state transition from a supercritically accreting slim disk state to a subcritically accreting Shakura-Sunyaev disk starts from the inner disk and propagates outward in a timescale of a day. The sudden drop of the X-ray flux observed in Swift J1644+57 in August 2012 can be explained by this transition. As long as ${\dot M}_{\rm supply}$ exceeds the threshold for the existence of a radiation pressure dominant disk, accumulation of the accreting gas in the subcritically accreting region triggers the transition from a gas pressure dominant Shakura-Sunyaev disk to a slim disk. This transition takes place at $t {\sim}~50/(α/0.1)$ days after the X-ray darkening. We expect that if $α> 0.01$, X-ray emission with luminosity $\gtrsim 10^{44}$ ${\rm erg}{\cdot}{\rm s}^{-1}$ and jet ejection will revive in Swift J1644+57 in 2013--2014.

preprint2012arXiv

Explicit-Implicit Scheme for Relativistic Radiation Hydrodynamics

We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion (FLD) approximation. For an hyperbolic term, of which the time scale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

preprint2012arXiv

Spectroscopy along Multiple, Lensed Sightlines through Outflowing Winds in the Quasar SDSS J1029+2623

We study the origin of absorption features on the blue side of the C IV broad emission line of the large-separation lensed quasar SDSS J1029+2623 at z_em ~ 2.197. The quasar images, produced by a foreground cluster of galaxies, have a maximum separation angle of ~ 22".5. The large angular separation suggests that the sight-lines to the quasar central source can go through different regions of outflowing winds from the accretion disk of the quasar, providing a unique opportunity to study the structure of outflows from the accretion disk, a key ingredient for the evolution of quasars as well as for galaxy formation and evolution. Based on medium- and high-resolution spectroscopy of the two brightest images conducted at the Subaru telescope, we find that each image has different intrinsic levels of absorptions, which can be attributed either to variability of absorption features over the time delay between the lensed images, ~ 774 days, or to the fine structure of quasar outflows probed by the multiple sight-lines toward the quasar. While both these scenarios are consistent with the current data, we argue that they can be distinguished with additional spectroscopic monitoring observations.

preprint2011arXiv

Global Structure of Three Distinct Accretion Flows and Outflows around Black Holes through Two-Dimensional Radiation-Magnetohydrodynamic Simulations

We present the detailed global structure of black hole accretion flows and outflows through newly performed two-dimensional radiation-magnetohydrodynamic simulations. By starting from a torus threaded with weak toroidal magnetic fields and by controlling the central density of the initial torus, rho_0, we can reproduce three distinct modes of accretion flow. In model A with the highest central density, an optically and geometrically thick supercritical accretion disk is created. The radiation force greatly exceeds the gravity above the disk surface, thereby driving a strong outflow (or jet). Because of the mild beaming, the apparent (isotropic) photon luminosity is ~22L_E (where L_E is the Eddington luminosity) in the face-on view. Even higher apparent luminosity is feasible if we increase the flow density. In model B with a moderate density, radiative cooling of the accretion flow is so efficient that a standard-type, cold, and geometrically thin disk is formed at radii greater than ~7R_S (where R_S is the Schwarzschild radius), while the flow is radiatively inefficient otherwise. The magnetic-pressure-driven disk wind appears in this model. In model C the density is too low for the flow to be radiatively efficient. The flow thus becomes radiatively inefficient accretion flow, which is geometrically thick and optically thin. The magnetic-pressure force, in cooperation with the gas-pressure force, drives outflows from the disk surface, and the flow releases its energy via jets rather than via radiation. Observational implications are briefly discussed.

preprint2011arXiv

New Method for Exploring Super-Eddington AGNs by Near-infrared Observations

We propose a new method to explore the candidate super-Eddington active galactic nuclei (AGNs). We examine the properties of infrared (IR) emission from the inner edge of the dusty torus in AGNs, which are powered by super- or sub-Eddington accretion flows around black holes, by considering the dependence of the polar angle on the radiation flux of accretion flows (Watarai et al. 2005). We find that for super-Eddington AGNs, of which the mass accretion rate is more than 10^2 times larger than the Eddington rate, the ratio of the AGN IR luminosity and the disc bolometric luminosity is less than 10^-2, unless the half opening angle of the torus (theta_torus) is small (theta_torus <65 degree). This is due to the self-occultation effect, whereby the self-absorption at the outer region of the super-Eddington flow dilutes the illumination of the torus. Such a small luminosity ratio is not observed in sub-Eddington AGNs, whose mass accretion rate is comparable to or no more than 10 times larger than the Eddington mass accretion rate, except for extremely thin tori (theta_torus >85 degree). We also consider the properties of the near-IR (NIR) emission radiated from hot dust >1000 K. We find that super-Eddington AGNs have a ratio of the NIR luminosity to the bolometric luminosity, L_NIR,AGN/L_bol,disc, at least one order of magnitude smaller than for sub-Eddington AGNs for a wide range of half opening angle (theta_torus > 65 degree), for various types of dusty torus model. Thus, a relatively low L_NIR,AGN/L_bol,disc is a property that allows identification of candidate super-Eddington AGNs. Lastly, we discuss the possibility that NIR-faint quasars at redshift z=6 discovered by a recent deep SDSS survey may be young quasars whose black holes grow via super-Eddington accretion.

preprint2010arXiv

Radiation Magnetohydrodynamics Simulation of Proto-Stellar Collapse: Two-Component Molecular Outflow

We perform a three-dimensional nested-grid radiation magneto-hydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is handled with the flux-limited diffusion approximation, operator-splitting and implicit time-integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and the size of first core is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array (ALMA).

preprint2009arXiv

Global Radiation-Magnetohydrodynamic Simulations of Black Hole Accretion Flow and Outflow: Unified Model of Three States

Black-hole accretion systems are known to possess several distinct modes (or spectral states), such as low/hard state, high/soft state, and so on. Since the dynamics of the corresponding flows is distinct, theoretical models were separately discussed for each state. We here propose a unified model based on our new, global, two-dimensional radiation-magnetohydrodynamic simulations. By controlling a density normalization we could for the first time reproduce three distinct modes of accretion flow and outflow with one numerical code. When the density is large (model A), a geometrically thick, very luminous disk forms, in which photon trapping takes place. When the density is moderate (model B), the accreting gas can effectively cool by emitting radiation, thus generating a thin disk, i.e., the soft-state disk. When the density is too low for radiative cooling to be important (model C), a disk becomes hot, thick, and faint; i.e., the hard-state disk. The magnetic energy is amplified within the disk up to about twice, 30%, and 20% of the gas energy in models A, B, and C, respectively. Notably, the disk outflows with helical magnetic fields, which are driven either by radiation pressure force or magnetic pressure force, are ubiquitous in any accretion modes. Finally, our simulations are consistent with the phenomenological alpha-viscosity prescription, that is, the disk viscosity is proportional to the pressure.

preprint2009arXiv

Modified Slim-Disk Model Based on Radiation-Hydrodynamic Simulation Data: The Conflict Between Outflow and Photon Trapping

Photon trapping and outflow are two key physics associated with the supercritical accretion flow. We investigate the conflict between these two processes based on two-dimensional radiation-hydrodynamic (RHD) simulation data and construct a simplified (radially) one-dimensional model. Mass loss due to outflow, which is not considered in the slim-disk model, will reduce surface density of the flow, and if very significant, it will totally suppress photon trapping effects. If the photon trapping is very significant, conversely, outflow will be suppressed because radiation pressure force will be reduced. To see what actually occurs, we examine the RHD simulation data and evaluate the accretion rate and outflow rate as functions of radius. We find that the former monotonically decreases, while the latter increases, as the radius decreases. However, the former is kept constant at small radii, inside several Schwarzschild radii, since the outflow is suppressed by the photon trapping effects. To understand the conflict between the photon trapping and outflow in a simpler way, we model the radial distribution of the accretion rate from the simulation data and build up a new (radially) one-dimensional model, which is similar to the slim-disk model but incorporates the mass loss effects due to the outflow. We find that the surface density (and, hence, the optical depth) is much reduced even inside the trapping radius, compared with the case without outflow, whereas the effective temperature distribution hardly changes. That is, the emergent spectra do not sensitively depend on the amount of mass outflow. We conclude that the slim-disk approach is valid for interpreting observations, even if the outflow is taken into account.

preprint2009arXiv

New Spectral State of Supercritical Accretion Flow with Comptonizing Outflow

Supercritical accretion flows inevitably produce radiation-pressure driven outflows, which will Compton up-scatter soft photons from the underlying accretion flow, thereby making hard emission. We perform two dimensional radiation hydrodynamic simulations of supercritical accretion flows and outflows, incorporating such Compton scattering effects, and demonstrate that there appears a new hard spectral state at higher photon luminosities than that of the slim-disk state. In this state, as the photon luminosity increases, the photon index decreases and the fraction of the hard emission increases. The Compton $y$-parameter is of the order of unity (and thus the photon index will be $\sim 2$) when the apparent photon luminosity is ${\sim}30L_{\rm E}$ (with $L_{\rm E}$ being the Eddington luminosity) for nearly face-on sources. This explains the observed spectral hardening of the ULX NGC1313 X-2 in its brightening phase and thus supports the model of supercritical accretion onto stellar mass black holes in this ULX.

preprint2007arXiv

Instability of Population III Black Hole Accretion Disks

We investigate the stability of black hole accretion disks in a primordial environment (POP III disks for short), by solving the vertical structure of optically thick disks, including convective energy transport, and by employing a one-zone model for optically thin isothermal disks. Because of the absence of metals in POP III disks, we find significant differences in stability associated with ionization between POP III disks and the disks of solar metallicity. An unstable branch in S-shaped equilibrium curves on the Mdot-Sigma (mass accretion rate - surface density) plane extends to a larger surface density compared with the case of disks of solar metallicity. The resulting equilibrium loci indicate that quasi-periodic oscillations in luminosity can also be driven in POP III disks, and their maximal luminosity is typically by an order of magnitude larger than that of the disks of solar metallicity. Such a strong outburst of POP III disks can be observed by future huge telescopes, in case that the mass is supplied onto the disks at the Bondi accretion rates in typical virialized small dark halos.

preprint2007arXiv

Two-dimensional Radiation-hydrodynamic Model for Supercritical Disk Accretion Flows onto Neutron Stars

We performed two-dimensional radiation hydrodynamic simulations of supercritical accretion flows around neutron stars (NSs). In contrast with the accretion flows onto black holes (BHs), we find that the shell-shaped high-density regions form around the NSs, since the radiation force is enhanced in the innermost regions. The enhanced radiation force drives strong outflows above and below the disk. The mass-accretion rate onto the NS exceeds the critical rate, $L_{\rm E}/c^2$, with $L_{\rm E}$ being the Eddington luminosity. However it is about $20-30%$ of that onto the BH, under the condition that we employ the same mass-input rate, $\dot{M}_{\rm input}$, which is mass injected from the outer disk boundary per unit time. The mass-outflow rate is a few-times larger in flows around NSs than in flows around BHs. The supercritical NS accretion flows mainly release the accretion energy as the kinetic energy of the outflows, though the disk luminosity is predominant over the kinetic energy output rate in the BH accretion flows. The resulting velocity and mass-outflow rate of the outflows are $0.2-0.3c$ and $150-700L_{\rm E}/c^2$, respectively, for the mass-input rate of $3\times 10^2\lsim \dot{M}_{\rm input}/(L_{\rm E}/c^2)\lsim 3\times 10^3$. This implies that the SS433 jets can be roughly explained by the supercritical accretion onto a NS. However, the collimation angle of the outflows in our simulations ($\sim 20^\circ$) is larger than that of the SS433 jets (a few degrees).