Structure and evolution of magnetohydrodynamic solitary waves with Hall and finite Larmor radius effects.

Nonlinear and low-frequency solitary waves are investigated in the framework of the one-dimensional Hall-magnetohydrodynamic model with finite Larmor effects and two different closure models for the pressures. For a double adiabatic pressure model, the organization of these localized structures in terms of the propagation angle with respect to the ambient magnetic field θ and the propagation velocity C is discussed. There are three types of regions in the θ-C plane that correspond to domains where either solitary waves cannot exist, are organized in branches, or have a continuous spectrum. A numerical method valid for the two latter cases, which rigorously proves the existence of the waves, is presented and used to locate many waves, including bright and dark structures. Some of them belong to parametric domains where solitary waves were not found in previous works. The stability of the structures has been investigated by performing a linear analysis of the background plasma state and by means of numerical simulations. They show that the cores of some waves can be robust, but, for the parameters considered in the analysis, the tails are unstable. The substitution of the double adiabatic model by evolution equations for the plasma pressures appears to suppress the instability in some cases and to allow the propagation of the solitary waves during long times.

Publisher's Note: Relativistic breather-type solitary waves with linear polarization in cold plasmas [Phys. Rev. E 91, 033102 (2015)].

This corrects the article DOI: 10.1103/PhysRevE.91.033102.

Relativistic breather-type solitary waves with linear polarization in cold plasmas.

Linearly polarized solitary waves, arising from the interaction of an intense laser pulse with a plasma, are investigated. Localized structures, in the form of exact numerical nonlinear solutions of the one-dimensional Maxwell-fluid model for a cold plasma with fixed ions, are presented. Unlike stationary circularly polarized solitary waves, the linear polarization gives rise to a breather-type behavior and a periodic exchange of electromagnetic energy and electron kinetic energy at twice the frequency of the wave. A numerical method based on a finite-differences scheme allows us to compute a branch of solutions within the frequency range Ωmin<Ω<ωpe, where ωpe and Ωmin are the electron plasma frequency and the frequency value for which the plasma density vanishes locally, respectively. A detailed description of the spatiotemporal structure of the waves and their main properties as a function of Ω is presented. Small-amplitude oscillations appearing in the tail of the solitary waves, a consequence of the linear polarization and harmonic excitation, are explained with the aid of the Akhiezer-Polovin system. Direct numerical simulations of the Maxwell-fluid model show that these solitary waves propagate without change for a long time.

Modeling relativistic soliton interactions in overdense plasmas: a perturbed nonlinear Schrödinger equation framework.

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrödinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrödinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude.

Anticorrelation between ion acceleration and nonlinear coherent structures from laser-underdense plasma interaction.

In laser-plasma experiments, we observed that ion acceleration from the Coulomb explosion of the plasma channel bored by the laser is prevented when multiple plasma instabilities, such as filamentation and hosing, and nonlinear coherent structures (vortices or postsolitons) appear in the wake of an ultrashort laser pulse. The tailoring of the longitudinal plasma density ramp allows us to control the onset of these instabilities. We deduced that the laser pulse is depleted into these structures in our conditions, when a plasma at about 10% of the critical density exhibits a gradient on the order of 250 µm (Gaussian fit), thus hindering the acceleration. A promising experimental setup with a long pulse is demonstrated enabling the excitation of an isolated coherent structure for polarimetric measurements and, in further perspectives, parametric studies of ion plasma acceleration efficiency.

Two-dimensional s-polarized solitary waves in relativistic plasmas. I. The fluid plasma model.

The properties of two-dimensional linearly s-polarized solitary waves are investigated by fluid-Maxwell equations and particle-in-cell (PIC) simulations. These self-trapped electromagnetic waves appear during laser-plasma interactions, and they have a dominant electric field component E(z), normal to the plane of the wave, that oscillates at a frequency below the electron plasma frequency ω(pe). A set of equations that describe the waves are derived from the plasma fluid model in the case of cold or warm plasma and then solved numerically. The main features, including the maximum value of the vector potential amplitude, the total energy, the width, and the cavitation radius are presented as a function of the frequency. The amplitude of the vector potential increases monotonically as the frequency of the wave decreases, whereas the width reaches a minimum value at a frequency of the order of 0.82 ω(pe). The results are compared with a set of PIC simulations where the solitary waves are excited by a high-intensity laser pulse.

Two-dimensional s-polarized solitary waves in plasmas. II. Stability, collisions, electromagnetic bursts, and post-soliton evolution.

The dynamics of two-dimensional s-polarized solitary waves is investigated with the aid of particle-in-cell (PIC) simulations. Instead of the usual excitation of the waves with a laser pulse, the PIC code was directly initialized with the numerical solutions from the fluid plasma model. This technique allows the analysis of different scenarios including the theoretical problems of the solitary wave stability and their collision as well as features already measured during laser-plasma experiments such as the emission of electromagnetic bursts when the waves reach the plasma-vacuum interface, or their expansion on the ion time scale, usually named post-soliton evolution. Waves with a single density depression are stable whereas multihump solutions decay to several waves. Contrary to solitons, two waves always interact through a force that depends on their relative phases, their amplitudes, and the distance between them. On the other hand, the radiation pattern at the plasma-vacuum interface was characterized, and the evolution of the diameter of different waves was computed and compared with the "snow plow" model.

Quasicollapse of oblique solitons of the weakly dissipative derivative nonlinear Schrödinger equation.

Numerical integrations of the derivative nonlinear Schrödinger equation for Alfvén waves, supplemented by a weak dissipative term (originating from diffusion or Landau damping), with initial conditions in the form of a bright soliton with nonvanishing conditions at infinity (oblique soliton), reveal an interesting phenomenon of "quasicollapse": as the dissipation parameter is reduced, larger amplitudes are reached and smaller scales are created, but on an increasing time scale. This process involves an early bifurcation of the initial soliton toward a breather that is analyzed by means of a numerical inverse scattering technique. This evolution leads to the formation of persistent dark solitons that are only weakly affected when crossed by the decaying breather which has the form of either a localized structure or an extended wave packet.