Diversity of Multi-wavelength Behavior of Relativistic Jet in 3C 279 Discovered During the Fermi Era

Transcription

1 Diversity of Multi-wavelength Behavior of Relativistic Jet in 3C 279 Discovered During the Fermi Era Rapid Variability of Blazar 3C 279 during Flaring States in with Joint Fermi-LAT, NuSTAR, Swift, and Ground- Based Multi-wavelength Observations Hayashida, Nalewajko, Madejski, Sikora+,2015, ApJ, in press (arxiv: ) M.Hayashida et al. 2015, ApJ, in press (arxiv: ) Relativistic Jets: Creation, Dynamics and Internal Physics in Krakow, Poland, 23 April 2015 Masaaki Hayashida (Institute for Cosmic-Ray Research, the University of Tokyo) 1

2 Emission from Jets (FSRQ) 3C 279 (MH, Madejski, Nalewajko, Sikora+12, ApJ, 754, 114)! synchrotron! (external) inverse- Compton! 43 GHz! 230 GHz! gamma ray! 23 April

3 Outline MWL observations in for 3C 279 during the flaring states! Fermi-LAT, NuSTAR, Swift and optical, radio! The First NuSTAR observations for the source! origin of the X-ray emission! orphan γ-ray flare in 2013 December.! where is the gamma-ray emission site?! what is the dominant component in jet?! what is the acceleration mechanism?! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 3

4 where is gamma-ray emission site? (HESS 07)! accre?on disk BH Dusty torus magnetic dominated?! pc (~10 3 R g ) Broad emission line region (BEL) e + e - UV photon γ γ particle dominated?! a few pc (~10 5 R g ) IR photon e + e - γ γ (fast variability in BL Lacs)! (TeV detection from FSRQs)! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 4

5 FSRQ 3C 279 (z=0.536) MAGIC detection (>100 GeV) in 2006 (MAGIC coll. 2008, Science)! Fermi+MWL in (Abdo+ 2010, Nature; Hayashida+ 2012, ApJ)! γ X-ray optical short time (~day) flare! PA <bent jet model: (pc scale)>! Nature, vol436, 887 (2010) A.Young PD radio Γ jet Γ jet no simultaneous HE (LAT) + VHE (IACT) spectra, 5 yet!

6 3C 279 activity for 6 years 2008 August 2014 August measured by Fermi-LAT! 2013/11/ /04/28 (MJD )! Gamma-ray flare activity reported in Atel! 2013/12/21: #5680 Fermi LAT detection of a GeV flare from the FSRQ 3C 279! 2014/04/01: #6036 Fermi LAT detection of renewed GeV activity from blazar 3C 279! 6

7 Fermi-LAT light curve reached! 10-5 ph/cm 2 /s! level in flux (> 100 MeV)!! only a several FSRQs show! such a high flux!! c.f.! Crab:! 2.5x10-6 ph/cm 2 /s!!! 23 April

8 Flare profile asymmetric profile! hourly scale variability: (2hr à R [emission region size] < ~ 4x10 15 (Γ/20) cm )! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 8

9 short time variability in FSRQs Fermi-LAT (>100 MeV)! 3C454.3 (Abdo+11, ApJL)! F peak ~8x10-5! 3/6 h bin! PKS (Saito+13 ApJL)! F peak ~4.5x10-5! 3 h bin! (see more in Saito s talk tomorrow)! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 9

10 LAT Spectrum Very hard index (1.71±0.10)! peaked at a few GeV! 10

11 Hard spectra in FSRQs 3C 279 (this work)! 1.71±0.10! Fermi 3rd AGN catalog (3LAC)! (= average state)! preliminary! (Pacciani+14, ApJ, flaring FSRQs)! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 11

12 Multi-band light curve Period (B):! no flare in other bands orphan γ-ray flare! (A), (C)! NuSTAR observations! Kanata, SMART! Kanata! SMA! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 12

13 X-ray bands (Swift+NuSTAR) Two NuSTAR observations (1: Dec , 2: Jan ) : ~ 40 ks exposure for each! Swift(XRT)+NuSTAR! Swift-XRT! photon index (Swift-XRT: > 0.5 kev) X-ray flux (0.5-5 kev: 10 erg cm -2 s -1 ) 2nd! 1st! 1st: simple power law: 1.74±0.01! 2nd : double broken power law:! 1.37±0.3, 1.72±0.02, 1.81±0.02! 23 April

