Kinetic Plasma Simulations in Astrophysics. Lorenzo Sironi

Size: px
Start display at page:

Download "Kinetic Plasma Simulations in Astrophysics. Lorenzo Sironi"

Transcription

1 Kinetic Plasma Simulations in Astrophysics Lorenzo Sironi

2 Outline Plasma physics in astrophysics. The Vlasov-Maxwell system. Fully-kinetic particle-in-cell codes. 1. Electrostatic codes. 2. Electromagnetic codes. Hybrid particle-in-cell codes. Astrophysical applications.

3 Time and length scales Characteristic time and length scales " skin = c /# p skin depth " c = eb mc Larmor Full kinetic models Hybrid models Fluid models MHD

4 Time and length scales Characteristic time and length scales " skin = c /# p skin depth " c = eb mc Larmor Full kinetic models Hybrid models Fluid models MHD [credit: A. Spitkovsky]

5 Binary vs collective effects Number of particles in Debye sphere If ND 1, binary interactions ( collisions ) are important. If ND 1, collective effects dominate. Most astrophysical plasmas have ND 1.

6 Kinetic plasma simulations

7 Lagrangian vs Eulerian LAGRANGIAN: observe/follow the motion of an individual element (fluid parcel, particle) as it moves through space and time. EULERIAN: focus on specific locations in space through which the various elements (fluid parcels, particles) pass in time.

8 The Vlasov equation LAGRANGIAN EULERIAN Collisionless Vlasov equation: i.e., conservation of the phase space density along characteristics (Liouville s theorem) The fluid equations follow from Vlasov equation (but closure problem!)

