2009/09/24 by Denis A. Areshkin, Branislav K. Nikolic, Branislav K. Nikolić · 2 citations
Engineering · Materials Science · Mathematics · Physics and Astronomy · #Algorithm #Complex plane #Condensed matter physics #Density functional theory #Electron #Electron density #Fermi level #Function (biology) #Graphene #Graphene research and applications #Mathematical analysis #Mathematics #Molecular Junctions and Nanostructures #Non-equilibrium thermodynamics #Physics #Quantum and electron transport phenomena #Quantum mechanics #cond-mat.mes-hall #cond-mat.stat-mech
paper · pdf · doi:10.1103/physrevb.81.155450
published as Phys. Rev. B 81, 155450 (2010). · 19 pages, 8 PDF figures, PDFLaTeX
arxiv created 2009/09/24 · openalex publication_date 2010/04/26 · arxiv updated 2010/04/28 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
The recent fabrication of graphene nanoribbon (GNR) field-effect transistors poses a challenge for first-principles modeling of carbon nanoelectronics due to many thousand atoms present in the device. The state of the art quantum transport algorithms, based on the nonequilibrium Green function formalism combined with the density-functional theory (NEGF-DFT), were originally developed to calculate self-consistent electron density in equilibrium and at finite bias voltage (as a prerequisite to obtain conductance or current-voltage characteristics, respectively) for small molecules attached to metallic electrodes where only a few hundred atoms are typically simulated. Here we introduce combination of two numerically efficient algorithms which make it possible to extend the NEGF-DFT framework to device simulations involving large number of atoms. Our first algorithm offers an alternative to the usual evaluation of the equilibrium part of electron density via numerical contour integration of the retarded Green function in the upper complex half-plane. It is based on the replacement of the Fermi function f(E) with an analytic function \stackrel\ifmmode \else \~\fif(E) coinciding with f(E) inside the integration range along the real axis, but decaying exponentially in the upper complex half-plane. Although \stackrel\ifmmode \else \~\fif(E) has infinite number of poles, whose positions and residues are determined analytically, only a finite number of those poles have non-negligible residues. We also discuss how this algorithm can be extended to compute the nonequilibrium contribution to electron density, thereby evading cumbersome real-axis integration (within the bias voltage window) of NEGFs which is very difficult to converge for systems with large number of atoms while maintaining current conservation. Our second algorithm combines the recursive formulas with the geometrical partitioning of an arbitrary multiterminal device into nonuniform segments in order to reduce the computational complexity of the retarded Green function evaluation by extracting only its submatrices required for electron density and transmission function. We illustrate fusion of these two algorithms into the NEGF-DFT-type code by computing charge transfer, charge redistribution and conductance in zigzag-GNR\ensuremath|variable-width-armchair-GNR\ensuremath|zigzag-GNR two-terminal device covered with a gate electrode made of graphene layer as well. The total number of carbon and edge-passivating hydrogen atoms within the simulated central region of this device is \ensuremath≃7000. Our self-consistent modeling of the gate voltage effect suggests that rather large gate voltage \ensuremath≃3 eV might be required to shift the band gap of the proposed AGNR interconnect and switch the transport from insulating into the regime of a single open conducting channel.