scholarly journals Particle acceleration in a reconnecting current sheet: PIC simulation

2009 ◽  
Vol 75 (5) ◽  
pp. 619-636 ◽  
Author(s):  
TARAS V. SIVERSKY ◽  
VALENTINA V. ZHARKOVA

AbstractThe acceleration of protons and electrons in a reconnecting current sheet (RCS) is simulated with a particle-in-cell (PIC) 2D3V (two-dimensional in space and three-dimensional in velocity space) code for the proton-to-electron mass ratio of 100. The electromagnetic configuration forming the RCS incorporates all three components of the magnetic field (including the guiding field) and a drifted electric field. PIC simulations reveal that there is a polarization electric field that appears during acceleration owing to a separation of electrons from protons towards the midplane of the RCS. If the plasma density is low, the polarization field is weak and the particle trajectories in the PIC simulations are similar to those in the test particle (TP) approach. For the higher plasma density the polarization field is stronger and it affects the trajectories of protons by increasing their orbits during acceleration. This field also leads to a less asymmetrical abundance of ejected protons towards the midplane in comparison with the TP approach. For a given magnetic topology electrons in PIC simulations are ejected to the same semispace as protons, in contrast to the TP results. This happens because the polarization field extends far beyond the thickness of a current sheet. This field decelerates the electrons, which are initially ejected into the semispace opposite to the protons, returns them back to the RCS, and, eventually, leads to the electron ejection into the same semispace as protons. The energy distribution of the ejected electrons is rather wide and single-peaked, in contrast to the two-peak narrow-energy distribution obtained in the TP approach. In the case of a strong guiding field, the mean energy of the ejected electrons is found to be smaller than it is predicted analytically and by the TP simulations. The beam of accelerated electrons is also found to generate turbulent electric field in the form of Langmuir waves.

2009 ◽  
Vol 75 (2) ◽  
pp. 159-181 ◽  
Author(s):  
V. V. ZHARKOVA ◽  
O. V. AGAPITOV

AbstractElectron and proton acceleration by a drifted super-Dreicer electric field is investigated in a strongly compressed non-neutral reconnecting current sheet (NRCS). The guiding field is assumed to be constant within a reconnecting current sheet (RCS) and parallel to the direction of the drifted electric field. The other two magnetic field components, transverse and tangential, are considered to vary exponentially and linearly with distances from the X-nullpoint. The proton and electron energy spectra are calculated numerically in a model RCS with different magnetic field topologies by solving an equation of motion in the test-particle approach with some test with a particle-in-cell (PIC) approach. Three kinds electric field generated inside a RCS are considered: a drifted electric field caused by the plasma inflows formed during a magnetic reconnection process; a polarization electric field induced by the accelerated protons and electrons; and a turbulent electric field induced by instabilities generated by accelerated particles. Electron and proton densities, and energy spectra inside a RCS and at ejection are found to be strongly affected by the magnetic field topology: for stronger magnetic fields the spectra are softer having a small higher-energy cutoff while for weaker magnetic fields the spectra are harder with much larger upper cutoff energies. Depending on the magnetic component ratios and drifted electric field magnitude, particles are found to be ejected either as quasi-thermal flows with very high temperatures or as focused power-law beams. A polarization field is found to reduce the acceleration time inside a RCS and to increase the energy gained by particles at acceleration by a pure drifted electric field by a few orders of magnitude. The turbulent electric field induced by the two beam instabilities of the same kind of particles leads to a significant increase in the number of particles with higher energies resulting in a flattening of their energy spectra.


2008 ◽  
Vol 26 (8) ◽  
pp. 2471-2483 ◽  
Author(s):  
K. G. Tanaka ◽  
A. Retinò ◽  
Y. Asano ◽  
M. Fujimoto ◽  
I. Shinohara ◽  
...  

Abstract. The magnetopause (MP) reconnection is characterized by a density asymmetry across the current sheet. The asymmetry is expected to produce characteristic features in the reconnection layer. Here we present a comparison between the Cluster MP crossing reported by Retinò et al. (2006) and virtual observations in two-dimensional particle-in-cell simulation results. The simulation, which includes the density asymmetry but has zero guide field in the initial condition, has reproduced well the observed features as follows: (1) The prominent density dip region is detected at the separatrix region (SR) on the magnetospheric (MSP) side of the MP. (2) The intense electric field normal to the MP is pointing to the center of the MP at the location where the density dip is detected. (3) The ion bulk outflow due to the magnetic reconnection is seen to be biased towards the MSP side. (4) The out-of-plane magnetic field (the Hall magnetic field) has bipolar rather than quadrupolar structure, the latter of which is seen for a density symmetric case. The simulation also showed rich electron dynamics (formation of field-aligned beams) in the proximity of the separatrices, which was not fully resolved in the observations. Stepping beyond the simulation-observation comparison, we have also analyzed the electron acceleration and the field line structure in the simulation results. It is found that the bipolar Hall magnetic field structure is produced by the substantial drift of the reconnected field lines at the MSP SR due to the enhanced normal electric field. The field-aligned electrons at the same MSP SR are identified as the gun smokes of the electron acceleration in the close proximity of the X-line. We have also analyzed the X-line structure obtained in the simulation to find that the density asymmetry leads to a steep density gradient in the in-flow region, which may lead to a non-stationary behavior of the X-line when three-dimensional freedom is taken into account.