9 VLASOV { MAXWELL The Vlasov-Maxwell system

10 How to solve the Vlasov equation? Two options: Discretize the Vlasov equation on a grid in phase space: 1. computationally expensive to solve in 6+1 dimensions 2. how to determine the boundaries of the grid in momentum space? 3. what if f<0? Sample the phase space density with particles, and follow them as LAGRANGIAN tracers. PIC codes in a nutshell:

11 Super-particles vs real particles VLASOV - MAXWELL PIC codes follow the trajectories of super-particles (constancy of f constancy of particle strength/charge) A super-particle is a computational particle that represents many real particles. This is OK: the Lorentz force depends only on q/m. The number of super-particles should be chosen so that the phase space is well sampled (good statistics!).

12 Numerical noise in PIC The particle granularity gives short-scale fluctuations of the electromagnetic fields, whose mean amplitude scales (Poissonlike) as n, where n is the particle density. The fractional contribution of the fluctuations (over the slowly varying fields) scales as 1/ n. This is problematic because the number of super-particles in particle-in-cell codes is number of real particles. We need to control the level of the fluctuations such that they give negligible effects over the timespan of the simulations.

13 Particles of finite size Since effects from binary interactions are negligible in astrophysical plasmas (ND 1), we want to minimize the effect of short-range Coulomb forces. Yet, computational super-particles (with charge e) over-estimate the strength of the Coulomb force, as compared to real particles. If super-particles are not point particles, but rather have finite size, the strength of binary interactions is reduced.

14 The particle-in-cell (PIC) loop

15 Particle-particle vs particle-in-cell PARTICLE-PARTICLE PARTICLE-MESH Calculating inter-particle forces directly is extremely expensive. Warnings: - instantaneous forces (c=infinity) - in a plasma, no binary forces for Debye length Computationally cheaper. The finite light propagation time is self-consistently included. Only macroscopic e.m. fields.

16 The basic PIC loop

17 Electrostatic codes

18 Electrostatic codes Timescales of the system >> light crossing time; magnetic fields static. [credit: A. Spitkovsky]

19 [credit: A. Spitkovsky]

20 [credit: A. Spitkovsky]

21 [credit: A. Spitkovsky]

22 [credit: A. Spitkovsky]

23 [credit: A. Spitkovsky]

24 main lobe:

25 Electromagnetic codes

26 The basic PIC loop

27 1. The field solver In electromagnetic PIC codes, only two equations need to be solved. The other two are satisfied as initial conditions, and they continue to be satisfied for appropriate choices of the numerical scheme.

28 STAGGERING in time (leapfrog): second-order accurate in time 1. The field solver STAGGERING in space (Yee s mesh): electric fields on cell edges, magnetic fields on cell faces second-order accurate in space maintains divergence-free B

29 2. The field interpolation The fields obtained from Maxwell s equations are determined only at the grid points, they need to be interpolated to the particle positions. The interpolation is done by assuming a particle shape function. The shape function needs to be: 1. isotropic 2. zero outside some range 3. higher order B-splines are computationally more expensive, but more accurate and less collisional

30 2. The field interpolation Area-weighting scheme The same shape function should be used for both field interpolation and current deposition, to ensure: global momentum conservation no self-force

31 3. The particle pusher If the number of ppc is >>1, most of the computing time is spent in pushing the particles. The BORIS pusher (leapfrog method) advance the position advance the momentum (semi-implicit scheme)

32 3. The particle pusher (A) Advance the momentum: electric acceleration (B) Advance the momentum: magnetic rotation no ambiguity! Can overstep magnetic rotation without stability issues!!! More precisely:

33 4. The current deposition Charge conservation is required to satisfy Poisson s equation. The current deposition scheme needs to be charge-conserving. Or, a divergence-cleaning solver should be employed.

34 Assume a plane electromagnetic wave in vacuum: (E,B)=(E 0,B 0 ) exp(ik x-iωt) Dispersion relation on the grid coming from the finite difference scheme of the field solver. Small t is required for stability in explicit schemes. Large t can be achieved in implicit schemes (but still need to resolve the relevant physics!). Stability of PIC codes In addition, to correctly capture the plasma physics, we need: Debye length plasma frequency

35 Numerical Cerenkov instability For relativistic particles with v>vph, unwanted grid- Cerenkov radiation can occur. The particles lose energy and momentum to grid- Cerenkov radiation. (partial) Solution: 4th-order spatial stencil of field solver. (partial) Solution: filter electric currents.

36 Boundary conditions and parallelization

37 Boundary conditions BCs for the fields: periodic conducting (absorbing PML) BCs for the particles: periodic reflecting (absorbing)

38 Parallelization

39 Parallelization Distributed memory systems using MPI Field solver and particle pusher are local - Communication only with neighboring nodes Each computing node acts as a section of the global physical space We will show results obtained with the PIC code TRISTAN-MP (Spitkovsky 05) Strong Scaling (fixed total load) Weak Scaling (fixed load per CPU) 1 node 1 node

40 Summary (PIC codes) The PIC model is a first-principle kinetic model for plasma (astro)physics It efficiently tracks particle trajectories to sample the particle phase space Its structure is scalable to massively parallel computer systems

41 Hybrid codes

42 Physical scales Electron time and length scales are << proton time and length scales (credit: D. Burgess)

43 Hybrid-PIC codes Hybrid approximation: kinetic ions (as in fully-kinetic PIC), fluid electrons with eq. of state ion kinetic physics can be captured, electron kinetic physics cannot the displacement current is neglected (non-relativistic) electrons are massless (non-relativistic) quasi neutrality

44 Particles in hybrid codes Ions are leap-frogged via the Boris pusher (as in PIC codes) semi-implicit for the velocity Electrons follow the equation of motion of a massless fluid EQUATION FOR THE ELECTRIC FIELD

45 Electromagnetic fields in hybrid codes No time advance of E is needed Time advance of B Implementation: two-step Lax-Wendroff with a spatially staggered grid (Yee mesh); assume to have B, E at step n, ion positions at n & n+1 and ion velocities at n+1/2 Lax-Wendroff

46 Complications in the field solver Now we need to determine but we do not know Common approaches: Predictor-corrector: predict fields at n+1, advance particles up to n+3/2, compute fields at n+3/2, average fields at n+1/2 and n+3/2 to obtain fields at n+1. extrapolate E fields extrapolate ion fluid speed CAM-CL (Current Advance Method - Cyclic Leapfrog)

47 Low-density regions in hybrid codes In regions of low density, the Alfven velocity increases, and the electric field explodes. Possible solutions: in the empty regions, set the density and the electric current to zero, and solve Laplace s equation. increase the number of ppc add a density floor (but unphysical waves will grow!)

48 Parallelization and scalings Distributed memory systems using MPI Field solver and particle pusher are local Each computing node acts as a section of the global physical space We will show results obtained with the hybrid code dhybrid (Gargate et al. 07) Strong Scaling (fixed total load) Weak Scaling (fixed load per CPU)

49 Summary (hybrid codes) The hybrid method treats ions kinetically (as PIC) and electrons as a fluid. It is faster than PIC codes, since it neglects the kinetic scales of the electrons. References: * D. Winske & K. B. Quest, 1986 and 1988 * A. P. Matthews, JCP 1994 * A. S. Lipatov, The hybrid multiscale simulation technology, Springer

50 Astrophysical applications

51 The collisionless life of galaxy clusters merger shocks in the outskirts of galaxy clusters

52

53 The collisionless life of relativistic jets Reconnection and shocks in astrophysical jets

54

55 The collisionless life of accretion disks Anisotropy-driven instabilities regulate electron heating and angular momentum transport in collisionless accretion flows, like Sgr A* at our Galactic center.

56 Blazars Blazars: jets from Active Galactic Nuclei pointing along our line of sight 90º 60º 30º 0º 3C 279 Blazars broadband spectrum, from radio to γ- rays (and even TeV energies) low-energy synchrotron + high-energy inverse Compton (IC) from non-thermal particles high degree of radio and optical polarization magnetic fields (Boettcher 07)

57 Internal dissipation in blazar jets 3C 120 Internal shocks in blazars and gamma-ray burst jets: trans-relativistic (γ 0 ~a few) magnetized (σ>0.01) δ (mas) = B n 0 m p c 2 toroidal field around the jet field to the shock normal RA (mas) Internal Dissipation: Shocks or Reconnection? B0 θ~90 γ 0

58 γ dn/dγ Shocks: no turbulence no acceleration σ=0.1 θ=90 γ 0 =15 e - -e + shock Density The Fermi process Maxwellian ε B γβ x Momentum space No returning particles Shock γ No returning particles No self-generated turbulence No self-generated turbulence No particle acceleration (LS+ 13, LS & Spitkovsky 09,11) Strongly magnetized (σ>10-3 ) quasi-perp γ 0 1 shocks are poor particle accelerators: B0 θ σ is large particles slide along field lines θ is large particles cannot outrun the shock unless v>c ( superluminal shock) Fermi acceleration is generally suppressed

59 Are relativistic shocks always inefficient? = B n 0 m p c 2 Internal shocks in blazars and gamma-ray burst jets: γ 0 ~a few quasi-perpendicular shocks σ~ internal shocks (Gehrels et al 02) external shocks Gamma-ray burst external shocks: γ 0 ~a few hundreds perpendicular shocks σ~10-9 B0~0

60 High-σ vs low-σ shocks B0 High-σ shocks: no returning particles no turbulence Density σ=0.1 perp shock γ 0 =15 e - -e + ε B γβ x (LS & Spitkovsky 11) Low-σ shocks: returning particles oblique & filamentation instabilities σ=0 γ 0 =15 e - -e + Density ε B returning stream γβ x incoming stream (LS et al 13)

61 B0 Low-σ shocks are filamentary Mediated by the filamentation (Weibel) instability, which generates small-scale sub-equipartition magnetic fields. σ=0 γ0=15 e--e+ shock y [c/ωp] εb εb x-xsh [c/ωp] Density z (LS et al. 13)

62 The Fermi process in low-σ shocks B0 σ=0 γ 0 =15 e - -e + shock ε B γ (LS et al 13) Particle acceleration via the Fermi process in self-generated turbulence, for initially unmagnetized (i.e., σ=0) or weakly magnetized flows.

63 Relativistic reconnection within jets reconnecting field v X v reconnection electric field reconnecting field Relativistic Reconnection = B n 0 m p c 2 1 v A c

64 Hierarchical reconnection σ=10 electron-positron Density box expanding at c y B0 x ε B y x (LS & Spitkovsky 14) Reconnection is a hierarchical process of island formation and merging. The field energy is transferred to the particles at the X-points, in between the magnetic islands.

65 Hierarchical reconnection σ=10 electron-positron Density ε B ε B -ε E The current sheet breaks into a series of secondary islands (e.g., Loureiro+ 07, Bhattacharjee+ 09, Uzdensky+ 10, Huang & Bhattacharjee 12, Takamoto 13). The field energy is transferred to the particles at the X-points, in between the magnetic islands. Localized regions exist at the X-points where E>B. (LS & Spitkovsky 14)

66 Inflows and outflows σ=10 electron-positron Density Inflow speed vin/c Outflow speed Inflow into the X-line is non-relativistic, at vin ~ 0.1 c (Lyutikov & Uzdensky 03, Lyubarsky 05) Outflow from the X-points is ultra-relativistic, reaching the Alfven speed r v A = c 1+

67 The particle energy spectrum σ=10 electron-positron MB p=2 Time At late times, the particle spectrum in the current sheet approaches a power law dn/ dγ γ - p of slope p~2. The normalization increases, as more and more particles enter the current sheet. The mean particle energy in the current sheet reaches ~σ/4 rough energy equipartition (LS & Spitkovsky 14) γmax γ max t Time [ω p -1] The max energy grows as γ max t (compare to γ max t 1/2 in shocks).

68 tearing drift-kink (LS & Spitkovsky 14) In 3D, the in-plane tearing mode and the out-of-plane drift-kink mode coexist. The drift-kink mode is the fastest to grow, but the physics at late times is governed by the tearing mode, as in 2D.

69 3D: particle spectrum σ=10 electron-positron 2D out-plane 2D in-plane Time At late times, the particle spectrum approaches a power-law tail of slope p~2, extending in time to higher and higher energies. The same as in 2D. γmax γ max t The maximum energy grows as γ max t (compare to γ max t 1/2 in shocks). The reconnection rate is v in /c~ 0.02 in 3D (compare to v in /c ~ 0.1 in 2D). (LS & Spitkovsky 14) Time [ω p -1]

70 The highest energy particles γ Density y x (LS & Spitkovsky 14) Two acceleration phases: (1) at the X-point; (2) in between merging islands

71 Plasmoids in reconnection layers Density B0 2L Magnetic energy Kinetic energy Outflow momentum the plasmoids are over-dense, by a factor of a few. they are in rough equipartition of magnetic and kinetic energy. they are moving outwards at ultra-relativistic speeds.

72 Plasmoid space-time tracks We can follow the plasmoids in space and time. First they grow, then they go: First, they grow in the center at non-relativistic speeds. Time [L/c] Island Size [L] Then, they accelerate outwards approaching the speed of light. (LS, Petropoulou & Giannios, in prep)

73 First they grow. Size [L] σ=3 σ=10 σ=50 The plasmoid size grows at a rate of ~0.1 c (dashed lines), with weak dependence on magnetization. Time [L/c] Time [L/c] Time [L/c] (LS, Petropoulou & Giannios, in prep) Outflow Momentum Then they go. σ=3 σ=10 σ=50 Distance/Size Distance/Size Distance/Size (LS, Petropoulou & Giannios, in prep) Size [L] After growing, they accelerate up to the terminal four-velocity ~ σ with a universal acceleration profile.

74 Dependence on the system size Max Size [L] Let us measure the system size in units of the Larmor radius of post-reconnection particles: r 0,hot = mc2 eb 0 Max Larmor radius [L] Larger systems contain larger plasmoids, and the size of the largest plasmoids is always ~ L. Larger systems contain higher energy particles, and the Larmor radius of the highest energy particles is always ~ L (Hillas criterion for reconnection). (LS, Petropoulou & Giannios, in prep)

75 Plasmoids can give bright flares Let us assume a jet with bulk Lorentz factor 10 and half-opening angle 1/10. We vary the plasmoid direction of motion inside the jet (θ =0 along the jet, θ =π/2 perp to the jet). We vary the plasmoid size at the end of the layer (and correspondingly, its final velocity). small and fast large and slow

76 Plasmoids can give fast flares Let us assume a jet with bulk Lorentz factor 10 and half-opening angle 1/10. We vary the plasmoid direction of motion inside the jet (θ =0 along the jet, θ =π/2 perp to the jet). We vary the plasmoid size at the end of the layer (and correspondingly, its final velocity). small and fast large and slow

77 Spectra and light curves small and fast large and slow

PLASMOIDS IN RELATIVISTIC RECONNECTION: THE BLOBS OF BLAZAR EMISSION? Maria Petropoulou Purdue University

PLASMOIDS IN RELATIVISTIC RECONNECTION: THE BLOBS OF BLAZAR EMISSION? Maria Petropoulou Purdue University PLASMOIDS IN RELATIVISTIC RECONNECTION: THE BLOBS OF BLAZAR EMISSION? Maria Petropoulou Purdue University in collaboration with Dimitrios Giannios (Purdue) Lorenzo Sironi(Columbia) October 19, 2016 Einstein

More information

Hybrid Simulations: Numerical Details and Current Applications

Hybrid Simulations: Numerical Details and Current Applications Hybrid Simulations: Numerical Details and Current Applications Dietmar Krauss-Varban and numerous collaborators Space Sciences Laboratory, UC Berkeley, USA Boulder, 07/25/2008 Content 1. Heliospheric/Space

More information

Hybrid Simulation Method ISSS-10 Banff 2011

Hybrid Simulation Method ISSS-10 Banff 2011 Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)

