Katsuaki Asano ICRR U Tokyo Collaborators Yuto Teraki YITP Kyoto Masaaki Hayashida ICRR Fumio Takahara Osaka etc Nonthermal Emission Supernova Remnant SN1006 ID: 797439
Download The PPT/PDF document "Turbulence Acceleration Model for the Br..." is the property of its rightful owner. Permission is granted to download and print the materials on this web site for personal, non-commercial use only, and to display it on your personal computer provided you do not modify the materials and that you retain all copyright notices contained in the materials. By downloading content from our website, you accept the terms of this agreement.
Slide1
Turbulence Acceleration Model for the Broad Band Blazar Spectra
Katsuaki
Asano
(
ICRR, U. Tokyo
)
Collaborators:
Yuto
Teraki
(YITP, Kyoto), Masaaki
Hayashida
(ICRR),
Fumio
Takahara
(Osaka) etc.
Slide2Non-thermal Emission
Supernova Remnant SN1006
Consistent with the diffusive shock acceleration (Fermi I)+Bohm Limit!
All non-thermal phenomena are due to shocks?
Acero
+ 2010
Electron energy distribution:
Single power-law with the index of 2.1,
cut off at 10 TeV and B=30μG.
Bohm Limit (
) &
Open Problems in Blazar Emission
Kino,
Takahara
,
Kusunose
2002
1.
Index harder than
2 for the electron injection spectrum.2. Lower maximum energyMrk 421:B=38mGIf even for .But actual Max. Energy
3. Too sharp break in the electronspectrum.Compared to the cooling break,the difference in the indicesseems large.Inoue & Takahara 2002
The Bohm factor should be
Electron Index in Fermi BL Lacs
Yan+
Slide5Many Parameters in Broken-power-law fit
Need
double break
and
low-energy cut-off
Radio should be different component in this case
Slide6Particle Acceleration in Blazar Jets
Shock Acceleration may be hard in highly magnetized plasma.
Kirk &
Heavens 1989, Lemoine
& Pelletier 2010, Sironi+ 2013 etc.In blazar
Sironi+ 2015
The parameter range is very limited.Acceleration via turbulence
Park & Petrosian 1995
Becker & Dermer 2006Stawarz & Petrosian 2008
Bottcher+ 1999Schlickeiser & Dermer 2000Lefa+ 2011Kakuwa+ 2015In BlazarsIn GRBsBykov & Meszaros 1996Dermer & Humi 2001Asano & Terasawa 2009, 2015Murase+ 2012
Slide7Phenomenological Treatment
Kolmogorov+Alfvenic
q=5/3, Compressible q=2 (hard sphere)
Diffusion
Acceleration
Cooling
Injection
No cooling, continuous injection, time evolution
harder than the shock case
Stochastic Acceleration due to turbulence
(2
nd
order Fermi acc.) leads to a hard spectrum.
Slide8Interaction with Alfven Wave
Energy gain per scattering
Alfvenic
Wave (transverse/incompressible)
pitch angle diffusion
→
mean free path
resonance condition
However, the ratio
is typically 0.1 in blazars.
It seems that the Alfven waves cannot be the
enrgy
source.
Steady outflow
Continuous shell ejection with a width of R
0
/Γ in commoving frame
Elecrton injection from R=R
0
to 2R
0
with stochastic accelerationTurbulence Index: Kolmogorov q=5/3Both injection and acceleration stop at R=2R0Our ModelElectron injectionStochastic accelerationSynchrotron emission and coolingInverse Compton emission and coolingAdiabatic cooling (V∝R2)Photon escapeNo electron escape!
Physical Processes
Slide10Simple case
Constant injection & constant diffusion coefficient
Electron spectrum
Energy density ratio
Slide111ES 1101-232
The simple model can reproduce the spectrum of 1ES 1101-232
Slide12But, in Mrk 421…
Slide13But, in Mrk 421…
To produce a soft spectrum, increasing injection rate
with radius is assumed.
Hard to reconcile GeV
data.
Magnetic Energy is subdominant.
Alfven Wave is actually responsible?
Slide14Alternative: Compressible Waves
Fast
Slow
x5 speed
Acceleration by compressible wave
Transit Time Damping (Scatter by mirror force)
Formulae of Energy Diffusion Coefficient:
Ptuskin
1988
Cho &
Lazarian
2006
Lynn+ 2014
Non-resonant (need spatial diffusion
coeff
.)
Resonant (
)
The largest scale eddy dominates
The amplitude
is dominated by the
smallest scale for q<2.
is written with the smallest eddy?
Projection
Test particle
simulation in pure linear
Fast Waves
Fast
, 10 mode
Kolmogorov
High-beta plasma
Teraki
& Asano in prep.
Slide17Results of the Test Particle Simulation
Inject
(weak
turb
.)
While only a small fraction of particles of
can be accelerated for
,
almost all of relativistic particles are accelerated in
Preliminary
Mirror force:
Interaction Time scale:
Electric Field:
Energy Gain Rate:
Hard Sphere
Slide18Hard-Sphere Model
Parameters are
Required for any model
Only
one
peculiar parameters
Cho &
Vishniac
2000
Supposing the compressional waves
are responsible, we model as following.
Parameters are constant during the dynamical
time scale
, and injection and acceleration
suddenly shutdown after that.
Mrk421
Low maximum energy and curved electron spectrum are naturally reproduced
with temporal evolution effects.
See also
Kakuwa
+ 2015
Slide20Diffusion coefficient in Mrk 421
If we assume
Adopting our formula for the index of 5/3
Then,
While the required value:
Reasonable
Slide21Very Hard Spectrum in 3C 279
Broken
power-law
Parameters
Slide223C 279 Steady Model
Hayashida+2012
Model for the Steady State
SSC Component constrains
the emission radius
Asano &
Hayashida
2015
Slide23Flare Model
The same radius, Lorentz factor, and UV field as those in the steady model.
Light curve
Slide25Energy Density
Drenkhahn
2002
Assumed Values
Very Weak Magnetic Field.
Negative for
the magnetic reconnection,and the jet acceleration by magnetic dissipation.
Extreme Hard Blazar 1ES 1101-232
Electron spectrum
Photon
spectrum
q=5/3
Asano+ 2014
Expanding jet
→
No
steady state
→
Temporal evolution is essential
Slide27FSRQ 3C 279
Hayashida+2012
Steady State
SSC Component constrains
the emission radius
Asano &
Hayashida
2015
Flare State
RX J1136.5+6737
Preliminary
MAGIC team
1ES 1959+650
Preliminary
PKS 1510−089
Preliminary
Slide31Summary
In relativistic jets, acceleration by fast wave turbulence may be expected.
Energy diffusion timescale becomes time-independent=Hard Sphere.
Simple model agrees with Blazar spectra, which supports the Hard Sphere acceleration.
Mizuno+ 2007
Slide32予備スライド
Slide33Test Particle Simulation
Particle injection in
MHD simulations,
where the s
low-mode is dominant.
Lynn+ 2014
Slide343C 279 Flare
Lightcurve
Evolution of energy density
Slide35Big Flare in FSRQ 3C 279
Hayashida
+ 2015
Slide36Spectrum
Extremely hard.
Stochastic Acc.
seems preferable.
Broken
power-lawModel parameters
Slide37Magnetic
Pressure Dominant Case
Fast Mode
Transversal to longitudinal
Nearly Perpendicular to B
Slow Mode
Longitudinal to Transversal
Nearly Parallel to B
Slide38RM不安定性による磁場の増幅・乱流スペクトル
Kolmogorov