[Wien] Units in Stoner criteria

[Alex]

Dear all,

I was wondering that what the correct units of Stoner parameter (*I*) and
DOS(*E*_F) (density of state at Fermi level) used in the Stoner criteria of
magnetism (*I**DOS(*E*_F) >‌ 1) is?
Does DOS(*E*_F) should have a unit of (atom/energy) or (atom/volume/energy)
... and then similar change for the unit of *I*?
Different defined units for them in different literatures cause confusion
(for me), or for example when Coulomb repulsion (*U*) is used instead of
*I,* the Stoner criteria change to *U**D*(E_*F*) > *(1/v) where v is the
unit cell volume*.*
Or for instance inspection of the formulas of 14--27 of below reference [1]
convinces me that the unit of volume is missing somewhere.

[1].
http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.92.9357=rep1=pdf

In my specific case I consider an simple FCC Ni, and then in non-magnetic
phase I calculated the DOS(*E*_F) for one of the Ni atoms such that the
integration of the corresponding DOS curve until the Fermi level gives
around 10 electron.

Thanks for you help.

Regards,
Alex
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] problem while running "x w2w" command

[pachineela rambabu]

Dear Gavin,
 thank you very much for the kind reply. According to your
suggestion, I installed 17.1 version and it is working properly at the
moment.

thanking you,


On Thu, Dec 7, 2017 at 6:59 PM, Gavin Abo  wrote:

> SIGSEGV errors can have many possible causes [
> https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14476.html
> ].
>
> In your case, I see W_14.2 below.  So perhaps it is due to a bug in an the
> older WIEN2k 14.2 or wien2wannier version that you are using.  You might
> try the newer wien2wannier version in WIEN2k 17.1.  However, you have to
> patch it [ https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.
> at/msg16099.html ].
>
> What ifort compiler version are you using?  If it is a 2016 and 2017,
> maybe it is caused by the IO problem with large files [
> https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg16802.html
> ].
> On 12/6/2017 3:58 AM, pachineela rambabu wrote:
>
> While running "x w2w" command, I encountered following error:
>
> forrtl: severe (174): SIGSEGV, segmentation fault occurred
> Image  PCRoutineLine
> Source
> w2wc   004EC771  Unknown   Unknown  Unknown
> w2wc   004EA8AB  Unknown   Unknown  Unknown
> w2wc   004935B4  Unknown   Unknown  Unknown
> w2wc   004933C6  Unknown   Unknown  Unknown
> w2wc   0043AFE9  Unknown   Unknown  Unknown
> w2wc   0043EBF6  Unknown   Unknown  Unknown
> libpthread-2.23.s  2AD859B82390  Unknown   Unknown  Unknown
> w2wc   00465AF7  Unknown   Unknown  Unknown
> w2wc   00426F83  read_vec_  52
> read_vec.f
> w2wc   00420F70  MAIN__159  main.f
> w2wc   00402FAE  Unknown   Unknown  Unknown
> libc-2.23.so   2AD85A65D830  __libc_start_main Unknown
> Unknown
> w2wc   00402EA9  Unknown   Unknown  Unknown
> 0.3u 0.2s 0:01.03 56.3% 0+0k 88016+16io 0pf+0w
> error: command   /home/kanchana/Soft/W_14.2/w2wc w2w.def   failed
>
>
> Please tell me how to rectify this problem.  For cubic systems, I didn't
> get any problem. For hexagonal system, this problem occurred.
>
> thanks in advance,
> --
> *P. Rambabu*
> PhD Scholor
> Department of Physics, IIT Hyderabad, India
> Mobile: 9074508220.
>
>
> ___
> Wien mailing list
> Wien@zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> SEARCH the MAILING-LIST at:  http://www.mail-archive.com/
> wien@zeus.theochem.tuwien.ac.at/index.html
>
>


-- 
*P. Rambabu*
PhD Scholor
Department of Physics, IIT Hyderabad, India
Mobile: 9074508220.
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Pavel Ondračka]