More information

Electron heating in shocks and reconnection

Electron heating in shocks and reconnection Electron heating in shocks and reconnection 11th Plasma Kinetics Working Meeting, Vienna, July 23rd 2018 Lorenzo Sironi (Columbia) with: Xinyi Guo & Michael Rowan (Harvard), Aaron Tran (Columbia) Electron

More information

Relativistic reconnection at the origin of the Crab gamma-ray flares

Relativistic reconnection at the origin of the Crab gamma-ray flares Relativistic reconnection at the origin of the Crab gamma-ray flares Benoît Cerutti Center for Integrated Plasma Studies University of Colorado, Boulder, USA Collaborators: Gregory Werner (CIPS), Dmitri

More information

Crab flares - explosive Reconnection Events in the Nebula

Crab flares - explosive Reconnection Events in the Nebula Crab flares - explosive Reconnection Events in the Nebula Maxim Lyutikov (Purdue) in collaboration with Sergey Komissarov (Leeds) Lorenzo Sironi (Columbia) Oliver Porth (Frankfurt) - ApJ 2017; - JPP, 2017abc

More information

Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays

Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays DESY Summer Student Programme, 2016 Olga Lebiga Taras Shevchenko National University of Kyiv, Ukraine Supervisors Reinaldo

More information

Tsvi Piran The Hebrew University

Tsvi Piran The Hebrew University Some new (old) ideas about particle acceleration and other topics Tsvi Piran The Hebrew University Evgeny Derishev, Daniel Kagan, Ehud Nakar, Glennys Farrar Paris Sept 13-16, 2016 Outline Shock Acceleration