14 X-ray bands (Swift+NuSTAR) Two NuSTAR observations (1: Dec , 2: Jan ) : ~ 40 ks exposure for each! Swift(XRT)+NuSTAR! <Light curve in the 2nd obs.>! NuSTAR! Swift-XRT! 2nd! 1st! EIC by low energy electrons? (should not be variable in a day)! SSC? (too hard of Γ x =1.37 in the Swift-XRT band)! another region? (slower (sheath) part? )! 14

15 Broad band SED (3 days)! (0.2 days)! (3 days)! (0.27 days)! 23 April

16 constraints of emission model parameter (see details in Nalewajko s talk tomorrow)! r BLR r IR 10 3 R g 10 4 R g 10 5 R g j B1 j j = 1 B2 L SSC = erg s -1 (u /u B = 200) L SSC = erg s -1 (u /u B = 2000) SSA = 300 µm L j,min = erg s -1 L j,min = 3x10 45 erg s E max,obs [TeV] 10 E cool = 100 MeV j j = 0.1 j j = C 279 MJD L = 6.0e+48 erg s -1 t var,obs = 2.0 h 0.01 pc 0.1 pc 1 pc 10 pc r r (distance from the center)! 16

17 -8 emission model for Period B mm µm kev MeV GeV one-zone leptonic model: BLAZAR (Moderski+2003)! 49 log 10 F [erg s -1 cm -2 ] SYN SSC ERC log 10 L [erg s -1 ] log 10 [Hz] B1 B2 D1 D2 1. Gamma-ray emission site should be inside BLR (< 0.1 pc)! 2. very matter dominated jet: L B /L jet ~ 10-4! 3. hard index (γ-ray band) in the fast cooling regime! 45 required very hard index for electron injection spectrum: p=1! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 17

18 -8 emission model for Period B mm µm kev MeV GeV one-zone leptonic model: BLAZAR (Moderski+2003)! 49 log 10 F [erg s -1 cm -2 ] SYN SSC ERC log 10 L [erg s -1 ] log 10 [Hz] B1 B2 D1 D2 1. Gamma-ray emission site should be inside BLR (< 0.1 pc)! 2. very matter dominated jet: L B /L jet ~ 10-4! 3. hard index (γ-ray band) in the fast cooling regime! 45 required very hard index for electron injection spectrum: p=1! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 18

19 where is the emission region? (this work)! 2009 (Hayashida+12)! -8 mm µm kev MeV GeV 49 F [erg s -1 cm -2 ] 1e-09 1e-10 1e-11 Period D Period E ev kev MeV GeV log 10 F [erg s -1 cm -2 ] SYN SSC ERC log 10 L [erg s -1 ] 1e-12 1e-13 1e+09 1e+12 1e+15 1e+18 1e+21 1e+24 [Hz] a few pc (outside BLR)! log 10 [Hz] B1 B2 D1 D pc (inside BLR)! 45 emission site is not unique!! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 19

20 Multi-components for gamma-ray emission Just examples! PKS (z=0.31)! 2011 (Nalewajko+12)! PKS (z=1.522)! (Buson+14)! See also Fink&Demer+10 for 3C454.3! and many other works! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 20

21 -8 emission model for Period B mm µm kev MeV GeV one-zone leptonic model: BLAZAR (Moderski+2003)! 49 log 10 F [erg s -1 cm -2 ] SYN SSC ERC log 10 L [erg s -1 ] log 10 [Hz] B1 B2 D1 D2 1. Gamma-ray emission site should be inside BLR (< 0.1 pc)! 2. very matter (kinetic power) dominated jet: L B /L jet ~ 10-4! 3. hard index (γ-ray band) in the fast cooling regime! 45 required very hard index for electron injection spectrum: p=1! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 21

22 Regions of AGN Jet Propagation Jet Launching Region Jet Collimation Region ( Launching Region) Modified from Graphic courtesy David Meier Alfven Point Modified Fast Point Sheath Slow MS Point Poynting Flux Dominated CD Unstable Magnetic Helicity Driven Region Fast MS Point High speed spine Combined CD/KH Unstable Region Collimation Shock Kinetic Energy Flux Dominated with Tangled (?) Field KH Unstable Velocity Shear Driven Region 23 April 2015 slide by Y.Mizuno! Masaaki Hayashida (ICRR) in Jets 2015 Krakow 22

