PBE0 cannot be used in vasp5.2?
Moderators: Global Moderator, Moderator
PBE0 cannot be used in vasp5.2?
DEAR ALL:
why the label "LHFCALC" can not be used in vasp5.2?
The calculation seems to get into endless loop after five interaction without any output any more.
Thanks!
why the label "LHFCALC" can not be used in vasp5.2?
The calculation seems to get into endless loop after five interaction without any output any more.
Thanks!
Last edited by yuanxun on Mon Aug 24, 2009 6:45 am, edited 1 time in total.
-
- Hero Member
- Posts: 584
- Joined: Tue Nov 16, 2004 2:21 pm
- License Nr.: 5-67
- Location: Germany
PBE0 cannot be used in vasp5.2?
The first 5 iterations are DFT only, then exact exchange is calculated. This is _very_ time consuming. Try a small system.
Cheers
Alex
Cheers
Alex
Last edited by alex on Sun Aug 30, 2009 9:57 am, edited 1 time in total.
- VASP001
- Jr. Member
- Posts: 72
- Joined: Tue Oct 27, 2009 2:56 pm
- License Nr.: 176, upgraded to paw
- Location: Hefei National Laboratory for Physical Sciences at the Microscale
PBE0 cannot be used in vasp5.2?
[quote="yuanxun"]DEAR ALL:
why the label "LHFCALC" can not be used in vasp5.2?
The calculation seems to get into endless loop after five interaction without any output any more.
Thanks! [/quote]
You can check the file of OUTCAR ,to find the LHFCALC=.TRUE. or .FALSE.
why the label "LHFCALC" can not be used in vasp5.2?
The calculation seems to get into endless loop after five interaction without any output any more.
Thanks! [/quote]
You can check the file of OUTCAR ,to find the LHFCALC=.TRUE. or .FALSE.
Last edited by VASP001 on Sat Jan 09, 2010 6:28 pm, edited 1 time in total.
-
- Full Member
- Posts: 201
- Joined: Thu Nov 02, 2006 4:35 pm
- License Nr.: 5-532
- Location: Ghent, Belgium
- Contact:
PBE0 cannot be used in vasp5.2?
the first 5 iterations are the normal speed (as alex mentioned) after this it becomes very slow, we are talking about orders of magnitude:
A small example of a single atom fcc system with 5x5x5 k-points and HSE06 calculation...there are only 10 kpoints in the IBZand the the steps become 100x slower than the first 5:
LOOP: cpu time 8.38: real time 8.38
LOOP: cpu time 2.46: real time 2.46
LOOP: cpu time 2.13: real time 2.14
LOOP: cpu time 2.08: real time 2.08
LOOP: cpu time 3.09: real time 3.09
LOOP: cpu time 311.79: real time 312.30
LOOP: cpu time 288.99: real time 289.54
LOOP: cpu time 285.76: real time 286.52
LOOP: cpu time 284.03: real time 284.25
LOOP: cpu time 282.32: real time 282.33
LOOP: cpu time 285.06: real time 285.34
LOOP: cpu time 281.36: real time 281.36
LOOP: cpu time 279.80: real time 279.81
LOOP: cpu time 283.62: real time 283.74
LOOP: cpu time 288.30: real time 289.55
LOOP: cpu time 287.12: real time 287.44
LOOP: cpu time 285.60: real time 286.15
LOOP: cpu time 286.33: real time 286.51
LOOP: cpu time 286.84: real time 287.50
Cheers Danny
A small example of a single atom fcc system with 5x5x5 k-points and HSE06 calculation...there are only 10 kpoints in the IBZand the the steps become 100x slower than the first 5:
LOOP: cpu time 8.38: real time 8.38
LOOP: cpu time 2.46: real time 2.46
LOOP: cpu time 2.13: real time 2.14
LOOP: cpu time 2.08: real time 2.08
LOOP: cpu time 3.09: real time 3.09
LOOP: cpu time 311.79: real time 312.30
LOOP: cpu time 288.99: real time 289.54
LOOP: cpu time 285.76: real time 286.52
LOOP: cpu time 284.03: real time 284.25
LOOP: cpu time 282.32: real time 282.33
LOOP: cpu time 285.06: real time 285.34
LOOP: cpu time 281.36: real time 281.36
LOOP: cpu time 279.80: real time 279.81
LOOP: cpu time 283.62: real time 283.74
LOOP: cpu time 288.30: real time 289.55
LOOP: cpu time 287.12: real time 287.44
LOOP: cpu time 285.60: real time 286.15
LOOP: cpu time 286.33: real time 286.51
LOOP: cpu time 286.84: real time 287.50
Cheers Danny
Last edited by Danny on Thu Mar 04, 2010 1:45 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
PBE0 cannot be used in vasp5.2?
it is the exact exchange part which makes the calculation slow down.
whenever possible, use a k-mesh that cen be reduced for the exact-exchange contributions so you can make use of NKRED
whenever possible, use a k-mesh that cen be reduced for the exact-exchange contributions so you can make use of NKRED
Last edited by admin on Thu Mar 04, 2010 4:30 pm, edited 1 time in total.
-
- Full Member
- Posts: 201
- Joined: Thu Nov 02, 2006 4:35 pm
- License Nr.: 5-532
- Location: Ghent, Belgium
- Contact:
PBE0 cannot be used in vasp5.2?
What are the properties of a k-mesh that can be reduced?
eg: k-point mesh Nx*Ny*Nz
with Nx a integer multiple of NKREDX, Ny of NKREDY and Nz of NKREDZ?
Danny
eg: k-point mesh Nx*Ny*Nz
with Nx a integer multiple of NKREDX, Ny of NKREDY and Nz of NKREDZ?
Danny
Last edited by Danny on Fri Mar 05, 2010 8:57 am, edited 1 time in total.
-
- Newbie
- Posts: 24
- Joined: Wed Feb 17, 2010 11:34 pm
- License Nr.: 1118
PBE0 cannot be used in vasp5.2?
I dont really undrestand this manual entry for NKRED.
http://cms.mpi.univie.ac.at/vasp/vasp/W ... plied.html
I cant get this to work at all with various entries. =>
k-point set generating an fcc lattice -> 4 k-points in BZ
0
diret
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
0 0 0
Is it intended that ISEMAR=-5 is used with this set??
Do i need to get the Wavecar first?
Is it suitable for a loopscript?
<span class='smallblacktext'>[ Edited Wed Oct 20 2010, 02:59PM ]</span>
http://cms.mpi.univie.ac.at/vasp/vasp/W ... plied.html
I cant get this to work at all with various entries. =>
k-point set generating an fcc lattice -> 4 k-points in BZ
0
diret
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
0 0 0
Is it intended that ISEMAR=-5 is used with this set??
Do i need to get the Wavecar first?
Is it suitable for a loopscript?
<span class='smallblacktext'>[ Edited Wed Oct 20 2010, 02:59PM ]</span>
Last edited by Sonny on Wed Oct 13, 2010 7:38 am, edited 1 time in total.