Laurence Marks píše v Čt 07. 12. 2017 v 07:37 -0600:
> Probably overlapping response from Peter...
> 
> Based upon the Check-mixing output, the density is very well
> converged and did not change with the jump in energy. The most likely
> way for :DEN alone to change, and nothing else, is that there was a
> change in how EXC is calculated, for instance using a metaGGA or
> changing case.in0, without a change in VXC.
> 
> Do "grep :POT *scf" to check this.
> 
This is a standard PBE calculation:
:POT  : POTENTIAL OPTION EX_PBE EC_PBE VX_PBE VC_PBEfor all the cycles.
In fact the surprising thing was that this started out of the blue
after at least 15 iterations when the calculation was already looking
quite converged.

BTW as mentioned in the beginning this might have to do something with
the atomic spheres, I converged this structure with the almost touching
spheres (0% reduction) and later I did reduce_rmt_lapw -r 3.5 to
produce this case, hence it was already pre-converged and I only wanted
to fine-converge it, however the energy oscillations showed...
Dunno of this might play a role.

Best regards
Pavel
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

[Wien] Large oscillations in total energy

[Pavel Ondračka]

Dear Wien2k mailing list,

I'm seeing some strange energy oscillations in one calculation near the
end when everything seems almost converged.

:ENERGY convergence:  0 0.0001 .00390097
:ENERGY convergence:  0 0.0001 .016890645000
:ENERGY convergence:  0 0.0001 .00965980
:ENERGY convergence:  0 0.0001 .016205335000
:ENERGY convergence:  0 0.0001 .00896832
:ENERGY convergence:  0 0.0001 55.639853345000
:ENERGY convergence:  0 0.0001 27.821044505000
:ENERGY convergence:  0 0.0001 27.82352526

the charge is converged already:

:CHARGE convergence:  1 0.001 -.000920
:CHARGE convergence:  1 0.001 -.000977
:CHARGE convergence:  1 0.001 -.000933
:CHARGE convergence:  1 0.001 -.000972
:CHARGE convergence:  1 0.001 -.000984
:CHARGE convergence:  1 0.001 -.000988
:CHARGE convergence:  1 0.001 -.000997
:CHARGE convergence:  1 0.001 -.000998

This is quite large calculation of TiSiO amorphous-like cell (96 atoms,
 struct file attached, compressed in order to pass the moderation
queue). The numerical parameters are quite standard,
2x2x2 shifted k-grid (can also reproduce this with 3x3x3), Rkmax 7.5,
PBE, rest of stuff at default. I have done a lot of calculation with
similar cells and such parameters without any problems (albeit mostly
some time ago with older Wien2k versions)...

I'm trying to calculate the binding energy of the O 1s for the first
oxygen using the Slaters transition state approach. Surprisingly when I
use almost touching spheres it converges (O:1.48, Ti: 1.72, Si:1.48),
however when I reduce them by 3% (as in the attached struct file) this
problems surfaces. However I need this reduced spheres for comparison
with another calculation. BTW the spheres around 1.5 for oxygen
should be reasonable anyway?

No errors or warnings anywhere, except for the expected 
WARN  :CHARGED CELL with  -0.500

This is with Wien2k 17.1, compiled with MKL 2017.0.098 and ifort.

Any ideas what might be causing this problem?
Best regards
Pavel

TiSiO.tar.xz
Description: application/xz-compressed-tar
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Laurence Marks]

I suggest using the "Check" command, which shows more useful information
and post the output. Without that it is hard to guess.

_
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu

On Dec 7, 2017 4:20 AM, "Pavel Ondračka"  wrote:

