KPOINTS related problem with Vasp 5.3.3
Posted: Fri Jan 11, 2013 5:36 pm
I am calculating for a cerium oxide fcc slab system with 108 ions (DFT+U). Its a 3 x 3 cell with 4 layers of O-Ce-O. Here are my POSCAR and INCAR files.
With 7x7x1 Gamma Points, the optimizations do not start and the standard log file just prints upto ..."WAVECAR not read". However, when I reduce the no of KPOINTS to 3x3x1 Gamma Points, The optimizations start but are extremely slow.
Please let me know why should it behave this way ??
POSCAR
---------------------------------------------------------------------------------
4TLCeO2
1.00000
11.478500366 0.000000000 0.000000000
-5.739250183 9.940672915 0.000000000
0.000000000 0.000000000 40.934101105
O Ce
72 36
Selective Dynamics
Direct
0.000000010 0.000000000 0.267120004 T T T
0.000000010 0.666670024 0.267120004 T T T
0.666670024 0.000000000 0.267120004 T T T
0.666670024 0.666670024 0.267120004 T T T
0.333330005 0.000000000 0.267120004 T T T
0.333330005 0.666670024 0.267120004 T T T
0.000000000 0.333330005 0.267120004 T T T
0.333330005 0.333330005 0.267120004 T T T
0.666670024 0.333330005 0.267120004 T T T
0.111110002 0.555559993 0.228960007 T T T
0.777779996 0.555559993 0.228960007 T T T
0.444440007 0.222220004 0.228960007 T T T
0.444440007 0.888890028 0.228960007 T T T
0.111110002 0.222220004 0.228960007 T T T
0.111110002 0.888890028 0.228960007 T T T
0.777779996 0.222220004 0.228960007 T T T
0.777779996 0.888890028 0.228960007 T T T
0.444440007 0.555559993 0.228960007 T T T
0.555559993 0.111110002 0.190799996 T T T
0.555559993 0.777779996 0.190799996 T T T
0.555559993 0.444440007 0.190799996 T T T
0.222220004 0.444440007 0.190799996 T T T
0.888890028 0.444440007 0.190799996 T T T
0.222220004 0.111110002 0.190799996 T T T
0.222220004 0.777779996 0.190799996 T T T
0.888890028 0.111110002 0.190799996 T T T
0.888890028 0.777779996 0.190799996 T T T
0.000000000 0.333330005 0.152640000 T T T
0.666670024 0.333330005 0.152640000 T T T
0.000000010 0.000000010 0.152640000 T T T
0.000000010 0.666670024 0.152640000 T T T
0.666670024 0.000000000 0.152640000 T T T
0.666670024 0.666670024 0.152640000 T T T
0.333330005 0.000000000 0.152640000 T T T
0.333330005 0.666670024 0.152640000 T T T
0.333330005 0.333330005 0.152640000 T T T
0.444440007 0.555559993 0.114480004 T T T
0.444440007 0.222220004 0.114480004 T T T
0.444440007 0.888890028 0.114480004 T T T
0.111110002 0.222220004 0.114480004 T T T
0.111110002 0.888890028 0.114480004 T T T
0.777779996 0.222220004 0.114480004 T T T
0.777779996 0.888890028 0.114480004 T T T
0.111110002 0.555559993 0.114480004 T T T
0.777779996 0.555559993 0.114480004 T T T
0.888890028 0.444440007 0.076320000 T T T
0.222220004 0.444440007 0.076320000 T T T
0.555559993 0.444440007 0.076320000 T T T
0.555559993 0.111110002 0.076320000 T T T
0.555559993 0.777779996 0.076320000 T T T
0.222220004 0.111110002 0.076320000 T T T
0.222220004 0.777779996 0.076320000 T T T
0.888890028 0.111110002 0.076320000 T T T
0.888890028 0.777779996 0.076320000 T T T
0.333330005 0.000000000 0.038160000 T T T
0.333330005 0.666670024 0.038160000 T T T
0.000000010 0.000000100 0.038160000 T T T
0.000000010 0.666670024 0.038160000 T T T
0.666670024 0.000000000 0.038160000 T T T
0.666670024 0.666670024 0.038160000 T T T
0.000000000 0.333330005 0.038160000 T T T
0.666670024 0.333330005 0.038160000 T T T
0.333330005 0.333330005 0.038160000 T T T
0.111110002 0.555559993 0.000000000 F F F
0.777779996 0.555559993 0.000000000 F F F
0.444440007 0.222220004 0.000000000 F F F
0.444440007 0.888890028 0.000000000 F F F
0.444440007 0.555559993 0.000000000 F F F
0.111110002 0.222220004 0.000000000 F F F
0.111110002 0.888890028 0.000000000 F F F
0.777779996 0.222220004 0.000000000 F F F
0.777779996 0.888890028 0.000000000 F F F
0.222220004 0.444440007 0.248040006 T T T
0.888890028 0.444440007 0.248040006 T T T
0.555559993 0.111110002 0.248040006 T T T
0.555559993 0.777779996 0.248040006 T T T
0.555559993 0.444440007 0.248040006 T T T
0.222220004 0.111110002 0.248040006 T T T
0.222220004 0.777779996 0.248040006 T T T
0.888890028 0.111110002 0.248040006 T T T
0.888890028 0.777779996 0.248040006 T T T
0.444440007 0.222220004 0.171719998 T T T
0.444440007 0.888890028 0.171719998 T T T
0.444440007 0.555559993 0.171719998 T T T
0.111110002 0.222220004 0.171719998 T T T
0.111110002 0.888890028 0.171719998 T T T
0.777779996 0.222220004 0.171719998 T T T
0.777779996 0.888890028 0.171719998 T T T
0.111110002 0.555559993 0.171719998 T T T
0.777779996 0.555559993 0.171719998 T T T
0.222220004 0.444440007 0.019080000 T T T
0.888890028 0.444440007 0.019080000 T T T
0.555559993 0.111110002 0.019080000 T T T
0.555559993 0.777779996 0.019080000 T T T
0.888890028 0.111110002 0.019080000 T T T
0.888890028 0.777779996 0.019080000 T T T
0.555559993 0.444440007 0.019080000 T T T
0.222220004 0.111110002 0.019080000 T T T
0.222220004 0.777779996 0.019080000 T T T
0.000000000 0.333330005 0.095399998 T T T
0.666670024 0.333330005 0.095399998 T T T
0.333330005 0.333330005 0.095399998 T T T
0.333330005 0.000000000 0.095399998 T T T
0.333330005 0.666670024 0.095399998 T T T
0.000000010 0.000000010 0.095399998 T T T
0.000000010 0.666670024 0.095399998 T T T
0.666670024 0.000000000 0.095399998 T T T
0.666670024 0.666670024 0.095399998 T T T
---------------------------------------------------------------------------------
INCAR
---------------------------------------------------------------------------------
SYSTEM = Ceria(111)
ENCUT=500
ISIF=2
IBRION=2
NSW=200
ISPIN=2
LWAVE=.FALSE. # change for re-runs of converged calculations
LCHARG=.FALSE.
SYMPREC=1E-8
#LVTOT=.TRUE. # use for getting the local absolute potential
LREAL=AUTO
#LAECHG=.TRUE. # use for bader analysis
ADDGRID=.TRUE.
#IDIPOL=3 # do use once the code really works
#LDIPOL=.TRUE.
#LORBIT=11 # use when DOS is needed
#NBANDS=300 # use if vasp suggests
ISMEAR=1
SIGMA=0.1
LASPH=.TRUE.
LDAU=.TRUE.
LDAUTYPE= 3
LDAUL=-1 3 #for O,Ce
LDAUU=0 5 #for O,Ce
LDAUJ=0 1 #for O,Ce
LMAXMIX=6
---------------------------------------------------------------------------------
With 7x7x1 Gamma Points, the optimizations do not start and the standard log file just prints upto ..."WAVECAR not read". However, when I reduce the no of KPOINTS to 3x3x1 Gamma Points, The optimizations start but are extremely slow.
Please let me know why should it behave this way ??
POSCAR
---------------------------------------------------------------------------------
4TLCeO2
1.00000
11.478500366 0.000000000 0.000000000
-5.739250183 9.940672915 0.000000000
0.000000000 0.000000000 40.934101105
O Ce
72 36
Selective Dynamics
Direct
0.000000010 0.000000000 0.267120004 T T T
0.000000010 0.666670024 0.267120004 T T T
0.666670024 0.000000000 0.267120004 T T T
0.666670024 0.666670024 0.267120004 T T T
0.333330005 0.000000000 0.267120004 T T T
0.333330005 0.666670024 0.267120004 T T T
0.000000000 0.333330005 0.267120004 T T T
0.333330005 0.333330005 0.267120004 T T T
0.666670024 0.333330005 0.267120004 T T T
0.111110002 0.555559993 0.228960007 T T T
0.777779996 0.555559993 0.228960007 T T T
0.444440007 0.222220004 0.228960007 T T T
0.444440007 0.888890028 0.228960007 T T T
0.111110002 0.222220004 0.228960007 T T T
0.111110002 0.888890028 0.228960007 T T T
0.777779996 0.222220004 0.228960007 T T T
0.777779996 0.888890028 0.228960007 T T T
0.444440007 0.555559993 0.228960007 T T T
0.555559993 0.111110002 0.190799996 T T T
0.555559993 0.777779996 0.190799996 T T T
0.555559993 0.444440007 0.190799996 T T T
0.222220004 0.444440007 0.190799996 T T T
0.888890028 0.444440007 0.190799996 T T T
0.222220004 0.111110002 0.190799996 T T T
0.222220004 0.777779996 0.190799996 T T T
0.888890028 0.111110002 0.190799996 T T T
0.888890028 0.777779996 0.190799996 T T T
0.000000000 0.333330005 0.152640000 T T T
0.666670024 0.333330005 0.152640000 T T T
0.000000010 0.000000010 0.152640000 T T T
0.000000010 0.666670024 0.152640000 T T T
0.666670024 0.000000000 0.152640000 T T T
0.666670024 0.666670024 0.152640000 T T T
0.333330005 0.000000000 0.152640000 T T T
0.333330005 0.666670024 0.152640000 T T T
0.333330005 0.333330005 0.152640000 T T T
0.444440007 0.555559993 0.114480004 T T T
0.444440007 0.222220004 0.114480004 T T T
0.444440007 0.888890028 0.114480004 T T T
0.111110002 0.222220004 0.114480004 T T T
0.111110002 0.888890028 0.114480004 T T T
0.777779996 0.222220004 0.114480004 T T T
0.777779996 0.888890028 0.114480004 T T T
0.111110002 0.555559993 0.114480004 T T T
0.777779996 0.555559993 0.114480004 T T T
0.888890028 0.444440007 0.076320000 T T T
0.222220004 0.444440007 0.076320000 T T T
0.555559993 0.444440007 0.076320000 T T T
0.555559993 0.111110002 0.076320000 T T T
0.555559993 0.777779996 0.076320000 T T T
0.222220004 0.111110002 0.076320000 T T T
0.222220004 0.777779996 0.076320000 T T T
0.888890028 0.111110002 0.076320000 T T T
0.888890028 0.777779996 0.076320000 T T T
0.333330005 0.000000000 0.038160000 T T T
0.333330005 0.666670024 0.038160000 T T T
0.000000010 0.000000100 0.038160000 T T T
0.000000010 0.666670024 0.038160000 T T T
0.666670024 0.000000000 0.038160000 T T T
0.666670024 0.666670024 0.038160000 T T T
0.000000000 0.333330005 0.038160000 T T T
0.666670024 0.333330005 0.038160000 T T T
0.333330005 0.333330005 0.038160000 T T T
0.111110002 0.555559993 0.000000000 F F F
0.777779996 0.555559993 0.000000000 F F F
0.444440007 0.222220004 0.000000000 F F F
0.444440007 0.888890028 0.000000000 F F F
0.444440007 0.555559993 0.000000000 F F F
0.111110002 0.222220004 0.000000000 F F F
0.111110002 0.888890028 0.000000000 F F F
0.777779996 0.222220004 0.000000000 F F F
0.777779996 0.888890028 0.000000000 F F F
0.222220004 0.444440007 0.248040006 T T T
0.888890028 0.444440007 0.248040006 T T T
0.555559993 0.111110002 0.248040006 T T T
0.555559993 0.777779996 0.248040006 T T T
0.555559993 0.444440007 0.248040006 T T T
0.222220004 0.111110002 0.248040006 T T T
0.222220004 0.777779996 0.248040006 T T T
0.888890028 0.111110002 0.248040006 T T T
0.888890028 0.777779996 0.248040006 T T T
0.444440007 0.222220004 0.171719998 T T T
0.444440007 0.888890028 0.171719998 T T T
0.444440007 0.555559993 0.171719998 T T T
0.111110002 0.222220004 0.171719998 T T T
0.111110002 0.888890028 0.171719998 T T T
0.777779996 0.222220004 0.171719998 T T T
0.777779996 0.888890028 0.171719998 T T T
0.111110002 0.555559993 0.171719998 T T T
0.777779996 0.555559993 0.171719998 T T T
0.222220004 0.444440007 0.019080000 T T T
0.888890028 0.444440007 0.019080000 T T T
0.555559993 0.111110002 0.019080000 T T T
0.555559993 0.777779996 0.019080000 T T T
0.888890028 0.111110002 0.019080000 T T T
0.888890028 0.777779996 0.019080000 T T T
0.555559993 0.444440007 0.019080000 T T T
0.222220004 0.111110002 0.019080000 T T T
0.222220004 0.777779996 0.019080000 T T T
0.000000000 0.333330005 0.095399998 T T T
0.666670024 0.333330005 0.095399998 T T T
0.333330005 0.333330005 0.095399998 T T T
0.333330005 0.000000000 0.095399998 T T T
0.333330005 0.666670024 0.095399998 T T T
0.000000010 0.000000010 0.095399998 T T T
0.000000010 0.666670024 0.095399998 T T T
0.666670024 0.000000000 0.095399998 T T T
0.666670024 0.666670024 0.095399998 T T T
---------------------------------------------------------------------------------
INCAR
---------------------------------------------------------------------------------
SYSTEM = Ceria(111)
ENCUT=500
ISIF=2
IBRION=2
NSW=200
ISPIN=2
LWAVE=.FALSE. # change for re-runs of converged calculations
LCHARG=.FALSE.
SYMPREC=1E-8
#LVTOT=.TRUE. # use for getting the local absolute potential
LREAL=AUTO
#LAECHG=.TRUE. # use for bader analysis
ADDGRID=.TRUE.
#IDIPOL=3 # do use once the code really works
#LDIPOL=.TRUE.
#LORBIT=11 # use when DOS is needed
#NBANDS=300 # use if vasp suggests
ISMEAR=1
SIGMA=0.1
LASPH=.TRUE.
LDAU=.TRUE.
LDAUTYPE= 3
LDAUL=-1 3 #for O,Ce
LDAUU=0 5 #for O,Ce
LDAUJ=0 1 #for O,Ce
LMAXMIX=6
---------------------------------------------------------------------------------