A closed-form algorithm for computing subsurface Green functions-the blocks of a material's Green function between the surface and the bulk-is presented, where we assume the system satisfies a common principal-layer approximation. By exploiting the block tridiagonal and nearly block Toeplitz structure of the Hamiltonian and overlap matrices, this method scales independently of the system size (constant scaling), allowing studies of large systems. As a proof-of-concept example, we investigate the decay of surface effects in an armchair graphene nanoribbon, demonstrating the persistence of surface effects hundreds of atomic layers (~0.5 μm) away from a surface. We finally compare the surface-to-bulk transitions of finite and semi-infinite systems, finding that finite systems exhibit amplified surface effects.
|Journal||Physical Review B - Condensed Matter and Materials Physics|
|Publication status||Published - Feb 15 2011|
ASJC Scopus subject areas
- Condensed Matter Physics
- Electronic, Optical and Magnetic Materials