In this work we observe for the first time the evolution of the electronic charge density in a collisional system. The model used is system S 3+ + H over the range of relative ion kinetic energies from [1–10] eV. To this aim, ab initio potential energy surfaces and non-adiabatic couplings together with wavepacket propagations have been used.