More information

arxiv: v1 [astro-ph.he] 1 Aug 2018

arxiv: v1 [astro-ph.he] 1 Aug 2018 arxiv:188.478v1 [astro-ph.he] 1 Aug 218 Relativistic magnetic reconnection in application to gamma-ray astrophysics Krzysztof Nalewajko 1 1. Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences

More information

Explosive X-point reconnection & Crab flares. Maxim Lyutikov (Purdue U.)

Explosive X-point reconnection & Crab flares. Maxim Lyutikov (Purdue U.) Explosive X-point reconnection & Crab flares Maxim Lyutikov (Purdue U.) Spectra of Crab nebula & flares Tavani et al. 20 Beuhler et al., 2011 ] E 2.F [erg cm 2 s 1 9 11 Break at ~ 0 MeV Fermi CGRO COMPTEL

More information

Understanding the pulsar magnetosphere through first-principle simulations

Understanding the pulsar magnetosphere through first-principle simulations Understanding the pulsar magnetosphere through first-principle simulations Alexander Y. Chen In collaboration with: Andrei Beloborodov Rui Hu The Many Faces of Neutron Stars August 25, 2015 Pulsars: Rotating

More information

Kinetic Plasma Simulations. Anatoly Spitkovsky (Princeton)

Kinetic Plasma Simulations. Anatoly Spitkovsky (Princeton) Kinetic Plasma Simulations Anatoly Spitkovsky (Princeton) Contents Plasma physics on computers How PIC works Electrostatic codes Charge assignment and shape factors Discretization effects Electromagnetic

More information

Space Plasma Physics Thomas Wiegelmann, 2012

Space Plasma Physics Thomas Wiegelmann, 2012 Space Plasma Physics Thomas Wiegelmann, 2012 1. Basic Plasma Physics concepts 2. Overview about solar system plasmas Plasma Models 3. Single particle motion, Test particle model 4. Statistic description

More information

Plasmoids in relativistic reconnection, from birth to adulthood: first they grow, then they go

Plasmoids in relativistic reconnection, from birth to adulthood: first they grow, then they go Advance Access publication 2016 July 6 doi:10.1093/mnras/stw1620 Plasmoids in relativistic reconnection, from birth to adulthood: first they grow, then they go Lorenzo Sironi, 1 Dimitrios Giannios 2 and

More information

Particle acceleration at relativistic shock waves and gamma-ray bursts

Particle acceleration at relativistic shock waves and gamma-ray bursts Particle acceleration at relativistic shock waves and gamma-ray bursts Martin Lemoine Institut d Astrophysique de Paris CNRS, Université Pierre & Marie Curie Outline: 1. Particle acceleration and relativistic

More information

Mikhail V. Medvedev (KU)

Mikhail V. Medvedev (KU) Students (at KU): Sarah Reynolds, Sriharsha Pothapragada Mikhail V. Medvedev (KU) Collaborators: Anatoly Spitkovsky (Princeton) Luis Silva and the Plasma Simulation Group (Portugal) Ken-Ichi Nishikawa

More information

Mechanisms for particle heating in flares

Mechanisms for particle heating in flares Mechanisms for particle heating in flares J. F. Drake University of Maryland J. T. Dahlin University of Maryland M. Swisdak University of Maryland C. Haggerty University of Delaware M. A. Shay University

More information

Simulations of relativistic reconnection in pulsar wind nebulae and pulsar winds

Simulations of relativistic reconnection in pulsar wind nebulae and pulsar winds Simulations of relativistic reconnection in pulsar wind nebulae and pulsar winds Benoît Cerutti Lyman Spitzer Jr. Fellow Princeton University, Dept. of Astrophysical Sciences Collaborators @ Colorado :

More information

arxiv: v2 [astro-ph.he] 18 May 2016

arxiv: v2 [astro-ph.he] 18 May 2016 Mon. Not. R. Astron. Soc. 000, 1?? (2015) Printed 19 May 2016 (MN LATEX style file v2.2) Plasmoids in relativistic reconnection, from birth to adulthood: first they grow, then they go arxiv:1605.02071v2

More information

Magnetic Reconnection: explosions in space and astrophysical plasma. J. F. Drake University of Maryland

Magnetic Reconnection: explosions in space and astrophysical plasma. J. F. Drake University of Maryland Magnetic Reconnection: explosions in space and astrophysical plasma J. F. Drake University of Maryland Magnetic Energy Dissipation in the Universe The conversion of magnetic energy to heat and high speed

More information

Cosmic Ray Acceleration by Magnetic Reconnection

Cosmic Ray Acceleration by Magnetic Reconnection Cosmic Ray Acceleration by Magnetic Reconnection Black Hole sources are Cosmic Ray (CR) accelerators and Very High Energy (VHE) emitters AGNs (blazars, radio-galaxies, seyferts) Black Hole Binaries (Microquasars)

More information

The MRI in a Collisionless Plasma

The MRI in a Collisionless Plasma The MRI in a Collisionless Plasma Eliot Quataert (UC Berkeley) Collaborators: Prateek Sharma, Greg Hammett, Jim Stone Modes of Accretion thin disk: energy radiated away (relevant to star & planet formation,

More information

Particle acceleration during 2D and 3D magnetic reconnection

Particle acceleration during 2D and 3D magnetic reconnection Particle acceleration during 2D and 3D magnetic reconnection J. Dahlin University of Maryland J. F. Drake University of Maryland M. Swisdak University of Maryland Astrophysical reconnection Solar and stellar

More information

Particle Acceleration by Reconnection and VHE emission Around Black Holes and Relativistic Jets

Particle Acceleration by Reconnection and VHE emission Around Black Holes and Relativistic Jets Particle Acceleration by Reconnection and VHE emission Around Black Holes and Relativistic Jets Deciphering the Violent Universe, Playa del Carmen, December 11-15, 2017 Accretion disk coronae Star Formation

More information

Simulation of Relativistic Jet-Plasma Interactions

Simulation of Relativistic Jet-Plasma Interactions Simulation of Relativistic Jet-Plasma Interactions Robert Noble and Johnny Ng Stanford Linear Accelerator Center SABER Workshop, Laboratory Astrophysics WG SLAC, March 15-16, 2006 Motivations High energy

More information

Dissipation Mechanism in 3D Magnetic Reconnection

Dissipation Mechanism in 3D Magnetic Reconnection Dissipation Mechanism in 3D Magnetic Reconnection Keizo Fujimoto Computational Astrophysics Laboratory, RIKEN Reconnection (in the Earth Magnetosphere) Coroniti [1985] 10 km 10 5 km 10 3 km Can induce

More information

2 Relativistic shocks and magnetic fields 17

2 Relativistic shocks and magnetic fields 17 Contents I Introduction 13 1 Overview 15 2 Relativistic shocks and magnetic fields 17 3 Basic Physics 19 3.1 Shock waves.............................. 20 3.2 The theory of relativity........................

More information

The Physics of Collisionless Accretion Flows. Eliot Quataert (UC Berkeley)

The Physics of Collisionless Accretion Flows. Eliot Quataert (UC Berkeley) The Physics of Collisionless Accretion Flows Eliot Quataert (UC Berkeley) Accretion Disks: Physical Picture Simple Consequences of Mass, Momentum, & Energy Conservation Matter Inspirals on Approximately

More information

Explosive reconnection of the double tearing mode in relativistic plasmas

Explosive reconnection of the double tearing mode in relativistic plasmas Explosive reconnection of the double tearing mode in relativistic plasmas Application to the Crab Jérôme Pétri 1 Hubert Baty 1 Makoto Takamoto 2, Seiji Zenitani 3 1 Observatoire astronomique de Strasbourg,

More information

Cosmic Pevatrons in the Galaxy

Cosmic Pevatrons in the Galaxy Cosmic Pevatrons in the Galaxy Jonathan Arons UC Berkeley Cosmic Rays Acceleration in Supernova Remnants Pulsar Wind Nebulae Cosmic rays Cronin, 1999, RMP, 71, S165 J(E) = AE! p, p " 2.7,1GeV < E

More information

8.2.2 Rudiments of the acceleration of particles

8.2.2 Rudiments of the acceleration of particles 430 The solar wind in the Universe intergalactic magnetic fields that these fields should not perturb them. Their arrival directions should thus point back to their sources in the sky, which does not appear

More information

Cosmic Accelerators. 2. Pulsars, Black Holes and Shock Waves. Roger Blandford KIPAC Stanford

Cosmic Accelerators. 2. Pulsars, Black Holes and Shock Waves. Roger Blandford KIPAC Stanford Cosmic Accelerators 2. Pulsars, Black Holes and Shock Waves Roger Blandford KIPAC Stanford Particle Acceleration Unipolar Induction Stochastic Acceleration V ~ Ω Φ I ~ V / Z 0 Z 0 ~100Ω P ~ V I ~ V 2 /Z

More information

Magnetic dissipation in pulsar winds

Magnetic dissipation in pulsar winds Magnetic dissipation in pulsar winds Benoît Cerutti, CNRS & Univ. Grenoble Alpes, France. In collaboration with Sasha Philippov, UC Berkeley, USA. Cerutti & Philippov, A&A (2017) Third Purdue Workshop

More information

The Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley)

The Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley) The Effects of Anisotropic Transport on Dilute Astrophysical Plasmas Eliot Quataert (UC Berkeley) in collaboration with Ian Parrish, Prateek Sharma, Jim Stone, Greg Hammett Hydra A w/ Chandra Galactic

