Numerical growth of emittance in simulations of laser-wakefield acceleration

Transverse emittance is a crucial feature of laser-wakefield accelerators, yet accurately reproducing its value in numerical simulations remains challenging. It is shown here that, when the charge of the bunch exceeds a few tens of picocoulombs, particle-in-cell (PIC) simulations erroneously overestimate the emittance. This is mostly due the interaction of spurious Cherenkov radiation with the bunch, which leads to a steady growth of emittance during the simulation. A new computational scheme is proposed, which is free of spurious Cherenkov radiation. It can be easily implemented in existing PIC codes and leads to a substantial reduction of the emittance growth.