Propagating stationary profile excitations (solitons) are ubiquitous in nature. Arising in areas as diverse as telecommunications, hydrodynamics and Space physics, pulse-shaped soliton forms are modeled via generic nonlinear differential equations. In the physics of charged matter (plasmas), electrostatic solitons are modelled by means of two theoretical approaches, namely the Sagdeev pseudopotential (SPP) and the reductive perturbation (RPM) method. Both methods come with serious restrictions in terms of the parameter values, which makes their use qualitatively attractive, but rather dubious in terms of quantitative predictions. This project focuses on a novel approach in modeling nonlinear plasma waves, based on a hybrid simulation method combining the SPP method with an original kinetic-Vlasov type simulation technique. The results of the simulation can be used to verify a variety of theoretical suggestions on modelling multi-species plasmas. The focus of this project will be on new types of solitons suggested recently, such as super-solitons and the possibility of existence of negative solitons. Moreover, the study will examine the properties of solitons during mutual collisions in order to establish their robustness.