More information

Acceleration of Particles in Gamma-Ray Bursts

Acceleration of Particles in Gamma-Ray Bursts Acceleration of Particles in Gamma-Ray Bursts Bing Zhang Department of Physics and Astronomy University of Nevada, Las Vegas Sep. 29, 2009 In Nonlinear Processes in Astrophysical Plasma: Particle Acceleration,

More information

Strong collisionless shocks are important sources of TeV particles. Evidence for TeV ions is less direct but very strong.

Strong collisionless shocks are important sources of TeV particles. Evidence for TeV ions is less direct but very strong. Collisionless Shocks in 12 minutes or less Don Ellison, North Carolina State Univ. Andrei Bykov, Ioffe Institute, St. Petersburg Don Warren, RIKEN, Tokyo Strong collisionless shocks are important sources

More information

Radia%ve Magne%c Reconnec%on. in Astrophysical Plasmas. Dmitri Uzdensky. (University of Colorado, Boulder) collaborators:

Radia%ve Magne%c Reconnec%on. in Astrophysical Plasmas. Dmitri Uzdensky. (University of Colorado, Boulder) collaborators: Radia%ve Magne%c Reconnec%on collaborators: in Astrophysical Plasmas Dmitri Uzdensky (University of Colorado, Boulder) - B. CeruF *, G. Werner, K. Nalewajko, M. Begelman (Univ. Colorado) - A. Spitkovsky

More information

Magnetically-dominated relativistic jets.

Magnetically-dominated relativistic jets. Magnetically-dominated relativistic jets. Serguei Komissarov University of Leeds UK N.Vlahakis, Y.Granot, A.Konigl, A.Spitkovsky, M.Barkov, J.McKinney, Y.Lyubarsky, M.Lyutikov, N.Bucciantini Plan 1. Astrophysical

More information

Plasma Physics for Astrophysics

Plasma Physics for Astrophysics - ' ' * ' Plasma Physics for Astrophysics RUSSELL M. KULSRUD PRINCETON UNIVERSITY E;RESS '. ' PRINCETON AND OXFORD,, ', V. List of Figures Foreword by John N. Bahcall Preface Chapter 1. Introduction 1

More information

Mikhail V. Medvedev (KU)