> Dear Wien2k mailing list,
>
> I'm seeing some strange energy oscillations in one calculation near the
> end when everything seems almost converged.
>
> :ENERGY convergence:  0 0.0001 .00390097
> :ENERGY convergence:  0 0.0001 .016890645000
> :ENERGY convergence:  0 0.0001 .00965980
> :ENERGY convergence:  0 0.0001 .016205335000
> :ENERGY convergence:  0 0.0001 .00896832
> :ENERGY convergence:  0 0.0001 55.639853345000
> :ENERGY convergence:  0 0.0001 27.821044505000
> :ENERGY convergence:  0 0.0001 27.82352526
>
> the charge is converged already:
>
> :CHARGE convergence:  1 0.001 -.000920
> :CHARGE convergence:  1 0.001 -.000977
> :CHARGE convergence:  1 0.001 -.000933
> :CHARGE convergence:  1 0.001 -.000972
> :CHARGE convergence:  1 0.001 -.000984
> :CHARGE convergence:  1 0.001 -.000988
> :CHARGE convergence:  1 0.001 -.000997
> :CHARGE convergence:  1 0.001 -.000998
>
> This is quite large calculation of TiSiO amorphous-like cell (96 atoms,
>  struct file attached, compressed in order to pass the moderation
> queue). The numerical parameters are quite standard,
> 2x2x2 shifted k-grid (can also reproduce this with 3x3x3), Rkmax 7.5,
> PBE, rest of stuff at default. I have done a lot of calculation with
> similar cells and such parameters without any problems (albeit mostly
> some time ago with older Wien2k versions)...
>
> I'm trying to calculate the binding energy of the O 1s for the first
> oxygen using the Slaters transition state approach. Surprisingly when I
> use almost touching spheres it converges (O:1.48, Ti: 1.72, Si:1.48),
> however when I reduce them by 3% (as in the attached struct file) this
> problems surfaces. However I need this reduced spheres for comparison
> with another calculation. BTW the spheres around 1.5 for oxygen
> should be reasonable anyway?
>
> No errors or warnings anywhere, except for the expected
> WARN  :CHARGED CELL with  -0.500
>
> This is with Wien2k 17.1, compiled with MKL 2017.0.098 and ifort.
>
> Any ideas what might be causing this problem?
> Best regards
> Pavel
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Laurence Marks]

I've posted it to the list several times. Peter added it to 17.1 as the
command "Check-mixing".

_
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu

On Dec 7, 2017 5:50 AM, "Pavel Ondračka"  wrote:

> Laurence Marks píše v Čt 07. 12. 2017 v 05:25 -0600:
> > I suggest using the "Check" command, which shows more useful
> > information and post the output. Without that it is hard to guess.
> >
> Dear prof. Marks,
>
> I'm not sure if I follow. Are you referring to the check (check_lapw)
> script? In my understanding it checks for running Wien2k instances,
> which are none in my case. Did you mean maybe analyse_lapw or some
> other check script which is not part of the default Wien2k
> installation? I can extract some data from the scf file by hand, if you
> let me know what info would be helpful in this case.
>
> Best regards
> Pavel
> ___
> Wien mailing list
> Wien@zeus.theochem.tuwien.ac.at
> https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.
> theochem.tuwien.ac.at_mailman_listinfo_wien=DwIGaQ=
> yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws=U_
> T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0=wmNIgGC9279w-
> ypk9mxq7YkSo3PDiJjxzomeLzPBKgI=Z5WBRqGtloTIO8Gojn6KlHimo55Fnm
> jIuuPgU2aG1LI=
> SEARCH the MAILING-LIST at:  https://urldefense.proofpoint.
> com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-
> 40zeus.theochem.tuwien.ac.at_index.html=DwIGaQ=
> yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws=U_
> T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0=wmNIgGC9279w-
> ypk9mxq7YkSo3PDiJjxzomeLzPBKgI=eEQ6FIg4tayKydo7w7T4rGnOdWsW43
> _MvuzYjcJn6ZU=
>
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Pavel Ondračka]

Laurence Marks píše v Čt 07. 12. 2017 v 05:25 -0600:
> I suggest using the "Check" command, which shows more useful
> information and post the output. Without that it is hard to guess.
> 
Dear prof. Marks,

I'm not sure if I follow. Are you referring to the check (check_lapw)
script? In my understanding it checks for running Wien2k instances,
which are none in my case. Did you mean maybe analyse_lapw or some
other check script which is not part of the default Wien2k
installation? I can extract some data from the scf file by hand, if you
let me know what info would be helpful in this case.

Best regards
Pavel
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] problem while running "x w2w" command

[Gavin Abo]

SIGSEGV errors can have many possible causes [ 
https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14476.html 
].


In your case, I see W_14.2 below.  So perhaps it is due to a bug in an 
the older WIEN2k 14.2 or wien2wannier version that you are using.  You 
might try the newer wien2wannier version in WIEN2k 17.1.  However, you 
have to patch it [ 
https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg16099.html 
].


