The Physics of Pulsar Magnetospheres June 6th 2016 - NASA GSFC - - PowerPoint PPT Presentation

the physics of pulsar magnetospheres
SMART_READER_LITE
LIVE PREVIEW

The Physics of Pulsar Magnetospheres June 6th 2016 - NASA GSFC - - PowerPoint PPT Presentation

The Physics of Pulsar Magnetospheres June 6th 2016 - NASA GSFC Developing of a specialized Particle-In-Cell code to study Pulsar Magnetospheres Gabriele Brambilla NASA GSFC, Universit degli Studi di Milano Im a grad student and I work


slide-1
SLIDE 1

Developing of a specialized Particle-In-Cell code to study Pulsar Magnetospheres The Physics of Pulsar Magnetospheres

June 6th 2016 - NASA GSFC

Gabriele Brambilla

NASA GSFC, Università degli Studi di Milano

I’m a grad student and I work with: Alice Harding Andrey Timokhin Constantinos Kalapotharakos Demosthenes Kazanas

slide-2
SLIDE 2

In PIC codes, particles moved by the fields form the currents that act on the fields themselves

Birdsall & Langdon 1985 Plasma Physics via Computer Simulation (New York: McGraw-Hill)

3D Cartesian

Electromagnetic Relativistic pusher Parallelized

PML

perfectly matched layer

Domain size

physical 3.0 RLC with PML 3.6 RLC

Resolution (max)

7203 cells 1.5 109 particles All the simulation run on Discover - NASA NCCS Most of the plots and videos are done with VisIt (produced by LLNL)

slide-3
SLIDE 3

We tested our code with typical plasma problems, like the two stream instability

slide-4
SLIDE 4

We tested our code with typical plasma problems, like the two stream instability

slide-5
SLIDE 5

We inject neutral plasma everywhere until we reach a lower magnetization limit inside a fixed radius

slide-6
SLIDE 6

[arbitrary units]

We inject neutral plasma everywhere until we reach a lower magnetization limit inside a fixed radius

slide-7
SLIDE 7

We inject neutral plasma everywhere until we reach a lower magnetization limit inside a fixed radius

slide-8
SLIDE 8

[arbitrary units]

We inject neutral plasma everywhere until we reach a lower magnetization limit inside a fixed radius

slide-9
SLIDE 9

We reached a configuration similar to the force free gradually changing the injection parameters

slide-10
SLIDE 10

In the force free solution the parallel electric field is screened except for the current sheet

slide-11
SLIDE 11

The currents flow on the current sheet and from the polar cap

[arbitrary units]

slide-12
SLIDE 12

It is necessary to keep the magnetization high to resemble the real pulsar behavior

Magnetization

slide-13
SLIDE 13

In PIC simulations the multiplicity is lower than in a real pulsar

Multiplicity

slide-14
SLIDE 14
slide-15
SLIDE 15

We can look at the different contribution of electrons and positrons to the current

Positrons Electrons

[arbitrary units]

slide-16
SLIDE 16

… And we can see where electrons mostly contribute to the current and where positrons do

[arbitrary units]

slide-17
SLIDE 17

We can also look at the direction of the flows

[arbitrary units]

slide-18
SLIDE 18

There are zones in which the flows counter stream and others in which they flow in the same direction

[arbitrary units]

J

Positrons Electrons

slide-19
SLIDE 19

Looking at the energy of the particles we see the most energetic flow in the branches of the current sheet

[arbitrary units] [arbitrary units] [arbitrary units]

charge

slide-20
SLIDE 20

The knowledge of the particle energy distribution is limited by the noise and the magnetic field

maximum noise

  • pen field voltage limit
slide-21
SLIDE 21

In the whole 3D structure the most energetic particles flow out on the current sheet

corotating frame inertial frame

slide-22
SLIDE 22