Researcher profile

Hiroyuki R. Takahashi

Hiroyuki R. Takahashi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
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

13 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

Directed flow in relativistic resistive magneto-hydrodynamic expansion for symmetric and asymmetric collision systems

We construct a dynamical model for high-energy heavy-ion collision based on the relativistic resistive magneto-hydrodynamic framework. Using our newly developed (3+1)-dimensional relativistic resistive magneto-hydrodynamics code, we investigate magneto-hydrodynamic expansion in symmetric and asymmetric collision systems as a first application to high-energy heavy-ion collisions. As a realistic initial condition for electromagnetic fields, we consider the solutions of the Maxwell equations with the source term of point charged particles moving in the direction of the beam axis, including finite constant electrical conductivity of the medium. We evaluate the directed flow in the symmetric and asymmetric collisions at RHIC energy. We find a significant effect of finite electrical conductivity on the directed flow in the asymmetric collision system. We confirm that a certain amount of energy transfer by dissipation associated with Ohmic conduction occurs in the asymmetric collision system because of asymmetry of the electric field produced by two different colliding nuclei. Because this energy transfer makes the pressure gradient of the medium flatter, the growth of directed flow decreases.

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

Magnetic support for neutrino-driven explosion of 3D non-rotating core-collapse supernova models

The impact of the magnetic field on postbounce supernova dynamics of non-rotating stellar cores is studied by performing three-dimensional magnetohydrodynamics simulations with spectral neutrino transport. The explodability of strongly and weakly magnetized models of $20$ and $27$ $M_{\odot}$ pre-supernova progenitors are compared. We find that although the efficiency for the conversion of the neutrino heating into turbulent energy including magnetic fields in the gain region is not significantly different between the strong and weak field models, the amplified magnetic field due to the neutrino-driven convection on large hot bubbles just behind stalled shock results in a faster and more energetic explosion in the strongly magnetized models. In addition, by comparing the difference between the 2nd- and 5th-order spatial accuracy of the simulation in the strong field model for $27$ $M_{\odot}$ progenitor, we also find that the higher order accuracy in space is beneficial to the explosion because it enhances the growth of neutrino-driven convection in the gain region. Based on our results of core-collapse supernova simulations for the non-rotating model, a new possibility for the origin of the magnetic field of the protoneutron star (PNS) is proposed. The magnetic field is accumulated and amplified to magnetar level, that is, $\mathcal{O}(10^{14})$ G, in the convectively stable shell near the PNS surface.

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

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

Linear Growth of the Kelvin-Helmholtz Instability with an Adiabatic Cosmic-ray Gas

We investigate effects of cosmic-rays on the linear growth of the Kelvin-Helmholtz instability. Cosmic-rays are treated as an adiabatic gas and allowed to diffuse along magnetic field lines. We calculated the dispersion relation of the instability for various sets of two free parameters, the ratio of the cosmic-ray pressure to the thermal gas pressure and the diffusion coefficient. Including cosmic-ray effects, a shear layer is more destabilized and the growth rates can be enhanced in comparison with the ideal magnetohydrodynamical case. Whether the growth rate is effectively enhanced or not depends on the diffusion coefficient of cosmic-rays. We obtain the criterion for effective enhancement by comparing the growing time scale of the instability with the diffusion time scale of cosmic-rays. These results can be applied to various astrophysical phenomena where a velocity shear is present, such as outflows from star-forming galaxies, AGN jet, channel flows resulting from the nonlinear development of the magnetorotational instability, and galactic disks.

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.

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.

preprint2011arXiv

Relativistic Expansion of Magnetic Loops at the Self-similar Stage II: Magnetized outflows interacting with the ambient plasma

We obtained self-similar solutions of relativistically expanding magnetic loops by assuming axisymmetry and a purely radial flow. The stellar rotation and the magnetic fields in the ambient plasma are neglected. We include the Newtonian gravity of the central star. These solutions are extended from those in our previous work (Takahashi, Asano, & Matsumoto 2009) by taking into account discontinuities such as the contact discontinuity and the shock. The global plasma flow consists of three regions, the outflowing region, the post shocked region, and the ambient plasma. They are divided by two discontinuities. The solutions are characterized by the radial velocity, which plays a role of the self-similar parameter in our solutions. The shock Lorentz factor gradually increases with radius. It can be approximately represented by the power of radius with the power law index of 0.25. We also carried out magnetohydrodynamic simulations of the evolution of magnetic loops to study the stability and the generality of our analytical solutions. We used the analytical solutions as the initial condition and the inner boundary conditions. We confirmed that our solutions are stable over the simulation time and that numerical results nicely recover the analytical solutions. We then carried out numerical simulations to study the generality of our solutions by changing the power law index δof the ambient plasma density ρ_0 \propto r^{-δ}. We alter the power law index δfrom 3.5 in the analytical solutions. The analytical solutions are used as the initial conditions inside the shock in all simulations. We observed that the shock Lorentz factor increases with time when δis larger than 3, while it decreases with time when δis smaller than 3. The shock Lorentz factor is proportional to t^{(δ-3)/2}. These results are consistent with the analytical studies by Shapiro (1979).

