Re: Initial Charge Densities ( No.1 ) 
 Date: 2016/11/08 19:36
 Name: Artem Pulkin
 Hi,
Do not worry: these are only initial occupation numbers. They are only important for magnetic systems where several alternative magnetic configurations are possible.
Now, your particular question is about the pseudopotential approximation. There, we define a frozen core and valence electrons for each atom. The particular Tin pseudopotential you are talking about has 36 frozen electrons and 14 valence. You may increase the balance towards 46 frozen electrons and 4 valence. This will simplify your calculations considerably, however, since only 4 electrons participate in bonding the properties of your solid will be described poorly.
So, in short, valence electrons are not defined strictly. When you look into the periodic table they just correspond to the outermost atomic shell, however, given a pseudopotential approximation one typically includes several more "semicore" electrons into the valence shell.
Keywords: pseudopotential, valence, semicore
Artem

Re: Initial Charge Densities ( No.2 ) 
 Date: 2016/11/09 17:45
 Name: Clarke
 Dear Artem,
Thank you for your reply and giving proper suggestion. I have completed my calculations successfully. I plotted transmission vs energychemical potential of the left lead (Emu_L). It is observed that the Fermi level is not at zero but it is shifted to around 0.75 eV. I have gone through forum posts, it is found that Openmx does not calculate Fermi level instead it calculates the chemical potential. The calculation are done at default value of temperature i.e. at T = 300 K.
So my question is:
1) the obtained results are correct or incorrect, and how to defend it if the Fermi level is not at zero
2) Is there any way to set the Fermi level at zero.
I have taken a look on the results reported in the literature using Openmx that almost everybody plotting transmission vs Emu, and they are showing Fermi level at zero, but I did not understand the setting of Fermi level at zero.
Your help would be highly appreciated.
Thanking you.
Clarke

Re: Initial Charge Densities ( No.3 ) 
 Date: 2016/11/09 19:40
 Name: Artem Pulkin
 1) I do not know. In earlier version of OpenMX everything related to transmissions was relative to the Fermi level. You may always run examples provided with the distribution and check if it is the case.
2) Maybe just subtract its value?

Re: Initial Charge Densities ( No.4 ) 
 Date: 2016/11/10 00:08
 Name: Clarke
 Dear Artem,
Thank you for your reply. Prof. Ozaki mentioned the following:
What OpenMX gives is not Fermi level (EF), but chemical potential(mu) for a given temperature (T). The mu is found in OpenMX so that the total number of electrons can be conserved at T. Also, we have a relation:
lim T>0 mu = EF
at the forum discussion http://www.openmxsquare.org/forum/patio.cgi?mode=view&no=975
Acording to the above discussion,
If the calculations are done at the finite temperature the Fermi level may not be at Zero. In this case can we say that the Fermi level is not as zero because the calculations are done at finite temperature that may shift the Fermi level from zero.
2) May be just subtract its value?
All the researchers are just subtracting its value to match the Fermi level at zero.
Your help would be highly appreciated because all my results depends on this statement. Please help me in this regard.
Thanking you.
Clarke

Re: Initial Charge Densities ( No.5 ) 
 Date: 2016/11/11 13:09
 Name: Clarke
 Dear Artem,
I checked with graphene nanoribbon. The Fermi level is coming near at zero but exactly at zero. My question is that if I am dealing with 'Tin'then why the Fermi level is not at zero.
Please help in this regard.
Thanking you.

Re: Initial Charge Densities ( No.6 ) 
 Date: 2016/11/11 19:16
 Name: Artem Pulkin
 "My question is that if I am dealing with 'Tin'then why the Fermi level is not at zero."
Look, you may set your energy origin at arbitrary value without losing physical sense. The Fermi level is set with respect to THE ORIGIN, thus, generally, the Fermi level is nonzero. When you plot a band structure for the publication you typically want to set your Fermi level to zero since it is the only "special" energy on your plot. So people take the calculated band structure energies and subtract the Fermi level from them setting THE ORIGIN to the Fermi level. Then, they state that the Fermi level on the plot at zero.
For example, when simulating 2D structures, your energy origin is close to the vacuum level. The Fermi level is always below it, thus, you always find it negative (for example, E_F = 5 eV). In 3D structures it is different (there is a reason for that): it can be both positive and negative. I.e. ARBITRARY.
Finally, a few words about the chemical potential. While there is an exact relation between the Fermi level and the chemical potential they are almost equal. As far as I remember the difference between them is less than the order of magnitude of temperature 300K = 25 meV. However, the error in the calculated Fermi level or chemical potential (whatever you call it) is typically larger: your chosen kpoint grid is just too small to capture the 25 meV difference. The chemical potential of band gap materials, for example, should always be in the conduction energy region. However one always finds it at some arbitrary levels inside the band gap. So while it is the chemical potential by concept, it smells, tastes and looks like the Fermi level.
Hope I was clear enough. Whether you have more questions or not I suggest you to follow some solid state physics courses.
Artem

