Abstract

Simulations of quantum transport in coherent conductors have evolved into mature techniques that are used in fields of physics ranging from electrical engineering to quantum nanoelectronics and material science. The most efficient general-purpose algorithms have a computational cost that scales as $L^{6 \dots 7}$ in 3D, which on the one hand is a substantial improvement over older algorithms, but on the other hand still severely restricts the size of the simulation domain, limiting the usefulness of simulations through strong finite-size effects. Here, we present a novel class of algorithms that, for certain systems, allows to directly access the thermodynamic limit. Our approach, based on the Green's function formalism for discrete models, targets systems which are mostly invariant by translation, i.e. invariant by translation up to a finite number of orbitals and/or quasi-1D electrodes and/or the presence of edges or surfaces. Our approach is based on an automatic calculation of the poles and residues of series expansions of the Green's function in momentum space. This expansion allows to integrate analytically in one momentum variable. We illustrate our algorithms with several applications: devices with graphene electrodes that consist of half an infinite sheet; Friedel oscillation calculations of infinite 2D systems in presence of an impurity; quantum spin Hall physics in presence of an edge; the surface of a Weyl semi-metal in presence of impurities and electrodes connected to the surface. In this last example, we study the conduction through the Fermi arcs of the topological material and its resilience to the presence of disorder. Our approach provides a practical route for simulating 3D bulk systems or surfaces as well as other setups that have so far remained elusive.

Highlights

  • Quantum nanoelectronics is changing from a domain of fundamental research into one of the main platforms for development of quantum technologies

  • Translationally invariant systems allow us to describe a great variety of quantum systems of high interest that are beyond the scope of traditional simulation techniques

  • Our conductance calculation provides direct evidence for the role in transport of Fermi arcs which are present at the surface of Weyl semimetals, a topic that has attracted a great deal of attention recently [26,27,28]

Read more

Summary

INTRODUCTION

Quantum nanoelectronics is changing from a domain of fundamental research into one of the main platforms for development of quantum technologies. There are many situations where one must simulate large three-dimensional (3D) systems It is the case in particular for any realistic geometry, in the presence of different characteristic length scales (such as Fermi wavelength and mean level spacing), or in the case of topological materials such as 3D topological insulators or Weyl semimetals [9] where well-separated surfaces are important to avoid surface states mixing. The need to simulate 3D quantum systems has provoked the development of new methods that scale linearly with the system size These techniques include the Kernel polynomial expansion [10,11] and variants of the Lanczos [12] method such as the Lanczos recursion method [13,14]. The main technical challenge lies in performing a momentum integral of matrices that contain Dirac and principal part distributions

MOSTLY TRANSLATIONALLY INVARIANT SYSTEMS
Problem formulation
Principle of the technique
APPLICATIONS
Friedel oscillations in a two-dimensional electron gas
Quantum spin Hall effect
Graphene nanoribbon
Weyl semimetal three-terminal device
RESIDUE PROBLEM
Formulation of the residue problem
Solution to the residue problem
NUMERICAL FOURIER TRANSFORM
GLUING SEQUENCE
BOUND-STATE PROBLEM
Findings
VIII. CONCLUSION
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call