What ifort compiler version are you using?  If it is a 2016 and 2017, 
maybe it is caused by the IO problem with large files [ 
https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg16802.html 
].


On 12/6/2017 3:58 AM, pachineela rambabu wrote:

While running "x w2w" command, I encountered following error:

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image  PC    Routine Line    Source
w2wc   004EC771  Unknown Unknown  Unknown
w2wc   004EA8AB  Unknown Unknown  Unknown
w2wc   004935B4  Unknown Unknown  Unknown
w2wc   004933C6  Unknown Unknown  Unknown
w2wc   0043AFE9  Unknown Unknown  Unknown
w2wc   0043EBF6  Unknown Unknown  Unknown
libpthread-2.23.s  2AD859B82390  Unknown Unknown  Unknown
w2wc   00465AF7  Unknown Unknown  Unknown
w2wc   00426F83 read_vec_  52  
read_vec.f

w2wc   00420F70 MAIN__    159  main.f
w2wc   00402FAE  Unknown Unknown  Unknown
libc-2.23.so  2AD85A65D830  
__libc_start_main Unknown  Unknown

w2wc   00402EA9  Unknown Unknown  Unknown
0.3u 0.2s 0:01.03 56.3% 0+0k 88016+16io 0pf+0w
error: command   /home/kanchana/Soft/W_14.2/w2wc w2w.def   failed


Please tell me how to rectify this problem.  For cubic systems, I 
didn't get any problem. For hexagonal system, this problem occurred.


thanks in advance,
--
*P. Rambabu*
PhD Scholor
Department of Physics, IIT Hyderabad, India
Mobile: 9074508220.
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Peter Blaha]

As the charge distance is not going bad, I don't expect it has to do 
with mixing, nor with any scf problem.


I'd guess, only E-tot is calculated wrongly in mixer or the run script 
has a problem .
It could be that you have 2 times the same iteration number in the scf 
file (:ITE) or that for some reasons once contribution to :ENE is not 
listed/listed twice, ...


Please check the scf file for :ITE, but also manually for the different 
E-tot contributions like (grep :ENE case.scf)


:ENE
:DEN
:SUM
and the that the core eigenvalues are listed only once.

ALSO: save the calculation, and continue with another run_lapw
Eventually it should stop after 3 cycles, because everything is converged.


On 12/07/2017 11:20 AM, Pavel Ondračka wrote:

Dear Wien2k mailing list,

I'm seeing some strange energy oscillations in one calculation near the
end when everything seems almost converged.

:ENERGY convergence:  0 0.0001 .00390097
:ENERGY convergence:  0 0.0001 .016890645000
:ENERGY convergence:  0 0.0001 .00965980
:ENERGY convergence:  0 0.0001 .016205335000
:ENERGY convergence:  0 0.0001 .00896832
:ENERGY convergence:  0 0.0001 55.639853345000
:ENERGY convergence:  0 0.0001 27.821044505000
:ENERGY convergence:  0 0.0001 27.82352526

the charge is converged already:

:CHARGE convergence:  1 0.001 -.000920
:CHARGE convergence:  1 0.001 -.000977
:CHARGE convergence:  1 0.001 -.000933
:CHARGE convergence:  1 0.001 -.000972
:CHARGE convergence:  1 0.001 -.000984
:CHARGE convergence:  1 0.001 -.000988
:CHARGE convergence:  1 0.001 -.000997
:CHARGE convergence:  1 0.001 -.000998

This is quite large calculation of TiSiO amorphous-like cell (96 atoms,
 struct file attached, compressed in order to pass the moderation
queue). The numerical parameters are quite standard,
2x2x2 shifted k-grid (can also reproduce this with 3x3x3), Rkmax 7.5,
PBE, rest of stuff at default. I have done a lot of calculation with
similar cells and such parameters without any problems (albeit mostly
some time ago with older Wien2k versions)...

I'm trying to calculate the binding energy of the O 1s for the first
oxygen using the Slaters transition state approach. Surprisingly when I
use almost touching spheres it converges (O:1.48, Ti: 1.72, Si:1.48),
however when I reduce them by 3% (as in the attached struct file) this
problems surfaces. However I need this reduced spheres for comparison
with another calculation. BTW the spheres around 1.5 for oxygen
should be reasonable anyway?

