Supernova remnant and heliopause termination shock plasmas may contain significant populations of minority heavy ions, with relative number densities n α /n i up to 50%. Preliminary kinetic simulations of collisionless shocks in these environments showed that the reformation cycle and acceleration mechanisms at quasi-perpendicular shocks can depend on the value of n α /n i . Shock reformation unfolds on ion spatio-temporal scales, requiring fully kinetic simulations of particle dynamics, together with the self-consistent electric and magnetic fields. This paper presents the first set of particle-in-cell simulations for two ion species, protons (n p ) and α-particles (n α ), with differing mass and charge-to-mass ratios, that spans the entire range of n α /n i from 0% to 100%. The interplay between the differing gyro length scales and timescales of the ion species is crucial to the time-evolving phenomenology of the shocks, the downstream turbulence, and the particle acceleration at different n α /n i . We show how the overall energization changes with n α /n i , and relate this to the processes individual ions undergo in the shock region and in the downstream turbulence, and to the power spectra of magnetic field fluctuations. The crossover between shocks dominated by the respective ion species happens when n α /n i = 25%, and minority ion energization is strongest in this regime. Energization of the majority ion species scales with injection energy. The power spectrum of the downstream turbulence includes peaks at sequential ion cyclotron harmonics, suggestive of ion ring-beam collective instability.

With reference to laboratory Q-machine studies we analyze the dynamics of a plasma diode under external forcing. Assuming a strong axial magnetic field, the problem is analyzed in one spatial dimension by a particle-in-cell code. The cathode is assumed to be operated in electron rich conditions, supplying an abundance of electrons. We compare different forcing schemes with the results obtained by solving the van der Pol equation. In one method of forcing we apply an oscillation in addition to the DC end plate bias and consider both amplitude and frequency variations. An alternative method of perturbation consists of modelling an absorbing grid at some internal position. Also in this case we can have a constant frequency with varying amplitude or alternatively an oscillation with chirped frequency but constant amplitude. We find that the overall features of the forced van der Pol equation are recovered, but the details in the plasma response need more attention to the harmonic responses, requiring extensions of the model equation. The analysis is extended by introducing collisional effects, where we emphasize charge exchange collisions of ions, since these processes usually have the largest cross sections and give significant modifications of the diode performance. In particular we find a reduction in oscillator frequency, although a linear scaling of the oscillation time with the system length remains also in this case. V C 2012 American Institute of Physics. [http://dx.

Small solid metallic objects in relative motion to thermal plasmas are studied by numerical simulations. We analyze supersonic motions, where a distinctive ion wake is formed behind obstacles. At these plasma drift velocities, ions enter the wake predominantly due to deflections by the electric field in the sheath around the obstacle. By irradiating the back side of the object by ultraviolet (UV) light, we can induce also an enhanced photo-electron population there. The resulting charge distribution gives rise to a pronounced local potential and plasma density well behind the object. This potential variation has the form of a three-dimensional ion acoustic double layer, containing also an ion phase space vortex. The analysis is supported also by one-dimensional numerical simulations to illustrate the importance of boundary conditions, Dirichlet and von Neumann conditions in particular.

The space-time evolution of an initial step-like plasma density variation is studied. We give particular attention to formulate the problem in a way that opens for the possibility of realizing the conditions experimentally. After a short transient time interval of the order of the electron plasma period, the solution is self-similar as illustrated by a video where the space-time evolution is reduced to be a function of the ratio x/t. Solutions of this form are usually found for problems without characteristic length and time scales, in our case the quasi-neutral limit. By introducing ion collisions with neutrals into the numerical analysis, we introduce a length scale, the collisional mean free path. We study the breakdown of the self-similarity of the solution as the mean free path is made shorter than the system length. Analytical results are presented for charge exchange collisions, demonstrating a short time collisionless evolution with an ensuing long time diffusive relaxation of the initial perturbation. For large times, we find a diffusion equation as the limiting analytical form for a charge-exchange collisional plasma, with a diffusion coefficient defined as the square of the ion sound speed divided by the (constant) ion collision frequency. The ion-neutral collision frequency acts as a parameter that allows a collisionless result to be obtained in one limit, while the solution of a diffusion equation is recovered in the opposite limit of large collision frequencies. V C 2013 AIP Publishing LLC. [http://dx.

No abstract