-
Production of intense episodic Alfvén pulses: GRMHD simulation of black hole accretion disks
Authors:
Akira Mizuta,
Toshikazu Ebisuzaki,
Toshiki Tajima,
Shigehiro Nagataki
Abstract:
The episodic dynamics of the magnetic eruption of a spinning black hole (BH) accretion disks and its associated intense shapeup of their jets is studied via three-dimensional general-relativistic magnetohydrodynamics (GRMHD). The embedded magnetic fields in the disk get amplified by the magnetorotational instability (MRI) so large as to cause an eruption of magnetic field (recconection) and large…
▽ More
The episodic dynamics of the magnetic eruption of a spinning black hole (BH) accretion disks and its associated intense shapeup of their jets is studied via three-dimensional general-relativistic magnetohydrodynamics (GRMHD). The embedded magnetic fields in the disk get amplified by the magnetorotational instability (MRI) so large as to cause an eruption of magnetic field (recconection) and large chunks of matter episodically accrete toward the roots of the jets upon such an event. We also find that the eruption events produce intensive Alfvén pulses, which propagate through the jets. After the eruption, the disk backs to the weakly magnetic states. Such disk activities cause short time variabilities in mass accretion rate at the event horizon as well as electromagnetic luminosity inside the jet. Since the dimensionless strength parameter $a_0=eE/m_e ωc$ of these Alfvén wave pulses is extremely high for a substantial fraction of Eddington accretion rate accretion flow onto a supermassive black hole, the Alfvén shocks turn into ultrarelativistic $(a_0\gg 1)$ bow wake acceleration, manifesting into the ultra-high energy cosmic rays and electrons which finally emit gamma-rays. Since our GRMHD model has universality in its spatial and temporal scales, it is applicable to a wide range of astrophysical objects ranging from those of AGN (which is the primary target of this research), to micro-quasars. Properties such as time variabilities of blazar gamma-ray flares and spectrum observed by {\it Fermi} Gamma-ray Observatory are well explained by linear acceleration of electrons by the bow wake.
△ Less
Submitted 4 June, 2018; v1 submitted 27 July, 2017;
originally announced July 2017.
-
Spectral and Polarization Properties of Photospheric Emission From Stratified Jets
Authors:
Hirotaka Ito,
Shigehiro Nagataki,
Jin Matsumoto,
Shiu-Hang Lee,
Alexey Tolstov,
Jirong Mao,
Maria Dainotti,
Akira Mizuta
Abstract:
We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the je…
▽ More
We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($θ$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.
△ Less
Submitted 24 May, 2014;
originally announced May 2014.
-
Photospheric emission from stratified jets
Authors:
Hirotaka Ito,
Shigehiro Nagataki,
Masaomi Ono,
Shiu-Hang Lee,
Jirong Mao,
Shoichi Yamada,
Asaf Pe'er,
Akira Mizuta,
Seiji Harikae
Abstract:
We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a…
▽ More
We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.
△ Less
Submitted 26 August, 2013; v1 submitted 20 June, 2013;
originally announced June 2013.
-
Opening Angles of Collapsar Jets
Authors:
Akira Mizuta,
Kunihito Ioka
Abstract:
We investigate the jet propagation and breakout from the stellar progenitor for gamma-ray burst (GRB) collapsars by performing two-dimensional relativistic hydrodynamic simulations and analytical modeling. We find that the jet opening angle is given by $θ_j \sim 1/5 Γ_{0}$, and infer the initial Lorentz factor of the jet at the central engine, $Γ_0$, is a few for existing observations of $θ_j$. Th…
▽ More
We investigate the jet propagation and breakout from the stellar progenitor for gamma-ray burst (GRB) collapsars by performing two-dimensional relativistic hydrodynamic simulations and analytical modeling. We find that the jet opening angle is given by $θ_j \sim 1/5 Γ_{0}$, and infer the initial Lorentz factor of the jet at the central engine, $Γ_0$, is a few for existing observations of $θ_j$. The jet keeps the Lorentz factor low inside the star by converging cylindrically via collimation shocks under the cocoon pressure, and accelerates at jet breakout before the free expansion to a hollow-cone structure. In this new picture the GRB duration is determined by the sound crossing time of the cocoon, after which the opening angle widens, reducing the apparent luminosity. Some bursts violating the maximum opening angle $θ_{j,\max}\sim 1/5 \sim 12^{\circ}$ imply the existence of a baryon-rich sheath or a long-acting jet. We can explain the slopes in both Amati and Yonetoku spectral relations using an off-centered photosphere model, if we make only one assumption that the total jet luminosity is proportional to the initial Lorentz factor of the jet. We also numerically calibrate the pre-breakout model (Bromberg et al.) for later use.
△ Less
Submitted 8 November, 2013; v1 submitted 31 March, 2013;
originally announced April 2013.
-
Gamma-Ray Burst without Baryonic and Magnetic Load?
Authors:
Kunihito Ioka,
Yutaka Ohira,
Norita Kawanaka,
Akira Mizuta
Abstract:
We show that, contrary to common belief, internal shocks can arise in an accelerating radiation-dominated jet if it is confined even weakly to a converging opening angle because the acceleration declines. The radiation-dominated internal shock (RDIS) enables a very efficient yet highly nonthermal emission by Fermi-like photon acceleration, keeping the electron-positron ($e^{\pm}$) pair photosphere…
▽ More
We show that, contrary to common belief, internal shocks can arise in an accelerating radiation-dominated jet if it is confined even weakly to a converging opening angle because the acceleration declines. The radiation-dominated internal shock (RDIS) enables a very efficient yet highly nonthermal emission by Fermi-like photon acceleration, keeping the electron-positron ($e^{\pm}$) pair photosphere and inertia up to a high Lorentz factor >1000. In gamma-ray bursts (GRBs), a weak confinement would persist beyond the progenitor star or surrounding matter because of the fast cocoon accompanying the breakout jet. The simplest model predicts few high-energy cosmic rays and neutrinos, and a correlation between the early afterglow and the GeV-TeV prompt emission. The central engine allows a less fine-tuned baryon load than previously thought, even including pure-leptonic unmagnetized outflows.
△ Less
Submitted 14 September, 2011; v1 submitted 29 March, 2011;
originally announced March 2011.
-
Thermal Radiation from GRB Jets
Authors:
Akira Mizuta,
Shigehiro Nagataki,
Junichi Aoi
Abstract:
In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around the head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-a…
▽ More
In this study, the light curves and spectrum of the photospheric thermal radiation from ultrarelativistic gamma-ray burst (GRB) jets are calculated using 2D relativistic hydrodynamic simulations of jets from a collapsar. As the jet advances, the density around the head of the jet decreases, and its Lorentz factor reaches as high as 200 at the photosphere and 400 inside the photosphere. For an on-axis observer, the photosphere appears concave shaped due to the low density and high beaming factor of the jet. The luminosity varies because of the abrupt change in the position of the photosphere due to the internal structure of the jet. Comparing our results with GRB090902B, the flux level of the thermal-like component is similar to our model, although the peak energy looks a little bit higher (but still within a factor of 2). From the comparison, we estimate that the bulk Lorentz factor of GRB090902B is $Γ\sim 2.4 \times 10^2 (r/10^{12} \rm cm$) where $r$ is the radius of the photosphere. The spectrum for an on-axis observer is harder than that for an off-axis observer. There is a time lag of a few seconds for high energy bands in the light curve. This may be the reason for the delayed onset of GeV emission seen in GRB080916C. The spectrum below the peak energy is a power law and the index is $2.3 \sim 2.6$ which is softer than that of single temperature plank distribution but still harder than that of typical value of observed one.
△ Less
Submitted 8 February, 2011; v1 submitted 12 June, 2010;
originally announced June 2010.
-
Hysteresis of Backflow Imprinted in Collimated Jets
Authors:
Akira Mizuta,
Motoki Kino,
Hiroki Nagakura
Abstract:
We report two different types of backflow from jets by performing 2D special relativistic hydrodynamical simulations. One is anti-parallel and quasi-straight to the main jet (quasi-straight backflow), and the other is bent path of the backflow (bent backflow). We find that the former appears when the head advance speed is comparable to or higher than the local sound speed at the hotspot while th…
▽ More
We report two different types of backflow from jets by performing 2D special relativistic hydrodynamical simulations. One is anti-parallel and quasi-straight to the main jet (quasi-straight backflow), and the other is bent path of the backflow (bent backflow). We find that the former appears when the head advance speed is comparable to or higher than the local sound speed at the hotspot while the latter appears when the head advance speed is slower than the sound speed bat the hotspot. Bent backflow collides with the unshocked jet and laterally squeezes the jet. At the same time, a pair of new oblique shocks are formed at the tip of the jet and new bent fast backflows are generated via these oblique shocks. The hysteresis of backflow collisions is thus imprinted in the jet as a node and anti-node structure. This process also promotes broadening of the jet cross sectional area and it also causes a decrease in the head advance velocity. This hydrodynamic process may be tested by observations of compact young jets.
△ Less
Submitted 18 December, 2009;
originally announced December 2009.
-
Angular Energy Distribution of Collapsar-Jets
Authors:
Akira Mizuta,
Miguel A. Aloy
Abstract:
Collapsars are fast-spinning, massive stars, whose core collapse liberates an energy, that can be channeled in the form of ultrarelativistic jets. These jets transport the energy from the collapsed core to large distances, where it is dissipated in the form of long-duration gamma-ray bursts. In this paper we study the dynamics of ultrarelativistic jets produced in collapsars. Also we extrapolate…
▽ More
Collapsars are fast-spinning, massive stars, whose core collapse liberates an energy, that can be channeled in the form of ultrarelativistic jets. These jets transport the energy from the collapsed core to large distances, where it is dissipated in the form of long-duration gamma-ray bursts. In this paper we study the dynamics of ultrarelativistic jets produced in collapsars. Also we extrapolate our results to infer the angular energy distribution of the produced outflows in the afterglow phase. Our main focus is to look for global energetical properties which can be imprinted by the different structure of different progenitor stars. Thus, we employ a number of pre-supernova, stellar models (with distinct masses and metallicities), and inject in all of them jets with fixed initial conditions. We assume that at the injection nozzle, the jet is mildly relativistic (Lorentz factor $\sim 5$), has a finite half-opening angle ($5^\circ$), and carries a power of $10^{51} $erg s$^{-1}$. These jets arrive intact to the stellar surface and break out of it. A large Lorentz factor region $Γ\simmore 100$ develops well before the jet reaches the surface of the star, in the unshocked part of the beam, located between the injection nozzle and the first recollimation shock. These high values of $Γ$ are possible because the finite opening angle of the jet allows for free expansion towards the radial direction. We find a strong correlation between the angular energy distribution of the jet, after its eruption from the progenitor surface, and the mass of the progenitors. The angular energy distribution of the jets from light progenitor models is steeper than that of the jets injected in more massive progenitor stars. This trend is also imprinted in the angular distribution of isotropic equivalent energy.
△ Less
Submitted 11 May, 2009; v1 submitted 28 December, 2008;
originally announced December 2008.
-
Numerical Study on GRB-Jet Formation in Collapsars
Authors:
Shigehiro Nagataki,
Rohta Takahashi,
Akira Mizuta,
Tomoya Takiwaki
Abstract:
Two-dimensional magnetohydrodynamic simulations are performed using the ZEUS-2D code to investigate the dynamics of a collapsar that generates a GRB jet, taking account of realistic equation of state, neutrino cooling and heating processes, magnetic fields, and gravitational force from the central black hole and self gravity. It is found that neutrino heating processes are not so efficient to la…
▽ More
Two-dimensional magnetohydrodynamic simulations are performed using the ZEUS-2D code to investigate the dynamics of a collapsar that generates a GRB jet, taking account of realistic equation of state, neutrino cooling and heating processes, magnetic fields, and gravitational force from the central black hole and self gravity. It is found that neutrino heating processes are not so efficient to launch a jet in this study. It is also found that a jet is launched mainly by B_φfields that are amplified by the winding-up effect. However, since the ratio of total energy relative to the rest mass energy in the jet is not so high as several hundred, we conclude that the jets seen in this study are not be a GRB jet. This result suggests that general relativistic effects, which are not included in this study, will be important to generate a GRB jet. Also, the accretion disk with magnetic fields may still play an important role to launch a GRB jet, although a simulation for much longer physical time (\sim 10-100 s) is required to confirm this effect. It is shown that considerable amount of 56Ni is synthesized in the accretion disk. Thus there will be a possibility for the accretion disk to supply sufficient amount of 56Ni required to explain the luminosity of a hypernova. Also, it is shown that neutron-rich matter due to electron captures with high entropy per baryon is ejected along the polar axis. Moreover, it is found that the electron fraction becomes larger than 0.5 around the polar axis near the black hole by ν_e capture at the region. Thus there will be a possibility that r-process and r/p-process nucleosynthesis occur at these regions. Finally, much neutrons will be ejected from the jet, which suggests that signals from the neutron decays may be observed as the delayed bump of afterglow or gamma-rays.
△ Less
Submitted 23 December, 2006; v1 submitted 11 August, 2006;
originally announced August 2006.
-
Collimated Jet or Expanding Outflow: Possible Origins of GRBs and X-Ray Flashes
Authors:
Akira Mizuta,
Tatsuya Yamasaki,
Shigehiro Nagataki,
Shin Mineshige
Abstract:
We investigate the dynamics of an injected outflow propagating in a progenitor in the context of the collapsar model for gamma-ray bursts (GRBs) through two dimensional axisymmetric relativistic hydrodynamic simulations. Initially, we locally inject an outflow near the center of a progenitor. We calculate 25 models, in total, by fixing its total input energy to be 10^{51} ergs s^{-1} and radius…
▽ More
We investigate the dynamics of an injected outflow propagating in a progenitor in the context of the collapsar model for gamma-ray bursts (GRBs) through two dimensional axisymmetric relativistic hydrodynamic simulations. Initially, we locally inject an outflow near the center of a progenitor. We calculate 25 models, in total, by fixing its total input energy to be 10^{51} ergs s^{-1} and radius of the injected outflow to be $7\times 10^7$ cm while varying its bulk Lorentz factor, $Γ_{0} = 1.05\sim 5$, and its specific internal energy, $ε_0/c^2 = 0.1\sim 30$. The injected outflow propagates in the progenitor and drives a large-scale outflow or jet. We find a smooth but dramatic transition from a collimated jet to an expanding outflow among calculated models. The maximum Lorentz factor is, on the other hand, sensitive to both of $Γ_0$ and $ε_0$; roughly $Γ_{\rm max} \sim Γ_0 (1+ε_0/c^2)$. Our finding will explain a smooth transition between the GRBs, X-ray rich GRBs (XRRs) and X-ray Flashes (XRFs) by the same model but with different $ε_0$ values.
△ Less
Submitted 24 July, 2006;
originally announced July 2006.
-
Formation of Pillars at the Boundaries between H II Regions and Molecular Clouds
Authors:
Akira Mizuta,
Jave O. Kane,
Marc W. Pound,
Bruce A. Remington,
Dmitri D. Ryutov,
Hideaki Takabe
Abstract:
We investigate numerically the hydrodynamic instability of an ionization front (IF) accelerating into a molecular cloud, with imposed initial perturbations of different amplitudes. When the initial amplitude is small, the imposed perturbation is completely stabilized and does not grow. When the initial perturbation amplitude is large enough, roughly the ratio of the initial amplitude to waveleng…
▽ More
We investigate numerically the hydrodynamic instability of an ionization front (IF) accelerating into a molecular cloud, with imposed initial perturbations of different amplitudes. When the initial amplitude is small, the imposed perturbation is completely stabilized and does not grow. When the initial perturbation amplitude is large enough, roughly the ratio of the initial amplitude to wavelength is greater than 0.02, portions of the IF temporarily separate from the molecular cloud surface, locally decreasing the ablation pressure. This causes the appearance of a large, warm HI region and triggers nonlinear dynamics of the IF. The local difference of the ablation pressure and acceleration enhances the appearance and growth of a multimode perturbation. The stabilization usually seen at the IF in the linear regimes does not work due to the mismatch of the modes of the perturbations at the cloud surface and in density in HII region above the cloud surface. Molecular pillars are observed in the late stages of the large amplitude perturbation case. The velocity gradient in the pillars is in reasonably good agreement with that observed in the Eagle Nebula. The initial perturbation is imposed in three different ways: in density, in incident photon number flux, and in the surface shape. All cases show both stabilization for a small initial perturbation and large growth of the second harmonic by increasing amplitude of the initial perturbation above a critical value.
△ Less
Submitted 26 April, 2006;
originally announced April 2006.
-
Explosive Nucleosynthesis in GRB Jets Accompanied by Hypernovae
Authors:
Shigehiro Nagataki,
Akira Mizuta,
Katsuhiko Sato
Abstract:
Two-dimensional hydrodynamic simulations are performed to investigate explosive nucleosynthesis in a collapsar using the model of MacFadyen and Woosley (1999). It is shown that 56Ni is not produced in the jet of the collapsar sufficiently to explain the observed amount of a hypernova when the duration of the explosion is \sim 10 sec, which is considered to be the typical timescale of explosion i…
▽ More
Two-dimensional hydrodynamic simulations are performed to investigate explosive nucleosynthesis in a collapsar using the model of MacFadyen and Woosley (1999). It is shown that 56Ni is not produced in the jet of the collapsar sufficiently to explain the observed amount of a hypernova when the duration of the explosion is \sim 10 sec, which is considered to be the typical timescale of explosion in the collapsar model. Even though a considerable amount of 56Ni is synthesized if all explosion energy is deposited initially, the opening angles of the jets become too wide to realize highly relativistic outflows and gamma-ray bursts in such a case. From these results, it is concluded that the origin of 56Ni in hypernovae associated with GRBs is not the explosive nucleosynthesis in the jet. We consider that the idea that the origin is the explosive nucleosynthesis in the accretion disk is more promising. We also show that the explosion becomes bi-polar naturally due to the effect of the deformed progenitor. This fact suggests that the 56Ni synthesized in the accretion disk and conveyed as outflows are blown along to the rotation axis, which will explain the line features of SN 1998bw and double peaked line features of SN 2003jd. Some fraction of the gamma-ray lines from 56Ni decays in the jet will appear without losing their energies because the jet becomes optically thin before a considerable amount of 56Ni decays as long as the jet is a relativistic flow. We show that abundance of nuclei whose mass number \sim 40 in the ejecta depends sensitively on the energy deposition rate. So it may be determined by observations of chemical composition in metal poor stars which model is the proper one as a model of a gamma-ray burst accompanied by a hypernova.
△ Less
Submitted 11 May, 2006; v1 submitted 5 January, 2006;
originally announced January 2006.
-
Hydrodynamic instability of ionization fronts in HII regions
Authors:
Akira Mizuta,
Jave O. Kane,
Marc W. Pound,
Bruce A. Remington,
Dmitri D. Ryutov,
Hideaki Takabe
Abstract:
We investigate hydrodynamic instability of accelerating ionization fronts using two dimensional hydrodynamic simulations that include detailed energy deposition and release due to the absorption of UV radiation, recombination of hydrogen, radiative molecular cooling, and magnetic pressure. We consider linear perturbation growth and find that the stabilization mechanism associated with non-accele…
▽ More
We investigate hydrodynamic instability of accelerating ionization fronts using two dimensional hydrodynamic simulations that include detailed energy deposition and release due to the absorption of UV radiation, recombination of hydrogen, radiative molecular cooling, and magnetic pressure. We consider linear perturbation growth and find that the stabilization mechanism associated with non-accelerated fronts remains a significant factor even when acceleration is present. Conversely, if recombination in the ionized region is turned off,Rayleigh-Taylor (RT) instability becomes effective, and the classical RT growth rate recovered.
△ Less
Submitted 20 November, 2004;
originally announced November 2004.
-
Hydrodynamics of Internal Shocks in Relativistic Outflows
Authors:
M. Kino,
A. Mizuta,
A. Celotti,
S. Yamada
Abstract:
We study the hydrodynamical effects of two colliding shells, adopted to model internal shocks in various relativistic outflows such as gamma-ray bursts and blazars. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of the two shells. The shell spreading with a few ten percent of the speed of…
▽ More
We study the hydrodynamical effects of two colliding shells, adopted to model internal shocks in various relativistic outflows such as gamma-ray bursts and blazars. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of the two shells. The shell spreading with a few ten percent of the speed of light is also shown to be a notable aspect. The conversion efficiency of the bulk kinetic energy to internal one shows deviations from the widely-used inelastic two-point-mass-collision model. Observational implications are also shortly discussed.
△ Less
Submitted 15 October, 2004;
originally announced October 2004.
-
Hydrodynamical effects in internal shock of relativistic outflows
Authors:
Motoki Kino,
Akira Mizuta,
Shoichi Yamada
Abstract:
We study both analytically and numerically hydrodynamical effects of two colliding shells, the simplified models of the internal shock in various relativistic outflows such as gamma-ray bursts and blazars. We pay particular attention to three interesting cases: a pair of shells with the same rest mass density (``{\it equal rest mass density}''), a pair of shells with the same rest mass (``{\it e…
▽ More
We study both analytically and numerically hydrodynamical effects of two colliding shells, the simplified models of the internal shock in various relativistic outflows such as gamma-ray bursts and blazars. We pay particular attention to three interesting cases: a pair of shells with the same rest mass density (``{\it equal rest mass density}''), a pair of shells with the same rest mass (``{\it equal mass}''), and a pair of shells with the same bulk kinetic energy (``{\it equal energy}'') measured in the intersteller medium (ISM) frame. We find that the density profiles are significantly affected by the propagation of rarefaction waves. A split-feature appears at the contact discontinuity of two shells for the ``equal mass'' case, while no significant split appears for the ``equal energy'' and ``equal rest mass density'' cases. The shell spreading with a few ten percent of the speed of light is also shown as a notable aspect caused by rarefaction waves. The conversion efficiency of the bulk kinetic energy to internal one is numerically evaluated. The time evolutions of the efficiency show deviations from the widely-used inellastic two-point-mass-collision model.
△ Less
Submitted 28 April, 2004;
originally announced April 2004.
-
Propagation and Dynamics of Relativistic Jets
Authors:
Akira Mizuta,
Shoichi Yamada,
Hideaki Takabe
Abstract:
We investigate the dynamics and morphology of jets propagating into the interstellar medium using 2D relativistic hydrodynamic simulations. The calculations are performed assuming axisymmetric geometry and trace a long propagation of jets. The jets are assumed to be 'light' with the density ratio between the beam to the ambient gas much less than unity. We examine the mechanism for the appearanc…
▽ More
We investigate the dynamics and morphology of jets propagating into the interstellar medium using 2D relativistic hydrodynamic simulations. The calculations are performed assuming axisymmetric geometry and trace a long propagation of jets. The jets are assumed to be 'light' with the density ratio between the beam to the ambient gas much less than unity. We examine the mechanism for the appearance of vortices at the head of jets in the hot spot. Such vortices are known as a trigger of a deceleration phase which appears after a short phase in which one dimensional analysis is dominant. We find that an oblique shock at the boundary rim near the end of the beam strongly affects the flow structure in and around the hot spot. Weakly shocked gas passed through this oblique shock and becomes a trigger for the generation of vortices. We also find the parameter dependence of these effects for the propagation and dynamics of the jets. The jet with slower propagation velocity is weakly pinched, has large vortices and shows very complex structure at the head of the jets and extended synchrotron emissivity.
△ Less
Submitted 15 February, 2004;
originally announced February 2004.
-
Explosive Nucleosynthesis Associated with Formation of Jet-induced GRBs in Massive Stars
Authors:
S. Nagataki,
A. Mizuta,
S. Yamada,
H. Takabe,
K. Sato
Abstract:
We perform 2-dimensional relativistic hydrodynamical simulations in the context of collapsar model. Calculations of explosive nucleosynthesis are also accomplished. We investigate the influence of the structure of the progenitor and energy deposition rate on the resulting explosive nucleosynthesis, assuming that 56Ni is mainly synthesized in the jet launched by the neutrino heating. We show the…
▽ More
We perform 2-dimensional relativistic hydrodynamical simulations in the context of collapsar model. Calculations of explosive nucleosynthesis are also accomplished. We investigate the influence of the structure of the progenitor and energy deposition rate on the resulting explosive nucleosynthesis, assuming that 56Ni is mainly synthesized in the jet launched by the neutrino heating. We show the amount of 56Ni is very sensitive to the energy deposition rate. Thus we conclude that it is quite natural not to detect an underlying supernova in some X-ray afterglows as in GRB 010921. We also point out the possibility that the relative abundance of the elements with intermediate mass number such as Si and S in the X-ray afterglow of GRB 011211 may be naturally explained if the energy deposition rate at the central engine is relatively long because little amount of 56Ni should be synthesized under such an environment. If this discussion is true, there should be correlation between the line features in the X-ray afterglow and duration of the GRB. It should be noted that the duration of GRB 011211 is 270 seconds, making it the longest burst ever observed by Beppo-SAX although it suffers from the effect of red-shift (z_host=2.14), which supports our conclusion. Our results also suggest that the type I collapsar model in which the energy deposition rate is relatively low (\dot{E} \sim 10^{51} erg/s) might have difficulty in reproducing the observed amount of 56Ni in a hypernova such as SN 1998bw. This means that the mechanism of the central engine of a hypernova accompanying GRB may be constrained by the discussion of explosive nucleosynthesis.
△ Less
Submitted 20 June, 2003;
originally announced June 2003.
-
Numerical Analysis of Jets Produced by Intense Laser
Authors:
Akira Mizuta,
Shoich Yamada,
Hideaki Takabe
Abstract:
In this paper we present a numerical study of plasma jets produced by intense laser matter interactions. Through this study we hope to better understand astrophysical jets and their recent experimental simulations in the laboratory. We paid special attention to radiation cooling and the interaction of the jet with ambient gas. Four cases are presented in this paper; two of them deal with the pro…
▽ More
In this paper we present a numerical study of plasma jets produced by intense laser matter interactions. Through this study we hope to better understand astrophysical jets and their recent experimental simulations in the laboratory. We paid special attention to radiation cooling and the interaction of the jet with ambient gas. Four cases are presented in this paper; two of them deal with the propagation of jets in vacuum, while in the other two the propagation takes place in the ambient gas. Available experimental results are reproduced to good accuracy in the vacuum case. For jets in ambient gas, we find that the existence of the surrounding gas confines the jet into a narrow cylindrical shape so that both the density and temperature of the jet remain high enough for effective radiation cooling. As a result, a collimated plasma jet is formed in these cases. The dimensionless parameters characterizing the laboratory jets and protostellar jets have overlapping domains. We also discuss the cooling lengths for our model and compare them with the corresponding values in the astrophysical jets. A plasma jet in the ambient gas experiment is proposed which is within the reach of present day technology, and can be relevant to astrophysical phenomena.
△ Less
Submitted 5 November, 2001; v1 submitted 18 October, 2001;
originally announced October 2001.