No errors or warnings anywhere, except for the expected
WARN  :CHARGED CELL with  -0.500

This is with Wien2k 17.1, compiled with MKL 2017.0.098 and ifort.

Any ideas what might be causing this problem?
Best regards
Pavel



___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html



--

  P.Blaha
--
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982
Email: bl...@theochem.tuwien.ac.atWIEN2k: http://www.wien2k.at
WWW:   http://www.imc.tuwien.ac.at/TC_Blaha
--
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Pavel Ondračka]

Peter Blaha píše v Čt 07. 12. 2017 v 13:34 +0100:
> As the charge distance is not going bad, I don't expect it has to do 
> with mixing, nor with any scf problem.
> 
> I'd guess, only E-tot is calculated wrongly in mixer or the run
> script 
> has a problem .
> It could be that you have 2 times the same iteration number in the
> scf 
> file (:ITE) or that for some reasons once contribution to :ENE is
> not 
> listed/listed twice, ...
> 
> Please check the scf file for :ITE, but also manually for the
> different 
> E-tot contributions like (grep :ENE case.scf)
> 
> :ENE
> :DEN
> :SUM
> and the that the core eigenvalues are listed only once.
> 
> ALSO: save the calculation, and continue with another run_lapw
> Eventually it should stop after 3 cycles, because everything is
> converged.
> 
Thanks for the advice,
so the :ITE check did not spot any problems, the :SUM also stays the
same, however the :DEN show the same fluctuations as the total energy,
like this...

:DEN  : DENSITY INTEGRAL  =-15670.17810081   (Ry)
:DEN  : DENSITY INTEGRAL  =-15614.53403108   (Ry)
:DEN  : DENSITY INTEGRAL  =-15614.55064266   (Ry)
:DEN  : DENSITY INTEGRAL  =-15614.55963292   (Ry)
:DEN  : DENSITY INTEGRAL  =-15670.17573346   (Ry)
:DEN  : DENSITY INTEGRAL  =-15670.17527628   (Ry)
:DEN  : DENSITY INTEGRAL  =-15670.17574871   (Ry)


I'm not sure how to check that "core eigenvalues are listed only once"
but there is nothing suspicious at the first glance.

Also attached is the output from the Check-mixed as requested by prof.
Marks.

