Example: Benzenedithiol junction: Wide-Band-Limit

Download green_WBL.run

In the wide-band limit (WBL) the coupling to the leads is assumed to be independent of energy. Therefore one does not need to calculate any self-energies. This also means that the eigenspace of the Green’s function is independent of energy. It can therefore be diagonalized in advance, greatly speeding up the calculation of the DOS and the transmission.

In the example $ADFHOME/examples/adf/green_Al/green_WBL/green_WBL.run of green, the transmission of benzenedithiol junction in the wide-band limit (WBL) is calculated. In order to model the molecule-metal interface, we do need to include a few gold layers in the calculation. However, unlike before, only a single atomic layer as the principal layer is used.

Compare green_BDT.run

FRAGMENTS
    left        layer.t21
    molecule    molecule.t21
    right       layer.t21
END

with green_WBL/green_WBL.run, where each layer of gold atoms has its own fragment:

FRAGMENTS
    left        layer.t21
    left.2      layer.t21
    left.3      layer.t21
    molecule    molecule.t21
    right.3     layer.t21
    right.2     layer.t21
    right       layer.t21
END

Because a single atomic layer is an unnatural configuration for gold, a minor amount of smearing is necessary to make the calculation converge. The molecule is sandwiched in between the electrodes just like before (see Fig. 2 in the example for benzenedithiol). However, this time each atomic layer of gold gets its own fragment. The reason for this configuration is that if the WBL is used on the entire gold contact the result is an an unphysical coupling to the leads; even the gold atoms contacting the molecule would have a direct coupling to the environment. A much better result can be obtained by only using the WBL on the back-most atomic layer and letting the electrons propagate naturally through the rest of the contact. Because the WBL is computationally so inexpensive, we can easily calculate the DOS and transmission for 10,000 points instead of 1000.

$ADFBIN/green << eor
DOS fock.t21
TRANS fock.t21
EPS -0.5 0 10000
ETA 1e-6
LEFT
    FRAGMENT left
    ETA 2e-2
END
RIGHT
    FRAGMENT right
    ETA 2e-2
END
NOSAVE DOS_B, TRANS_B
eor

A comparison of the resulting transmission with the calculation with self-energies is shown in the following figure:

../_images/green_WBL.png

The WBL shows good agreement with the non-WBL transmission around the Fermi energy (-0.195 Hartree or -5.306 eV). Note that the quality of the WBL depends on the choice of the coupling (ETA). For this particular contact geometry we obtain good agreement for ETA = 0.02 Hartree, but a better value may be found for other electrodes. Finally, the WBL can be incrementally improved by adding more gold layers to the extended molecule. For many layers it converges to the calculation with full self-energies.