Search

The evolution of a deformed subcritical fast magnetosonic shock front is compared between two two-dimensional PIC simulations with different orientations of the magnetic field relative to the simulation box. All other initial and simulation conditions are kept identical. Shock boundary oscillations are observed in the simulation where the magnetic field direction is resolved. This oscillation is caused by the reformation of the shock front. One part of the front acts as a shock, while the other functions as a magnetic piston, with both halves changing their states in antiphase. The oscillation period corresponds to the time required for one shock wave to grow as the other collapses. In contrast, the corrugated fast magnetosonic shock does not oscillate in the second simulation, where the magnetic field is oriented out of the simulation plane. This dependence on magnetic field orientation suggests that the shock oscillation is induced by magnetic tension, which is only effective in the first simulation. In both simulations, the shock perturbation does not grow over time, indicating that the shocks are stable. The potential relevance of these findings for the Alfv\'enic oscillations of the supercritical Earth's bow shock, detected by the MMS multi-spacecraft mission, is also discussed., Comment: 16 pages, 7 figures, accepted for publication in Plasma Physics and Controlled Fusion

We study with particle-in-cell (PIC) simulations the stability of fast magnetosonic shocks. They expand across a collisionless plasma and an orthogonal magnetic field that is aligned with one of the directions resolved by the 2D simulations. The shock speed is 1.6 times the fast magnetosonic speed when it enters a layer with a reduced density of mobile ions, which decreases the shock speed by up to 15\% in 1D simulations. In the 2D simulations, the density of mobile ions in the layer varies sinusoidally perpendicularly to the shock normal. We resolve one sine period. This variation only leads to small changes in the shock speed evidencing a restoring force that opposes a shock deformation. As the shock propagates through the layer, the ion density becomes increasingly spatially modulated along the shock front and the magnetic field bulges out where the mobile ion density is lowest. The perturbed shock eventually reaches a steady state. Once it leaves the layer, the perturbations of the ion density and magnetic field oscillate along its front at a frequency close to the lower-hybrid frequency; the shock is mediated by a standing wave composed of obliquely propagating lower-hybrid waves. We perform three 2D simulations with different box lengths along the shock front. The shock front oscillations are aperiodically damped in the smallest box with the fastest variation of the ion density, strongly damped in the intermediate one, and weakly damped in the largest box. The shock front oscillations perturb the magnetic field in a spatial interval that extends by several electron skin depths upstream and downstream of the shock front and could give rise to Whistler waves that propagate along the shock's magnetic field overshoot. Similar waves were observed in hybrid and PIC simulations and by the MMS satellite mission., Comment: 25 pages, 12 figures, accepted for publication in Physica Scripta

We study with a 3D PIC simulation discontinuities between an electron-positron pair plasma and magnetized electrons and protons. A pair plasma is injected at one simulation boundary with a speed 0.6$c$ along its normal. It expands into an electron-proton plasma and a magnetic field that points orthogonally to the injection direction. Diamagnetic currents expel the magnetic field from within the pair plasma and pile it up in front of it. It pushes electrons, which induces an electric field pulse ahead of the magnetic one. This initial electromagnetic pulse (EMP) confines the pair plasma magnetically and accelerates protons electrically. The fast flow of the injected pair plasma across the protons behind the initial EMP triggers the filamentation instability. Some electrons and positrons cross the injection boundary and build up a second EMP. Electron-cyclotron drift instabilities perturb the plasma ahead of both EMPs seeding a Rayleigh-Taylor-type instability. Despite equally strong perturbations ahead of both EMPs, the second EMP is much more stable than the initial one. We attribute the rapid collapse of the initial EMP to the filamentation instability, which perturbed the plasma behind it. The Rayleigh-Taylor-type instability transforms the planar EMPs into transition layers, in which magnetic flux ropes and electrostatic forces due to uneven numbers of electrons and positrons slow down and compress the pair plasma and accelerate protons. In our simulation, the expansion speed of the pair cloud decreased by about an order of magnitude and its density increased by the same factor. Its small thickness implies that it is capable of separating a relativistic pair outflow from an electron-proton plasma, which is essential for collimating relativistic jets of pair plasma in collisionless astrophysical plasma., Comment: 25 pages, 12 figures, provisionally accepted for publication by the New Journal of Physics

By means of a particle-in-cell (PIC) simulation, we study the interaction between a uniform magnetized ambient electron-proton plasma at rest and an unmagnetized pair plasma, which we inject at one simulation boundary with a mildly relativistic mean speed and temperature. The magnetic field points out of the simulation plane. The injected pair plasma expels the magnetic field and piles it up at its front. It traps ambient electrons and drags them across the protons. An electric field grows, which accelerates protons into the pair cloud's expansion direction. This electromagnetic pulse separates the pair cloud from the ambient plasma. Electrons and positrons, which drift in the pulse's nonuniform field, trigger an instability that disrupts the current sheet ahead of the pulse. The wave vector of the growing perturbation is orthogonal to the magnetic field direction and magnetic tension cannot stabilize it. The electromagnetic pulse becomes permeable for pair plasma, which forms new electromagnetic pulses ahead of the initial one. A transition layer develops with a thickness of a few proton skin depths, in which protons and positrons are accelerated by strong electromagnetic fields. Protons form dense clumps surrounded by a strong magnetic field. The thickness of the transition layer grows less rapidly than we would expect from the typical speeds of the pair plasma particles and the latter transfer momentum to protons; hence, the transition layer acts as a discontinuity, separating the pair plasma from the ambient plasma. Such a discontinuity is an important building block for astrophysical pair plasma jets., Comment: Accepted for publication in Physics of Plasmas, 13 pages, 15 figures