non-collinear spin vasp 4.6.28 vs 5.2.11
Moderators: Global Moderator, Moderator
non-collinear spin vasp 4.6.28 vs 5.2.11
Dear all,
I'm investigating non-collinear magnetization in BiFeO3.
I've tested my calculation with vasp version 4.6.28, 4.6.35 and 5.2.11 and I find that vasp 5.2.11 gives different value for the magnetization w.r.t. previous vasp version!
In my input BiFeO3 is in the rhombohebral cell, SAXIS 100 and a magnetization is expected along 'y' (for this choice of coordinates).
To summarize the creepy results:
vasp 4.6.28
magnetization 0.0000198 -0.0363467 -0.0000001
vasp 4.6.35
magnetization 0.0000198 -0.0363468 -0.0000001
vasp 5.2.11
magnetization 0.0000119 -0.0762098 -0.0000001
=> the magnetization along y is doubled in vasp 5.2.11!
vasp 4.6.35 and 5.2.11 have been compiled with the same options and run on the same machine.
I've started each calculation from the CHGCAR of a previous collinear calculation, so everything is the same, just the vasp version is changed.
I've also checked against vasp 4.6.28 compiled and run on another machine, and I've got confirmation of the 4.6.35 results.
The question is:
which vasp version should I trust?
Hopefully the latest one should be the better one, maybe some bug fixes have occurred. Can anyone comment on this?
Thanks.
All the best,
Zeila
I'm investigating non-collinear magnetization in BiFeO3.
I've tested my calculation with vasp version 4.6.28, 4.6.35 and 5.2.11 and I find that vasp 5.2.11 gives different value for the magnetization w.r.t. previous vasp version!
In my input BiFeO3 is in the rhombohebral cell, SAXIS 100 and a magnetization is expected along 'y' (for this choice of coordinates).
To summarize the creepy results:
vasp 4.6.28
magnetization 0.0000198 -0.0363467 -0.0000001
vasp 4.6.35
magnetization 0.0000198 -0.0363468 -0.0000001
vasp 5.2.11
magnetization 0.0000119 -0.0762098 -0.0000001
=> the magnetization along y is doubled in vasp 5.2.11!
vasp 4.6.35 and 5.2.11 have been compiled with the same options and run on the same machine.
I've started each calculation from the CHGCAR of a previous collinear calculation, so everything is the same, just the vasp version is changed.
I've also checked against vasp 4.6.28 compiled and run on another machine, and I've got confirmation of the 4.6.35 results.
The question is:
which vasp version should I trust?
Hopefully the latest one should be the better one, maybe some bug fixes have occurred. Can anyone comment on this?
Thanks.
All the best,
Zeila
Last edited by Zeila on Wed Jun 01, 2011 6:30 pm, edited 1 time in total.
-
- Full Member
- Posts: 106
- Joined: Wed Jan 16, 2008 1:27 pm
- License Nr.: 567
- Location: France
non-collinear spin vasp 4.6.28 vs 5.2.11
Hi
I observed the same problem using the LDA+U method. There's a difference between the two versions and I couldn't figure out why.
It seems something was changed and not mentioned in the manual, as usual...
Maybe your problem is related to mine, so I will wait for an answer as well.
Regards
<span class='smallblacktext'>[ Edited Wed Jun 01 2011, 07:49PM ]</span>
I observed the same problem using the LDA+U method. There's a difference between the two versions and I couldn't figure out why.
It seems something was changed and not mentioned in the manual, as usual...
Maybe your problem is related to mine, so I will wait for an answer as well.
Regards
<span class='smallblacktext'>[ Edited Wed Jun 01 2011, 07:49PM ]</span>
Last edited by boris on Wed Jun 01, 2011 7:43 pm, edited 1 time in total.
non-collinear spin vasp 4.6.28 vs 5.2.11
I'm running the calculation with LSDA+U !
(sorry, I've forgot to mention it...)
so, yes my problem is related to your.
thanks for giving a feedback,
Zeila
(sorry, I've forgot to mention it...)
so, yes my problem is related to your.
thanks for giving a feedback,
Zeila
Last edited by Zeila on Thu Jun 02, 2011 3:20 pm, edited 1 time in total.
-
- Full Member
- Posts: 106
- Joined: Wed Jan 16, 2008 1:27 pm
- License Nr.: 567
- Location: France
non-collinear spin vasp 4.6.28 vs 5.2.11
Ok so this is a problem I've been working on for a while now.
When using the LDA+U, there's a difference in the total energy between the two versions. It seems that it's because the occupation matrices are calculated differently in the two versions, therefore yielding different magnetic moments.
It might be the cause of the discrepancy.
So maybe they changed the way to calculate the occupancies, or they changed something in the PAW subroutine, I have no idea....
Maybe we'll get an answer. But they should be a bit concerned by such a glitch in the code.
Boris
<span class='smallblacktext'>[ Edited Fri Jun 03 2011, 12:17AM ]</span>
When using the LDA+U, there's a difference in the total energy between the two versions. It seems that it's because the occupation matrices are calculated differently in the two versions, therefore yielding different magnetic moments.
It might be the cause of the discrepancy.
So maybe they changed the way to calculate the occupancies, or they changed something in the PAW subroutine, I have no idea....
Maybe we'll get an answer. But they should be a bit concerned by such a glitch in the code.
Boris
<span class='smallblacktext'>[ Edited Fri Jun 03 2011, 12:17AM ]</span>
Last edited by boris on Fri Jun 03, 2011 12:15 am, edited 1 time in total.
non-collinear spin vasp 4.6.28 vs 5.2.11
Indeed, I've repeated my calculation without "U" and now I'm getting the same results for both code versions.
Who can we ask to know which version is reliable?
Zeila
Who can we ask to know which version is reliable?
Zeila
Last edited by Zeila on Fri Jun 03, 2011 3:18 pm, edited 1 time in total.
non-collinear spin vasp 4.6.28 vs 5.2.11
[quote="boris"]Ok so this is a problem I've been working on for a while now.
When using the LDA+U, there's a difference in the total energy between the two versions. It seems that it's because the occupation matrices are calculated differently in the two versions, therefore yielding different magnetic moments.
It might be the cause of the discrepancy.
So maybe they changed the way to calculate the occupancies, or they changed something in the PAW subroutine, I have no idea....
Maybe we'll get an answer. But they should be a bit concerned by such a glitch in the code.
Boris
[/quote]
Dear Dr. Boris Dorado,
Sorry this is not directly related to this thread, but I just have a question about your methods for avoiding/locating local minima by controlling occupation matrices for DFT+U calculations using the VASP code.
You mentioned "As mentioned above, we implemented the possibility of defining initial occupation matrices and impose them during the calculation of the DFT+U potential in the VASP code. We
thus precondition the calculation of the potential which is then applied as a correction to the standard DFT potential. Occupation matrices are imposed during the first ten electronic steps. After these ten steps, the constraint is lifted and the calculation is left to converge self-consistently on its own." (Phys. Rev. B. 82, 035114 (2010)).
Could you please elaborate on how you define the initial occupation matrix and impose them during the first ten electronic steps in VASP? For example, maybe you do it by setting FERWE or FERDO tags in INCAR? Or you modified part of the source codes (e.g., the LDApU.F subroutine) to realize this implementation?
We are much interested in your invention and am trying to use your methods in our calculations using LDA+U. If you could kindly provide more details, we will appreciate much!
Cheers,
CMG
When using the LDA+U, there's a difference in the total energy between the two versions. It seems that it's because the occupation matrices are calculated differently in the two versions, therefore yielding different magnetic moments.
It might be the cause of the discrepancy.
So maybe they changed the way to calculate the occupancies, or they changed something in the PAW subroutine, I have no idea....
Maybe we'll get an answer. But they should be a bit concerned by such a glitch in the code.
Boris
[/quote]
Dear Dr. Boris Dorado,
Sorry this is not directly related to this thread, but I just have a question about your methods for avoiding/locating local minima by controlling occupation matrices for DFT+U calculations using the VASP code.
You mentioned "As mentioned above, we implemented the possibility of defining initial occupation matrices and impose them during the calculation of the DFT+U potential in the VASP code. We
thus precondition the calculation of the potential which is then applied as a correction to the standard DFT potential. Occupation matrices are imposed during the first ten electronic steps. After these ten steps, the constraint is lifted and the calculation is left to converge self-consistently on its own." (Phys. Rev. B. 82, 035114 (2010)).
Could you please elaborate on how you define the initial occupation matrix and impose them during the first ten electronic steps in VASP? For example, maybe you do it by setting FERWE or FERDO tags in INCAR? Or you modified part of the source codes (e.g., the LDApU.F subroutine) to realize this implementation?
We are much interested in your invention and am trying to use your methods in our calculations using LDA+U. If you could kindly provide more details, we will appreciate much!
Cheers,
CMG
Last edited by CMG on Sat Jun 04, 2011 2:52 am, edited 1 time in total.
non-collinear spin vasp 4.6.28 vs 5.2.11
Hi everybody!
I've clarified the problem. I was using LSDA+U with the Durarev approach, which is buggy in vasp4.6.
I've re-did my calculations with the Liechtenstein approach and now vasp4.6 and vasp5.2 give the same results which, by the way, are the same as Dudarev in vasp5.2 (since, dudarev = liechestein with J=0).
Here are the results:
vasp 4.6.28, Dudarev: magnetization 0.0000198 -0.0363467 -0.0000001
vasp 4.6.35, Dudarev: magnetization 0.0000198 -0.0363468 -0.0000001
vasp 5.2.11, Dudarev: magnetization 0.0000119 -0.0762098 -0.0000001
Liechtenstein:
vasp4.6, U=4 eV, J=0: magnetization??? 0.0000119? -0.0762081? -0.0000001
vasp5.2, U=4 eV, J=0: magnetization??? 0.0000119? -0.0762096? -0.0000001
So, don't use vasp4.6 with dudarev LDA+U!!!
Cheers!
Zeila
I've clarified the problem. I was using LSDA+U with the Durarev approach, which is buggy in vasp4.6.
I've re-did my calculations with the Liechtenstein approach and now vasp4.6 and vasp5.2 give the same results which, by the way, are the same as Dudarev in vasp5.2 (since, dudarev = liechestein with J=0).
Here are the results:
vasp 4.6.28, Dudarev: magnetization 0.0000198 -0.0363467 -0.0000001
vasp 4.6.35, Dudarev: magnetization 0.0000198 -0.0363468 -0.0000001
vasp 5.2.11, Dudarev: magnetization 0.0000119 -0.0762098 -0.0000001
Liechtenstein:
vasp4.6, U=4 eV, J=0: magnetization??? 0.0000119? -0.0762081? -0.0000001
vasp5.2, U=4 eV, J=0: magnetization??? 0.0000119? -0.0762096? -0.0000001
So, don't use vasp4.6 with dudarev LDA+U!!!
Cheers!
Zeila
Last edited by Zeila on Sun Jun 05, 2011 4:53 pm, edited 1 time in total.
-
- Full Member
- Posts: 106
- Joined: Wed Jan 16, 2008 1:27 pm
- License Nr.: 567
- Location: France
non-collinear spin vasp 4.6.28 vs 5.2.11
@CMG: I will gladly help you do that. I guess we should talk about that via email. My email is: bdorado [at] lanl.gov
Feel free to contact me.
@Zeila: This is very interesting. I will try to run the calculation with the Liechtenstein approach to see if there is any difference. Thanks for noticing that
Boris
Edit : I still have the total energy difference even with the Liechtenstein approach:
VASP 4.6:
free energy TOTEN = 1542.041737 eV
VASP 5.2:
free energy TOTEN = 1542.038787 eV
However, I get the same results if I consider a d electron system instead of a f electron system.
Therefore, it seems that there is a difference in the two versions only with f electrons...
<span class='smallblacktext'>[ Edited Mon Jun 06 2011, 08:04PM ]</span>
Feel free to contact me.
@Zeila: This is very interesting. I will try to run the calculation with the Liechtenstein approach to see if there is any difference. Thanks for noticing that
Boris
Edit : I still have the total energy difference even with the Liechtenstein approach:
VASP 4.6:
free energy TOTEN = 1542.041737 eV
VASP 5.2:
free energy TOTEN = 1542.038787 eV
However, I get the same results if I consider a d electron system instead of a f electron system.
Therefore, it seems that there is a difference in the two versions only with f electrons...
<span class='smallblacktext'>[ Edited Mon Jun 06 2011, 08:04PM ]</span>
Last edited by boris on Mon Jun 06, 2011 7:23 pm, edited 1 time in total.
-
- Newbie
- Posts: 6
- Joined: Thu Jan 25, 2007 1:26 pm
non-collinear spin vasp 4.6.28 vs 5.2.11
Hi everybody,
actually, besides all the discrepancies mentioned already above, I have difference in total energy even without U, if I calculate system with f-electrons with different VASP versions:
system UO22+, charged, so cubic unit cell, multipole corrections are switched on
XC = PW91
Total energy v 4.6 -4.72588
Total energy v 5.2 -4.74441
Waiting for the answer,
regards
actually, besides all the discrepancies mentioned already above, I have difference in total energy even without U, if I calculate system with f-electrons with different VASP versions:
system UO22+, charged, so cubic unit cell, multipole corrections are switched on
XC = PW91
Total energy v 4.6 -4.72588
Total energy v 5.2 -4.74441
Waiting for the answer,
regards
Last edited by Alenushka on Tue Jun 21, 2011 8:30 pm, edited 1 time in total.
-
- Full Member
- Posts: 106
- Joined: Wed Jan 16, 2008 1:27 pm
- License Nr.: 567
- Location: France
non-collinear spin vasp 4.6.28 vs 5.2.11
Hi
I also have a difference in the total energy without the U. So it seems that as soon as I consider f electrons, there is a difference between 4.6 and 5.2.
Something has changed, and I don't know what it is...
I also have a difference in the total energy without the U. So it seems that as soon as I consider f electrons, there is a difference between 4.6 and 5.2.
Something has changed, and I don't know what it is...
Last edited by boris on Fri Jul 08, 2011 2:10 am, edited 1 time in total.
-
- Full Member
- Posts: 106
- Joined: Wed Jan 16, 2008 1:27 pm
- License Nr.: 567
- Location: France
non-collinear spin vasp 4.6.28 vs 5.2.11
Hi everyone,
Here's the answer to our problems: the 4.6 version does not take into account the one-center multipoles up to L=6 (only up to L=4).
This has been corrected in the new 4.6 version available on the vasp FTP server.
Regards
Here's the answer to our problems: the 4.6 version does not take into account the one-center multipoles up to L=6 (only up to L=4).
This has been corrected in the new 4.6 version available on the vasp FTP server.
Regards
Last edited by boris on Tue Jul 12, 2011 11:29 pm, edited 1 time in total.