The relativistic wind in PWNe Niccolo Bucciantini Astronomy Group, - - PowerPoint PPT Presentation

the relativistic wind in pwne niccolo bucciantini
SMART_READER_LITE
LIVE PREVIEW

The relativistic wind in PWNe Niccolo Bucciantini Astronomy Group, - - PowerPoint PPT Presentation

The relativistic wind in PWNe Niccolo Bucciantini Astronomy Group, NORDITA, Albanova University http://www.nordita.edu L. Del Zanna, E. Amato, D. Volpi, J.Arons, S. Komissarov, N. Camus 28/06/2011 N. Bucciantini: HEPRO III 2011 1 Tuesday,


slide-1
SLIDE 1

28/06/2011

  • N. Bucciantini: HEPRO III 2011

1

The relativistic wind in PWNe

Astronomy Group, NORDITA, Albanova University http://www.nordita.edu

  • L. Del Zanna, E. Amato, D. Volpi, J.Arons, S. Komissarov, N. Camus

Niccolo’ Bucciantini

Tuesday, June 28, 2011

slide-2
SLIDE 2

28/06/2011

  • N. Bucciantini: HEPRO III 2011

2

PWNe

  • PWNe are hot bubbles (plerions)
  • f relativistic particles and

magnetic field emitting non- thermal radiation (synchrotron - IC) from Radio to γ-ray.

  • Originated by the interaction of

the ultra-relativistic magnetized pulsar wind with the expanding SNR (or with the ISM)

  • Crab Nebula in optical: central

amorphous mass (continuum) + external filaments (lines)

Tuesday, June 28, 2011

slide-3
SLIDE 3

28/06/2011

  • N. Bucciantini: HEPRO III 2011

3

PWNe analytical MHD theory

  • Theoretical model for PWNe - 1-D steady-state (Rees

& Gunn 1974; Kennel & Coroniti, 1984) and self-similar

(Emmering & Chevalier, 1987) - free expansion phase. Basic assumptions:

  • The wind terminates with a strong MHD shock
  • Particles are accelerated at TS
  • Relativistic MHD flow in the PWN region
  • Synchrotron losses inside the nebula
  • Wind parameters derived by comparison with observations:

Tuesday, June 28, 2011

slide-4
SLIDE 4

28/06/2011

  • N. Bucciantini: HEPRO III 2011

4

Sketch of PWN / SNR interaction

  • The SNR consists of a blast

wave expanding in the ISM and by ejecta in free expansion

  • The ultrarelativistic magnetized

pulsar wind inflate a plerion inside the SNR Pulsar wind Relativistic bubble (PWN, plerion) SNR ejecta SNR shell TS CD

Tuesday, June 28, 2011

slide-5
SLIDE 5

28/06/2011

  • N. Bucciantini: HEPRO III 2011

5

Crab Nebula at various energies

RADIO X-RAYS OPTICAL γ−rays (<100 MeV)

?

Tuesday, June 28, 2011

slide-6
SLIDE 6

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Polarization radio

6

Vela, Dodson et al 03 G106.6+29, Kothes et al 06

Old nebulae interacting with ejecta - distorted stretched field in the back direction

Tuesday, June 28, 2011

slide-7
SLIDE 7

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Polarization radio

6

Vela, Dodson et al 03 G106.6+29, Kothes et al 06

Old nebulae interacting with ejecta - distorted stretched field in the back direction

Toroidal field

Tuesday, June 28, 2011

slide-8
SLIDE 8

28/06/2011

  • N. Bucciantini: HEPRO III 2011

7

Fine structures

  • Crab nebula (Weisskopf et al., 2000; Hester et al., 2002)
  • Vela pulsar (Helfand et al., 2001; Pavlov et al., 2003)

Crab Vela

Tuesday, June 28, 2011

slide-9
SLIDE 9

28/06/2011

  • N. Bucciantini: HEPRO III 2011

8

Wind models

Force-free (Contopulos et al 1999, Gruzinov 2005, Spitkovsky 2006) RMHD (Bogovalov 2001, Komissarov 2006, Bucciantini et al. 2006)

Lorentz factor ~ sin(θ) Energy flux ~ sin2(θ)

Tuesday, June 28, 2011

slide-10
SLIDE 10

28/06/2011

  • N. Bucciantini: HEPRO III 2011

9

TS structure and flow pattern

  • The wind anisotropy shapes

the TS structure. Downstream flow - equatorial collimation due to the TS shape:

  • A: ultrarelativistic pulsar wind
  • B: subsonic equatorial outflow
  • C: supersonic equatorial funnel
  • D: super-fastmagnetosonic flow
  • a: termination shock front
  • b: rim shock
  • c: fastmagnetosonic surface

σ=0.03

Tuesday, June 28, 2011

slide-11
SLIDE 11

28/06/2011

  • N. Bucciantini: HEPRO III 2011

10

Formation of polar jets by hoop stresses

σ=0.003 σ=0.03 σ=0.01

  • The global nebular flow changes with σ
  • Flow is diverted to the axis when equipartition is reached
  • For high magnetization (σ > 0.01) a supersonic jet is formed
  • Equipartition must be reached inside the PWN