preprint2011arXiv

Scaling Law of Relativistic Sweet-Parker Type Magnetic Reconnection

Relativistic Sweet-Parker type magnetic reconnection is investigated by relativistic resistive magnetohydrodynamic (RRMHD) simulations. As an initial setting, we assume anti-parallel magnetic fields and a spatially uniform resistivity. A perturbation imposed on the magnetic fields triggers magnetic reconnection around a current sheet, and the plasma inflows into the reconnection region. The inflows are then heated due to ohmic dissipation in the diffusion region, and finally become relativistically hot outflows. The outflows are not accelerated to ultra-relativistic speeds (i.e., Lorentz factor ~ 1), even when the magnetic energy dominates the thermal and rest mass energies in the inflow region. Most of the magnetic energy in the inflow region is converted into the thermal energy of the outflow during the reconnection process. The energy conversion from magnetic to thermal energy in the diffusion region results in an increase in the plasma inertia. This prevents the outflows from being accelerated to ultra-relativistic speeds. We find that the reconnection rate R obeys the scaling relation R S^{-0.5}, where S is the Lundquist number. This feature is the same as that of non-relativistic reconnection. Our results are consistent with the theoretical predictions of Lyubarsky (2005) for Sweet-Parker type magnetic reconnection.

preprint2010arXiv

Magnetic Field Decay Due to the Wave-Particle Resonances in the Outer Crust of the Neutron Star

Bearing in mind the application to the outer crust of the neutron stars (NSs), we investigate the magnetic field decay by means of the fully relativistic Particle-In-Cell simulations. Numerical computations are carried out in 2-dimensions, in which the initial magnetic fields are set to be composed both of the uniform magnetic fields that model the global fields penetrating the NS and of the turbulent magnetic fields that would be originated from the Hall cascade of the large-scale turbulence. Our results show that the whistler cascade of the turbulence transports the magnetic energy preferentially in the direction perpendicular to the uniform magnetic fields. It is also found that the distribution function of electrons becomes anisotropic because electrons with lower energies are predominantly heated in the direction parallel to the uniform magnetic fields due to the Landau resonance, while electrons with higher energies are heated mainly by the cyclotron resonance that makes the distribution function isotropic for the high energy tails. Furthermore we point out that the degree of anisotropy takes maximum as a function of the initial turbulent magnetic energy. As an alternative to the conventional ohmic dissipation, we propose that the magnetic fields in the outer crust of NSs, cascading down to the electron inertial scale via the whistler turbulence, would decay predominantly by the dissipation processes through the Landau damping and the cyclotron resonance.

preprint2010arXiv

Stability of MRI turbulent Accretion Disks

Based on the characteristics of the magnetorotational instability (MRI) and the MRI-driven turbulence, we construct a steady model for a geometrically thin disk using "non-standard" $α$-prescription. The efficiency of the angular momentum transport depends on the magnetic Prandtl number, $Pm = ν/η$, where $ν$ and $η$ are the microscopic viscous and magnetic diffusivities. In our disk model, Shakura-Sunyaev's $α$-parameter has a power-law dependence on the magnetic Prandtl number, that is $α\propto Pm^δ$ where $δ$ is the constant power-law index. Adopting Spitzer's microscopic diffusivities, the magnetic Prandtl number becomes a decreasing function of the disk radius when $δ> 0$. The transport efficiency of the angular momentum and the viscous heating rate are thus smaller in the outer part of the disk, while these are impacted by the size of index $δ$. We find that the disk becomes more unstable to the gravitational instability for a larger value of index $δ$. The most remarkable feature of our disk model is that the thermal and secular instabilities can grow in its middle part even if the radiation pressure is negligibly small in the condition $δ> 2/3$. In the realistic disk systems, it would be difficult to maintain the steady mass accretion state unless the $Pm$-dependence of MRI-driven turbulence is relatively weak.