-
Relaxation to universal non-Maxwellian equilibria in a collisionless plasma
Authors:
Robert J. Ewart,
Michael L. Nastac,
Pablo J. Bilbao,
Thales Silva,
Luís O. Silva,
Alexander A. Schekochihin
Abstract:
Generic equilibria are derived for turbulent relaxing plasmas via an entropy-maximization procedure that accounts for the short-time conservation of certain collisionless invariants. The conservation of these collisionless invariants endows the system with a partial `memory' of its prior conditions, but is imperfect on long time scales due to the development of a turbulent cascade to small scales,…
▽ More
Generic equilibria are derived for turbulent relaxing plasmas via an entropy-maximization procedure that accounts for the short-time conservation of certain collisionless invariants. The conservation of these collisionless invariants endows the system with a partial `memory' of its prior conditions, but is imperfect on long time scales due to the development of a turbulent cascade to small scales, which breaks the precise conservation of phase volume, making this memory imprecise. The equilibria are still determined by the short-time collisionless invariants, but the invariants themselves are driven to a universal form by the nature of the turbulence. This is numerically confirmed for the case of beam instabilities in one-dimensional electrostatic plasmas, where sufficiently strong turbulence appears to cause the distribution function of particle energies to develop a universal power-law tail, with exponent -2.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
Ring momentum distributions as a general feature of Vlasov dynamics in the synchrotron dominated regime
Authors:
Pablo. J. Bilbao,
Robert J. Ewart,
Francisco Assunçao,
Thales Silva,
Luis O. Silva
Abstract:
We study how radiation reaction leads plasmas initially in kinetic equilibrium to develop features in momentum space, such as anisotropies and population inversion, resulting in a ring-shaped momentum distribution that can drive kinetic instabilities. We employ the Landau-Lifshiftz radiation reaction model for a plasma in a strong magnetic field, and we obtain the necessary condition for the devel…
▽ More
We study how radiation reaction leads plasmas initially in kinetic equilibrium to develop features in momentum space, such as anisotropies and population inversion, resulting in a ring-shaped momentum distribution that can drive kinetic instabilities. We employ the Landau-Lifshiftz radiation reaction model for a plasma in a strong magnetic field, and we obtain the necessary condition for the development of population inversion, we show that isotropic Maxwellian and Maxwell-Jüttner plasmas, with thermal temperature $T>m_e c^2/\sqrt{3}$, will develop a ring-like momentum distribution. The timescales and features for forming ring-shaped momentum distributions, the effect of collisions and non-uniform magnetic fields are disscussed, and compared with typical astrophysical and laboratory plasmas parameters. Our results show the pervasiveness of ring-like momentum distribution functions in synchrotron dominated plasma conditions.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
Laboratory study of magnetic reconnection in lunar-relevant mini-magnetospheres
Authors:
Lucas Rovige,
Filipe D. Cruz,
Robert S. Dorst,
Jessica J. Pilgram,
Carmen G. Constantin,
Stephen Vincena,
Fábio Cruz,
Luis O. Silva,
Christoph Niemann,
Derek B. Schaeffer
Abstract:
Mini-magnetospheres are small ion-scale structures that are well-suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres…
▽ More
Mini-magnetospheres are small ion-scale structures that are well-suited to studying kinetic-scale physics of collisionless space plasmas. Such ion-scale magnetospheres can be found on local regions of the Moon, associated with the lunar crustal magnetic field. In this paper, we report on the laboratory experimental study of magnetic reconnection in laser-driven, lunar-like ion-scale magnetospheres on the Large Plasma Device (LAPD) at the University of California - Los Angeles. In the experiment, a high-repetition rate (1 Hz), nanosecond laser is used to drive a fast moving, collisionless plasma that expands into the field generated by a pulsed magnetic dipole embedded into a background plasma and magnetic field. The high-repetition rate enables the acquisition of time-resolved volumetric data of the magnetic and electric fields to characterize magnetic reconnection and calculate the reconnection rate. We notably observe the formation of Hall fields associated with reconnection. Particle-in-cell simulations reproducing the experimental results were performed to study the micro-physics of the interaction. By analyzing the generalized Ohm's law terms, we find that the electron-only reconnection is driven by kinetic effects, through the electron pressure anisotropy. These results are compared to recent satellite measurements that found evidence of magnetic reconnection near the lunar surface.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
OSIRIS-GR: General relativistic activation of the polar cap of a compact neutron star
Authors:
R. Torres,
T. Grismayer,
F. Cruz,
R. A. Fonseca,
L. O. Silva
Abstract:
We present ab initio global general-relativistic Particle-in-cell (GR-PIC) simulations of compact millisecond neutron star magnetospheres in the axisymmetric aligned rotator configuration. We investigate the role of GR and plasma supply on the polar cap particle acceleration efficiency - the precursor of coherent radio emission - employing a new module for the PIC code OSIRIS, designed to model pl…
▽ More
We present ab initio global general-relativistic Particle-in-cell (GR-PIC) simulations of compact millisecond neutron star magnetospheres in the axisymmetric aligned rotator configuration. We investigate the role of GR and plasma supply on the polar cap particle acceleration efficiency - the precursor of coherent radio emission - employing a new module for the PIC code OSIRIS, designed to model plasma dynamics around compact objects with fully self-consistent GR effects. We provide a detailed description of the main sub-algorithms of the novel PIC algorithm, including a charge-conserving current deposit scheme for curvilinear coordinates. We demonstrate efficient particle acceleration in the polar caps of compact neutron stars with denser magnetospheres, numerically validating the spacelike current extension provided by force-free models. We show that GR relaxes the minimum required poloidal magnetospheric current for the transition of the polar cap to the accelerator regime, thus justifying the observation of weak pulsars beyond the expected death line. We denote that spin-down luminosity intermittency and radio pulse nullings for older pulsars might arise from the interplay between the polar and outer gaps. Also, narrower radio beams are expected for weaker low-obliquity pulsars.
△ Less
Submitted 1 June, 2024; v1 submitted 5 January, 2024;
originally announced January 2024.
-
Laboratory realization of relativistic pair-plasma beams
Authors:
C. D. Arrowsmith,
P. Simon,
P. Bilbao,
A. F. A. Bott,
S. Burger,
H. Chen,
F. D. Cruz,
T. Davenne,
I. Efthymiopoulos,
D. H. Froula,
A. M. Goillot,
J. T. Gudmundsson,
D. Haberberger,
J. Halliday,
T. Hodge,
B. T. Huffman,
S. Iaquinta,
F. Miniati,
B. Reville,
S. Sarkar,
A. A. Schekochihin,
L. O. Silva,
R. Simpson,
V. Stergiou,
R. M. G. M. Trines
, et al. (4 additional authors not shown)
Abstract:
Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dy…
▽ More
Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dynamics of such compact objects is believed to be fundamental. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies which are rather limited. We present first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN's Super Proton Synchrotron (SPS) accelerator. The produced pair beams have a volume that fills multiple Debye spheres and are thus able to sustain collective plasma oscillations. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations.
△ Less
Submitted 8 December, 2023;
originally announced December 2023.
-
Particle-in-cell simulations of pulsar magnetospheres: transition between electrosphere and force-free regimes
Authors:
Fábio Cruz,
Thomas Grismayer,
Alexander Y. Chen,
Anatoly Spitkovsky,
Ricardo A. Fonseca,
Luis O. Silva
Abstract:
Global particle-in-cell (PIC) simulations of pulsar magnetospheres are performed with a volume, surface and pair production-based plasma injection schemes to systematically investigate the transition between electrosphere and force-free pulsar magnetospheric regimes. A new extension of the PIC code OSIRIS to model pulsar magnetospheres using a two-dimensional axisymmetric spherical grid is present…
▽ More
Global particle-in-cell (PIC) simulations of pulsar magnetospheres are performed with a volume, surface and pair production-based plasma injection schemes to systematically investigate the transition between electrosphere and force-free pulsar magnetospheric regimes. A new extension of the PIC code OSIRIS to model pulsar magnetospheres using a two-dimensional axisymmetric spherical grid is presented. The sub-algorithms of the code and thorough benchmarks are presented in detail, including a new first-order current deposition scheme that conserves charge to machine precision. It is shown that all plasma injection schemes produce a range of magnetospheric regimes. Active solutions can be obtained with surface and volume injection schemes when using artificially large plasma injection rates, and with pair production-based plasma injection for sufficiently large separation between kinematic and pair production energy scales.
△ Less
Submitted 9 September, 2023;
originally announced September 2023.
-
High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations
Authors:
K. M. Schoeffler,
T. Grismayer,
D. Uzdensky,
L. O. Silva
Abstract:
The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiat…
▽ More
The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field $B$ and plasma density~$n$, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the $n\mbox{-}B$ space are developed and used to visualize and quantify the effects of compression. The $n\mbox{-}B$ histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic (QED) effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.
△ Less
Submitted 29 March, 2023;
originally announced March 2023.
-
Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star
Authors:
R. Torres,
T. Grismayer,
F. Cruz,
L. O. Silva
Abstract:
Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell's equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly-rotating magnetised neu…
▽ More
Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell's equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly-rotating magnetised neutron stars taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. $R_s \sim 0.5 [R_*]$, we show that the new terms lead to a percent order correction in the magnetic field orientation and strength compared to the case with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.
△ Less
Submitted 7 January, 2024; v1 submitted 20 February, 2023;
originally announced February 2023.
-
Radiation reaction cooling as a source of anisotropic momentum distributions with inverted populations
Authors:
Pablo J. Bilbao,
Luis O. Silva
Abstract:
Under the presence of strong electromagnetic fields and radiation reaction, plasmas develop anisotropic momentum distributions, characterized by a population inversion. This is a general property of collisionless plasmas when the radiation reaction force is taken into account. We study the case of a plasma in a strong magnetic field and demonstrate the development of ring momentum distributions. T…
▽ More
Under the presence of strong electromagnetic fields and radiation reaction, plasmas develop anisotropic momentum distributions, characterized by a population inversion. This is a general property of collisionless plasmas when the radiation reaction force is taken into account. We study the case of a plasma in a strong magnetic field and demonstrate the development of ring momentum distributions. The timescales for ring formation are derived for this configuration. The analytical results for the ring properties and the timescales for ring formation are confirmed with particle-in-cell simulations. The resulting momentum distributions are kinetically unstable and are known to lead to coherent radiation emission in astrophysical plasmas and laboratory setups.
△ Less
Submitted 29 March, 2023; v1 submitted 23 December, 2022;
originally announced December 2022.
-
Model of pulsar pair cascades in non uniform electric fields: growth rate, density profile and screening time
Authors:
Fábio Cruz,
Thomas Grismayer,
Simon Iteanu,
Paolo Tortone,
Luis O. Silva
Abstract:
Time-dependent cascades of electron-positron pairs are thought to be the main source of plasma in pulsar magnetospheres and a primary ingredient to explain the nature of pulsar radio emission, a longstanding open problem in high-energy astrophysics. During these cascades - positive feedback loops of gamma-ray photon emission, via curvature radiation by TeV electrons and positrons, and pair product…
▽ More
Time-dependent cascades of electron-positron pairs are thought to be the main source of plasma in pulsar magnetospheres and a primary ingredient to explain the nature of pulsar radio emission, a longstanding open problem in high-energy astrophysics. During these cascades - positive feedback loops of gamma-ray photon emission, via curvature radiation by TeV electrons and positrons, and pair production -, the plasma self-consistently develops inductive waves that couple to electromagnetic modes capable of escaping the pulsar dense plasma. In this work, we present an analytical description of pair cascades relevant in pulsars, including their onset, exponential growth and saturation stages. We study this problem in the case of a background linear electric field, relevant in pulsar polar caps, and using an heuristic model of the pair production process. The analytical results are confirmed with particle-in-cell simulations performed with OSIRIS including heuristic pair production.
△ Less
Submitted 7 April, 2022;
originally announced April 2022.
-
Coherent emission from QED cascades in pulsar polar caps
Authors:
Fábio Cruz,
Thomas Grismayer,
Alexander Y. Chen,
Anatoly Spitkovsky,
Luis O. Silva
Abstract:
Pulsar magnetospheres are thought to be filled with electron-positron plasma generated in pair cascades. The driving mechanism of these cascades is the emission of gamma-ray photons and their conversion into pairs via Quantum Electrodynamics (QED) processes. In this work, we present 2D particle-in-cell simulations of pair cascades in pulsar polar caps with realistic magnetic field geometry that in…
▽ More
Pulsar magnetospheres are thought to be filled with electron-positron plasma generated in pair cascades. The driving mechanism of these cascades is the emission of gamma-ray photons and their conversion into pairs via Quantum Electrodynamics (QED) processes. In this work, we present 2D particle-in-cell simulations of pair cascades in pulsar polar caps with realistic magnetic field geometry that include the relevant QED processes from first principles. Our results show that, due to variation of magnetic field curvature across the polar cap, pair production bursts self-consistently develop an inclination with respect to the local magnetic field that favors the generation of coherent electromagnetic modes with properties consistent with pulsar radio emission. We show that this emission is peaked along the magnetic axis and close to the polar cap edge and may thus offer an explanation for the core and conal components of pulsar radio emission.
△ Less
Submitted 26 August, 2021;
originally announced August 2021.
-
Neutrino-electron magnetohydrodynamics in an expanding Universe
Authors:
L. M. Perrone,
G. Gregori,
B. Reville,
L. O. Silva,
R. Bingham
Abstract:
We derive a new model for neutrino-plasma interactions in an expanding universe that incorporates the collective effects of the neutrinos on the plasma constituents. We start from the kinetic description of a multi-species plasma in the flat Friedmann-Robertson-Walker metric, where the particles are coupled to neutrinos through the charged- and neutral-current forms of the weak interaction. We the…
▽ More
We derive a new model for neutrino-plasma interactions in an expanding universe that incorporates the collective effects of the neutrinos on the plasma constituents. We start from the kinetic description of a multi-species plasma in the flat Friedmann-Robertson-Walker metric, where the particles are coupled to neutrinos through the charged- and neutral-current forms of the weak interaction. We then derive the fluid equations and specialize our model to (a) the lepton epoch, where we consider a pair electron-positron plasma interacting with electron (anti-)neutrinos, and (b) after the electron-positron annihilation, where we model an electron-proton plasma and take the limit of slow ions and inertia-less electrons to obtain a set of neutrino-electron magnetohydrodynamics (NEMHD) equations. In both models, the dynamics of the plasma is affected by the neutrino motion through a ponderomotive force and, as a result, new terms appear in the induction equation that can act as a source for magnetic field generation in the early universe. A brief discussion on the possible applications of our model is proposed.
△ Less
Submitted 14 October, 2021; v1 submitted 28 June, 2021;
originally announced June 2021.
-
Kinetic model of large-amplitude oscillations in neutron star pair cascades
Authors:
Fábio Cruz,
Thomas Grismayer,
Luis O. Silva
Abstract:
Electron-positron pair cascades developed in the extreme electromagnetic fields of neutron star polar caps are considered a key source of magnetospheric plasma in these objects. We use a simplified model that maps the Quantum Electrodynamics processes governing the pair cascades to analytically and numerically model the development of the pair cascade, and show that large-amplitude oscillations of…
▽ More
Electron-positron pair cascades developed in the extreme electromagnetic fields of neutron star polar caps are considered a key source of magnetospheric plasma in these objects. We use a simplified model that maps the Quantum Electrodynamics processes governing the pair cascades to analytically and numerically model the development of the pair cascade, and show that large-amplitude oscillations of the electric field are inductively driven by the resulting plasma. A plasma instability arises in these oscillations, and particles accelerated in growing electric field perturbations can drive secondary pair bursts that damp the large-amplitude oscillations. An analytical model is proposed to describe this interplay between the pair production and kinetic collective plasma processes. All analytical results are shown to be in excellent agreement with particle-in-cell simulations.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
Generating ultra-dense pair beams using 400 GeV/c protons
Authors:
C. D. Arrowsmith,
N. Shukla,
N. Charitonidis,
R. Boni,
H. Chen,
T. Davenne,
D. H. Froula,
B. T. Huffman,
Y. Kadi,
B. Reville,
S. Richardson,
S. Sarkar,
J. L. Shaw,
L. O. Silva,
R. M. G. M. Trines,
R. Bingham,
G. Gregori
Abstract:
A previously unexplored experimental scheme is presented for generating low-divergence, ultra-dense, relativistic, electron-positron beams using 400 GeV/c protons available at facilities such as HiRadMat and AWAKE at CERN. Preliminary Monte-Carlo and Particle-in-cell simulations demonstrate the possibility of generating beams containing $10^{13}-10^{14}$ electron-positron pairs at sufficiently hig…
▽ More
A previously unexplored experimental scheme is presented for generating low-divergence, ultra-dense, relativistic, electron-positron beams using 400 GeV/c protons available at facilities such as HiRadMat and AWAKE at CERN. Preliminary Monte-Carlo and Particle-in-cell simulations demonstrate the possibility of generating beams containing $10^{13}-10^{14}$ electron-positron pairs at sufficiently high densities to drive collisionless beam-plasma instabilities, which are expected to play an important role in magnetic field generation and the related radiation signatures of relativistic astrophysical phenomena. The pair beams are quasi-neutral, with size exceeding several skin-depths in all dimensions, allowing for the first time the examination of the effect of competition between transverse and longitudinal instability modes on the growth of magnetic fields. Furthermore, the presented scheme allows for the possibility of controlling the relative density of hadrons to electron-positron pairs in the beam, making it possible to explore the parameter spaces for different astrophysical environments.
△ Less
Submitted 9 November, 2020;
originally announced November 2020.
-
Magnetized Current Filaments as a Source of Circularly Polarized Light
Authors:
U. Sinha,
J. Martins,
J. Vieira,
K. M. Schoeffler,
R. A. Fonseca,
L. O. Silva
Abstract:
We show that the Weibel or currente filamentation instability can lead to the emission of circularly polarized radiation. Using particle-in-cell (PIC) simulations and a radiation post-processing numerical algorithm, we demonstrate that the level of circular polarization increases with the initial plasma magnetization, saturating at ~13% when the magnetization, given by the ratio of magnetic energy…
▽ More
We show that the Weibel or currente filamentation instability can lead to the emission of circularly polarized radiation. Using particle-in-cell (PIC) simulations and a radiation post-processing numerical algorithm, we demonstrate that the level of circular polarization increases with the initial plasma magnetization, saturating at ~13% when the magnetization, given by the ratio of magnetic energy density to the electron kinetic energy density, is larger than 0.05. Furthermore, we show that this effect requires an ion-electron mass ratio greater than unity. These findings, which could also be tested in currently available laboratory conditions, show that the recent observation of circular polarization in gamma ray burst afterglows could be attributed to the presence of magnetized current filaments driven by the Weibel or the current filamentation instability.
△ Less
Submitted 17 December, 2018;
originally announced December 2018.
-
Bright gamma-ray flares powered by magnetic reconnection in QED-strength magnetic fields
Authors:
K. M. Schoeffler,
T. Grismayer,
D. Uzdensky,
R. A. Fonseca,
L. O. Silva
Abstract:
Strong magnetic fields in magnetospheres of neutron stars (especially magnetars) and other astrophysical objects may release their energy in violent, intense episodes of magnetic reconnection. While reconnection has been studied extensively, the extreme field strength near neutron stars introduces new effects: synchrotron cooling and electron-positron pair production. Using massively parallel part…
▽ More
Strong magnetic fields in magnetospheres of neutron stars (especially magnetars) and other astrophysical objects may release their energy in violent, intense episodes of magnetic reconnection. While reconnection has been studied extensively, the extreme field strength near neutron stars introduces new effects: synchrotron cooling and electron-positron pair production. Using massively parallel particle-in-cell simulations that self-consistently incorporate these new quantum-electrodynamic effects, we investigate relativistic magnetic reconnection in the strong-field regime. We show that reconnection in this regime can efficiently convert magnetic energy to X-ray and gamma-ray radiation and thus power bright high-energy astrophysical flares. Rapid radiative cooling causes strong plasma and magnetic field compression in compact plasmoids. In the most extreme cases, the field can approach the critical quantum limit, leading to copious pair production.
△ Less
Submitted 29 October, 2018; v1 submitted 25 July, 2018;
originally announced July 2018.
-
Fully kinetic large scale simulations of the collisionless Magnetorotational instability
Authors:
Giannandrea Inchingolo,
Thomas Grismayer,
Nuno F. Loureiro,
Ricardo A. Fonseca,
Luis O. Silva
Abstract:
We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $β$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the…
▽ More
We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $β$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the underlying mechanisms for the development of this turbulent state is the drift-kink instability (DKI) of the current sheets resulting from the nonlinear evolution of the channel modes. The onset of the DKI can only be observed for simulation domain sizes exceeding several linear MRI wavelengths. The DKI, together with ensuing magnetic reconnection, activate the turbulent motion of the plasma in the late stage of the nonlinear evolution of the MRI. At steady state, the magnetic energy has an MHD-like spectrum with a slope of $k^{-5/3}$ for $kρ<1$ and $k^{-3}$ for sub-Larmor scale ($kρ>1$). We also examine the role of the collisionless MRI and associated magnetic reconnection in the development of pressure anisotropy. We study the stability of the system due to this pressure anisotropy, observing the development of mirror instability during the early-stage of the MRI. We further discuss the importance of magnetic reconnection for particle acceleration during the turbulence regime. In particular, consistent with reconnection studies, we show that at late times the kinetic energy presents a characteristic slope of $ε^{-2}$ in the high-energy region.
△ Less
Submitted 29 April, 2018; v1 submitted 25 January, 2018;
originally announced January 2018.
-
Physics of collisionless shocks - theory and simulation
Authors:
A. Stockem Novo,
A. Bret,
R. A. Fonseca,
L. O. Silva
Abstract:
Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical envir…
▽ More
Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations.
△ Less
Submitted 16 September, 2015;
originally announced September 2015.
-
Slow down of a globally neutral relativistic $e^-e^+$ beam shearing the vacuum
Authors:
E. P. Alves,
T. Grismayer,
M. G. Silveirinha,
R. A. Fonseca,
L. O. Silva
Abstract:
The microphysics of relativistic collisionless sheared flows is investigated in a configuration consisting of a globally neutral, relativistic $e^-e^+$ beam streaming through a hollow plasma/dielectric channel. We show through multidimensional PIC simulations that this scenario excites the Mushroom instability (MI), a transverse shear instability on the electron-scale, when there is no overlap (no…
▽ More
The microphysics of relativistic collisionless sheared flows is investigated in a configuration consisting of a globally neutral, relativistic $e^-e^+$ beam streaming through a hollow plasma/dielectric channel. We show through multidimensional PIC simulations that this scenario excites the Mushroom instability (MI), a transverse shear instability on the electron-scale, when there is no overlap (no contact) between the $e^-e^+$ beam and the walls of the hollow plasma channel. The onset of the MI leads to the conversion of the beam's kinetic energy into magnetic (and electric) field energy, effectively slowing down a globally neutral body in the absence of contact. The collisionless shear physics explored in this configuration may operate in astrophysical environments, particularly in highly relativistic and supersonic settings where macroscopic shear processes are stable.
△ Less
Submitted 29 June, 2015;
originally announced June 2015.
-
Formation of Lunar Swirls
Authors:
R. A. Bamford,
E. P. Alves,
F. Cruz,
B. J. Kellett,
R. A. Fonseca,
L. O Silva,
R. M. G. M. Trines,
J. S. Halekas,
G. Kramer,
E. Harnett,
R. A. Cairns,
R. Bingham
Abstract:
In this paper we show a plausible mechanism that could lead to the formation of the Dark Lanes in Lunar Swirls, and the electromagnetic shielding of the lunar surface that results in the preservation of the white colour of the lunar regolith.
We present the results of a fully self-consistent 2 and 3 dimensional particle-in-cell simulations of mini-magnetospheres that form above the lunar surface…
▽ More
In this paper we show a plausible mechanism that could lead to the formation of the Dark Lanes in Lunar Swirls, and the electromagnetic shielding of the lunar surface that results in the preservation of the white colour of the lunar regolith.
We present the results of a fully self-consistent 2 and 3 dimensional particle-in-cell simulations of mini-magnetospheres that form above the lunar surface and show that they are consistent with the formation of `lunar swirls' such as the archetypal formation Reiner Gamma. The simulations show how the microphysics of the deflection/shielding of plasma operates from a kinetic-scale cavity, and show that this interaction leads to a footprint with sharp features that could be the mechanism behind the generation of `dark lanes'. The physics of mini-magnetospheres is described and shown to be controlled by space-charge fields arising due to the magnetized electrons and unmagnetized ions. A comparison between model and observation is shown for a number of key plasma parameters.
△ Less
Submitted 23 May, 2015;
originally announced May 2015.
-
Transverse electron-scale instability in relativistic shear flows
Authors:
E. P. Alves,
T. Grismayer,
R. A. Fonseca,
L. O. Silva
Abstract:
Electron-scale surface waves are shown to be unstable in the transverse plane of a shear flow in an initially unmagnetized plasma, unlike in the (magneto)hydrodynamics case. It is found that these unstable modes have a higher growth rate than the closely related electron-scale Kelvin-Helmholtz instability in relativistic shears. Multidimensional particle-in-cell simulations verify the analytic res…
▽ More
Electron-scale surface waves are shown to be unstable in the transverse plane of a shear flow in an initially unmagnetized plasma, unlike in the (magneto)hydrodynamics case. It is found that these unstable modes have a higher growth rate than the closely related electron-scale Kelvin-Helmholtz instability in relativistic shears. Multidimensional particle-in-cell simulations verify the analytic results and further reveal the emergence of mushroom-like electron density structures in the nonlinear phase of the instability, similar to those observed in the Rayleigh Taylor instability despite the great disparity in scales and different underlying physics. Macroscopic ($\gg c/ω_{pe}$) fields are shown to be generated by these microscopic shear instabilities, which are relevant for particle acceleration, radiation emission and to seed MHD processes at long time-scales.
△ Less
Submitted 22 May, 2015;
originally announced May 2015.
-
Shock formation in electron-ion plasmas: mechanism and timing
Authors:
Anne Stockem Novo,
Antoine Bret,
Ricardo A. Fonseca,
Luis O. Silva
Abstract:
We analyse the full shock formation process in electron-ion plasmas in theory and simulations. It is accepted that electromagnetic shocks in initially unmagnetised relativistic plasmas are triggered by the filamentation instability. However, the transition from the first unstable phase to the quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, taking in…
▽ More
We analyse the full shock formation process in electron-ion plasmas in theory and simulations. It is accepted that electromagnetic shocks in initially unmagnetised relativistic plasmas are triggered by the filamentation instability. However, the transition from the first unstable phase to the quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, taking into account the filament merging in the non-linear phase of the filamentation instability. This process is much slower than in electron-positron pair shocks, so that the shock formation is longer by a factor proportional to sqrt(m_i/m_e) ln(m_i/m_e).
△ Less
Submitted 2 April, 2015;
originally announced April 2015.
-
Electromagnetic field generation in the downstream of electrostatic shocks due to electron trapping
Authors:
A. Stockem,
T. Grismayer,
R. A. Fonseca,
L. O. Silva
Abstract:
A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $ \tilde \, 10^4 \, ω_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instabil…
▽ More
A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $ \tilde \, 10^4 \, ω_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock.
△ Less
Submitted 20 August, 2014;
originally announced August 2014.
-
SEP acceleration in CME driven shocks using a hybrid code
Authors:
L. Gargate,
R. A. Fonseca,
L. O. Silva,
R. A. Bamford,
R. Bingham
Abstract:
We preform hybrid simulations of super Alfvenic quasi-parallel shock, driven by a Coronal Mass Ejection (CME), propagating in the Outer Coronal or Solar Wind at distances of between 3 to 6 solar radii. The hybrid treatment of the problem enable the study of the shock propagation on the ion time scale, preserving ion kinetics and allowing for a self consistent treatment of the shock propagation and…
▽ More
We preform hybrid simulations of super Alfvenic quasi-parallel shock, driven by a Coronal Mass Ejection (CME), propagating in the Outer Coronal or Solar Wind at distances of between 3 to 6 solar radii. The hybrid treatment of the problem enable the study of the shock propagation on the ion time scale, preserving ion kinetics and allowing for a self consistent treatment of the shock propagation and particle acceleration. The CME plasma drags the embedded magnetic field lines stretching from the sun, and propagates out into interplanetary space at a greater velocity than the in-situ solar wind, driving the shock, and producing very energetic particles. Our results show electromagnetic Alfven waves are generated at the shock front. The waves propagate upstream of the shock and are produced by the counter streaming ions of the solar wind plasma being reflected at the shock. A significant fraction of the particles are accelerated in two distinct phases: first, particles drift from the shock and are accelerated in the upstream region and, second, particles arriving at the shock get trapped, and are accelerated at the shock front. A fraction of the particles diffused back to the shock, which is consistent with the Fermi acceleration mechanism.
△ Less
Submitted 20 June, 2014;
originally announced June 2014.
-
Collisionless Weibel shocks: full formation mechanism and timing
Authors:
Antoine Bret,
Anne Stockem,
Ramesh Narayan,
Luis O. Silva
Abstract:
Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation…
▽ More
Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.
△ Less
Submitted 16 June, 2014;
originally announced June 2014.
-
Electron-scale shear instabilities: magnetic field generation and particle acceleration in astrophysical jets
Authors:
E. P. Alves,
T. Grismayer,
R. A. Fonseca,
L. O. Silva
Abstract:
Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields make shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the c…
▽ More
Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields make shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the collisionless, unmagnetized Kelvin-Helmholtz instability (KHI) and a large-scale dc magnetic field generation mechanism on the electron scales. We show that these processes are important candidates to generate magnetic fields in the presence of strong velocity shears, which may naturally originate in energetic matter outburst of active galactic nuclei and gamma-ray bursters. We show that the KHI is robust to density jumps between shearing flows, thus operating in various scenarios with different density contrasts. Multidimensional particle-in-cell (PIC) simulations of the KHI, performed with OSIRIS, reveal the emergence of a strong and large-scale dc magnetic field component, which is not captured by the standard linear fluid theory. This dc component arises from kinetic effects associated with the thermal expansion of electrons of one flow into the other across the shear layer, whilst ions remain unperturbed due to their inertia. The electron expansion forms dc current sheets, which induce a dc magnetic field. Our results indicate that most of the electromagnetic energy developed in the KHI is stored in the dc component, reaching values of equipartition on the order of $10^{-3}$ in the electron time-scale, and persists longer than the proton time-scale. Particle scattering/acceleration in the self generated fields of these shear flow instabilities is also analyzed.
△ Less
Submitted 2 April, 2014;
originally announced April 2014.
-
Interaction of Ultra Relativistic e- e+ Fireball Beam with Plasma
Authors:
P. Muggli,
S. F. Martins,
N. Shukla,
J. Vieira,
L. O. Silva
Abstract:
Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electron-positron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma wi…
▽ More
Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electron-positron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma with a density 10^17/cc. The different regimes of the instability, from the purely transverse to the mixed mode filamentation, can be accessed by varying the background plasma density. The instability generates large local plasma gradients, intense transverse magnetic fields, and enhanced emission of radiation. We suggest that these effects may be observed experimentally for the first time.
△ Less
Submitted 11 September, 2018; v1 submitted 18 June, 2013;
originally announced June 2013.
-
Collisionless shock formation, spontaneous electromagnetic fluctuations and streaming instabilities
Authors:
A. Bret,
A. Stockem,
F. Fiuza,
C. Ruyer,
L. Gremillet,
R. Narayan,
L. O. Silva
Abstract:
Collisionless shocks are ubiquitous in astrophysics and in the lab. Recent numerical simulations and experiments have shown how they can arise from the encounter of two collisionless plasma shells. When the shells interpenetrate, the overlapping region turns unstable, triggering the shock formation. As a first step towards a microscopic understanding of the process, we analyze here in detail the i…
▽ More
Collisionless shocks are ubiquitous in astrophysics and in the lab. Recent numerical simulations and experiments have shown how they can arise from the encounter of two collisionless plasma shells. When the shells interpenetrate, the overlapping region turns unstable, triggering the shock formation. As a first step towards a microscopic understanding of the process, we analyze here in detail the initial instability phase. On the one hand, 2D relativistic PIC simulations are performed where two symmetric initially cold pair plasmas collide. On the other hand, the instabilities at work are analyzed, as well as the field at saturation and the seed field which gets amplified. For mildly relativistic motions and onward, Weibel modes govern the linear phase. We derive an expression for the duration of the linear phase in good agreement with the simulations. This saturation time constitutes indeed a lower-bound for the shock formation time.
△ Less
Submitted 9 April, 2013; v1 submitted 17 March, 2013;
originally announced March 2013.
-
The impact of kinetic effects on the properties of relativistic electron-positron shocks
Authors:
A. Stockem,
F. Fiuza,
R. A. Fonseca,
L. O. Silva
Abstract:
We assess the impact of non-thermally shock-accelerated particles on the magnetohydrodynamic (MHD) jump conditions of relativistic shocks. The adiabatic constant is calculated directly from first principle particle-in-cell simulation data, enabling a semi-kinetic approach to improve the standard fluid model and allowing for an identification of the key parameters that define the shock structure. W…
▽ More
We assess the impact of non-thermally shock-accelerated particles on the magnetohydrodynamic (MHD) jump conditions of relativistic shocks. The adiabatic constant is calculated directly from first principle particle-in-cell simulation data, enabling a semi-kinetic approach to improve the standard fluid model and allowing for an identification of the key parameters that define the shock structure. We find that the evolving upstream parameters have a stronger impact than the corrections due to non-thermal particles. We find that the decrease of the upstream bulk speed yields deviations from the standard MHD model up to 10%. Furthermore, we obtain a quantitative definition of the shock transition region from our analysis. For Weibel-mediated shocks the inclusion of a magnetic field in the MHD conservation equations is addressed for the first time.
△ Less
Submitted 9 October, 2012;
originally announced October 2012.
-
Acceleration in perpendicular relativistic shocks for plasmas consisting of leptons and hadrons
Authors:
A. Stockem,
F. Fiuza,
R. A. Fonseca,
L. O. Silva
Abstract:
We investigate the acceleration of light particles in perpendicular shocks for plasmas consisting of a mixture of leptonic and hadronic particles. Starting from the full set of conservation equations for the mixed plasma constituents, we generalize the magneto-hydrodynamical jump conditions for a multi-component plasma, including information about the specific adiabatic constants for the different…
▽ More
We investigate the acceleration of light particles in perpendicular shocks for plasmas consisting of a mixture of leptonic and hadronic particles. Starting from the full set of conservation equations for the mixed plasma constituents, we generalize the magneto-hydrodynamical jump conditions for a multi-component plasma, including information about the specific adiabatic constants for the different species. The impact of deviations from the standard model of an ideal gas is compared in theory and particle-in-cell simulations, showing that the standard-MHD model is a good approximation. The simulations of shocks in electron-positron-ion plasmas are for the first time multi-dimensional, transverse effects are small in this configuration and 1D simulations are a good representation if the initial magnetization is chosen high. 1D runs with a mass ratio of 1836 are performed, which identify the Larmor frequency ω_{ci} as the dominant frequency that determines the shock physics in mixed component plasmas. The maximum energy in the non-thermal tail of the particle spectra evolves in time according to a power-law proportional to t^αwith αin the range 1/3 < α< 1, depending on the initial parameters. A connection is made with transport theoretical models by Drury (1983) and Gargate & Spitkovsky (2011), which predict an acceleration time proportional to γand the theory for small wavelength scattering by Kirk & Reville (2010), which predicts a behavior rather as proportional to γ^2. Furthermore, we compare different magnetic field orientations with B_0 inside and out of the plane, observing qualitatively different particle spectra than in pure electron-ion shocks.
△ Less
Submitted 12 June, 2012;
originally announced June 2012.
-
DC magnetic field generation in unmagnetized shear flows
Authors:
T. Grismayer,
E. P. Alves,
R. Fonseca,
L. O. Silva
Abstract:
The generation of DC magnetic fields in unmagnetized plasmas with velocity shear is predicted for non relativistic and relativistic scenarios either due to thermal effects or due to the onset of the Kelvin-Helmholtz instability (KHI). A kinetic model describes the growth and the saturation of the DC field. The predictions of the theory are confirmed by multidimensional particle-in-cell simulations…
▽ More
The generation of DC magnetic fields in unmagnetized plasmas with velocity shear is predicted for non relativistic and relativistic scenarios either due to thermal effects or due to the onset of the Kelvin-Helmholtz instability (KHI). A kinetic model describes the growth and the saturation of the DC field. The predictions of the theory are confirmed by multidimensional particle-in-cell simulations, demonstrating the formation of long lived magnetic fields ($t \sim 100s ω_{pi}^{-1}$) along the full longitudinal extent of the shear layer, with transverse width on the electron length scale ($\sqrt{γ_0}c/ω_{pe}$), reaching magnitudes $eB_{\mathrm{DC}}/m_ecω_{pe}\sim β_0\sqrt{γ_0}$.
△ Less
Submitted 14 November, 2012; v1 submitted 10 May, 2012;
originally announced May 2012.
-
Weibel-instability-mediated collisionless shocks in laboratory with ultraintense lasers
Authors:
F. Fiuza,
R. A. Fonseca,
J. Tonge,
W. B. Mori,
L. O. Silva
Abstract:
The formation of non-relativistic collisionless shocks in laboratory with ultrahigh intensity lasers is studied via \emph{ab initio} multi-dimensional particle-in-cell simulations. The microphysics behind shock formation and dissipation, and the detailed shock structure are analyzed, illustrating that the Weibel instability plays a crucial role in the generation of strong subequipartition magnetic…
▽ More
The formation of non-relativistic collisionless shocks in laboratory with ultrahigh intensity lasers is studied via \emph{ab initio} multi-dimensional particle-in-cell simulations. The microphysics behind shock formation and dissipation, and the detailed shock structure are analyzed, illustrating that the Weibel instability plays a crucial role in the generation of strong subequipartition magnetic fields that isotropize the incoming flow and lead to the formation of a collisionless shock, similarly to what occurs in astrophysical scenarios. The possibility of generating such collisionless shocks in laboratory opens the way to the direct study of the physics associated with astrophysical shocks.
△ Less
Submitted 26 April, 2012;
originally announced April 2012.
-
Large-scale magnetic field generation via the kinetic Kelvin-Helmholtz instability in unmagnetized scenarios
Authors:
E. P. Alves,
T. Grismayer,
S. F. Martins,
F. Fiúza,
R. A. Fonseca,
L. O. Silva
Abstract:
Collisionless plasma instabilities are fundamental in magnetic field generation in astrophysical scenarios, but their role has been addressed in scenarios where velocity shear is absent. In this work we show that velocity shears must be considered when studying realistic astrophysical scenarios, since these trigger the collisionless Kelvin-Helmholtz instability (KHI). We present the first self-con…
▽ More
Collisionless plasma instabilities are fundamental in magnetic field generation in astrophysical scenarios, but their role has been addressed in scenarios where velocity shear is absent. In this work we show that velocity shears must be considered when studying realistic astrophysical scenarios, since these trigger the collisionless Kelvin-Helmholtz instability (KHI). We present the first self-consistent three-dimensional (3D) particle-in-cell (PIC) simulations of the KHI in conditions relevant for unmagnetized relativistic outflows with velocity shear, such as active galactic nuclei (AGN) and gamma-ray bursts (GRBs). We show the generation of a strong large-scale DC magnetic field, which extends over the entire shear-surface, reaching thicknesses of a few tens of electron skin depths, and persisting on time-scales much longer than the electron time scale. This DC magnetic field is not captured by MHD models since it arises from intrinsically kinetic effects. Our results indicate that the KHI can generate intense magnetic fields yielding equipartition values up to ε_B/ε_p ~ 10^-3 in the electron time-scale. The KHI-induced magnetic fields have a characteristic structure that will lead to a distinct radiation signature, and can seed the turbulent dynamo amplification process. The dynamics of the KHI are relevant for non-thermal radiation modeling and can also have a strong impact on the formation of relativistic shocks in presence of velocity shears.
△ Less
Submitted 6 December, 2011; v1 submitted 29 July, 2011;
originally announced July 2011.
-
The nonlinear saturation of the non-resonant kinetically driven streaming instability
Authors:
L. Gargate,
R. A. Fonseca,
J. Niemiec,
M. Pohl,
R. Bingham,
L. O. Silva
Abstract:
A non-resonant instability for the amplification of the interstellar magnetic field in young Supernova Remnant (SNR) shocks was predicted by Bell (2004), and is thought to be relevant for the acceleration of cosmic ray (CR) particles. For this instability, the CRs streaming ahead of SNR shock fronts drive electromagnetic waves with wavelengths much shorter than the typical CR Larmor radius, by i…
▽ More
A non-resonant instability for the amplification of the interstellar magnetic field in young Supernova Remnant (SNR) shocks was predicted by Bell (2004), and is thought to be relevant for the acceleration of cosmic ray (CR) particles. For this instability, the CRs streaming ahead of SNR shock fronts drive electromagnetic waves with wavelengths much shorter than the typical CR Larmor radius, by inducing a current parallel to the background magnetic field. We explore the nonlinear regime of the non-resonant mode using Particle-in-Cell (PIC) hybrid simulations, with kinetic ions and fluid electrons, and analyze the saturation mechanism for realistic CR and background plasma parameters. In the linear regime, the observed growth rates and wavelengths match the theoretical predictions; the nonlinear stage of the instability shows a strong reaction of both the background plasma and the CR particles, with the saturation level of the magnetic field varying with the CR parameters. The simulations with CR-to-background density ratios of n_CR/n_b=10^(-5) reveal the highest magnetic field saturation levels, with energy also being transferred to the background plasma and to the perpendicular velocity components of the CR particles. The results show that amplification factors >10 for the magnetic field can be achieved, and suggest that this instability is important for the generation of magnetic field turbulence, and for the acceleration of CR particles.
△ Less
Submitted 8 February, 2010;
originally announced February 2010.
-
Ion dynamics and acceleration in relativistic shocks
Authors:
S. F. Martins,
R. A. Fonseca,
L. O. Silva,
W. B. Mori
Abstract:
Ab-initio numerical study of collisionless shocks in electron-ion unmagnetized plasmas is performed with fully relativistic particle in cell simulations. The main properties of the shock are shown, focusing on the implications for particle acceleration. Results from previous works with a distinct numerical framework are recovered, including the shock structure and the overall acceleration featur…
▽ More
Ab-initio numerical study of collisionless shocks in electron-ion unmagnetized plasmas is performed with fully relativistic particle in cell simulations. The main properties of the shock are shown, focusing on the implications for particle acceleration. Results from previous works with a distinct numerical framework are recovered, including the shock structure and the overall acceleration features. Particle tracking is then used to analyze in detail the particle dynamics and the acceleration process. We observe an energy growth in time that can be reproduced by a Fermi-like mechanism with a reduced number of scatterings, in which the time between collisions increases as the particle gains energy, and the average acceleration efficiency is not ideal. The in depth analysis of the underlying physics is relevant to understand the generation of high energy cosmic rays, the impact on the astrophysical shock dynamics, and the consequent emission of radiation.
△ Less
Submitted 3 April, 2009; v1 submitted 20 March, 2009;
originally announced March 2009.
-
Physical problems (microphysics) in relativistic plasma flows
Authors:
L. O. Silva
Abstract:
Many problems in astrophysics involve relativistic outflows. The plasma dynamics in these scenarios is critical to determine the conditions for the self-consistent evolution of the fields and particle acceleration. Advances in computer power now allow for kinetic plasma simulations, based on the particle-in-cell (PIC) paradigm, capable of providing information about the role of plasma instabilit…
▽ More
Many problems in astrophysics involve relativistic outflows. The plasma dynamics in these scenarios is critical to determine the conditions for the self-consistent evolution of the fields and particle acceleration. Advances in computer power now allow for kinetic plasma simulations, based on the particle-in-cell (PIC) paradigm, capable of providing information about the role of plasma instabilities in relativistic outflows. A discussion of the key issues associated with PIC simulations is presented, along with some the most important results and open questions, with a particular emphasis on the long time evolution of the filamentation, or Weibel, instability, and on the possible collisionless mechanisms for particle acceleration arising in the collision of relativistic plasma shells.
△ Less
Submitted 11 October, 2006;
originally announced October 2006.
-
Baryon loading and the Weibel instability in gamma-ray bursts
Authors:
M. Fiore,
L. O. Silva,
C. Ren,
M. A. Tzoufras,
W. B. Mori
Abstract:
The dynamics of two counter-streaming electron-positron-ion unmagnetized plasma shells with zero net charge is analyzed in the context of magnetic field generation in GRB internal shocks due to the Weibel instability. The effects of large thermal motion of plasma particles, arbitrary mixture of plasma species and space charge effects are taken into account. We show that, although thermal effects…
▽ More
The dynamics of two counter-streaming electron-positron-ion unmagnetized plasma shells with zero net charge is analyzed in the context of magnetic field generation in GRB internal shocks due to the Weibel instability. The effects of large thermal motion of plasma particles, arbitrary mixture of plasma species and space charge effects are taken into account. We show that, although thermal effects slow down the instability, baryon loading leads to a non-negligible growth rate even for large temperatures and different shell velocities, thus guaranteeing the robustness and the occurrence of the Weibel instability for a wide range of scenarios.
△ Less
Submitted 24 August, 2006;
originally announced August 2006.
-
Very High Mach Number Electrostatic Shocks in Collisionless Plasmas
Authors:
G. Sorasio,
M. Marti,
R. Fonseca,
L. O. Silva
Abstract:
The kinetic theory of collisionless electrostatic shocks resulting from the collision of plasma slabs with different temperatures and densities is presented. The theoretical results are confirmed by self-consistent particle-in-cell simulations, revealing the formation and stable propagation of electrostatic shocks with very high Mach numbers ($M \gg 10$), well above the predictions of the classi…
▽ More
The kinetic theory of collisionless electrostatic shocks resulting from the collision of plasma slabs with different temperatures and densities is presented. The theoretical results are confirmed by self-consistent particle-in-cell simulations, revealing the formation and stable propagation of electrostatic shocks with very high Mach numbers ($M \gg 10$), well above the predictions of the classical theories for electrostatic shocks.
△ Less
Submitted 23 December, 2005;
originally announced December 2005.
-
Cluster magnetic fields from large-scale-structure and galaxy-cluster shocks
Authors:
Mikhail V. Medvedev,
Luis O. Silva,
Marc Kamionkowski
Abstract:
The origin of the micro-Gauss magnetic fields in galaxy clusters is one of the outstanding problem of modern cosmology. We have performed three-dimensional particle-in-cell simulations of the nonrelativistic Weibel instability in an electron-proton plasma, in conditions typical of cosmological shocks. These simulations indicate that cluster fields could have been produced by shocks propagating t…
▽ More
The origin of the micro-Gauss magnetic fields in galaxy clusters is one of the outstanding problem of modern cosmology. We have performed three-dimensional particle-in-cell simulations of the nonrelativistic Weibel instability in an electron-proton plasma, in conditions typical of cosmological shocks. These simulations indicate that cluster fields could have been produced by shocks propagating through the intergalactic medium during the formation of large-scale structure or by shocks within the cluster. The strengths of the shock-generated fields range from tens of nano-Gauss in the intercluster medium to a few micro-Gauss inside galaxy clusters.
△ Less
Submitted 3 December, 2005;
originally announced December 2005.
-
Long-time evolution of magnetic fields in relativistic GRB shocks
Authors:
Mikhail V. Medvedev,
Massimiliano Fiore,
Ricardo A. Fonseca,
Luis O. Silva,
Warren B. Mori
Abstract:
We investigate the long-time evolution of magnetic fields generated by the two-stream instability at ultra- and sub-relativistic astrophysical collisionless shocks. Based on 3D PIC simulation results, we introduce a 2D toy model of interacting current filaments. Within the framework of this model, we demonstrate that the field correlation scale in the region far downstream the shock grows nearly…
▽ More
We investigate the long-time evolution of magnetic fields generated by the two-stream instability at ultra- and sub-relativistic astrophysical collisionless shocks. Based on 3D PIC simulation results, we introduce a 2D toy model of interacting current filaments. Within the framework of this model, we demonstrate that the field correlation scale in the region far downstream the shock grows nearly as the light crossing time, lambda(t) ~ ct, thus making the diffusive field dissipation inefficient. The obtained theoretical scaling is tested using numerical PIC simulations. This result extends our understanding of the structure of collisionless shocks in gamma-ray bursts and other astrophysical objects.
△ Less
Submitted 15 September, 2004;
originally announced September 2004.
-
Kinetics of Particles in Relativistic Collisionless Shocks
Authors:
Mikhail V. Medvedev,
Luis O. Silva,
Ricardo A. Fonseca,
J. W. Tonge,
Warren B. Mori
Abstract:
Charged plasma particles form an anisotropic counter-streaming distribution at the front of a shock. In the near- and ultra-relativistic regimes, Weibel (or two-stream) instability produces near-equipartition, chaotic, small-scale magnetic fields. The fields introduce effective collisions and "thermalize" the plasma particles via pitch-angle scattering. The properties of jitter radiation emitted…
▽ More
Charged plasma particles form an anisotropic counter-streaming distribution at the front of a shock. In the near- and ultra-relativistic regimes, Weibel (or two-stream) instability produces near-equipartition, chaotic, small-scale magnetic fields. The fields introduce effective collisions and "thermalize" the plasma particles via pitch-angle scattering. The properties of jitter radiation emitted by accelerated electrons from these small-scale fields are markedly different from synchrotron spectra. Here we theoretical summary and the results of recent numerical 3D PIC plasma kinetic simulations. The relation of the obtained results to the theory of gamma-ray bursts is outlined.
△ Less
Submitted 14 October, 2003;
originally announced October 2003.
-
Interpenetrating plasma shells: near-equipartition magnetic field generation and non-thermal particle acceleration
Authors:
L. O. Silva,
R. A. Fonseca,
J. Tonge,
J. M. Dawson,
W. B. Mori,
M. V. Medvedev
Abstract:
We present the first three-dimensional fully kinetic electromagnetic relativistic particle-in-cell simulations of the collision of two interpenetrating plasma shells. The highly accurate plasma-kinetic "particle-in-cell" (with the total of $10^8$ particles) parallel code OSIRIS has been used. Our simulations show: (i) the generation of long-lived near-equipartition (electro)magnetic fields, (ii)…
▽ More
We present the first three-dimensional fully kinetic electromagnetic relativistic particle-in-cell simulations of the collision of two interpenetrating plasma shells. The highly accurate plasma-kinetic "particle-in-cell" (with the total of $10^8$ particles) parallel code OSIRIS has been used. Our simulations show: (i) the generation of long-lived near-equipartition (electro)magnetic fields, (ii) non-thermal particle acceleration, and (iii) short-scale to long-scale magnetic field evolution, in the collision region. Our results provide new insights into the magnetic field generation and particle acceleration in relativistic and sub-relativistic colliding streams of particles, which are present in gamma-ray bursters, supernova remnants, relativistic jets, pulsar winds, etc..
△ Less
Submitted 29 July, 2003;
originally announced July 2003.
-
Nuclear reaction rates and energy in stellar plasmas : The effect of highly damped modes
Authors:
Merav Opher,
Luis O. Silva,
Dean E. Dauger,
Viktor K. Decyk,
John M. Dawson
Abstract:
The effects of the highly damped modes in the energy and reaction rates in a plasma are discussed. These modes, with wavenumbers $k \gg k_{D}$, even being only weakly excited, with less than $k_{B}T$ per mode, make a significant contribution to the energy and screening in a plasma. When the de Broglie wavelength is much less than the distance of closest approach of thermal electrons, a classical…
▽ More
The effects of the highly damped modes in the energy and reaction rates in a plasma are discussed. These modes, with wavenumbers $k \gg k_{D}$, even being only weakly excited, with less than $k_{B}T$ per mode, make a significant contribution to the energy and screening in a plasma. When the de Broglie wavelength is much less than the distance of closest approach of thermal electrons, a classical analysis of the plasma can
△ Less
Submitted 9 May, 2001;
originally announced May 2001.
-
Comment on "Ponderomotive force due to neutrinos"
Authors:
L. O. Silva,
R. Bingham,
J. M. Dawson,
P. K. Shukla,
N. L. Tsintsadze,
J. T. Mendonca
Abstract:
The derivation of the ponderomotive force due to neutrinos by Hardy and Melrose, Phys.Rev.D v.54, p.6491 (1996) is based on a flawed analysis of the ponderomotive force concept. Their conclusions also contain an erroneous physical assumption related to the neutrino emission in supernovae. A correct analysis shows the importance of the ponderomotive force due to neutrinos in type II supernovae ex…
▽ More
The derivation of the ponderomotive force due to neutrinos by Hardy and Melrose, Phys.Rev.D v.54, p.6491 (1996) is based on a flawed analysis of the ponderomotive force concept. Their conclusions also contain an erroneous physical assumption related to the neutrino emission in supernovae. A correct analysis shows the importance of the ponderomotive force due to neutrinos in type II supernovae explosions.
△ Less
Submitted 10 December, 1998; v1 submitted 28 July, 1998;
originally announced July 1998.
