Exact solutions for topological surface states of three-dimensional lattice models

TL;DR

This paper introduces a generalized transfer matrix method that provides exact solutions for topological surface states in 3D lattice models, improving accuracy and efficiency over previous approaches.

Contribution

The authors develop a transfer matrix formalism that avoids singular matrices and applies it to derive exact surface states in topological models, surpassing prior numerical methods.

Findings

Exact solutions for surface states in 3D topological models

Method accurately reproduces numerical diagonalization results

Applicable to models beyond low-energy approximations

Abstract

In this work, we establish a generalized transfer matrix method that provides exact analytical and numerical solutions for lattice versions of topological models with surface termination in one direction. We construct a generalized eigenvalue equation, equivalent to the conventional transfer matrix, which neither suffers from nor requires singular (non-invertible) inter-layer hopping matrices, in contrast to previous works. We then apply this formalism to derive, with exactness, the topological surface states and Fermi arc states in two prototypical topological models: the 3D Bernevig-Hughes-Zhang model and a lattice model exhibiting Weyl semimetal behavior. Our results show that the surface states and bulk bands, across the projected 2D Brillouin zone, agree perfectly with those obtained through direct numerical diagonalization of the corresponding Hamiltonians in a slab geometry. This highlights that the generalized transfer matrix method is not only a powerful tool but also a highly efficient alternative to fully numerical methods for investigating surface physics and interfaces in topological systems, particularly when it is required to go beyond low-energy effective descriptions.