Abstract:
Recent fluid simulations show that the ponderomotive potentials and ponderomotive frequencies of
electrons and ions can be used as proxies to identify steepening and breaking of the ion acoustic
solitary waves (IASWs) in plasmas. However, the behavior of these proxies may deviate in the
presence of kinetic effects such as particle trapping. We performed one-dimensional particle-incell
(PIC) simulations to examine the effects of kinetic processes on the behavior of these proxies
at the breaking of IASWs in plasmas. The electron and ion equilibrium densities were superimposed
by a long-wavelength Gaussian type perturbation, which initially evolves into two IASWs
observed as two phase space vortices due to the trapping of electrons in the ion acoustic (IA) potentials.
These IASW structures grow due to the steepening of their trailing edges, and later they break
into a chain of IA phase space vortices. Each of these vortices is associated with a bipolar electric
field resulting in a positive potential structure. We examined the amplitude, width, and phase velocity
of the IASWs at their breaking process to clarify their link with the trapping velocity. In addition,
we estimated electron and ion ponderomotive potentials and frequencies from the PIC
simulations to verify their applicability in identifying wave breaking limit under the kinetic regime.
The present study shows that the behavior of the ponderomotive potential during the IA wave
breaking process is similar to the one, which is proposed through fluid simulations. We find that IA
wave breaking occurs when the maximum trapping velocity of the electron (Vtrap þ Vs) exceeds its
thermal velocity. The present simulation study shows that both maximum electron trapping velocity
and ponderomotive potential can be used to identify the IA wave breaking processes in plasmas.