I saved the calculations as advised and restarted. However I've done
this already several times (I've also tried clean_lapw) without any
effect.

Best regards
Pavel:DIS  :  CHARGE DISTANCE   (  0.02 for atom   68 spin 1)   0.00
:PLANE:  INTERSTITIAL TOTAL 10.52178 RMS  6.327E-04 , 3.154E-05 %
:CHARG:  CLM CHARGE   /ATOM  4.74789 RMS  1.320E-06 , 3.267E-06 %
:RANK :  ACTIVE   6.00/7  =  85.76 %
:DIRM :  MEMORY  7/8  RESCALE3.98 RED  0.53 PRED  0.92 NEXT  0.81 
:DIRP :  |MSR1|= 8.401E-06 |PRATT|= 1.321E-05 ANGLE=   7.6 DEGREES
:DIRQ :  |MSR1|= 7.540E-06 |PRATT|= 1.489E-05 ANGLE=  23.7 DEGREES
:DIRT :  |MSR1|= 1.129E-05 |PRATT|= 1.990E-05 ANGLE=  18.7 DEGREES
:MIX  :   MSD1   REGULARIZATION: 4.15E-04  GREED: 0.3000  Newton 1.00 0.567 
   
:ENE  : *WARNING** TOTAL ENERGY IN Ry =   -34887.62081542
:DIS  :  CHARGE DISTANCE   (  0.01 for atom4 spin 1)   0.00
:PLANE:  INTERSTITIAL TOTAL 10.52178 RMS  2.231E-04 , 1.112E-05 %
:CHARG:  CLM CHARGE   /ATOM  4.74789 RMS  5.877E-07 , 1.455E-06 %
:DIRM :  MEMORY  0/ 8 RESCALE1.00
:MIX  :   PRATT  REGULARIZATION: 1.00E-04  GREED: 0.1500
:ENE  : *WARNING** TOTAL ENERGY IN Ry =   -34887.63749053
:DIS  :  CHARGE DISTANCE   (  0.00 for atom4 spin 1)   0.00
:PLANE:  INTERSTITIAL TOTAL 10.52178 RMS  1.854E-04 , 9.245E-06 %
:CHARG:  CLM CHARGE   /ATOM  4.74789 RMS  3.431E-07 , 8.492E-07 %
:DIRM :  MEMORY  1/ 8 RESCALE4.75
:MIX  :   PRATT  REGULARIZATION: 1.00E-04  GREED: 0.2508  Expand 1.67 0.251 
   
:ENE  : *WARNING** TOTAL ENERGY IN Ry =   -34887.64647620
:DIS  :  CHARGE DISTANCE   (  0.00 for atom   65 spin 1)   0.00
:PLANE:  INTERSTITIAL TOTAL 10.52178 RMS  1.643E-04 , 8.194E-06 %
:CHARG:  CLM CHARGE   /ATOM  4.74789 RMS  3.088E-07 , 7.644E-07 %
:DIRM :  MEMORY  0/ 8 RESCALE1.00
:MIX  :   PRATT  REGULARIZATION: 1.00E-04  GREED: 0.1500
:ENE  : *WARNING** TOTAL ENERGY IN Ry =   -34943.26257364
:DIS  :  CHARGE DISTANCE   (  0.00 for atom   66 spin 1)   0.00
:PLANE:  INTERSTITIAL TOTAL 10.52178 RMS  1.376E-04 , 6.860E-06 %
:CHARG:  CLM CHARGE   /ATOM  4.74789 RMS  2.105E-07 , 5.210E-07 %
:DIRM :  MEMORY  1/ 8 RESCALE3.65
:MIX  :   PRATT  REGULARIZATION: 1.00E-04  GREED: 0.2000  Expand 1.33 0.200 
   
:ENE  : *WARNING** TOTAL ENERGY IN Ry =   -34943.26210909
:DIS  :  CHARGE DISTANCE   (  0.00 for atom   65 spin 1)   0.00
:PLANE:  INTERSTITIAL TOTAL 10.52178 RMS  1.107E-04 , 5.520E-06 %
:CHARG:  CLM CHARGE   /ATOM  4.74789 RMS  1.898E-07 , 4.697E-07 %
:RANK :  ACTIVE   1.05/2  =  52.36 %
:DIRM :  MEMORY  2/8  RESCALE3.66 RED  0.85 PRED  1.00 NEXT  1.00 
:DIRP :  |MSR1|= 5.317E-07 |PRATT|= 2.126E-06 ANGLE=   0.0 DEGREES
:DIRQ :  |MSR1|= 0.000E+00 |PRATT|= 0.000E+00 ANGLE=   0.0 DEGREES
:DIRT :  |MSR1|= 7.544E-07 |PRATT|= 3.017E-06 ANGLE=   0.0 DEGREES
:MIX  :   MSD1   REGULARIZATION: 2.00E-04  GREED: 0.2493  Newton 1.00 0.250 
   
:ENE  : *WARNING** TOTAL ENERGY IN Ry =   -34943.26257358
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Laurence Marks]

Probably overlapping response from Peter...

Based upon the Check-mixing output, the density is very well converged and
did not change with the jump in energy. The most likely way for :DEN alone
to change, and nothing else, is that there was a change in how EXC is
calculated, for instance using a metaGGA or changing case.in0, without a
change in VXC.

Do "grep :POT *scf" to check this.

_
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu

On Dec 7, 2017 7:19 AM, "Pavel Ondračka"  wrote:

