Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Gavin Abo
Just a thought, maybe it doesn't help, but have you checked the :log file for which programs are involved in your calculation. Then, used that to check if there are differences in the code for those programs that could affect your calculation using the updates list on WIEN2k website:

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Laurence Marks
1) Do "grep :MVORB *.scf", it is very possible to get fake convergence where the density is converged but the orbital potential is not. 2) Take the *converged* case run with 16.1, without changing anything except removing the case.bro* files, and run it with 17.1 3) Repeat this the other way

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Kefeng wang
Dear Dr. Tran, Thanks. I checked the input files carefully and repeated both of them twice. I still obtained the same results. For all of them, the SCF convergence has been reached and there is no error observed. :DIS at the last iteration for wien17.1 is 0.38 and for wien16.1 is 0.000190.

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Kefeng wang
Dear Prof. Blaha, Yes. I know the difference of 10 eV is unbelievable. I used the same struct file, k list, and coupling U. All the input files are the same. Only the difference is the different wien2k version. For wien2k16.1, the total energies for FM+U and AFM+U are -29235.40982268 and

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread tran
Hi, As suggested by others, maybe you made a mistake with one of the calculations. Repeat the two FM+U calculations and check that all input file are the same and that the SCF convergence has been reached. Maybe one of the calculations was not converged at all? What about :DIS in case.scf? FT

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Kefeng wang
Dear Dr. Tran, Thanks a lot for your reply. The magnetic moments for two calculations are also different. the one for wien16.1 is 0.32uB/Fe and the one for wien17.1 is 0.39uB/Fe. Best, K. Wang On Wed, Apr 25, 2018 at 12:05 PM, Kefeng wang wrote: > Dear all, > > I am using

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Lyudmila Dobysheva
25.04.2018 20:05, Kefeng wang wrote: the FM+U calculations for CoAs. It turns out that the difference of the total energy is about 10 eV. ... However,  the difference of the total energy of theAFM+U calculations for CoAs Wien16.1 and  Wien17.1 is almost zero. A bit unclear... Nevertheless,

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-26 Thread Peter Blaha
Something must have been gone very wrong in the FM case. 10 eV is an enormous difference. Even the difference between magnetic states (FM,AFM) should be on the mRy level. Compare the energies of the AFM and FM calculations. They should be "similar" and then you know which of the FM

Re: [Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-25 Thread tran
Hi, What about the magnetic moments of the two FM+U calculations? Are they the same or also different? FT On Wednesday 2018-04-25 18:05, Kefeng wang wrote: Date: Wed, 25 Apr 2018 18:05:28 From: Kefeng wang Reply-To: A Mailing list for WIEN2k users

[Wien] the difference of the total energy is about 10 eV using wien16.1 and wien17.1

2018-04-25 Thread Kefeng wang
Dear all, I am using Wien16.1 and Wien17.1 to perform the FM+U calculations for CoAs. It turns out that the difference of the total energy is about 10 eV. In both calculations, I used the same struct file as following: CoAs P LATTICE,NONEQUIV.ATOMS: 2 62 Pnma MODE OF CALC=RELA unit=ang