Ion Acoustic Soliton and the Lambert Function

The Sagdeev potential method is employed to compute the width of (Ion-acoustic) soliton propagated in a cold plasma. The computation indicates that the soliton width is a continuous function (of the Mach number M), which is expressed in terms of the Lambert Function. Despite the (fairly) complex form of the function, the numerical plotting makes sense about its changes.


Introduction
Over the last decades, there has been a great deal of interest and significant progress in the study of nonlinear plasma theories and many of works and researches in plasma physics devote much attention to these theories.The nonlinear theories include a large number of effects and phenomena such as the nonlinear coherent structures as shock waves, solitary waves (solitons), vortices, etc.
The collective electrical and magnetic properties of plasmas could produce interactions that take the place of collisions and permit shocks and solitons to form.A shock wave is a sudden transition (a type of propagating disturbance) in the properties of a fluid medium (liquid or gas), involving a difference in flow velocity across a narrow (ideally, abrupt) transition.In high-energy density physics, nearly any experiment involves at least one shock wave.Such shock waves may be also produced by applying pressure to a surface or by creating a collision between two materials.
Nonlinear effects in plasmas occur when a large amplitude wave is excited by an external means.Soliton (shock) waves are formed as a result of a balance between the nonlinearity and dispersion (dissipation) of the medium.For example, the dispersion in the ion acoustic wave can be counter-balanced by nonlinearity and an ion acoustic soliton can propagate.From the mathematical point of view, solitons are the stationary solutions to the Korteweg-de Vries equation: where ( ) = is a function of position x and time t.Some historical events and discoveries led to the soliton theory and Equation ( 1) is not the only equation that has solitonic solution.The sine-Gordon and the nonlinear Schrödinger equations are among them.The solutions to the Equation ( 1) are traveling wave with constant profile in time and they describe the different types of models in various branches of physics and natural science (for instance, Equation ( 1) is modeled to describe the waves on shallow water surfaces).
As mentioned, Equation (1) may arise in the study of such diverse physical systems as fluids and plasmas.In plasma ambient, Equation (1) may be due to the balancing between nonlinearity and dispersion and for this reason its study is of special interest.
One of the possible approaches to study of the nonlinear effects is the so-called reductive perturbation technique [1]- [5].This technique is widely employed to investigate the asymptotic behavior of nonlinear excitations and is more convenient to study of small-amplitude nonlinear perturbations, or to treat of plasma waves in a state very close to thermodynamic equilibrium.
Another successful approach to study of the electrostatic solitons and shock waves has been the Sagdeev Potential (SP) or Pseudo Potential (PP) method.The SP is one particular notion that has become immensely important in soliton and shock research [6].The main advantage of this method over reductive perturbation technique is that, it is appropriate for arbitrary amplitude waves and one can derive all the soliton results of perturbation methods and compare it with the exact results obtained by the SP method [7].
In the present work, our aim is to compute the soliton width (SW) for an ion acoustic solitary wave propagated in a cold plasma.The computation is based on analogy between pseudo-particle in PP well and real particle in conservative potential well.Knowing the angular frequency of (small) oscillation of real particle about its equilibrium position and comparing this with corresponding quantity in the plasma system, one attains a formula to compute the SW.To understand briefly, the SW is defined as spatial length corresponding to a (complete) spatial oscillation of psuedo-particle in the PP well; comparing this with the definition of period of temporal oscillation for real particle in the well, the computation is straightforward.Also, for better understanding of the changes, the graph of the width function is visualized, with the help of numerical plotting (computer algebra).
The work is organized as follows: The model and the main calculation are presented in the next section, and conclusions are given in the last section.

