Dirac fermions on a two-dimensional lattice with disorder are considered. The Dirac mass, which controls the gap between the two bands of the fermions, is subject to random fluctuations. Another type of disorder is discussed presented by a random vector potential. It is shown that the imaginary part of the one-particle Green's function can be written as the imaginary part of another Green's function, which has only poles on the lower half plane. Therefore, it is possible to perform a Cauchy integration for a Lorentzian distribution in analogy with the Lloyd model. The results are compared with calculations performed in the continuum limit based on renormalization group and bosonization methods. \textcopyright{} 1996 The American Physical Society.
We investigate the simulations of the the Schrödinger equation using the onedimensional quantum lattice Boltzmann (QLB) scheme and the irregular behaviour of solution. We isolate error due to approximation of the Schrödinger solution with the non-relativistic limit of the Dirac equation and numerical error in solving the Dirac equation. Detailed analysis of the original scheme showed it to be first order accurate. By discretizing the Dirac equation consistently on both sides we derive a second order accurate QLB scheme with the same evolution algorithm as the original and requiring only a one-time unitary transformation of the initial conditions and final output. We show that initializing the scheme in a way that is consistent with the non-relativistic limit supresses the oscillations around the Schrödinger solution. However, we find the QLB scheme better suited to simulation of relativistic quantum systems governed by the Dirac equation and apply it to the Klein paradox. We reproduce the quantum tunnelling results of previous research and show second order convergence to the theoretical wave packet transmission probability. After identifying and correcting the error in the multidimensional extension of the original QLB scheme that produced asymmetric solutions, we expand our second order QLB scheme to multiple dimensions. Next we use the QLB scheme to simulate Klein tunnelling of massless charge carriers in graphene, compare with theoretical solutions and study the dependence of charge transmission on the incidence angle, wave packet and potential barrier shape. To do this we derive a representation of the Dirac-like equation governing charge carriers in graphene for the one-dimensional QLB scheme, and derive a two-dimensional second order graphene QLB scheme for more accurate simulation of wave packets. We demonstrate charge confinement in a graphene device using a configuration of multiple smooth potential barriers, thereby achieving a high ratio of on/off current with potential application in graphene field effect transistors for logic devices. To allow simulation in magnetic or pseudo-magnetic fields created by deformation of graphene, we expand the scheme to include vector potentials. In addition, we derive QLB schemes for bilayer graphene and the non-linear Dirac equation governing Bose-Einstein condensates in hexagonal optical lattices.
Two-dimensional (2D) disordered superconductor (SC) in class D exhibits a disorder-induced quantum multicritical phenomenon among diffusive thermal metal (DTM), topological superconductor (TS), and conventional localized (AI) phases. To characterize the quantum tricritical point where these three phases meet, we carry out a two-loop renormalization group (RG) analysis for 2D Dirac fermion with random mass in terms of the $\ensuremath{\epsilon}$-expansion in the spatial dimension $d=2\ensuremath{-}\ensuremath{\epsilon}$. In two dimensions $(\ensuremath{\epsilon}=0)$, the random mass is marginally irrelevant around a clean-limit fixed point of the gapless Dirac fermion, while there exists an IR unstable fixed point at finite disorder strength that corresponds to the tricritical point. The critical exponent, dynamical exponent, and scaling dimension of the (uniform) mass term are evaluated around the tricritical point by the two-loop RG analysis. Using a mapping between an effective theory for the 2D random-mass Dirac fermion and the
(1+1)-dimensional Gross-Neveu model, we further deduce the four-loop evaluation of the critical exponent, and the scaling dimension of the uniform mass around the tricritical point. Both the two-loop and four-loop results suggest that criticalities of a AI-DTM transition line as well as TS-DTM transition line are controlled by other saddle-point fixed point(s) at finite uniform mass.