Source author record

Tomohisa Kawashima

Tomohisa Kawashima 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

12works
2topics
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

12 published item(s)

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

Mushroom-instability-driven Magnetic Reconnections in Collisionless Relativistic Jets

We study the kinetic plasma dynamics in collisionless relativistic jets with velocity shear, by carrying out particle-in-cell simulations in the transverse plane of a jet. It is discovered that intermittent magnetic reconnections (MRs) are driven by Mushroom instability (MI), which is an important kinetic-scale plasma instability in the plasma shear-flows with relativistic bulk speed. We refer to this sequence of kinetic plasma phenomena as "MI-driven MR". The MI-driven MRs intermittently occur with moving the location of the reconnection points from the vicinity of the initial velocity-shear surface towards the center of the jet. As a consequence, the number density of high energy electrons, which are accelerated by MI-driven MRs, increases with time in the region inside the initial velocity-shear surface with accompanying the generation and subsequent amplification of magnetic fields by MI. The maximum Lorentz factor of electrons increases with initial bulk Lorentz factor of the jet. A possible relation of MI-driven MR to the bright synchrotron emission in jet-spine of active galactic nucleus jets is also discussed.

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.

preprint2021arXiv

The intrinsic structure of Sagittarius A* at 1.3 cm and 7 mm

Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets to resolve the innermost region of SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multi-wavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ~704$\pm$102 $μ$as (axial ratio $\sim$1.19$^{+0.24}_{-0.19}$) and $\sim$300$\pm$25 $μ$as (axial ratio $\sim$1.28$\pm$0.2) at 1.349 cm and 6.95 mm respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power-law, with an index $\sim$1.2$\pm$0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes.

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.

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.

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.

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.

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.