2003 ◽  
Vol 10 (1/2) ◽  
pp. 53-63 ◽  
Author(s):  
N. Singh

Abstract. We report here further results from the three-dimensional particle-in-cell simulations of the electron-beam driven electron holes. We focus here on (i) the transformation of oscillatory waves driven by the electron-beam instability into electron holes, (ii) the continued evolution and propagation of electron holes after their formation, including merging of electron holes, and (iii) the effects of the evolution on the plasma density and ion velocity distribution function. We find that initially electron-beam modes with perpendicular wave numbers k^ = 0 and as well as k^ ≠ 0 are driven resonantly below the electron plasma frequency of the target plasma. The modes interact nonlinearly and modulate each other both in space and time, producing wave structures with finite perpendicular scale lengths. Nonlinear evolution of such wave structures generates the electron holes in the simulations. Initially, a large number of electron holes form in the plasma. Their merging yields continuously a decreasing number of electron holes. The propagation velocity of the electron holes evolves dynamically and is affected by their merging. At late times only a few electron holes are left in the simulation and they decay by emitting low-frequency electrostatic whistler waves just above the lower hybrid (LH) frequency vlh . These waves, which are long structures parallel to the ambient magnetic field B0 and quite short transverse to B0, are associated with similar structures in the plasma density, producing density filaments. It turns out that electron-beam driven plasmas, in general, develop such filaments at some stage of the evolution of the beam-driven waves. In view of the excitation of the LH waves near vlh, which could resonate with the ions, an analysis shows that it is possible to heat transversely the ions in a time scale of a few seconds in the auroral return current plasma, in which electron holes and transversely heated ions have been simultaneously observed.


2007 ◽  
Vol 15 (2) ◽  
pp. 83-94 ◽  
Author(s):  
Joseph Wang ◽  
Yong Cao ◽  
Raed Kafafy ◽  
Viktor Decyk

A parallel, three-dimensional electrostatic PIC code is developed for large-scale electric propulsion simulations using parallel supercomputers. This code uses a newly developed immersed-finite-element particle-in-cell (IFE-PIC) algorithm designed to handle complex boundary conditions accurately while maintaining the computational speed of the standard PIC code. Domain decomposition is used in both field solve and particle push to divide the computation among processors. Two simulations studies are presented to demonstrate the capability of the code. The first is a full particle simulation of near-thruster plume using real ion to electron mass ratio. The second is a high-resolution simulation of multiple ion thruster plume interactions for a realistic spacecraft using a domain enclosing the entire solar array panel. Performance benchmarks show that the IFE-PIC achieves a high parallel efficiency of ≥ 90%


2011 ◽  
Vol 18 (5) ◽  
pp. 727-733 ◽  
Author(s):  
C. Huang ◽  
Q. Lu ◽  
Z. Yang ◽  
M. Wu ◽  
Q. Dong ◽  
...  

Abstract. Two-dimensional (2-D) particle-in-cell (PIC) simulations are performed to investigate the evolution of the electron current sheet (ECS) in guide field reconnection. The ECS is formed by electrons accelerated by the inductive electric field in the vicinity of the X line, which is then extended along the x direction due to the imbalance between the electric field force and Ampere force. The tearing instability is unstable when the ECS becomes sufficiently long and thin, and several seed islands are formed in the ECS. These tiny islands may coalesce and form a larger secondary island in the center of the diffusion region.


2014 ◽  
Vol 981 ◽  
pp. 859-862
Author(s):  
Hui Hui Wang ◽  
Da Gang Liu ◽  
La Qun Liu ◽  
Lin Meng

In this paper, the code of Particle-In-Cell/Monte Carlo Collision (PIC/MCC) for oxygen microwave breakdown is developed. This code is based on the three dimensional particle-in-cell platform CHIPIC, and with a module for increasing the charge of each super-particle. With this PIC/MCC code, the multiplication rate of the electron density and the delay time in oxygen breakdown at atmospheric conditions are researched. The results show: the multiplication rate of the electron density is periodic, and its period is the half of the electric field period; the breakdown delay time in the gas breakdown increases while the frequency of electric field or the gas pressure increases.


Sign in / Sign up

Export Citation Format

Share Document