A path towards modeling 10–100 GeV plasma accelerator stages is described. Modeling such stages including the self‐consistent evolution of the driver, the plasma dynamics, and the beam loading of the trailing particles necessitates the use of particle methods. In a previous proceedings, I reviewed the status of particle based methods and stated that new methods needed to be developed if one hoped to routinely model the full scale in three‐dimensions of 10+ GeV plasma accelerator stages. In this article, I describe the development of a new fully parallelized algorithm that reproduces the results from standard particle‐in‐cell methods with at least a 100 times savings in CPU time. I also describe how standard methods are being used to discover new physics. © 2004 American Institute of Physics