No hybrid functional approximation is implemented at the current version. The computational cost for the hybrid functional calculation of a bulk or molecular system is quite expensive. Since the number of atoms in a two-probe device system is much larger than a bulk or molecule system, it becomes quite challenging or un-feasible to perform the hybrid functional calculation for a two-probe device system. It may be more desirable to try the low-cost but moderately accurate functional, e.g., LDA+U, self-interaction correction, and so on.