Source author record

Sanshiro Shibata

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

3works
2topics
3close collaborators

Published work

3 publications

preprint2015arXiv

Nuclear Composition of Magnetized GRB Jets

We investigate the fraction of metal nuclei in the relativistic jets of gamma-ray bursts associated with core-collapse supernovae. We simulate the fallback in jet-induced explosions with two-dimensional relativistic hydrodynamics calculations and the jet acceleration with steady, radial, relativistic magnetohydrodynamics calculations, and derive detail nuclear composition of the jet by postprocessing calculation. We found that if the temperature at the jet launch site is above 4.7×1094.7\times 10^9K, quasi-statistical equilibrium (QSE) is established and heavy nuclei are dissociated to light particles such as 4^4He during the acceleration of the jets. The criterion for the survival of metal nuclei is written in terms of the isotropic jet luminosity as Ljiso≲3.9×1050(Ri/107cm)2(1+σi) erg s−1L_{\rm j}^{\rm iso} \lesssim 3.9\times 10^{50}(R_{\rm i}/10^7{\rm cm})^2 (1+σ_{\rm i})~{\rm erg~s^{-1}}, where RiR_{\rm i} and σiσ_{\rm i} are the initial radius of the jets and the initial magnetization parameter, respectively. If the jet is initially dominated by radiation field (i.e., σi≪1σ_{\rm i} \ll 1) and the isotropic luminosity is relatively high (Ljiso≳4×1052 ergs−1L_{\rm j}^{\rm iso}\gtrsim 4\times 10^{52}~{\rm erg s^{-1}}), the metal nuclei cannot survive in the jet. On the other hand, if the jet is mainly accelerated by magnetic field (i.e., σi≫1σ_{\rm i} \gg1), metal nuclei initially contained in the jet can survive without serious dissociation even for the case of high luminosity jet. If the jet contains metal nuclei, the dominant nuclei are 28^{28}Si, 16^{16}O, and 32^{32}S and the mean mass number can be ⟨\langleA⟩∼25\rangle\sim25.

preprint2015arXiv

Time-dependent Multi-group Multidimensional Relativistic Radiative Transfer Code Based On Spherical Harmonic Discrete Ordinate Method

We develop a time-dependent multi-group multidimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) that evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with a ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed frame approach; the source function is evaluated in the comoving frame whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated with various test problems and comparisons with results of a relativistic Monte Carlo code. These validations confirm that the code correctly calculates intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations.

preprint2014arXiv

Random walks and effective optical depth in relativistic flow

We investigate the random walk process in relativistic flow. In the relativistic flow, photon propagation is concentrated in the directions of the flow velocity due to relativistic beaming effect. We show that, in the pure scattering case, the number of scatterings is proportional to the size parameter ξ≡L/l0ξ\equiv L/l_0 if the flow velocity β≡v/cβ\equiv v/c satisfies β/Γ≫ξ−1β/Γ\gg ξ^{-1}, while it is proportional to ξ2ξ^2 if β/Γ≪ξ−1β/Γ\ll ξ^{-1} where LL and l0l_0 are the size of the system in the observer frame and the mean free path in the comoving frame, respectively. We also examine the photon propagation in the scattering and absorptive medium. We find that, if the optical depth for absorption τaτ_{\rm a} is considerably smaller than the optical depth for scattering τsτ_{\rm s} (τa/τs≪1τ_{\rm a}/τ_{\rm s} \ll 1) and the flow velocity satisfies β≫2τa/τsβ\gg \sqrt{2τ_{\rm a}/τ_{\rm s}}, the effective optical depth is approximated by τ∗≃τa(1+β)/βτ_*\simeqτ_{\rm a}(1+β)/β. Furthermore, we perform Monte Carlo simulations of radiative transfer and compare the results with the analytic expression for the number of scattering. The analytic expression is consistent with the results of the numerical simulations. The expression derived in this Letter can be used to estimate the photon production site in relativistic phenomena, e.g., gamma-ray burst and active galactic nuclei.

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.