  Three dimensional, particle-in-cell, fully electromagnetic simulations of
electron plasma wake field acceleration in the blow out regime are presented.
Earlier results are extended by (i) studying the effect of longitudinal density
gradient; (ii) avoiding use of co-moving simulation box; (iii) inclusion of ion
motion; and (iv) studying fully electromagnetic plasma wake fields. It is
established that injecting driving and trailing electron bunches into a
positive density gradient of ten-fold increasing density over 10 cm long
Lithium vapor plasma, results in spatially more compact and three times larger,
compared to the uniform density case, electric fields ($-6.4 \times 10^{10}$
V/m), leading to acceleration of the trailing bunch up to 24.4 GeV (starting
from initial 20.4 GeV), with an energy transfer efficiencies from leading to
trailing bunch of 75 percent. In the uniform density case $-2.5 \times 10^{10}$
V/m wake is created leading to acceleration of the trailing bunch up to 22.4
GeV, with an energy transfer efficiencies of 65 percent. It is also established
that injecting the electron bunches into a negative density gradient of
ten-fold decreasing density over 10 cm long plasma, results in spatially more
spread and two-and-half smaller electric fields ($-1.0 \times 10^{10}$ V/m),
leading to a weaker acceleration of the trailing bunch up to 21.4 GeV, with an
energy transfer efficiencies of 45 percent. Inclusion of ion motions into
consideration shows that in the plasma wake ion number density can increase
over few times the background value. It is also shown that transverse
electromagnetic fields in plasma wake are of the same order as the longitudinal
(electrostatic) ones.
