Differences

This shows you the differences between two versions of the page.

smeagol [2012/12/04 11:49]
prokop
smeagol [2012/12/04 15:59] (current)
prokop
Line 14: Line 14:
====== Nonequlibrim Khon-Sham  ====== ====== Nonequlibrim Khon-Sham  ======
-Aktualne to vypada ze to konverguje alespon pro H-Chain pri nulovem napeti +Aktualne to konverguje (pri vysoke teplote tempfe=1000) a dava "nejakou" vodivost. Ale vysledny potencial (fftpot.xsf) nevypada vubec smysluplne - je symetricky.
- +
-nefungje pro non-zero VBias. Musis se poresit nejak jinak. Melo by se nastavovat VBias vstupujici do smeaglu, nikoli ho umele pricitat k hamiltonianu jako se to deje v Fir2smeagol_eq v diag_k_NET_eq.f90+
=== To Do === === To Do ===
-    * spocitat  Au_1D_chain, Au_bulk, Graphene BULK, SiBULK 
-    * vykreslit electorstatic potential 
-    * podivat se na Delta, sigma (imaginarni cast).... jsou dve, kde je ktera, zjistit jejich vliv na konvergenci, a na "spike" v proudu  
=== Problemy === === Problemy ===
-    * Au1D nekonverguje pro vyzsi napeti (> 0.2 eV), prodlouzeni screening regionu nepomaha. Normalne KS-SCF taky konverguje blbe, ale pomuze zvysit tempfe=1000. to ale ve smeaglu nejde.+    * Hlavni problem ze potencial pro nenulovy VBias nekonverguje k smysluplnemu reseni (asymetricky prohnute rampe), ale k necemu kde jsou LEADy +/- symetricky posunute oproti scattenring regionu. :-
 +    * Smeagol z nejakeho duvodu nezintegruje spravne hustotu stavu leadu. To se pokousi resit workarround v kspace_NET_KS.f90 v sekci //CORRECTION TO EDGE//, ale nejsem si jisty jetli by to tak melo byt, asi ne.
    * zda se mi ze to nenacita/neuklada sprave denmat.dat. protoze kdyz spustim smeagol znovu, jen s trochu pozmenenym napetim, tak je dost velika sigma=1E-1. => __Podivat se jak se nacita a uklada s denmat.__     * zda se mi ze to nenacita/neuklada sprave denmat.dat. protoze kdyz spustim smeagol znovu, jen s trochu pozmenenym napetim, tak je dost velika sigma=1E-1. => __Podivat se jak se nacita a uklada s denmat.__
-    * Je lepsi pouzivat fermi level ze smeagol-SCF nez z LEADu. S Ef z LEADu to nekonverguje vubec. Mozna by si mel smeagol __urcit fermi-level sam pomoci integrace__?????   
    * FIXME denmat se ve smeaglu nacita/uklada spatne, :-(  kdyz spustim 2x beh programu (1x uz zkonvergovany) je stejne v prvnim kromu sigma=1E-1     * FIXME denmat se ve smeaglu nacita/uklada spatne, :-(  kdyz spustim 2x beh programu (1x uz zkonvergovany) je stejne v prvnim kromu sigma=1E-1
 +
=== Kosmeticke vady === === Kosmeticke vady ===
Line 34: Line 30:
=== FAQ === === FAQ ===
-    * denmat.dat se nacita. Vypise to  "Read denmat.dat". deje se to v initbasics.f90 -> initdenmat.f90 +    * denmat.dat se nacita. Vypise to  "Read denmat.dat". deje se to v initbasics.f90 -> initdenmat.f90. ale funguje to nejak spatne asi, nechapu zatim proc, ale nevypada to ze by uz jednou zkonvergovany vypocet konvergoval vyrazne rychleji.
    * "tempfe" ovlivnuje "Tf" ktera se pouziva v negfk.F a selfenergy.F/CHARGEBIAS_EQ, prenasobi se kB=6.33363D-6 [Ry/K]     * "tempfe" ovlivnuje "Tf" ktera se pouziva v negfk.F a selfenergy.F/CHARGEBIAS_EQ, prenasobi se kB=6.33363D-6 [Ry/K]
-    * denmat se ve smeaglu nacita/uklada spatne +    * Delta resp Sigma (proste imaginarni cast, rozsireni peaku stavu) se nastavuje podle teploty - zavisi na tempfe 
 +    * Je lepsi pouzivat fermi level ze smeagol-SCF nez z LEADu. S Ef z LEADu to nekonverguje vubec. Mozna by si mel smeagol __urcit fermi-level sam pomoci integrace__????? 
 
smeagol.txt · Last modified: 2012/12/04 15:59 by prokop
 
Except where otherwise noted, content on this wiki is licensed under the following license:CC Attribution-Noncommercial-Share Alike 3.0 Unported
Recent changes RSS feed Donate Powered by PHP Valid XHTML 1.0 Valid CSS Driven by DokuWiki