The Soliton Width as a Function of the Mach Number
We will consider a one dimensional model for propagation of ion acoustic wave in a cold plasma.The solitary wave is generated by spatial oscillation of electrostatic potential ( ) which is the so-called psuedo-particle.The wave is travelling to the left in the (to say) x direction with a constant speed 0 u .Therefore, in the wave- frame position, the SP is given by [8]: This PP is subjected to the boundary condition ( ) where e T to be electron temperature, D λ is the ion Debye length and M is the Mach number.
The form of the PP would determine whether soliton like solutions may exist or not.The conditions for the existence of solitary waves are: where m χ is the maximum value of χ beyond which the PP becomes imaginary1 and it is often called am- plitude of the soliton (or shock).By virtue of above conditions, it is easily to show that the values of the Mach number are confined to the interval It may be useful to plot the graph of the PP (2) versus its argument χ .The graph is illustrated in Figure 1 for three values of M. As the figure shows, the depth of the potential well and amplitude of the soliton increase with increasing Mach number.
The PP (2) satisfies the energy condition: ( ) which is analogous to the principle of energy conservation ( ) for a real particle of mass m moving in a (conservative) potential ( ) U x .In view of this analogy, Equation ( 4) can be regarded as an energy integral of moving pseudo-particle of unite mass, pseudo-position χ , pseudo- velocity d d χ ξ and pseudo-potential ( ) V χ .Hence, by the following replacements: 1, , and , one can obtain similar formulas and results in the plasma system, as discussed below.
We know that for a real particle moving in the potential ( ) U x , the angular frequency of small oscillations ω , about the equilibrium position (a minimum at 0 x ) is given by ( ) and as a result, the period of oscillation motion T becomes ( ) The natural question, then, is what conclusions can be drawn from this for our plasma system, more clearly, what quantities are given by similar formulas which are obtained by the above replacements in ( 5) and (6), that is and ( ) where we also use 2π and the two replacement 2 .It is easily to check that s λ has the dimensions of the Length, and in below, we see that s λ is SW.
We first remind the following definitions: T: a time interval corresponding to one complete temporal oscillation of real particle (in the well).
On the other hand the width of the solitary wave is equal to the length that disturbance of the solitary wave takes place, or in SP terminology, SW: a distance interval corresponding to one complete spatial oscillation of pseudo-particle (in the PP well), thus, it is natural to drawing the result To use the above formula, we need to know the point 0 χ , that is minimum of PP.This is given by vanishing of the first derivative of PP (  the smaller one, then, the solitons corresponding to larger M have shorter width. In order to have a quantitative measure of the changes of the width function, let us compute the ratio ( ) ( ) This number is the average decreasing rate of the width per mach number, that is the change in the width divided by the change in Mach number.The negative sign indicates the decreasing nature of the function.As it is evident the amount of decreasing is relatively large.
In the above discussion, the width changes were described in terms of the Mach number itself, but, due to definition of the Mach number , the following easy corollary can be deduced: 1) A positive (negative) change in the soliton wave velocity 0 u leads to a negative (positive) change in the SW, and 2) A positive (negative) change in the ratio ( e T m ) leads to a positive (negative) change in the SW.In the end of this section, it is necessary to say that our ability to plot the width function s λ (as continues function) is because of the analytic form of the SP (2).Indeed, the method presented here may be employed for any SP of analytic form [9]. When the SP has not analytic form, there is no such ability, but we may still calculate numerically s λ at each allowable equilibrium (a minimum) position.This is based on estimation of equili- brium position and substituting into the second derivative of SP.In this case, one has a (discrete) point diagram, instead of a continues curve.One such computation has been represented in [10] for an electron-positron plasma including stationary ions in the background.

Conclusions and Results
The shock and soliton waves occur most likely because of the nonlinear disturbances, namely discontinuities in

2 .
(9), we obtain the final expression for the SW as function of the Mach Number, that is presence of Lambert function, we have no clear idea of the behavior of the width function (13), then it is instructive to graph the function.Because of the transcendental form of the function, we have to use the numerical plotting and the corresponding graph is illustrated in Figure It represents a monotonically decreasing function of M. The function takes on arbitrarily large values when 1 M → , and decreases rapidly over the (allowablementioned, the soliton width is the length of a complete sweep of the pseudo-particle in the potential well and as is clear from the Figure1, the potential well corresponding to the larger values of M is more narrow than

Figure 2 .
Figure 2. The graph of the soliton width λ s versus M: rapidly monotonic decreasing function.