Abstract:An efficient implementation of the nonequilibrium Green function (NEGF) method combined with the density functional theory (DFT) using localized pseudo-atomic orbitals (PAOs) is presented for electronic transport calculations of a system connected with two leads under a finite bias voltage. In the implementation, accurate and efficient methods are developed especially for evaluation of the density matrix and treatment of boundaries between the scattering region and the leads. Equilibrium and nonequilibrium contributions in the density matrix are evaluated with very high precision by a contour integration with a continued fraction representation of the Fermi-Dirac function and by a simple quadratureon the real axis with a small imaginary part, respectively. The Hartree potential is computed efficiently by a combination of the two dimensional fast Fourier transform (FFT) and a finite difference method, and the charge density near the boundaries is constructed with a careful treatment to avoid the spurious scattering at the boundaries. The efficiency of the implementation is demonstrated by rapid convergence properties of the density matrix. In addition, as an illustration, our method is applied for zigzag graphene nanoribbons, a Fe/MgO/Fe tunneling junction, and a LaMnO$_3/$SrMnO$_3$ superlattice, demonstrating its applicability to a wide variety of systems.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is how to achieve efficient and accurate calculations by combining the non - equilibrium Green's function (NEGF) method with density - functional theory (DFT) to study the electron - transport properties of a two - electrode - connected system under a finite bias voltage. Specifically, the paper focuses on the following two key issues:
1. **Expanding the application range of the NEGF method**: In order to apply the NEGF method to larger - scale and more complex systems, an efficient and accurate implementation method must be developed. The main challenges faced by existing NEGF methods in dealing with large - scale systems lie in the efficiency and accuracy issues when calculating the density matrix. Therefore, this paper proposes a new efficient method for evaluating the density matrix, which can better handle large - scale systems.
2. **Avoiding spurious scattering**: When the NEGF method is extended to include many - body effects (such as electron - phonon interactions, electron - electron interactions, etc.), spurious scattering may obscure these many - body effects. Spurious scattering usually occurs when the boundary conditions are not properly handled, especially at the boundary between the source and drain regions. Therefore, this paper also pays special attention to how to avoid spurious scattering by accurately handling the boundary conditions to ensure the accuracy of the calculation results.
To solve these problems, the author proposes an efficient implementation scheme based on the contour - integral method and the continued - fraction representation of the Fermi - Dirac function, and combines two - dimensional fast Fourier transform (FFT) and the finite - difference method to calculate the Hartree potential, ensuring the correct handling of the boundary conditions. In addition, the author also discusses in detail how to construct the charge density near the boundary and shows the application effects of this method in various material systems, including zigzag graphene nanoribbons, Fe/MgO/Fe tunnel junctions, and LaMnO₃/SrMnO₃ superlattices, etc.
### Summary of Key Formulas
1. **Expression of the Green's function**:
\[
G(k)_{\sigma,C}(Z) = \left(ZS(k)_C - H(k)_{\sigma,C} - \Sigma(k)_{\sigma,L}(Z) - \Sigma(k)_{\sigma,R}(Z)\right)^{-1}
\]
where $\Sigma(k)_{\sigma,L}(Z)$ and $\Sigma(k)_{\sigma,R}(Z)$ are self - energy terms.
2. **Iterative calculation of the surface Green's function**:
\[
Gs,11 = \lim_{i \to \infty} \epsilon^{-1}_{s,i}
\]
3. **Calculation of the equilibrium - state density matrix**:
\[
\rho^{(eq)}_{\sigma,Rn} = \frac{1}{V_c} \int_{BZ} dk^3 \left(\rho(k)_{\sigma,+} - \rho(k)_{\sigma,-}\right) e^{-ik \cdot R_n}
\]
where
\[
\rho(k)_{\sigma,\pm} = \frac{i}{2\pi} \int_{-\infty}^{\infty} dE G(k)_{\sigma,C}(E \pm i0^+) f(E - \mu)
\]
4. **Continued - fraction representation of the Fermi - Dirac function**:
\[
\frac{1}{1 + \exp(x)} = \frac{1}{2} + \sum_{p = 1}^{\infty} \left( \frac{R_p}{x - iz_p} + \frac{R_p}{x + iz_p} \right)
\]
where $x=\beta(z - \mu)$, $\beta = \frac{1}{k_B T}$.
5. **Solution of the Poisson equation**:
\[
\nabla^2 \Delta V_H(r) = - 4\pi \Delta n(r)
\]
where $\Delta n(r) =$