We present particle-in-cell simulations, using the VORPAL framework, of 10 GeV laser plasma wakefield accelerator stages. Scaling of the physical parameters with the plasma density allows us to perform these simulations at reasonable cost and to design high performance stages. In particular we show that, by choosing to operate in the quasi-linear regime, we can use higher order laser modes to tailor the focusing forces. This makes it possible to increase the matched electron beam radius and hence the total charge in the bunch while preserving the low bunch emittance required for applications.