Mikhail V. Medvedev (KU) KIAA GRB school Beijing, China, 9 June 2009 Mikhail V. Medvedev (KU) Students: Sarah Reynolds, Sriharsha Pothapragada Collaborators: Silva, Fonseca, Mori & others (Inst Superior Tecniko/GOLP, Portugal

More information

Summer College on Plasma Physics. 30 July - 24 August, The forming of a relativistic partially electromagnetic planar plasma shock

Summer College on Plasma Physics. 30 July - 24 August, The forming of a relativistic partially electromagnetic planar plasma shock 1856-31 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 The forming of a M. E. Dieckmann Institut fuer Theoretische Physik IV, Ruhr-Universitaet, Bochum, Germany The forming of a The forming

More information

Electrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley

Electrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley Electrodynamics of Magnetized Rotators Anatoly Spitkovsky,, UC Berkeley Magnetized rotators are ubiquitous: pulsars, AGN, GRBs (?) Rotation very efficient at long-term energy storage Extraction of rotational

More information

Macroscopic plasma description

Macroscopic plasma description Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion

More information

Plasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University

Plasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative

More information

Magnetic Fields in Blazar Jets

Magnetic Fields in Blazar Jets Magnetic Fields in Blazar Jets Bidzina Z. Kapanadze Ilia State University, Tbilisi, Georgia MFPO 2010- Cracow, May 17-21 Blazars are defined as a AGN class with following features: featureless spectra;

More information

PLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES: (references therein)

PLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES:  (references therein) PLASMA ASTROPHYSICS ElisaBete M. de Gouveia Dal Pino IAG-USP NOTES:http://www.astro.iag.usp.br/~dalpino (references therein) ICTP-SAIFR, October 7-18, 2013 Contents What is plasma? Why plasmas in astrophysics?

More information

Reconnection and the Formation of Magnetic Islands in MHD Models

Reconnection and the Formation of Magnetic Islands in MHD Models Reconnection and the Formation of Magnetic Islands in MHD Models N. F. Loureiro, D. A. Uzdensky, A. A. Schekochihin, R. Samtaney and S. C. Cowley Yosemite 2010 Reconnection Workshop Introduction (I) In

More information

Scaling of Magnetic Reconnection in Collisional and Kinetic Regimes

Scaling of Magnetic Reconnection in Collisional and Kinetic Regimes Scaling of Magnetic Reconnection in Collisional and Kinetic Regimes William Daughton Los Alamos National Laboratory Collaborators: Vadim Roytershteyn, Brian Albright H. Karimabadi, Lin Yin & Kevin Bowers

More information

PULSAR WIND NEBULAE AS COSMIC ACCELERATORS. Elena Amato INAF-Osservatorio Astrofisico di Arcetri

PULSAR WIND NEBULAE AS COSMIC ACCELERATORS. Elena Amato INAF-Osservatorio Astrofisico di Arcetri PULSAR WIND NEBULAE AS COSMIC ACCELERATORS Elena Amato INAF-Osservatorio Astrofisico di Arcetri WHY PWNe ARE INTERESTING! PULSAR PHYSICS: THEY ENCLOSE MOST OF THE PULSAR SPIN- DOWN ENERGY ( L,, L PWN 0.1E

More information

3D hybrid-kinetic turbulence and phase-space cascades

3D hybrid-kinetic turbulence and phase-space cascades 3D hybrid-kinetic turbulence and phase-space cascades ( in a β = 1 plasma ) Silvio Sergio Cerri Department of Astrophysical Sciences, Princeton University, USA 11th Plasma Kinetics Working Meeting WPI

More information

Particle-in-Cell Codes for plasma-based particle acceleration

Particle-in-Cell Codes for plasma-based particle acceleration A. Pukhov Institute for Theoretical Physics I University of Dusseldorf, Germany Particle-in-Cell Codes for plasma-based particle acceleration Outline Relativistic plasmas, acceleration and the simulation

More information

Simulations of External Shocks in. Gamma-Ray Bursts

Simulations of External Shocks in. Gamma-Ray Bursts Simulations of External Shocks in Gamma-Ray Bursts Sarah Wellons Advisor: Anatoly Spitkovsky Submitted in partial fulfillment of the requirements for the degree of Bachelor of Arts Department of Astrophysical

More information

Outline. 1 Why PIC Simulations Make Sense. 2 The Algorithm. 3 Examples Surface High Harmonics Generation. 4 Extensions Of The PIC Algorithm

Outline. 1 Why PIC Simulations Make Sense. 2 The Algorithm. 3 Examples Surface High Harmonics Generation. 4 Extensions Of The PIC Algorithm PIC Simulations an Introduction GRK 1203 Meeting February 12-15, 2008, Oelde Outline 1 Simulations Make Sense 2 3 Surface High Harmonics Generation 4 Of PIC Plasma Physics Is Complex Experiment real thing

More information

Pulsar Wind INAF

Pulsar Wind INAF Pulsar Wind Nebulae @ INAF Niccolo Bucciantini INAF - Osservatorio di Arcetri INFN - Sezione di Firenze UniFi - Dipartimento di Fisica & Astronomia 1 Pulsar Wind Nebulae PWN PWNe are hot bubbles of relativistic

More information

A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1

A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,

More information

Magnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas

Magnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas Magnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas Nick Murphy Harvard-Smithsonian Center for Astrophysics namurphy@cfa.harvard.edu http://www.cfa.harvard.edu/ namurphy/ November 18,

More information

Particle acceleration at relativistic shock waves

Particle acceleration at relativistic shock waves Particle acceleration at relativistic shock waves Martin Lemoine Institut d Astrophysique de Paris CNRS, Université Pierre & Marie Curie Introduction Why relativistic Fermi acceleration? Relativistic outflows

More information

Chapter 1. Introduction to Nonlinear Space Plasma Physics

Chapter 1. Introduction to Nonlinear Space Plasma Physics Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large

More information

Can blazar flares be triggered by the VHE gamma-rays from the surrounding of a supermassive black hole?

Can blazar flares be triggered by the VHE gamma-rays from the surrounding of a supermassive black hole? Can blazar flares be triggered by the VHE gamma-rays from the surrounding of a supermassive black hole? Department of Astrophysics, University of Lodz, Lodz, Poland E-mail: p.banasinski@uni.lodz.pl Wlodek

More information

The Physics of Fluids and Plasmas

The Physics of Fluids and Plasmas The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the

More information

The Weibel Instability in Collisionless Relativistic Shocks

The Weibel Instability in Collisionless Relativistic Shocks The Weibel Instability in Collisionless Relativistic Shocks Tanim Islam University of Virginia The Weibel instability[1], purely electromagnetic in nature, has been used to explain the strong magnetic

More information

On (shock. shock) acceleration. Martin Lemoine. Institut d Astrophysique d. CNRS, Université Pierre & Marie Curie

On (shock. shock) acceleration. Martin Lemoine. Institut d Astrophysique d. CNRS, Université Pierre & Marie Curie On (shock ( shock) acceleration of ultrahigh energy cosmic rays Martin Lemoine Institut d Astrophysique d de Paris CNRS, Université Pierre & Marie Curie 1 Acceleration Hillas criterion log 10 (B/1 G) 15

More information

Solar Flare. A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona)

Solar Flare. A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona) Solar Flares Solar Flare A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona) Flares release 1027-1032 ergs energy in tens of minutes. (Note: one H-bomb: 10

More information

Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction

Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction 16/05/2017, CTU in Prague Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej

More information

Balmer-Dominated Supernova Remnants and the Physics of Collisionless Shocks

Balmer-Dominated Supernova Remnants and the Physics of Collisionless Shocks Balmer-Dominated Supernova Remnants and the Physics of Collisionless Shocks Parviz Ghavamian SNR 0509-67.5 HST ACS Hα (F657N) Supernova Remnants Heat and Enrich the ISM and Accelerate Cosmic Rays reverse-shocked

More information

Waves in plasma. Denis Gialis

Waves in plasma. Denis Gialis Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.

More information

NUMERICAL STUDY OF NON-THERMAL EMISSION FROM LAGRANGIAN PARTICLES IN AGN ENVIRONMENTS.

NUMERICAL STUDY OF NON-THERMAL EMISSION FROM LAGRANGIAN PARTICLES IN AGN ENVIRONMENTS. 13th AGILE Workshop, ASI, Rome May 25, 2015 NUMERICAL STUDY OF NON-THERMAL EMISSION FROM LAGRANGIAN PARTICLES IN AGN ENVIRONMENTS. Dr. Bhargav Vaidya Università degli Studi di Torino, Torino. Collaborators:

More information

CHAPTER 4. Basics of Fluid Dynamics

CHAPTER 4. Basics of Fluid Dynamics CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,

More information

AST 553. Plasma Waves and Instabilities. Course Outline. (Dated: December 4, 2018)

AST 553. Plasma Waves and Instabilities. Course Outline. (Dated: December 4, 2018) AST 553. Plasma Waves and Instabilities Course Outline (Dated: December 4, 2018) I. INTRODUCTION Basic concepts Waves in plasmas as EM field oscillations Maxwell s equations, Gauss s laws as initial conditions

More information

PIC modeling of particle acceleration and high-energy radiation in pulsars

PIC modeling of particle acceleration and high-energy radiation in pulsars PIC modeling of particle acceleration and high-energy radiation in pulsars Benoît Cerutti IPAG, CNRS, Université Grenoble Alpes In collaboration with : Sasha Philippov (Princeton), Anatoly Spitkovsky (Princeton),

More information

ブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学,

ブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学, RKK ( ) BHmag2012,, 2012.2.29 Outline Motivation and basis: Magnetic reconnection around astrophysical black holes Standard equations of resistive GRMHD Test calculations of resistive GRMHD A simulation

More information

Pulsar Winds. John Kirk. Max-Planck-Institut für Kernphysik Heidelberg, Germany. < > p.1/18

Pulsar Winds. John Kirk. Max-Planck-Institut für Kernphysik Heidelberg, Germany. < > p.1/18 Pulsar Winds John Kirk Max-Planck-Institut für Kernphysik Heidelberg, Germany < > p.1/18 About 50 years after... The Crab Nebula Central star is source of particles and magnetic field (Piddington 1957)

More information

Lecture 2 Relativistic Shocks in GRBs 2

Lecture 2 Relativistic Shocks in GRBs 2 Lecture 2 Relativistic Shocks in GRBs 2 Shiho Kobayashi (Liverpool JMU) We have discussed a blast wave. the dynamics: simple: single parameter E /" Blast wave model: applicable to any central engine model

More information

Rela%vis%c Nonthermal Par%cle Accelera%on in Magne%c Reconnec%on

Rela%vis%c Nonthermal Par%cle Accelera%on in Magne%c Reconnec%on Rela%vis%c Nonthermal Par%cle Accelera%on in Magne%c Reconnec%on Dmitri Uzdensky (University of Colorado Boulder) Thanks to Univ. Colorado Theore9cal Plasma Astrophysics Group: Greg Werner Vladimir Zhdankin

More information

Jacob Trier Frederiksen, Niels Bohr Institute, Copenhagen, DK

Jacob Trier Frederiksen, Niels Bohr Institute, Copenhagen, DK Jacob Trier Frederiksen, Niels Bohr Institute, Copenhagen, DK Computational Astrophysics group Niels Bohr Institute Computational power History Present state at the Niels Bohr Institute Prospects Plasma

More information

Kinetic modelling of pulsar magnetospheres

Kinetic modelling of pulsar magnetospheres Kinetic modelling of pulsar magnetospheres Benoît Cerutti IPAG, CNRS, Université Grenoble Alpes In collaboration with : Sasha Philippov (Princeton), Anatoly Spitkovsky (Princeton), Jérémy Mortier (U. Grenoble

More information

Structure of nuclei of extragalactic radio sources and the link with GAIA

Structure of nuclei of extragalactic radio sources and the link with GAIA Structure of nuclei of extragalactic radio sources and the link with GAIA J Roland, IAP & S Lambert, SYRTE I General properties of extragalactic radio sources Radio galaxies : associated with elliptical

More information

Research on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics

Research on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics Chapter 4 Epoch Making Simulation Research on Structure Formation of Plasmas Dominated by Multiple Hierarchical Dynamics Group Representative Yasuaki Kishimoto Naka Fusion Research Establishment, Japan

More information

Turbulence & particle acceleration in galaxy clusters. Gianfranco Brunetti

Turbulence & particle acceleration in galaxy clusters. Gianfranco Brunetti Turbulence & particle acceleration in galaxy clusters Gianfranco Brunetti Mergers & CR-acceleration Cluster-cluster mergers are the most energetic events in the present Universe (10 64 erg/gyr). They can

More information

- Potentials. - Liénard-Wiechart Potentials. - Larmor s Formula. - Dipole Approximation. - Beginning of Cyclotron & Synchrotron

- Potentials. - Liénard-Wiechart Potentials. - Larmor s Formula. - Dipole Approximation. - Beginning of Cyclotron & Synchrotron - Potentials - Liénard-Wiechart Potentials - Larmor s Formula - Dipole Approximation - Beginning of Cyclotron & Synchrotron Maxwell s equations in a vacuum become A basic feature of these eqns is the existence

More information

Instabilities of relativistic jets

Instabilities of relativistic jets Instabilities of relativistic jets G. Bodo INAF Osservatorio Astrofisico di Torino, Italy A. Mignone, P. Rossi, G. Mamatsashvili, S. Massaglia, A. Ferrari show f Universality of relativistic jet phenomenon

More information

Cosmic Rays in CMSO. Ellen Zweibel University of Wisconsin-Madison Santa Fe, 2014

Cosmic Rays in CMSO. Ellen Zweibel University of Wisconsin-Madison Santa Fe, 2014 Cosmic Rays in CMSO Ellen Zweibel University of Wisconsin-Madison Santa Fe, 2014 Galaxies are Pervaded by Magnetic Fields & Relativistic Particles Synchrotron radiation from M51 (MPIfR/NRAO) Galactic molecular

More information

Reduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014

Reduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014 Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton

More information

GAMMA-RAYS FROM MASSIVE BINARIES

GAMMA-RAYS FROM MASSIVE BINARIES GAMMA-RAYS FROM MASSIVE BINARIES W lodek Bednarek Department of Experimental Physics, University of Lódź, Poland 1. Sources of TeV gamma-rays PSR 1259+63/SS2883 - (HESS) LS 5039 - (HESS) LSI 303 +61 o

More information

Stability of strong waves and its implications for pulsar wind shocks

Stability of strong waves and its implications for pulsar wind shocks Stability of strong waves and its implications for pulsar wind shocks Iwona Mochol in collaboration with John Kirk Max-Planck-Institut für Kernphysik Heidelberg, Germany Madrid, 22/05/2013 Pulsar winds

More information

Magnetic fields generated by the Weibel Instability

Magnetic fields generated by the Weibel Instability Magnetic fields generated by the Weibel Instability C. M. Ryu POSTECH, KOREA FFP14 Marseille 14.7.15-7.18 Outline I. Why Weibel instability? II. Simulations III. Conclusion Why Weibel instability? The

More information

Numerical Simulations of the Jet in the Crab Nebula

Numerical Simulations of the Jet in the Crab Nebula Numerical Simulations of the Jet in the Crab Nebula A. Mignone 1, A. Ferrari 1, E. Striani 2, M. Tavani 2 1 Dipartimento di Fisica, Università di Torino 2 Iasf/iaps università di tor vergata (roma) 1.

More information

Particle in cell simulations

Particle in cell simulations Particle in cell simulations Part III: Boundary conditions and parallelization Benoît Cerutti IPAG, CNRS, Université Grenoble Alpes, Grenoble, France. 1 Astrosim, Lyon, June 26 July 7, 2017. Plan of the

More information

Asynchronous 3D HYPERS Simulations of Compact Toroids and Magnetoplasmas

Asynchronous 3D HYPERS Simulations of Compact Toroids and Magnetoplasmas Asynchronous 3D HYPERS Simulations of Compact Toroids and Magnetoplasmas Y.A. Omelchenko Trinum Research, Inc., San Diego, CA, USA US-Japan CT2016 Workshop on Compact Tori, Irvine, CA, August 22-24, 2016

More information

Radiative processes in GRB (prompt) emission. Asaf Pe er (STScI)

Radiative processes in GRB (prompt) emission. Asaf Pe er (STScI) Radiative processes in GRB (prompt) emission Asaf Pe er (STScI) May 2009 Outline Historical approach Synchrotron: pro s and co s Compton scattering in prompt emission (and why it is different than in afterglow)

More information

from Fermi (Higher Energy Astrophysics)

from Fermi (Higher Energy Astrophysics) Particle Acceleration Results from Fermi (Higher Energy Astrophysics) Roger Blandford KIPAC Stanford 3 viii 2011 SLAC SSI 1 Fermi Joint NASA-DOE-Italy- France-Japan- Sweden, Germany mission Launch June

More information

Lecture 4: The particle equations (1)

Lecture 4: The particle equations (1) Lecture 4: The particle equations (1) Presenter: Mark Eric Dieckmann Department of Science and Technology (ITN), Linköping University, Sweden July 17, 2014 Overview We have previously discussed the leapfrog

More information

Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas

Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas Frank Jenko Gyrokinetics an efficient framework for studying turbulence and reconnection in magnetized plasmas Max-Planck-Institut für Plasmaphysik, Garching Workshop on Vlasov-Maxwell Kinetics WPI, Vienna,

More information

Accretion onto the Massive Black Hole in the Galactic Center. Eliot Quataert (UC Berkeley)

Accretion onto the Massive Black Hole in the Galactic Center. Eliot Quataert (UC Berkeley) Accretion onto the Massive Black Hole in the Galactic Center Eliot Quataert (UC Berkeley) Why focus on the Galactic Center? GR! Best evidence for a BH (stellar orbits) M 4x10 6 M Largest BH on the sky

More information

The Mystery of Fast Radio Bursts and its possible resolution. Pawan Kumar

The Mystery of Fast Radio Bursts and its possible resolution. Pawan Kumar The Mystery of Fast Radio Bursts and its possible resolution Outline Pawan Kumar FRBs: summary of relevant observations Radiation mechanism and polarization FRB cosmology Wenbin Lu Niels Bohr Institute,

More information

Heating and current drive: Radio Frequency

Heating and current drive: Radio Frequency Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)

More information

SOLAR WIND ION AND ELECTRON DISTRIBUTION FUNCTIONS AND THE TRANSITION FROM FLUID TO KINETIC BEHAVIOR

SOLAR WIND ION AND ELECTRON DISTRIBUTION FUNCTIONS AND THE TRANSITION FROM FLUID TO KINETIC BEHAVIOR SOLAR WIND ION AND ELECTRON DISTRIBUTION FUNCTIONS AND THE TRANSITION FROM FLUID TO KINETIC BEHAVIOR JUSTIN C. KASPER HARVARD-SMITHSONIAN CENTER FOR ASTROPHYSICS GYPW01, Isaac Newton Institute, July 2010

More information

PIC Simulation of Magnetic Reconnection with Adaptive Mesh Refinement

PIC Simulation of Magnetic Reconnection with Adaptive Mesh Refinement PIC Simulation of Magnetic Reconnection with Adaptive Mesh Refinement Keizo Fujimoto National Astronomical Observatory of Japan Overview Introduction Particle-in-Cell (PIC) model with adaptive mesh refinement

More information

GTC Simulation of Turbulence and Transport in Tokamak Plasmas

GTC Simulation of Turbulence and Transport in Tokamak Plasmas GTC Simulation of Turbulence and Transport in Tokamak Plasmas Z. Lin University it of California, i Irvine, CA 92697, USA and GPS-TTBP Team Supported by SciDAC GPS-TTBP, GSEP & CPES Motivation First-principles

More information

Lectures on basic plasma physics: Introduction

Lectures on basic plasma physics: Introduction Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters

More information