Researcher profile

Ken Ohsuga

Ken Ohsuga contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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