I computed (G03/G09) orbital energies of molecules containing heavy atoms such as Ge. Now I'd like to compare the respective orbital energies with VASP calculations but the results seem to be rather different? I used the samePOTCAR-files as one would use for bulk calculations with VASP. What would be the best approach to calculate the orbital energies of a molecule?
Thanks a lot for your help!
Orbital energies of molecules
Moderators: Global Moderator, Moderator
-
...andRobin!
Orbital energies of molecules
Last edited by ...andRobin! on Fri Aug 10, 2012 4:01 am, edited 1 time in total.
-
alex
- Hero Member
- Posts: 586
- Joined: Tue Nov 16, 2004 2:21 pm
- License Nr.: 5-67
- Location: Germany
Orbital energies of molecules
Hi,
it's not that straightforward, because in VASP you calculate a 3D periodic image of you your molecule and hence there is no vacuum for a zero energy level. You would like to calc. the workfunction of the system to gain this value. Check the forum for an appropriate post.
Cheers,
alex
it's not that straightforward, because in VASP you calculate a 3D periodic image of you your molecule and hence there is no vacuum for a zero energy level. You would like to calc. the workfunction of the system to gain this value. Check the forum for an appropriate post.
Cheers,
alex
Last edited by alex on Mon Aug 13, 2012 9:13 am, edited 1 time in total.
-
admin
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
Orbital energies of molecules
I would proceed as follows:
-- in 3D periodic systems, the E=0 level is not defined to within a constant offset Evacuum, this is correct; HOWEVER, you can use the tool vtotav.f (to process the LOCPOT file
written if LVTOT= .True. is set) in order to estimate Evacuum in the molecule calculation (i.e. a unit cell with a sufficiently large vacuum region). The energies of the Kohn-Sham levels should then be referred to that E, or to the workfunction (which is the energy of the Fermi level with respect to Evacuum). Please choose which one is more appropriate in your case.
-- please check the character of a certain Kohn-Sham state by having a look at the orbital (charge) character as written in the file PROCAR (set LORBIT = 11 or 12 in INCAR). As you should use the Gamma point only for the free molecule, there is no band dispersion of the states)
[ Edited Thu Aug 30 2012, 03:53PM ]
-- in 3D periodic systems, the E=0 level is not defined to within a constant offset Evacuum, this is correct; HOWEVER, you can use the tool vtotav.f (to process the LOCPOT file
written if LVTOT= .True. is set) in order to estimate Evacuum in the molecule calculation (i.e. a unit cell with a sufficiently large vacuum region). The energies of the Kohn-Sham levels should then be referred to that E, or to the workfunction (which is the energy of the Fermi level with respect to Evacuum). Please choose which one is more appropriate in your case.
-- please check the character of a certain Kohn-Sham state by having a look at the orbital (charge) character as written in the file PROCAR (set LORBIT = 11 or 12 in INCAR). As you should use the Gamma point only for the free molecule, there is no band dispersion of the states)
[ Edited Thu Aug 30 2012, 03:53PM ]
Last edited by admin on Thu Aug 30, 2012 3:49 pm, edited 1 time in total.