1D and 2D electrostatic solitary waves (ESW) were observed by spacecraft in Earth’s plasma sheet boundary layer and auroral zone. Recently, the observation of magnetic perturbations caused by ESW has been found in Earth’s plasma sheet. Because the time scales of ESW waveforms are related to the dynamics of electrons, we study the nonlinear evolution of electron two-stream instability in one- and two-dimensional systems. We develop one- and two-dimensional electromagnetic particle-in-cell codes and perform simulations for various sets of electron beam drift velocity, initial electron thermal velocity and number of beam electrons. We discuss the effects of different parameters on the evolution of electron solitary waves, and analyze the results of simulation by linear theory. We also compare with the differences between one- and two-dimensional simulation results.