Tuesday, June 28, 2011

slide-12
SLIDE 12

28/06/2011

  • N. Bucciantini: HEPRO III 2011

11

Modeling a striped wind case

  • Initial magnetic field with a narrow equatorial neutral sheet
  • Dissipation in a striped wind

b=10

Tuesday, June 28, 2011

slide-13
SLIDE 13

28/06/2011

  • N. Bucciantini: HEPRO III 2011

12

Comparison with Observations

Main torus Inner ring (wisps structure) Knot Back side of the inner ring

No jet - Axisymmetric assumption

Knot Ring Torus

Hester et al. 1995

Komissarov & Lyubarky 2004

Tuesday, June 28, 2011

slide-14
SLIDE 14

28/06/2011

  • N. Bucciantini: HEPRO III 2011

13

Striped wind and jet properties

What effect has the striped wind region size on the appearance of the jet? Where the wind magnetization is maximum? Equatorial region or polar? Wind with higher polar magnetization proposed by Arons 1998. σ=0.025, b=10 σ=0.1, b=1

Tuesday, June 28, 2011

slide-15
SLIDE 15

28/06/2011

  • N. Bucciantini: HEPRO III 2011

14

Comparison with Observations

Tuesday, June 28, 2011

slide-16
SLIDE 16

28/06/2011

  • N. Bucciantini: HEPRO III 2011

14

Comparison with Observations

Tuesday, June 28, 2011

slide-17
SLIDE 17

28/06/2011

  • N. Bucciantini: HEPRO III 2011

15

Photon Index

Brighter inner ring

Mori 02

Higher index in the torus

(recompression and boosting)

Jet not correctly reproduced

(required reacceleration)

Tuesday, June 28, 2011

slide-18
SLIDE 18

28/06/2011

  • N. Bucciantini: HEPRO III 2011

16

Time variability - wisps

  • Wisp moving outward
  • Year long limit cycle
  • Variability in the knot
  • Bubble in the jet v~ 0.6 c

Slane 05, DeLaney 06

Variability in the knot structure Jet feature moving at 0.6 c Local instabilities or global modes?

Tuesday, June 28, 2011

slide-19
SLIDE 19

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - Flow

17 Instability of the shear layers creates eddies at the rim shock Eddies are advected outward and a toroidal pressure wave is launched There is no wave reflection from the boundary Waves reflected on the axis modulate the TS shape The equatorial channel is kink unstable

Tuesday, June 28, 2011

slide-20
SLIDE 20

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - Flow

17 Instability of the shear layers creates eddies at the rim shock Eddies are advected outward and a toroidal pressure wave is launched There is no wave reflection from the boundary Waves reflected on the axis modulate the TS shape The equatorial channel is kink unstable

Tuesday, June 28, 2011

slide-21
SLIDE 21

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - SASI

18 The stable TS configuration

Shear layer Shear layer

Tuesday, June 28, 2011

slide-22
SLIDE 22

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - SASI

18

Tuesday, June 28, 2011

slide-23
SLIDE 23

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - SASI

18 The SASI TS instability

Eddy forms A pressure wave is launched inside the nebula

Tuesday, June 28, 2011

slide-24
SLIDE 24

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - SASI

18 The SASI TS instability

Eddy forms A pressure wave is launched inside the nebula Eddy advected

  • utward

Pressure wave converge to axis Pressure wave moves outward

Tuesday, June 28, 2011

slide-25
SLIDE 25

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - SASI

18 The SASI TS instability

Eddy forms A pressure wave is launched inside the nebula Eddy advected

  • utward

Pressure wave converge to axis Pressure wave moves outward Pressure wave reflected by the axis Compression waves in the jet Pressure wave compresses the TS

Tuesday, June 28, 2011

slide-26
SLIDE 26

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - Emission

19 Outgoing wave pattern Large luminosity variations Features slow down as they move

  • utward

Variability observed both in the knot and in the sprite Pressure waves produce variability in the axial emissivity Large striped wind are favored to produce a bright torus

Tuesday, June 28, 2011

slide-27
SLIDE 27

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - Emission

19 Outgoing wave pattern Large luminosity variations Features slow down as they move

  • utward

Variability observed both in the knot and in the sprite Pressure waves produce variability in the axial emissivity Large striped wind are favored to produce a bright torus

Tuesday, June 28, 2011

slide-28
SLIDE 28

28/06/2011

  • N. Bucciantini: HEPRO III 2011

MHD variability - High Energy

20

Wilson-Hodge et al. 2010

~ 2 yr Timescale

Tuesday, June 28, 2011

slide-29
SLIDE 29

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Flares

21

Not from pulsar :

  • flares are not pulsed or

in phase

  • no variations in the

timing residual

Sep 2010 Flare Apr 2011 Flare! Feb 2009 Flare

Wilson-Hodge et al. 2011

Tuesday, June 28, 2011

slide-30
SLIDE 30

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Flares

22

