J. C. Adam

OSIRIS: A Three-Dimensional, Fully Relativistic Particle in Cell Code for Modeling Plasma Based Accelerators

We describe OSIRIS, a three-dimensional, relativistic, massively parallel, object oriented particle-in-cell code for modeling plasma based accelerators. Developed in Fortran 90, the code runs on multiple platforms (Cray T3E, IBM SP, Mac clusters) and can be easily ported to new ones. Details on the codes capabilities are given. We discuss the object-oriented design of the code, the encapsulation of system dependent code and the parallelization of the algorithms involved. We also discuss the implementation of communications as a boundary condition problem and other key characteristics of the code, such as the moving window, open-space and thermal bath boundaries, arbitrary domain decomposition, 2D (cartesian and cylindric) and 3D simulation modes, electron sub-cycling, energy conservation and particle and field diagnostics. Finally results from three-dimensional simulations of particle and laser wakefield accelerators are presented, in connection with the data analysis and visualization infrastructure developed to post-process the scalar and vector results from PIC simulations.

Study of stationary plasma thrusters using two-dimensional fully kinetic simulations

Stationary plasma thrusters are devices that use crossed electric and magnetic fields to accelerate ions to high velocities. Ions are created by collisional ionization of a propellant gas with electrons injected from a hollow cathode external to the thruster. A major issue is the electron transport through the magnetic field. It is known to exceed considerably the values predicted by the classical theory. Various 2D models have shown that wall collisions, which have often been invoked as the origin of this anomalous transport, are in fact insufficient. Anomalous turbulent transport has to be added to the model to recover an adequate conductivity. In the present paper the first 2D kinetic model that shows that, indeed, plasma turbulence can explain the observed conductivity is presented. Without any free parameter the model is able to reproduce numerous experimental features. At the end of the paper a preliminary theoretical analysis of the observed instability is provided.

Modulational and Raman instabilities in the relativistic regime

A large amplitude electromagnetic wave propagating in a plasma is known to be subject to severe modulational and Raman instabilities. Previous works were devoted to the weakly relativistic limit and applied mainly to a cold underdense plasma. One extends these works to include the fully relativistic limit for a circularly polarized light for which one derives the dispersion relation in a one‐dimensional plasma. The characteristics of the instabilities are also calculated in the case where the plasma is classically overdense, with 1<(ωp/ω0)2<γ, where ωp is the plasma frequency, ω0 is the laser frequency, and γ is the relativistic factor of an electron in the laser field. Particle‐in‐cell simulations confirm the results of the numerical solutions of the dispersion relation. For (ωp/ω0)2/γ=0.57 the growth rate can be as large as 0.52ω0. The nonlinear stage of the instability results in a strong heating of the electron distribution function. The theory is further extended to the case of an initially hot plasm...

Application of adaptive mesh refinement to particle-in-cell simulations of plasmas and beams

Plasma simulations are often rendered challenging by the disparity of scales in time and in space which must be resolved. When these disparities are in distinctive zones of the simulation domain, a method which has proven to be effective in other areas (e.g. fluid dynamics simulations) is the mesh refinement technique. We briefly discuss the challenges posed by coupling this technique with plasma Particle-In-Cell simulations, and present examples of application in Heavy Ion Fusion and related fields which illustrate the effectiveness of the approach. We also report on the status of a collaboration under way at Lawrence Berkeley National Laboratory between the Applied Numerical Algorithms Group (ANAG) and the Heavy Ion Fusion group to upgrade ANAGs mesh refinement library Chombo to include the tools needed by Particle-In-Cell simulation codes.

Propagation of ultraintense laser pulses through overdense plasma layers

Due to relativistic effects, a large amplitude electromagnetic wave can propagate in a classically overdense plasma with ω2p≳ω2≳ω2p/γ, where ωp is the plasma frequency, ω the laser frequency, and γ the relativistic factor of an electron in the laser field. Particle‐in‐cell simulations are used to study the interaction of an ultrahigh intensity laser pulse in normal incidence on a one‐dimensional preformed plasma layer. Both electrons and ions dynamics are included. The width of the layer is 10 to 30 μm and the plasma is characterized by (ωp/ω)2=1.5. During the penetration of the electromagnetic wave a large longitudinal electric field is generated. It results in a strong longitudinal heating of electrons which reach relativistic temperatures. This heating further lowers the effective plasma frequency ωp/γ so the layer becomes almost transparent after the plasma crossing by the wave front. Velocity of the wave front, reflection and transmission rates are studied as functions of the incident energy flux, th...

Formation of ion‐acoustic double layers

The formation of an ion‐acoustic double layer from a potential well is studied by computer simulation and explained by a modified Korteweg–de Vries (KdV) equation taking into account the reflections of the electrons.

Electron sub-cycling in particle simulation of plasma

Abstract A straightforward modification which reduces by half the computational cost of standard particle-in-cell algorithms for simulation of plasmas is described. The saving is obtained by integrating only the electrons through a number of time steps (sub-cycles) in order to resolve their evolution, while integrating the much slower ions only once per cycle, i.e., to match the time step of each species to their characteristic frequencies. A dispersion relation is derived which describes the numerical instabilities expected by sampling frequency arguments. Simulations support the broad features of the analytical results, viz., the maximum growth rate and domain of the instability, and its stabilization by the addition of weak damping. An implicit sub-cycling algorithm is suggested which may provide further saving while avoiding a limitation of implicit algorithms described elsewhere.

Efficiency of resonant absorption of electromagnetic waves in an inhomogeneous plasma

The resonant absorption of electromagnetic waves in an inhomogeneous plasma is studied for values of the resonant plasma wave field such that the ponderomotive effects modify wave propagation. The stationary solutions are obtained and found unstable with low dissipation in most cases. The absorption coefficient is then time dependent. Density cavities are emitted at the resonance and they can induce large frequency shifts for the reflected light.

Strongly enhanced laser absorption and electron acceleration via resonant excitation of surface plasma waves

Two-dimensional (2D) particle-in-cell numerical simulations of the interaction between a high-intensity short-pulse p-polarized laser beam and an overdense plasma are presented. It is shown that, under appropriate physical conditions, a surface plasma wave can be resonantly excited by a short-pulse laser wave, leading to strong relativistic electron acceleration together with a dramatic increase, up to 70%, of light absorption by the plasma. Purely 2D effects contribute to enhancement of electron acceleration. It is also found that the angular distribution of the hot electrons is drastically affected by the surface wave. The subsequent ion dynamics is shown to be significantly modified by the surface plasma wave excitation.

Subfemtosecond, coherent, relativistic, and ballistic electron bunches generated at ω0 and 2ω0 in high intensity laser-matter interaction

Harmonics of the laser light have been observed from the rear side of solid targets irradiated by a laser beam at relativistic intensities. This emission evidences the acceleration of subfemtosecond electron bunches by the laser pulse in front of the target. These bunches emit coherent transition radiation (CTR) when passing through the back surface of the target. The spectral features of the signal recorded for targets of thicknesses up to several hundred microns are consistent with the electrons being accelerated by both the laser electric field—via vacuum heating and/or resonance absorption,—and the v×B component of the Lorentz force. The spatial study of the radiation shows that the relativistic electrons causing the CTR radiation are coherent and propagate ballistically through the target, originating from a source with a size of the order of the laser focal spot.