General Relativistic MHD Simulations of Jet Production

1. General Relativistic MHD Simulations of Jet Production Y.Mizuno (NSSTC/MSFC/NRC) Collaborators K.-I. Nishikawa (NSSTC/UA), P. Hardee(UA), S. Koide (Toyama Univ.), K. Ghosh (NSSTC/USRA) G.J. Fishman (NASA/MSFC)

2. Introduction • Relativistic jets have been observed in AGNs, microquasars and GRBs, and it is believed that they originate in the region near the black hole. • Magnetic accretion model is one of the promising model for the jet formations. • Recent long-term MHD simulations for the evolution of accretion disk show the turbulent and complex structures in accretion disk by the non-axisymmetric growth of magneto-rotational instability

3. Global 3D MHD simulations of black hole accretion disk(Machida & Matsumoto 2003) density Magnetic field line • They have performed 3D non-relativistic MHD simulations of geometrically thick accretion disk with toroidal magnetic fields. • After the non-axisymmetric MRI grows, the disk region becomes turbulent. • Magnetic field lines are less turbulent and globally show a bisymmetric spiral shape. • One-armed spiral mode dominates in the density distribution. • The formation of non-axisymmetric structure is a general picture for the accretion disk.

4. 3D MHD simulations of jets from accretion disk (Kigure & Shibata 2005) Sinusoidal random • They have performed non-relativistic MHD simulations of a jet formation by the interaction between an accretion disk and a large scale magnetic field including non-axisymmetric modes (sinusoidal (m=4) and random perturbation into rotation velocity of disk) • They found the formed jet have non-axisymmetric structure with m=2in both of cases • The non-axisymmetric structure in the jet originates in the accretion disk, not in the jet itself.

5. Purpose of Present Study • We investigate the connection between the non-axisymmetric structure in the jet and that in the disk based on GRMHD simulations • We compare the characteristics of jets found in 3D non-relativistic MHD simulations.

6. 4D General Relativistic MHD Equation • General relativistic equation of conservation laws and Maxwell equations: 　　　　　　　∇n ( n U n) = 0(conservation law of particle-number) 　　　　　　　∇nT mn= 0 (conservation law of energy momentum) 　　　　　　　∂mFnl + ∂nFlm + ∂lF mn = 0 　　　　　　　∇mF mn= - J n • Frozen-in condition:FnmUn= 0 • metric：ds2 = gmn dxm dxn ; g00 = - h02 ; gii= - hi2 ; g0i= - hi2wi (i=1,2,3) ; gij= 0 (i≠j) (Maxwell equations) We neglect the evolution of metric and the essential micro physics (we use gamma-law EOS) n: proper particle number density. p : proper pressure. c: speed of light. e : proper total energy density, e=mnc2 + p / (G -1).; Γ=5/3 m : rest mass of particles. G: specific heat ratio. Umu : velocity four vector. Amu : potential four vector. Jmu : current density four vector. ∇mn : covariant derivative. gmn : metric. Tmn : energy momentum tensor, Tmn = pgmn + (e+p)Um Un+FmsFns -gmnFlkFlk/4. Fmn : field-strength tensor, Fmn =∂mAn -∂n Am.

7. Vector Form of General Relativistic MHD Equation(3+1 Formalism) Special relativistic mass density, gr (conservation law of particle-number) general relativistic effect Special relativistic total momentum density (equation of motion) special relativistic effect Special relativistic total energy density (energy equation) (Maxwell equations) D: density P: momentum density T: energy-momentum tensor e: energy density (ideal MHD condition) Where : (shift velocity) : (Lapse function) : (shift vector)

8. Simulation Model of Jet formation • Numerical code • 3D general relativistic MHD code (Koide et al. 2000; Koide 2003; Mizuno et al. 2005) • Initial condition • Geometrically thin accretion disk (rd/rc=100) rotates around the rotating black hole (a=0.95) • The back ground corona is free-falling to a black hole (Bondi solution) • The global vertical magnetic field (Wald solution; B0=0.1(r0c2)1/2) • Input the sinusoidal perturbation (m=5 mode, 15% of Keplerian velocity) into the rotation of the disk (dvf=0.15VKsin5f) • Numerical Region and Mesh points • 1.0 rS< r < 60 rS, 0< q < p/2, 0< f < 2pwith 150*64*150 mesh points Initial condition of simulations Color: denisty White lines: magnetic field lines

9. Time evolution of density(x-z plane (y=0)); Non-perturb case • The matter in the disk loses its angular momentum by magnetic field and falls to the black hole • A centrifugal barrier decelerates the falling matter and make a shock around r=2rS. • The matter near the shock region is accelerated by the J×B force and gas pressure and forms relativistic outflows. Color: density White lines: poloidal magnetic filed Vector: poloidal velocity

10. Time evolution of density(x-z plane (y=0)); Sinusoidal-perturb case • The general properties are almost same as that of non-perturb case (shock and jet formation). • Due to the initial perturbation of the disk, the disk structure is slightly different with that of non-perturbation case (thicker and less concentration of density near the black hole) Color: density White lines: poloidal magnetic filed Vector: poloidal velocity