> Peter Blaha píše v Čt 07. 12. 2017 v 13:34 +0100:
> > As the charge distance is not going bad, I don't expect it has to do
> > with mixing, nor with any scf problem.
> >
> > I'd guess, only E-tot is calculated wrongly in mixer or the run
> > script
> > has a problem .
> > It could be that you have 2 times the same iteration number in the
> > scf
> > file (:ITE) or that for some reasons once contribution to :ENE is
> > not
> > listed/listed twice, ...
> >
> > Please check the scf file for :ITE, but also manually for the
> > different
> > E-tot contributions like (grep :ENE case.scf)
> >
> > :ENE
> > :DEN
> > :SUM
> > and the that the core eigenvalues are listed only once.
> >
> > ALSO: save the calculation, and continue with another run_lapw
> > Eventually it should stop after 3 cycles, because everything is
> > converged.
> >
> Thanks for the advice,
> so the :ITE check did not spot any problems, the :SUM also stays the
> same, however the :DEN show the same fluctuations as the total energy,
> like this...
>
> :DEN  : DENSITY INTEGRAL  =-15670.17810081   (Ry)
> :DEN  : DENSITY INTEGRAL  =-15614.53403108   (Ry)
> :DEN  : DENSITY INTEGRAL  =-15614.55064266   (Ry)
> :DEN  : DENSITY INTEGRAL  =-15614.55963292   (Ry)
> :DEN  : DENSITY INTEGRAL  =-15670.17573346   (Ry)
> :DEN  : DENSITY INTEGRAL  =-15670.17527628   (Ry)
> :DEN  : DENSITY INTEGRAL  =-15670.17574871   (Ry)
>
>
> I'm not sure how to check that "core eigenvalues are listed only once"
> but there is nothing suspicious at the first glance.
>
> Also attached is the output from the Check-mixed as requested by prof.
> Marks.
>
> I saved the calculations as advised and restarted. However I've done
> this already several times (I've also tried clean_lapw) without any
> effect.
>
> Best regards
> Pavel
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Re: [Wien] Large oscillations in total energy

[Laurence Marks]

Did the :ENE/:DEN change occur only after you changed the RMT? If so, then
it is to be expected although I am very surprised that :DIS & :PLAN did not
change a little.

If it occurred at constant RMT, RKMAX etc then my only explanation is that
somehow EC/EX changed. The charge/CLM/PW values only depend upon the
potential, whereas :DEN depends upon the exchange-correlation energy and
the potential.

Beyond this, I have no clue.

On Thu, Dec 7, 2017 at 8:44 AM, Pavel Ondračka 
wrote:

> Laurence Marks píše v Čt 07. 12. 2017 v 07:37 -0600:
> > Probably overlapping response from Peter...
> >
> > Based upon the Check-mixing output, the density is very well
> > converged and did not change with the jump in energy. The most likely
> > way for :DEN alone to change, and nothing else, is that there was a
> > change in how EXC is calculated, for instance using a metaGGA or
> > changing case.in0, without a change in VXC.
> >
> > Do "grep :POT *scf" to check this.
> >
> This is a standard PBE calculation:
> :POT  : POTENTIAL OPTION EX_PBE EC_PBE VX_PBE VC_PBEfor all the cycles.
> In fact the surprising thing was that this started out of the blue
> after at least 15 iterations when the calculation was already looking
> quite converged.
>
> BTW as mentioned in the beginning this might have to do something with
> the atomic spheres, I converged this structure with the almost touching
> spheres (0% reduction) and later I did reduce_rmt_lapw -r 3.5 to
> produce this case, hence it was already pre-converged and I only wanted
> to fine-converge it, however the energy oscillations showed...
> Dunno of this might play a role.
>
> Best regards
> Pavel
> ___
> Wien mailing list
> Wien@zeus.theochem.tuwien.ac.at
> https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.
> theochem.tuwien.ac.at_mailman_listinfo_wien=DwIGaQ=
> yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws=U_
> T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0=p3o__
> T45kxKCNwRjQiKGvwqlO8WqCRN2C_c3f-a_GY8=57S8TTdUmZlCdeWGChaGOK-
> 5e5rVCXhwwJaycX--yPs=
> SEARCH the MAILING-LIST at:  https://urldefense.proofpoint.
> com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-
> 40zeus.theochem.tuwien.ac.at_index.html=DwIGaQ=
> yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws=U_
> T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0=p3o__
> T45kxKCNwRjQiKGvwqlO8WqCRN2C_c3f-a_GY8=tGAFPo0AFcaVSdbW8jdi9HkUAUF53A
> DdevPByqcdlak=
>



-- 
Professor Laurence Marks
"Research is to see what everybody else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu ; Corrosion in 4D: MURI4D.numis.northwestern.edu
Partner of the CFW 100% program for gender equity, www.cfw.org/100-percent
Co-Editor, Acta Cryst A
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html