Kinetic Plasma Simulations in Astrophysics. Lorenzo Sironi
|
|
- Austen Kelly
- 5 years ago
- Views:
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 in collaboration with Dimitrios Giannios (Purdue) Lorenzo Sironi(Columbia) October 19, 2016 Einstein
More informationHybrid 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 informationHybrid 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 informationElectron 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 informationRelativistic 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 informationCrab 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 informationLinear 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 informationTsvi 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 informationarxiv: 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 informationExplosive 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 informationUnderstanding 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 informationKinetic 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 informationSpace 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 informationPlasmoids 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 informationParticle 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 informationMikhail 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 informationMechanisms 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 informationSimulations 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 informationarxiv: 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 informationMagnetic 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 informationCosmic 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 informationThe 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 informationParticle 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 informationParticle 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 informationSimulation 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 informationDissipation 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 information2 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 informationThe 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 informationExplosive 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 informationCosmic 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 information8.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 informationCosmic 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 informationMagnetic 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 informationThe 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 informationAcceleration 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 informationStrong 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 informationRadia%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 informationMagnetically-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 informationPlasma 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 informationMikhail 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 informationSummer 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 informationElectrodynamics 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 informationMacroscopic 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 informationPlasma 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 informationMagnetic 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 informationPLASMA 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 informationReconnection 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 informationScaling 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 informationPULSAR 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 information3D 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 informationParticle-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 informationSimulations 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 informationOutline. 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 informationPulsar 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 informationA 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 informationMagnetic 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 informationParticle 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 informationChapter 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 informationCan 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 informationThe 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 informationThe 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 informationOn (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 informationSolar 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 informationTight-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 informationBalmer-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 informationWaves 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 informationNUMERICAL 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 informationCHAPTER 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 informationAST 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 informationPIC 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,, 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 informationPulsar 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 informationLecture 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 informationRela%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 informationJacob 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 informationKinetic 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 informationStructure 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 informationResearch 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 informationTurbulence & 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 Maxwell s equations in a vacuum become A basic feature of these eqns is the existence
More informationInstabilities 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 informationCosmic 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 informationReduced 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 informationGAMMA-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 informationStability 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 informationMagnetic 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 informationNumerical 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 informationParticle 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 informationAsynchronous 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 informationRadiative 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 informationfrom 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 informationLecture 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 informationGyrokinetics 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 informationAccretion 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 informationThe 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 informationHeating 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 informationSOLAR 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 informationPIC 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 informationGTC 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 informationLectures 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