Unlikely MHD origin like the slow variability of the wisps:

  • MHD effects are

achromatic

  • size of the accelerator is

very small (day-light)

  • unlikely high magnetic

field

Wilson-Hodge et al. 2011

Tuesday, June 28, 2011

slide-31
SLIDE 31

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Flares

22

Unlikely MHD origin like the slow variability of the wisps:

  • MHD effects are

achromatic

  • size of the accelerator is

very small (day-light)

  • unlikely high magnetic

field

Wilson-Hodge et al. 2011

Electrostatic acceleration?

Tuesday, June 28, 2011

slide-32
SLIDE 32

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Free Expansion into Ejecta

23

  • Continuous energy injection - High synchrotron luminosity
  • PWN expands supersonically, RPWN ∝ t 6/5
  • Pulsar at the center of PWN

SNR G21.5-0.9 (X-rays) Matheson & Safi-Harb 2005 SNR G11.2-0.3 (X-rays) Kaspi et al. 2001

Tuesday, June 28, 2011

slide-33
SLIDE 33

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Reverberation

24

  • Reverse interacts with PWN after time
  • Compression; synchrotron burn-off at high energies
  • Effects of inhomogeneous ISM
  • Offset pulsar; filamentary structure; mixing

t ~ 7M10M sun

5/ 6

E51

−1/ 2n0 −1/ 3 kyr

van der Swaluw et al. (2004)

forward shock PWN reverse shock

Tuesday, June 28, 2011

slide-34
SLIDE 34

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Reverberation

24

  • Reverse interacts with PWN after time
  • Compression; synchrotron burn-off at high energies
  • Effects of inhomogeneous ISM
  • Offset pulsar; filamentary structure; mixing

t ~ 7M10M sun

5/ 6

E51

−1/ 2n0 −1/ 3 kyr

Vela (radio) Duncan et al. 1996 van der Swaluw et al. (2004)

forward shock PWN reverse shock

Tuesday, June 28, 2011

slide-35
SLIDE 35

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Transition to bow shock

25

  • PWN expands into shocked ejecta
  • “Relic” radio PWN left behind
  • New PWN around pulsar (X-ray)

van der Swaluw (2004)

Tuesday, June 28, 2011

slide-36
SLIDE 36

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Transition to bow shock

25

  • PWN expands into shocked ejecta
  • “Relic” radio PWN left behind
  • New PWN around pulsar (X-ray)

SNR G327.1-1.1, Gaensler & van der Swaluw (2004)

Tuesday, June 28, 2011

slide-37
SLIDE 37

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Transition to bow shock

25

  • PWN expands into shocked ejecta
  • “Relic” radio PWN left behind
  • New PWN around pulsar (X-ray)

SNR G327.1-1.1, Gaensler & van der Swaluw (2004) SNR W44 (Frail et al. 1996, Giacani et al. 1997)

Tuesday, June 28, 2011

slide-38
SLIDE 38

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Bow Shock PWN

26

  • Most pulsars kick velocity is supersonic in ISM
  • Forward shock visible in Hα
  • PWN visible as a radio and X-rays tail

Hα X- rays

15”

PSR B1957+20 (Stappers et al. 2003) PSR B2224+65 (Chatterjee & Cordes 2002)

Tuesday, June 28, 2011

slide-39
SLIDE 39

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Are we over-biased by Crab ?

27

Not every PWN is Crab Not all Crab properties are universal

Tuesday, June 28, 2011

slide-40
SLIDE 40

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Are we over-biased by Crab ?

27

Not every PWN is Crab Not all Crab properties are universal

Tuesday, June 28, 2011

slide-41
SLIDE 41

28/06/2011

  • N. Bucciantini: HEPRO III 2011

Are we over-biased by Crab ?

27

Not every PWN is Crab Not all Crab properties are universal

Tuesday, June 28, 2011

slide-42
SLIDE 42

28/06/2011

  • N. Bucciantini: HEPRO III 2011

28

Summary and conclusions

  • PWNe are well defined system, where a base agreement on engine and

emission properties has beed reached

  • There is now a “canonical” model which can explain many of the
  • bserved features
  • There are still open areas of research for old objects
  • Tools for PWNe study are easy to use, and widely available now
  • Still hard to derive fully consistent models that can account for all of the
  • bserved properties
  • Holes remain in the entire picture that NEED to be filled
  • Not all PWNe are Crab “clones” - Nor should every model!

Tuesday, June 28, 2011

slide-43
SLIDE 43

28/06/2011

  • N. Bucciantini: HEPRO III 2011

28

Summary and conclusions

  • PWNe are well defined system, where a base agreement on engine and

emission properties has beed reached

  • There is now a “canonical” model which can explain many of the
  • bserved features
  • There are still open areas of research for old objects
  • Tools for PWNe study are easy to use, and widely available now
  • Still hard to derive fully consistent models that can account for all of the
  • bserved properties
  • Holes remain in the entire picture that NEED to be filled
  • Not all PWNe are Crab “clones” - Nor should every model!

Thank you

Tuesday, June 28, 2011