11. 2D (x-z plane) snapshot of density and Plasma beta • In both cases, the relativistic jet is formed near the black hole and propagate outward with the twisted magnetic field (torsional Alfven waves) • The accretion disk becomes geometrically thick near the black hole by the propagation of shock produced by the centrifugal barrier around r=2rS. • Due to the initial perturbation of the disk, the disk structure is slightly different with that of non-perturbation case (thicker and less concentration of density near the black hole) Contour Bf

12. Time evolution of density and magnetic field Sinusoidal perturbation case Time: 0 ~ 120 tS (step 5 tS) Color contour: density, slice (at equatorial plane and x=0). White curves: magnetic field lines

13. 3D snapshot of density and magnetic field structure Sinusoidal perturbation case (t/tS=114) Non-perturbation case (t/tS=120) • The overall structure of magnetic field is almost same in both of cases (magnetic field near the black hole is strongly twisted and propagate outward).

14. Time evolution of density and magnetic field into the disk (x-y plane) Sinusoidal perturbation case Time: 0 ~ 120 tS (step 5 tS) Color contor: density at z=1.0 rS. White curves: magnetic field lines • Initial perturbation (m=5) is clearly seen in early stage • Initial perturbation becomes weak by the coupling of accretion matter • High density region (shocked region) show the non-axisymmetric structure clearly (initial m=5 mode and other?). • In late stage, by the coupling them density becomes more complicated structure

15. 2D (x-y plane) plot of density into accretion disk Non-perturbation case (t/tS=120) Sinusoidal perturbation case (t/tS=114) • In the sinusoidal perturbation case, the density into the disk show the complex non-axisymmetric structure. • Some part of the non-axisymmetric structure reflects the initial sinusoidal mode (m=5) and some part of that mixes other (lower?) mode. • The initial perturbation of rotation velocity make homogeneity of accretion to the black hole. • Although the homogeneity of accretion becomes weak by the coupling of accretion matter, it still survives and produce other mode by the coupling.

16. 2D (x-y plane) plot of density into jet • On the other hand, in the jet some non-axisymmetric structure are seen in both of cases, but they are very faint. • The non-axisymmetric structure in the sinusoidal perturbation case reflects the initial sinusoidal perturbation of disk and may mix with other lower mode (is not seen clearly).

17. Outflow structure Total velocity density • The outflow has the complicated density and velocity structure. • Especially, the velocity is fast near the rotation axis and gradually decrease with distance from rotation axis and tend to inflow after the transition region. • These inflow/outflow structure become a strong theoretical support to understand the relativistic outflow and/or inflow properties seen in the absorption line spectra of some narrow emission line quasars.

18. Out flow structure 3D GRMHD, a=0.95, b0=0.1(r0c2)1/2, non-perturb 1D slice along jet (x=5rS, y=0rS) Vtot > Vkand V tot ~ VA In the part of jet, density is lower than coronal density

19. Relativistic radiation transfer Collaborative works with S.Fuerst and K. Wu (astro-ph/0509601) • We have calculated the optically thin thermal free-free emission from a disk-outflow system based on the results of our 2D GRMHD simulations. • We consider a covariant radiation transfer formulation (Fuerst & Wu 2004, A&A, 424, 733) and solve the transfer equation using a ray-tracing algorithm. • In this algorithm, we treat general relativistic effect (light bending, gravitational lensing, gravitational redshift, frame-dragging effect etc.). Emission, absorption & scattering

20. Relativistic radiation transfer Project image of thermal emission (<20 rs) 2D GRMHD simulations (a = 0.95, B = 0.1 (ρc2)-2) • The radiation image shows the front side of the accretion disk and the other side of the disk at the top and bottom regions. • It is because the general relativistic effects. • We can see the propagation of waves and the strong radiation from geometrically thick disk near the BHs. • The jet generated in GRMHD simulation is not visible in the radiation image. • This is because we assume the thermal free-free emission. It has a strong density dependence and the jet is less dense than the disk. • If we calculate the emission with weaker dependence on the density, such as non-thermal process or Compton scattering, the jet would be visible.

21. Summary and Discussion • We have performed 3D GRMHD simulations of the jet formation from the accretion disk with/without the initial perturbation (sinusoidal perturbation) around the rotating black hole. • Although the initial perturbation of accretion disk becomes weak by the interaction with accretion matter, it survive and affect the disk and jet structure. • Some other modes are happen in the accretion disk. • Newtonian MHD simulations of jet formation with non-axisymmetric mode show the growth of m=2 mode but our GRMHD simulations have not seen the clear growth of m=2 mode. • We have calculated the free-free emission from a disk/outflow near a rotating black hole based on 2D GRMHD simulations using a covariant radiation transfer formulation. • The calculation shows the propagation of waves and the strong radiation from geometrically thick disk near the black hole.