Numerical issues of the two-dimensional Dirac equation

The two-dimensional Dirac equation has been widely used in graphene physics, the surface of topological insulators, and especially quantum scarring. Although a numerical approach to tackling an arbitrary confining problem was proposed several years ago, several fundamental issues must be thoroughly understood and solved. In this work, we conceal and address these challenges and finally develop a complete method, validated by comparison with analytical results.