LOCAL TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE COLLISIONLESS MAGNETOROTATIONAL INSTABILITY

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the plasma in the disk is essentially collisionless. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell plasma simulations. We focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field, B, until the Alfven velocity, v(A), is comparable to the speed of light, c (independent of the initial value of v(A)/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the pressure perpendicular to B being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume-averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission-from the radio to the gamma-rays-of systems such as Sgr A*.

TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10(15.5) eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized nonpropagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times, we observe saturation of the instability in the simulation, but the mechanism responsible is an artifact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

Two-dimensional particle-in-cell simulation of the expansion of a plasma into a rarefied medium

The expansion of a dense plasma through a more rarefied ionized medium has been studied by means of two-dimensional particle-in-cell simulations. The initial conditions involve a density jump by a factor of 100, located in the middle of an otherwise equally dense electron-proton plasma with uniform proton and electron temperatures of 10 eV and 1 keV, respectively. Simulations show the creation of a purely electrostatic collisionless shock together with an ion-acoustic soliton tied to its downstream region. The shock front is seen to evolve in filamentary structures consistently with the onset of the ion-ion instability. Meanwhile, an un-magnetized drift instability is triggered in the core part of the dense plasma. Such results explain recent experimental laser-plasma experiments, carried out in similar conditions, and are of intrinsic relevance to non-relativistic shock scenarios in the solar and astrophysical systems.

Particle-in-cell simulation study of a lower-hybrid shock

The expansion of a magnetized high-pressure plasma into a low-pressure ambient medium is examined with particle-in-cell simulations. The magnetic field points perpendicular to the plasma's expansion direction and binary collisions between particles are absent. The expanding plasma steepens into a quasi-electrostatic shock that is sustained by the lower-hybrid (LH) wave. The ambipolar electric field points in the expansion direction and it induces together with the background magnetic field a fast E cross B drift of electrons. The drifting electrons modify the background magnetic field, resulting in its pile-up by the LH shock. The magnetic pressure gradient force accelerates the ambient ions ahead of the LH shock, reducing the relative velocity between the ambient plasma and the LH shock to about the phase speed of the shocked LH wave, transforming the LH shock into a nonlinear LH wave. The oscillations of the electrostatic potential have a larger amplitude and wavelength in the magnetized plasma than in an unmagnetized one with otherwise identical conditions. The energy loss to the drifting electrons leads to a noticeable slowdown of the LH shock compared to that in an unmagnetized plasma.

A Particle-in-cell scheme of the RFQ in the SSC-Linac

A 52 MHz Radio Frequency Quadrupole (RFQ) linear accelerator (linac) is designed to serve as an initial structure for the SSC-Linac system (injector into Separated Sector Cyclotron). The designed injection and output energy are 3.5 keV/u and 143 keV/u, respectively. The beam dynamics in this RFQ have been studied using a three-dimensional Particle-In-Cell (PIC) code BEAMPATH. Simulation results show that this R,FQ structure is characterized by stable values of beam transmission efficiency (at least 95%) for both zero-current mode and the space charge dominated regime. The beam accelerated in the RFQ has good quality in both transverse and longitudinal directions, and could easily be accepted by Drift Tube Linac (DTL). The effect of the vane error and that of the space charge on the beam parameters have been studied as well to define the engineering tolerance for RFQ vane machining and alignment.

Plasma Ionization in Low-Pressure Radio-Frequency Discharges—Part II: Particle-in-Cell Simulation

Plasma ionization in the low-pressure operation regime ( $«$ 5 Pa) of RF capacitively coupled plasmas (CCPs) is governed by a complex interplay of various mechanisms, such as field reversal, sheath expansion, and wave–particle interactions. In a previous paper, it was shown that experimental observations in a hydrogen CCP operated at 13.56 MHz are qualitatively well described in a 1-D symmetrical particle-in-cell (PIC) simulation. In this paper, a spherical asymmetrical PIC simulation that is closer to the conditions of the highly asymmetrical experimental device is used to simulate a low-pressure neon CCP operated at 2 MHz. The results show a similar behavior, with pronounced ionization through field reversal, sheath expansion, and wave–particle interactions, and can be exploited for more accurate quantitative comparisons with experimental observations.

Particle-in-cell simulation study of the interaction between a relativistically moving leptonic micro-cloud and ambient electrons

Context. The jets of compact accreting objects are composed of electrons and a mixture of positrons and ions. These outflows impinge on the interstellar or intergalactic medium and both plasmas interact via collisionless processes. Filamentation (beam-Weibel) instabilities give rise to the growth of strong electromagnetic fields. These fields thermalize the interpenetrating plasmas. 

Aims. Hitherto, the effects imposed by a spatial non-uniformity on filamentation instabilities have remained unexplored. We examine the interaction between spatially uniform background electrons and a minuscule cloud of electrons and positrons. The cloud size is comparable to that created in recent laboratory experiments and such clouds may exist close to internal and external shocks of leptonic jets. The purpose of our study is to determine the prevalent instabilities, their ability to generate electromagnetic fields and the mechanism, by which the lepton micro-cloud transfers energy to the background plasma. 

Methods. A square micro-cloud of equally dense electrons and positrons impinges in our particle-in-cell (PIC) simulation on a spatially uniform plasma at rest. The latter consists of electrons with a temperature of 1 keV and immobile ions. The initially charge- and current neutral micro-cloud has a temperature of 100 keV and a side length of 2.5 plasma skin depths of the micro-cloud. The side length is given in the reference frame of the background plasma. The mean speed of the micro-cloud corresponds to a relativistic factor of 15, which is relevant for laboratory experiments and for relativistic astrophysical outflows. The spatial distributions of the leptons and of the electromagnetic fields are examined at several times. 

Results. A filamentation instability develops between the magnetic field carried by the micro-cloud and the background electrons. The electromagnetic fields, which grow from noise levels, redistribute the electrons and positrons within the cloud, which boosts the peak magnetic field amplitude. The current density and the moduli of the electromagnetic fields grow aperiodically in time and steadily along the direction that is anti-parallel to the cloud's velocity vector. The micro-cloud remains conjoined during the simulation. The instability induces an electrostatic wakefield in the background plasma. 

Conclusions. Relativistic clouds of leptons can generate and amplify magnetic fields even if they have a microscopic size, which implies that the underlying processes can be studied in the laboratory. The interaction of the localized magnetic field and high-energy leptons will give rise to synchrotron jitter radiation. The wakefield in the background plasma dissipates the kinetic energy of the lepton cloud. Even the fastest lepton micro-clouds can be slowed down by this collisionless mechanism. Moderately fast charge- and current neutralized lepton micro-clouds will deposit their energy close to relativistic shocks and hence they do not constitute an energy loss mechanism for the shock.

Bare-tether sheath and current: comparison of asymptotic theory and kinetic simulations in stationary plasma

Analytical expressions for current to a cylindrical Langmuir probe at rest in unmagnetized plasma are compared with results from both steady-state Vlasov and particle-in-cell simulations. Probe bias potentials that are much greater than plasma temperature (assumed equal for ions and electrons), as of interest for bare conductive tethers, are considered. At a very high bias, both the electric potential and the attracted-species density exhibit complex radial profiles; in particular, the density exhibits a minimum well within the plasma sheath and a maximum closer to the probe. Excellent agreement is found between analytical and numerical results for values of the probe radiusR close to the maximum radius Rmax for orbital-motion-limited (OML) collection at a particular bias in the following number of profile features: the values and positions of density minimum and maximum, position of sheath boundary, and value of a radius characterizing the no-space-charge behavior of a potential near the high-bias probe. Good agreement between the theory and simulations is also found for parametric laws jointly covering the following three characteristic R ranges: sheath radius versus probe radius and bias for Rmax; density minimum versus probe bias for Rmax; and (weakly bias-dependent) current drop below the OML value versus the probe radius for R > Rmax.

Laser-driven coherent betatron oscillation in a laser-wakefield cavity

The origin of beam disparity in emittance and betatron oscillation orbits, in and out of the polarization plane of the drive laser of laser-plasma accelerators, is explained in terms of betatron oscillations driven by the laser field. As trapped electrons accelerate, they move forward and interact with the laser pulse. For the bubble regime, a simple model is presented to describe this interaction in terms of a harmonic oscillator with a driving force from the laser and a restoring force from the plasma wake field. The resulting beam oscillations in the polarization plane, with period approximately the wavelength of the driving laser, increase emittance in that plane and cause microbunching of the beam. These effects are observed directly in 3D particle-in-cell simulations.

Electron self-injection and acceleration driven by a tightly focused intense laser beam in an underdense plasma

A scheme for electron self-injection in the laser wakefield acceleration is proposed. In this scheme, the transverse wave breaking of the wakefield and the tightly focused geometry of the laser beam play important roles. A large number of the background electrons are self-injected into the acceleration phase of the wakefield during the defocusing of the tightly focused laser beam as it propagates through an underdense plasma. Particle-in-cell simulations performed using a 2D3V code have shown generation of a collimated electron bunch with a total number of 1.4 x 109 and energies up to 8 MeV. (C) 2005 American Institute of Physics.

Electron injection by a nanowire in the bubble regime

The triggering of wave-breaking in a three-dimensional laser plasma wake (bubble) is investigated. The Coulomb potential from a nanowire is used to disturb the wake field to initialize the wave-breaking. The electron acceleration becomes more stable and the laser power needed for self-trapping is lowered. Three-dimensional particle-in-cell simulations were performed. Electrons with a charge of about 100 pC can be accelerated stably to energy about 170 MeV with a laser energy of 460 mJ. The first step towards tailoring the electron beam properties such as the energy, energy spread, and charge is discussed. (C) 2007 American Institute of Physics.

Electron acceleration by a propagating laser pulse in vacuum

Electrons accelerated by a propagating laser pulse of linear or circular polarization in vacuum have been investigated by one-dimensional particle-in-cell simulations and analytical modeling. A stopping target is used to stop the laser pulse and extract the energetic electrons from the laser field. The effect of the reflected light is taken into account. The maximum electron energy depends on the laser intensity and initial electron energy. There is an optimal acceleration length for electrons to gain maximum energy where electrons meet the peak of the laser pulse. The optimal acceleration length depends strongly on the laser pulse duration and amplitude. (C) 2007 American Institute of Physics.

Bubble regime for ion acceleration in a laser-driven plasma

Proton trapping and acceleration by an electron bubble-channel structure in laser interaction with high-density plasma is investigated by using three-dimensional particle-in-cell simulations. It is shown that protons can be trapped, bunched, and efficiently accelerated for appropriate laser and plasma parameters, and the proton acceleration is enhanced if the plasma consists mainly of heavier ions such as tritium. The observed results are analyzed and discussed in terms of a one-dimensional analytical three-component-plasma wake model.

Macroscopic evidence of soliton formation in multiterawatt laser-plasma interaction

A novel physical phenomenon has been observed following the interaction of an intense (10(19) W/cm(2)) laser pulse with an underdense plasma. Long-lived, macroscopic bubblelike structures have been detected through the deflection that the associated electric charge separation causes in a proton probe beam. These structures are interpreted as the remnants of a cloud of relativistic solitons generated in the plasma by the ultraintense laser pulse. This interpretation is supported by an analytical study of the soliton cloud evolution, by particle-in-cell simulations, and by a reconstruction of the proton-beam deflection.

Plasma ion evolution in the wake of a high-intensity ultrashort laser pulse

Experimental investigations of the late-time ion structures formed in the wake of an ultrashort, intense laser pulse propagating in a tenuous plasma have been performed using the proton imaging technique. The pattern found in the wake of the laser pulse shows unexpectedly regular modulations inside a long, finite width channel. On the basis of extensive particle in cell simulations of the plasma evolution in the wake of the pulse, we interpret this pattern as due to ion modulations developed during a two-stream instability excited by the return electric current generated by the wakefield.