Graphene modeling by TLM approach

We introduce for the first time a condensed node scheme for solving the Dirac equation in 2D graphene. This scheme satisfies the standard charge conservation requirement and allows adopting boundary conditions for graphene circuits. The correlation between the graphene equations and its self-consistent symmetrical condensed node -transmission line matrix formulation is highlighted. This concept, in turn, is related to the generalized Huygens principle for the Dirac equations.