1 / 51

Extragalactic Magnetic Fields and Dynamo

Extragalactic Magnetic Fields and Dynamo. Extragalactic magnetic fields and Magnetized Universe Project Implications on Dynamo: a) “Kinematic” dynamo in AGN accretion disks b) Magnetic Relaxation (Flux conversion dynamo?): Dynamical magnetic relaxation in force-free plasmas

Download Presentation

Extragalactic Magnetic Fields and Dynamo

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Extragalactic Magnetic Fields and Dynamo • Extragalactic magnetic fields and Magnetized Universe Project • Implications on Dynamo: a) “Kinematic” dynamo in AGN accretion disks b) Magnetic Relaxation (Flux conversion dynamo?): • Dynamical magnetic relaxation in force-free plasmas • Helicity and energy transport and dissipation • Astrophysical Implications Hui Li with K. Bowers, X. Tang, and S. Colgate (Los Alamos National Lab)

  2. Radio Galaxies

  3. Estimating Magnetic Energy • Total magnetic energy: ~ 1060 ergs • Total volume: ~ 1072 cm3 • Typical size: 30 kpc wide, 300 kpc long • Electron energy: 10 GeV – 10 TeV (gmin unknown) • Magnetic fields: 0.5 – 5 m Gauss • Density: ~10-6 cm-3 (thermal), < 10-7 cm-3 (relativi) • Total current: I ~ 5 B R ~ 1018 – 1019 A. • drift velocity: ~10 mm/sec ! to relativistic ~ c • Radio luminosity, spectral index • Estimating the volume, filling factor (~0.1) • Use equipartition/minimum energy assumption • Estimate the total particle and magnetic field energy in the lobes.

  4. Magnetic Energy of Radio Lobes Kronberg, Dufton, Li, Colgate’02 Giants Cluster sources High z sources

  5. Faraday Rotation Measure(Taylor & Perley’93; Colgate & Li’00) Very high FRM, giving mean B fields ~ 30 mG, over size L ~ 50 kpc • implying total magnetic energy 4x1059 ergs, and coherent flux of 8x1041 G cm2. Only supermassive black holes can perhaps provide such energy and flux.

  6. Galaxy Cluster: Perseus A in X-ray 300 Kpc

  7. Perseus A: X-ray + Radio

  8. Ubiquity of Supermassive Black Holes(Kormendy et al. 2001) rSMBH = 5h2 x 105 Msun / Mpc3

  9. Rationale: Energy, Energy, and Energy Black Hole Mass Growth Magnetic Energy Growth

  10. The Magnetized Universe Project • Energy Transport: • (1) how do jets/helix collimate? • (2) how do radio lobes form? • Energy Production: • (1) how to form SMBH? • (2) Accretion disk physics? • Energy Conversion: • Gravitational  Magnetic, dynamo • Astrophys. implications: • (1) will lobes expand? • (2) how do they impact structure formation? • (3) how to prove the existence of B fields? • Energy Dissipation: • (1) how do magnetic • fields dissipate? • (2) how to accelerate • particles?

  11. Implications on Dynamo Magnetic energy and flux of radio lobes and their impact to the general IGM really emphasize the need of understanding: (1) “kinematic” dynamo in accretion disk around SMBHs:How to convert gravitational energy to magnetic energy? * seed field perhaps a non-issue due to large number of rotations * this dynamo seems to saturate at the limit of extracting a significant fraction of the available energy during the SMBH formation Colgate et al.: star-disk collision model for dynamo and liquid sodium experiment at NM-Tech (2) magnetic relaxation (flux-conversion dynamo):How would these lobes evolve in the IGM?---- Similar to Spheromak and RFP? * degree of magnetization of the IGM, impact on galaxy formation? * ultimate fate of the magnetic energy --- extra-galactic cosmic rays? Li et al.: kinetic simulation of collisionless force-free plasmas

  12. AGN Disk Dynamo:Star-Disk Collisions( Colgate et al; Pariev & Colgate’03 ) W-phase: disk rotation – toroidal fields a-phase (helicity injection): rotation of the rising plumes made by star-disk collisions

  13. Liquid Sodium Experiment(Colgate et al.)

  14. Magnetic Lobe Relaxation • Particles are continuously accelerated in-situ, implying continuous energy conversion. • Lobes made in relatively short time (107-108 yrs), in a finite volume, with a finite amount of energy and helicity. Since it is over-pressured compared to its surrounding, it should evolve (by relaxation?). • Kinetic physics should be included in reconnection in lobes: Kinetic scales: c/wpi ~ 1011 cm (n ~10-6 /cc) Sweet-Parker layer width: (Lh/v)1/2 ~ 109 cm (filaments: L ~10 kpc, eta ~ 103, v ~ 108 cm/s)

  15. An idealized Problem Sheet-Pinch: In astrophysical plasmas, the condition is often assumed and it is nearly force-free . Q: Is this sheet-pinch configuration unstable? Q: If so, how does it convert B2 into plasmas?

  16. Why would it evolve? Lz < Lx, Ly --- relaxation Lz is the longest length scale --- no relaxation

  17. Lx Flipping … Lz Lz Lx • Predicting final Bz flux: Bzf = B0 nx (Lz/Lx) • Predicting final magnetic Energy: B2(t=0) = By2 + Bx2 B2 (tf) = By2 + Bz2 DEB = 1 – (Lz/Lx)2

  18. Stability z • For a background B0 = (Bx,By,0), • consider a perturbed Bz, with modes k=(kx, ky), one gets: • gBz = i(k .B0)uz + non-ideal terms • At resonant layer, k .B0 = 0, energy flows into the layer: • for D’ > 0. Dissipation dominates in this layer: • resistivity: Furth et al’63 • collisionless: Drake & Lee’77, Bobrova et al.’01, Li et al.’03 2p MHD kinetic dissipation MHD Lz /2 kinetic dissipation 0 MHD x

  19. Resonant Layers in 3D • In 2D, two layers: az = p/2, 3p/2 • In 3D, large number of modes and layers! Layer-Layer Interaction in 3D is expected to play an important role.

  20. q-Profile of Resonant Surfaces q = By / Bx 2D z 3D

  21. Collisionless Tearing: Linear Growth Rate(Li et al’03, PoP) Short Wavelength Limit Long Wavelength Limit

  22. V4PICA Particle-in-Cell (PIC) Kinetic Code • First principles simulation of the relativistic Maxwell-Vlasov system in three dimensions • Does not need an equation-of-state for closure • Particles are advanced using fields interpolated from a mesh; fields are advanced using sources accumulated from particles • V4PIC was designed from the ground up for ultra-high performance on modern commodity processors

  23. V4PIC -- Under the Hood Lagrangian Eulerian

  24. V4PIC -- Current Status • Sustained 7.1M particles advanced per second per processor (0.14 ms per particle) in the common case was demonstrated on a Pentium 4 2.5GHz. • Memory subsystem is at theoretical limits. • Floating point subsystem is near theoretical limits (~60-80%). • Substantially faster (well over an order of magnitude in some cases) than other PIC codes. • A simple parallelization of V4PIC has been done. • Tens of processors on particle dominated simulations. • Routinely running ~1003 meshes with ~0.5B particles for ~50K time steps on 16 to 32 processors (ranging from overnight to a couple of days per run). • V4PIC has been ported to several x86 clusters and LANL’s Q machine and validated against simple test problems, magnetic reconnection and plasma instabilities simulations.

  25. PIC Simulation Parameters

  26. Total Energy Evolution (Nishimura et al’02,03; Li et al’03a; Li et al’03b) I: Linear Stage; II: Layer Interaction Stage; III: Saturation Stage 2D 3D I II III I II III

  27. Global Evolution (I): Tearing with Island Growth and Transition to Stochastic Field lines (1,0) (0,1) (1,-1) (1,1)

  28. Global Evolution (II-III):Multi-layer Interactions, Transition to Turbulence, Relaxation, and Re-Orientation

  29. Helicity and Energy Dissipation Run 9c4 Run 9c1

  30. 9c4 9c1

  31. Run 9c4 Total H H at k = a H (k < a) • Two Stage: • Total H & W conserved but with significant spectral transfer. • Net H & W dissipation. H (k > a)

  32. Run 9c1 Total H Total W H at k = a W at k = a

  33. Inertial Range ? Dissipation Range 2p/Lx 2p/Lz 2p/di 2p/de

  34. Relaxation with intermittency (Li et al’03) • Evolution of |J| in pdf. • Its mean is decreasing with time, i.e., relaxing. • But with a significant high |J| tail, i.e., localized high |J| filaments where reconnection is occurring. t=0-1.6 15.6 7.8 f(|J|) 14 11 9.35 |J|

  35. CurrentFilaments tWi = 8

  36. Current Filamentation

  37. Field-Aligned E Generalized Ohm’s Law:

  38. Summary • System evolves by constantly choosing more “relaxed states” within the constraints of the geometry. In so doing, converting the excess magnetic energy to particle heating/acceleration. • 3D simulation shows the existence of intermittent regions, with large local shear, current density and magnetic dissipation rates. • Needs more dynamic range to cleanly separate inertial and dissipation ranges, might recover the helicity conservation?

  39. New Experiments Expanding and relaxing magnetic bubble experiments without external driving.

More Related