Self-consistent iteration

The nextnano.NEGF tool is based on the non-equilibrium Green’s functions (NEGF) formalism (also known as the Keldysh, or Kadanoff-Baym formalism). This formalism allows to account for both quantum transport effects (i.e. coherent transport effects such as resonant tunneling), as well as scattering mechanisms.

In the NEGF formalism, scattering processes are described in terms of self-energies. Self-energies and Green’s functions are calculated in a self-consistent way, as both elastic and inelastic scattering processes are accounted within the self-consistent Born approximation (Figure 4.2.1). For the details of the Green’s functions, see Green’s functions.

../../../_images/nnNEGF_model_self_consistent_iteration.gif

Figure 4.2.1 Self-consistent solution algorithm below and above lasing threshold.

At the beginning of the calculation, the one-dimensional Schrödinger equation is solved in real space. The energy levels (minibands in the case of periodic system) are selected by the energy interval specified in the input file (see Axial energy cut-off). The energy levels are transformed into a localized basis of modes which we call Reduced Real Space basis. These wave functions are then used as the basis of the Green’s functions and self-energies.

Secondly, the scattering coupling terms are calculated for each of the accounted mechanism (optical and acoustic phonons, charged impurities, interface roughness, alloy disorder, and electron-electron scattering).

Then, the Green’s functions and self-energies are calculated iteratively. Simultaneously, the mean-field electrostatic potential is calculated self-consistently via Poisson’s equation. Such iterations continue until Convergence is reached for the Green’s functions as well as for the calculated current.

After the convergence, the carrier and current densities are calculated from the Green’s functions. To gain more insights into analysis, the populations, density matrix, oscillator strengths are output in different bases.

The Optical gain/absorption spectrum is calculated if requested in the input file (see Gain{ }). In this case, an additional self-consistent routine can be used to calculate the linear response of the Green’s functions to an a.c. electromagnetic field perturbation.

Convergence

The NEGF self-consistent loop ends successfully if the relative difference between two consecutive iterations has reached ConvergenceValueGF for the lesser Green’s function and ConvergenceValueCurrent for the current density. If the convergence is not achieved, the iteration is terminated after NMaxIterations.

Note

Small values for the convergence values, together with sufficiently large NMaxIterations, will give the most accurate results.

Note

At zero bias, when the current asymptotically approaches 0, the current convergence factor does not converge to zero. Is this ok?

When the current approaches 0, this behavior is normal. In this case, the convergence should be checked by the other convergence factors based on the lesser Green’s function.


Last update: 30/10/2024