23 Poynting flux dominated? Kinetic energy flux dominated? if jet is derived by the magnetic field (e.g., Blandford-Znajek process),,,,! à jet should be Poynting-flux dominated jet < 10 3 r g (= inside BRL)! Leptonic models can explain well the broad band SED inside BLR (0.03 pc < 10 3 r g for 5x10 8 M solar )! the emission model results suggest kinetic energy dominated jets (some models with equipartition see e.g., Dermer+14, ApJ, 782 for 3C 279)! Hadronic models require stronger magnetic fields ( G) than the Leptonic models ( G), but also requires very high power of relativistic protons, erg/s (e.g.,zdziarski & Boettcher 15)! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 23

24 -8 emission model for Period B mm µm kev MeV GeV one-zone leptonic model: BLAZAR (Moderski+2003)! 49 log 10 F [erg s -1 cm -2 ] SYN SSC ERC log 10 L [erg s -1 ] log 10 [Hz] B1 B2 D1 D2 1. Gamma-ray emission site should be inside BLR (< 0.1 pc)! 2. very matter dominated jet: L B /L jet ~ 10-4! 3. hard index (γ-ray band) in the fast cooling regime! 45 required very hard index for electron injection spectrum: p=1! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 24

25 hard (p<2) electron index p : injected electron index! p 2: normal standard shock (Fermi-I) acceleration too soft!!! magnetic reconnection! (Giannios+09)! Our result:! jet magnetization: σ < 10-3! (e.g., Guo+14, PRL)! the reconnection will efficiently work in this condition?! very localized acceleration sites?! can generate erg/s emission?! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 25

26 Stochastic acceleration (Fermi-II) (Model: Asano+2014, ApJ 784, 64)! Steady outflow Continuous shell ejection with a width of R 0 /Γ in commoving frame Electron injection from R=R 0 to 2R 0 with stochastic acceleration Turbulence Index: q=2 (hard-sphere scattering) Both injection and acceleration stop at R=2R 0 Physical Processes Electron injection Stochastic acceleration Synchrotron emission and cooling Inverse Compton emission and cooling Adiabatic cooling (V R 2 ) Photon escape No electron escape!

27 Steady (base line) model A high state in 2009 as reference! (Asano & Hayashida, in prep)! ε ε ε ε electrons in one shell! emission (photons) from (data from MH+12)! the all shells of R > R 0! ε ε R 0 = pc, Γ = 15, B 0 = 7 G! K (energy diffusion coefficient) = s 1! N e (electron injection rate) = s 1! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 27

28 application for the 2013 flare emission from a single shell! R 0 = pc, Γ = 15, B 0 : 7 G à 0.25 G! K: s 1 à s 1! N e : s 1 à s 1! ε ε low B in the dflare??! electron! ε ε photon! (Asano & Hayashida, in prep)! light curve (> 100 MeV)! ε ε the turbulence is generated by the hydrodynamical instability?!

29 Summary & Conclusion 3C 279 showed the highest γ-ray flux level in ! orphan γ-ray flare was detected! where is the gamma-ray emission site?! inside BRL (~0.03 pc < 10 3 r g ) for hourly scale variability at 100 MeV! (both inside and outside BLR ( to r g ) event by event)! what is the dominant component in jet?! emission model : kinetic energy dominated : L B /L jet ~ 10-4! jet simulation: Poynting-flux dominated (< 10 3 r g )! Any ideas for this issue?! what is the acceleration mechanism?! not only shock accelerations! stochastic acceleration (Fermi-II) can also work for rapid γ-ray flares! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 29

30 back up! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 30

31 energetics Lγ erg s -1! Lj Lγ/(ηΓ 2 ) (η=0.1)! Lj erg s -1! L B erg s -1,! Lj/L disk 25! L B /Lj 10-5! Lj > L Edd! L disk erg s -1! L Edd erg s -1! M BH M,! 23 April 2015 Masaaki Hayashida (ICRR) in Jets 2015 Krakow 31