Modeling Klein tunneling and caustics of electron waves in graphene
We employ the tight-binding propagation method to study Klein tunneling and quantum interference in large graphene systems. With this efficient numerical scheme, we model the propagation of a wave packet through a potential barrier and determine the tunneling probability for different incidence angles. We consider both sharp and smooth potential barriers in n-p-n and n-n' junctions and find good agreement with analytical and semiclassical predictions. When we go outside the Dirac regime, we observe that sharp n-p junctions no longer show Klein tunneling because of intervalley scattering. However, this effect can be suppressed by considering a smooth potential. Klein tunneling holds for potentials changing on the scale much larger than the interatomic distance. When the energies of both the electrons and holes are above the Van Hove singularity, we observe total reflection for both sharp and smooth potential barriers. Furthermore, we consider caustic formation by a two-dimensional Gaussian potential. For sufficiently broad potentials we find a good agreement between the simulated wave density and the classical electron trajectories.