Self-energies are calculated from a pristine, bulk calculation of only the electrode.
This is a requirement of the self-energy algorithm. This also forces the boundary conditions of the electrodes.
Since, the calculation is: $$ \mathbf G^{-1} = \begin{bmatrix} \mathbf S_1 e - \mathbf H_1 - \Sigma_1 & \cdots \\ \vdots \\ \end{bmatrix} $$ The $\Sigma_1$ is the self-energy which is derived from the bulk electrode named $1$. However, the above requires that the electrode in the device region is equivalent to the bulk electrode region.
The initiated reader will notice that one can choose $\mathbf H_1$ between two values:
- the pristine bulk Hamiltonian:
TS.Elec.<>.Bulk true
forces the coupled region in the device region to be very electrode bulk like - the device region Hamiltonion:
TS.Elec.<>.Bulk false
forces the electrode part of the device to couple to a bulk electrode with bulk properties. This can in some cases spare you some layers of electrode in the device region.
It is always the user's responsibility to ensure the boundary conditions are fulfilled! TranSIESTA will/can not check this.