Re: Initial Charge Densities ( No.7 ) 
 Date: 2016/11/14 00:19
 Name: Clarke
 Dear Artem,
Thanking you for your reply.
First two paragraphs are OK to me no problem. My main concern is
>>However, the error in the calculated Fermi level or chemical potential (whatever you call it) is typically larger: your chosen kpoint grid is just too small to capture the 25 meV difference.
In this, you talk about kpoint grid:
Is it: 1) scf kpoint grid of leads
or 2) scf.Kgrid (in NEGF file) or 3) NEGF.scf.Kgrid or 4) NEGF.tran.Kgrid
One more thing about NEGF.gate.voltage (Is it not working when we work with system like left leadcentral regionright lead) because I am dealing with this type of system and applying gate voltage under zero and finite bias voltages, but I am not seeing any difference between the results with and wtithout gate voltages. However, the gate voltage mainly effect the central region, in my case the central region is very thin, may be not more than 810 atoms.
Your suggestions would be highly appreciated.
Thanking you.

Re: Initial Charge Densities ( No.8 ) 
 Date: 2016/11/14 21:44
 Name: Artem Pulkin
 "In this, you talk about kpoint grid:"
The Fermi level is calculated prior to any Greens function calculation so you should be concerned about the first option: scf.kgrid of your leads.
Regarding the rest: if you want to have answers please ask questions.

Re: Initial Charge Densities ( No.9 ) 
 Date: 2016/11/15 15:29
 Name: Clarke
 Dear Artem,
Thank you very much for your help. Now I am doing some calculation on graphene nanoribbon with spinorbit coupling but the results are not differnt looks same. The results are not affected by SOC because it is weak in graphene, is it the case. For spin orbit coupling I used the following input file
# # File Name #
System.CurrrentDirectory ./ # default=./ System.Name leadl6zgnrsoc level.of.stdout 1 # default=1 (13) level.of.fileout 1 # default=1 (02)
# # Definition of Atomic Species #
Species.Number 2 <Definition.of.Atomic.Species C C5.0s2p1 C_CA13 H H5.0s2 H_CA13 Definition.of.Atomic.Species>
NEGF.output_hks true NEGF.filename.hks leadl6zgnrsoc.hks
# # Atoms #
Atoms.Number 28 Atoms.SpeciesAndCoordinates.Unit Ang # AngAU <Atoms.SpeciesAndCoordinates 1 C 0.000000 0.000000 0.000000 2.500000 1.500000 2 C 0.000000 2.836000 0.000000 2.500000 1.500000 3 C 2.456048 0.000000 0.000000 2.500000 1.500000 4 C 1.228024 0.709000 0.000000 2.500000 1.500000 5 C 1.228024 2.127000 0.000000 2.500000 1.500000 ........ 25 H 0.000000 1.090000 0.000000 0.500000 0.500000 26 H 2.456048 1.090000 0.000000 0.500000 0.500000 27 H 0.000000 12.434000 0.000000 0.500000 0.500000 28 H 2.456048 12.434000 0.000000 0.500000 0.500000 Atoms.SpeciesAndCoordinates> Atoms.UnitVectors.Unit Ang # AngAU <Atoms.UnitVectors 4.902196 0.000000 0.000000 0.000000 23.524000 0.000000 0.000000 0.000000 10.000000 Atoms.UnitVectors>
# # SCF or Electronic System #
scf.XcType LSDACA # LDALSDACALSDAPWGGAPBE scf.SpinPolarization NC # OnOffNC scf.SpinOrbit.Coupling on # OnOff scf.ElectronicTemperature 300.0 # default=300 (K) scf.energycutoff 200.0 # default=150 (Ry) #scf.Ngrid 140 140 100 # about=200 (Ry) scf.maxIter 1000 # default=40 scf.EigenvalueSolver Band # DCGDCClusterBand scf.lapack.dste dstevx # dstegrdstedcdstevx, default=dstegr scf.Kgrid 201 1 1 # means n1 x n2 x n3 scf.Mixing.Type rmmdiisk # SimpleRmmDiisGrPulayKerkerRmmDiisk scf.Kerker.factor 10.0 # default=1 scf.Init.Mixing.Weight 0.0001 # default=0.30 scf.Min.Mixing.Weight 0.0001 # default=0.001 scf.Max.Mixing.Weight 0.105 # default=0.40 scf.Mixing.History 150 # default=5 scf.Mixing.StartPulay 10 # default=6 scf.Mixing.EveryPulay 1 # default=6 scf.criterion 1.0e6 # default=1.0e6 (Hartree)
# # MD or Geometry Optimization #
