QuantumATK Forum

QuantumATK => General Questions and Answers => Topic started by: kaypu on October 7, 2012, 05:13

Title: Question about self-energy for scf and TE
Post by: kaypu on October 7, 2012, 05:13
Dear QuantumWise:
    i've got some problem about my system, there are some negetive Transmissionspectrum(TS) near Ef. i've found some way to deal with the problem in the forum. So, i change self-energy calculator to Direct for TS. but for SCF, i continue to use krylov for real, i got TS (named A). finally, the TS is unchanged, why?

If i use Direct for SCF and TS, and got TS(named B), is there some different between A and B?
Title: Re: Question about self-energy for scf and TE
Post by: Anders Blom on October 8, 2012, 21:27
The Direct (and Recursion) method is very stable, but a bit slow, whereas Krylov can give instabilities, but it's very fast. If no numerical issues occur, they should all give the same numbers in the end (within some precision), but as you noted if the Krylov method gives negative transmission, one should switch.

It should be said that we have improved the Krylov method in 12.2 (and some additional fixes in 12.8.x) such that if negative results are discovered, the engine switches internally to Direct by itself, in order not to return unphysical values, and also the Krylov method itself has been made more stable. Thus we hope the problems with negative transmission should be more or less eliminated in the newer versions of ATK.

To return to the question - it's a bit unchartered territory. If you want to be on the absolutely safe side, you should redo the scf with Recursion. Pragmatically, however, as long as you get reasonable transmission with the mix of methods, I think you are quite ok. The negative transmission only occur for specific energies, after all, and you are somewhat unlikely to hit those on the contour. But don't hold me to that answer - it's your own responsibility to check :)
Title: Re: Question about self-energy for scf and TE
Post by: kaypu on October 9, 2012, 02:09
 i  recalculate the scf with Direct, the problem is gone. thank you  very much professor Anders.