  A Korteweg-de Vries (KdV) equation including the effect of Landau damping is
derived to study the propagation of weakly nonlinear and weakly dispersive ion
acoustic waves in a collisionless unmagnetized plasma consisting of warm
adiabatic ions and two different species of electrons at different
temperatures. The hotter energetic electron species follows the nonthermal
velocity distribution of Cairns et al. [Geophys. Res. Lett. 22, 2709 (1995)]
whereas the cooler electron species obeys the Boltzmann distribution. It is
found that the coefficient of the nonlinear term of this KdV like evolution
equation vanishes along different family of curves in different parameter
planes. In this context, a modified KdV (MKdV) equation including the effect of
Landau damping effectively describes the nonlinear behaviour of ion acoustic
waves. It has also been observed that the coefficients of the nonlinear terms
of the KdV and MKdV like evolution equations including the effect of Landau
damping, are simultaneously equal to zero along a family of curves in the
parameter plane. In this situation, we have derived a further modified KdV
(FMKdV) equation including the effect of Landau damping to describe the
nonlinear behaviour of ion acoustic waves. In fact, different modified KdV like
evolution equations including the effect of Landau damping have been derived to
describe the nonlinear behaviour of ion acoustic waves in different region of
parameter space. The method of Ott & Sudan [Phys. Fluids 12, 2388 (1969)] has
been applied to obtain the solitary wave solution of the evolution equation
having the nonlinear term $(\phi^{(1)})^{r}\frac{\partial \phi^{(1)}}{\partial
\xi}$, where $\phi^{(1)}$ is the first order perturbed electrostatic potential
and $r =1,2,3$. We have found that the amplitude of the solitary wave solution
decreases with time for all $r =1,2,3$.
