Hi, Great you have tried to fully understand this issue.
Den man. 20. apr. 2020 kl. 22.01 skrev El-abed Haidar < ehai2...@uni.sydney.edu.au>: > Following up on that: In the manual it says that if we use* TBT kgrid > Monkhorst Pack (block):* > This means it would act as the same as * kgrid Monkhorst Pack*, however, > this is used to have a separate k-sampling in the tbtrans utility. If it > does not exist it will use kgrid Monkhorst Pack > But as i look into the output it seems to make it back to the default of 1 > 1 1. > Not sure why? > > El-abed Haidar | Doctor of Philosophy (Science) > Condensed Matter Theory (CMT) Group > | School of Physics > THE UNIVERSITY OF SYDNEY | NSW | 2006 > > ------------------------------ > *From:* El-abed Haidar > *Sent:* Friday, 17 April 2020 2:44 AM > *To:* siesta-l@uam.es <siesta-l@uam.es> > *Subject:* Still confused about K points in siesta, transiesta and > tbtrans in siesta 4.1b4 > > Good afternoon to whom it may concern, > i really want to understand the in depth behind the kpoints in electron > transport. I still feel after going through the siesta mail where I have > asked about it before in > https://www.mail-archive.com/siesta-l@uam.es/msg10916.html leading to > https://www.mail-archive.com/siesta-l@uam.es/msg09386.html . However, the > answers were not clear to be honest as i did not get the complete > understanding out of this. I went back and read more about periodicity > as suggested back here in > https://www.mail-archive.com/siesta-l@uam.es/msg08489.html and as Cubot > mentioned its not about periodicity that is usually defined in solid states.* > Any > resources on k points for transport we can rely on or tutorials? Because* > I am more confused on the k points when defined for electron transport > study. So im trying really and really would appreciate if there is a > tutorial in the future for the latest siesta version to how we use k points > BECAUSE > > So far i understand that we will go through 3 steps > 1-transiesta <electrod.fdf >electrode.out to calculate Hamiltonian and > overlap matrix of electrodes > 2-transiesta <scat.fdf >scattering.out to calculate density matrix of > scattering region ] > 3-tbtran <scat.fdf >.tbtrans.out to calculate current and transmission > When papers define the k points to be 1 1 100 along electron transport i > would assume we use in all siesta versions %block Monkhorst Pack and > through out the 3 steps i assume i am using 1 1 100 in my entire > calculation. *I found out recently that once siesta run is done in step > 2 and starts a transiesta run the k points change suddenly to 1 1 1. > WHY???? It was > mentioned https://answers.launchpad.net/siesta/+question/686603 > <https://answers.launchpad.net/siesta/+question/686603> but not enough > explanation on the reason why. * > > I am not sure why especially that i want my study to be infinite along x > and y ( or lattice vector a and b) and finite in z ( or c vector). > *Now if I am suppose to add a block in each step above mentioned, what > should I do to make sure all 3 steps have 1 1 100? IF NOT AGAIN WHY?* > *I fully understand that i should do convergence tests for k points but i > need to know also if that is always the case. Because i read that 1 1 100 > is more than enough (based on literature and reading around the siesta > mail)* > > *Would anyone be able to give their in sight thoughts as detailed as > possible. Thank you* > I would suggest you to first conceptually figure out what the NEGF code does. As per your 3 steps above 1. transiesta elec.fdf This calculates a _truly_ bulk electrode. Ask your selve this: a) how do you approach the truly bulk charge density in a regular calculation? Which parameter should be *tuned* to find the truly bulk electrode? b) You want the electrode to be added to a device region. In the device calculation you have a particular direction of the electrodes that are "semi-infinite". You want the two systems to be commensurate in terms of calculation parameters. How does this relate to a)? What does this imply for the semi-infinite direction, and for the transverse directions? 2. transiesta scat.fdf Here you calculate the device region where you place your _truly_ bulk electrodes. What does the NEGF implementation state about electrodes? What does an electrode mean in a device region system? And more importantly, if you have 2 electrodes both with the same semi-infinite direction. What does this mean for the k-point sampling along that direction? Why shouldn't it matter whether you use any number of k-points? (once you figure out/understand this it means you understand the NEGF concept ;) ) 3. tbtrans scat.fdf This calculates the spectroscopic quantities, such as DOS. There is no conceptual difference from TBtrans and from a regular PDOS calculation using Siesta. E.g. if you use a k-point sampling [nx=3 ny=3 nz=3] in siesta, would that be enough to get a good PDOS estimate? If not, why? And how do you resolve this? > > El-abed Haidar | Doctor of Philosophy (Science) > Condensed Matter Theory (CMT) Group > | School of Physics > THE UNIVERSITY OF SYDNEY | NSW | 2006 > > > -- > SIESTA is supported by the Spanish Research Agency (AEI) and by the > European H2020 MaX Centre of Excellence (http://www.max-centre.eu/) > -- Kind regards Nick
-- SIESTA is supported by the Spanish Research Agency (AEI) and by the European H2020 MaX Centre of Excellence (http://www.max-centre.eu/)