negative occupancies with tetrahedron method
Moderators: Global Moderator, Moderator
- eariel
- Newbie
- Posts: 16
- Joined: Sat Sep 13, 2008 6:45 pm
- License Nr.: 490
negative occupancies with tetrahedron method
I find several bands with occupancies different from 0 and 1 while using ISMEAR = -5. I have values smaller than zero and larger than 1. I have always understood that with ISMEAR = -5 one has occupantions 0 or (or 2 for non-spin-polarized calcualtions).
I read the articule by Bl?chl et al (PRB 49, 16223 (1994)) and there the occupations are 0 and 1.
In a different calculation, after obtaing the SCF state using ISMEAR=2, I ran a calculation using ALGO=None with ISMEAR = -5. Then I see the following line in the OUTCAR
None: do nothing, only one-electron occupancies are recalculated
and again I got strange occupations
grep ' 484 ' OUTCAR-none
484 3.1231 1.00000
484 3.6826 0.79753
484 3.6823 0.79795
484 4.3278 -0.30484
484 3.6901 0.78905
484 3.2231 1.00000
484 3.8080 0.60281
484 3.8083 0.60224
484 4.4551 -0.39050
484 3.8063 0.60579
484 3.1231 1.00000
484 3.6826 0.79753
484 3.6823 0.79794
484 4.3278 -0.30485
484 3.6901 0.78904
484 3.2231 1.00000
484 3.8080 0.60280
484 3.8083 0.60223
484 4.4551 -0.39051
484 3.8063 0.60579
Is it a bug? I am using VASP 5.2.12
Regards, Eduardo
I read the articule by Bl?chl et al (PRB 49, 16223 (1994)) and there the occupations are 0 and 1.
In a different calculation, after obtaing the SCF state using ISMEAR=2, I ran a calculation using ALGO=None with ISMEAR = -5. Then I see the following line in the OUTCAR
None: do nothing, only one-electron occupancies are recalculated
and again I got strange occupations
grep ' 484 ' OUTCAR-none
484 3.1231 1.00000
484 3.6826 0.79753
484 3.6823 0.79795
484 4.3278 -0.30484
484 3.6901 0.78905
484 3.2231 1.00000
484 3.8080 0.60281
484 3.8083 0.60224
484 4.4551 -0.39050
484 3.8063 0.60579
484 3.1231 1.00000
484 3.6826 0.79753
484 3.6823 0.79794
484 4.3278 -0.30485
484 3.6901 0.78904
484 3.2231 1.00000
484 3.8080 0.60280
484 3.8083 0.60223
484 4.4551 -0.39051
484 3.8063 0.60579
Is it a bug? I am using VASP 5.2.12
Regards, Eduardo
Last edited by eariel on Tue Feb 05, 2013 12:31 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
negative occupancies with tetrahedron method
This is not a bug. Bloechls correction to the tetrahedron method can lead to occupancies >2 or <0.
Last edited by admin on Tue Feb 05, 2013 2:56 pm, edited 1 time in total.
- eariel
- Newbie
- Posts: 16
- Joined: Sat Sep 13, 2008 6:45 pm
- License Nr.: 490
negative occupancies with tetrahedron method
Thank you. Where can I read the method to obtain the occupancies? In Bloechl's paper the occupancies are fixed.
Last edited by eariel on Tue Feb 05, 2013 5:42 pm, edited 1 time in total.
- eariel
- Newbie
- Posts: 16
- Joined: Sat Sep 13, 2008 6:45 pm
- License Nr.: 490
negative occupancies with tetrahedron method
So, having no answer, I wonder if the optimization of the occupancies is an unpublished recipe ? Is there any keyword in VASP to enforce the occupancies to their physical values at 0 K ?
Last edited by eariel on Thu Feb 07, 2013 1:46 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
negative occupancies with tetrahedron method
Modified tetrahedron method was published by Bloechl in PRB 49(1994)16223. The method itself allow only slight deviation of occupancies from phasical values. The values like 0.5 are caused by symmetry of energy levels.
Try another smearing, e.g. Gaussian (ISMEAR=0) and tune the SIGMA to small values until you will get occupancies only 0 and 1.
Try another smearing, e.g. Gaussian (ISMEAR=0) and tune the SIGMA to small values until you will get occupancies only 0 and 1.
Last edited by admin on Thu Feb 07, 2013 2:24 pm, edited 1 time in total.
- eariel
- Newbie
- Posts: 16
- Joined: Sat Sep 13, 2008 6:45 pm
- License Nr.: 490
negative occupancies with tetrahedron method
I began this thread mentioning that I had read Bloechl in PRB 49(1994)16223 (see my first post), and I found no mention to occupations within the tetrahedron method. I see a discussion around Eq. (27), on page 16231, but I understand that this concerns smearing methods and their difference with the tetrahedron method. Yes, I understand the other smearing methods, and the way to tune SIGMA. I also read in the manual that the tetrahedron method is the most precise method for total energies, even in metals if the forces are not need. There is mention to unphysical occupations if the Methfessel-Paxton method is used with non-metals, but no mention to unphysical occupations (i.e., negative) with the modified tetrahedron method.
Last edited by eariel on Thu Feb 07, 2013 2:49 pm, edited 1 time in total.
- eariel
- Newbie
- Posts: 16
- Joined: Sat Sep 13, 2008 6:45 pm
- License Nr.: 490
negative occupancies with tetrahedron method
In Bloechl's paper it is mentioned a possibility to make the method variational. It is juat a possibility mentioned in Ref. 21. Is this comment related with the optimization of occupation numbres in VASP?
Last edited by eariel on Thu Feb 07, 2013 2:54 pm, edited 1 time in total.