<div dir="ltr">Yes....and no.<div><br></div><div>In 14.2 Peter added a "Push" for the 1st iteration if :DIS is small, some lines starting with</div><div><div>#check :dis of first cycle and rerun mixer with increased mixing</div></div><div>after the execution of mixer</div><div><br></div><div>This was added to avoid stagnation where it just takes a few small steps and stops.</div><div><br></div><div>In your case this "disturbed" 14.2 and it is adjusting the greed down to try and stabilize. In 13.1 this push is missing, so it starts with a small greed and finds what it thinks might be a solution.</div><div><br></div><div>The next version does not have this push (it is done inside mixer), but still needs work.</div><div><br></div><div>Unfortunately both cases are really bad -- :FCheck should be close to 1.0 or less, and your values indicate something really wrong. For reasons I don't understand MSR1a does not behave well with non-centro symmetric structures. I don't know why this is.</div><div><br></div><div>N.B., MSR1a is not perfect, and there are cases where it does not want to behave well. Whenever the model is susceptible to ghostbands it is not ideal. I am currently trying a few things, but.....</div><div><br></div><div>N.N.B., Peter's explanation of why :FCHECK can be large is that the linearization of different atoms is unbalanced. Maybe he has ideas about how to cure this, it is beyond me.</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 30, 2015 at 9:00 AM, Elias Assmann <span dir="ltr"><<a href="mailto:elias.assmann@gmail.com" target="_blank">elias.assmann@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi List,<br>
<br>
I reply to this thread<br>
<<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg13276.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg13276.html</a>>.<br>
This is still about the same set of oxide calculations.<br>
<br>
Laurence Marks alerted me that :FCHECK (sum of forces) was large (see<br>
below).  This is the same over a set of similar calculations (somewhat<br>
different structures, GGA+U vs. plain GGA), and does not seem to get<br>
better even after hundreds of calculations.<br>
<br>
I have been using the latest Wien2k (14.2), but I did similar things<br>
before with 13.1 and I did not remember having this problem, at least<br>
to this extent.  So I compiled 13.1 (same options as 14.2) and<br>
switched one of my calculations to it for 10 iterations.<br>
<br>
Indeed, it seems that 13.1 is more â€œwell-behaved” in this case than<br>
14.2 â€” in any case they are different.  Below, I include some output<br>
of the two versions, starting from the same charge density in both<br>
cases.<br>
<br>
As you can see, the two versions are the more or less the same in the<br>
first iteration, but diverge from there; and 13.1 seems to converge,<br>
in contrast to 14.2.  Also, 13.1 seems to be already satisfied with<br>
the structure (it switches MSR1a off), while 14.2 continues<br>
minimizing.  The greed (on whose smallness Prof. Marks also commented)<br>
is pretty similar.<br>
<br>
To me, that suggests that I should go on with these calculations using<br>
the older Wien.  Is there any reason that doing so might be â€œunsafe”?<br>
<br>
<br>
  Â  Â  Â  Elias<br>
<br>
<br>
14.2 :FCHECK:  Â  Â 11.008001043  Â  12.164812176  Â -14.416571844<br>
  Â  Â :FCHECK:  Â  249.814119002  Â 264.664202568  -694.164686632<br>
  Â  Â :FCHECK:  Â  Â 97.292437031  Â 103.207655853  -256.561081883<br>
  Â  Â :FCHECK:  Â  126.184534662  Â 136.666423509  Â -39.224005265<br>
  Â  Â :FCHECK:  Â  176.349427337  Â 188.815368048  Â  36.914398561<br>
  Â  Â :FCHECK:  Â  Â 50.400659984  Â  67.046955772  Â  76.978472788<br>
  Â  Â :FCHECK:  Â  -74.640170442  Â  25.157628216  Â 121.498217602<br>
  Â  Â :FCHECK:  Â  -82.337983720  Â  28.388321933  Â 121.952351989<br>
  Â  Â :FCHECK:  Â -130.382154150  Â  24.859863620  Â 204.399962720<br>
  Â  Â :FCHECK:  Â  -39.381668995  Â -31.379111345  Â  91.507308943<br>
<br>
13.1 :FCHECK:  Â  Â 11.009473749  Â  12.153718710  Â -14.482142314<br>
  Â  Â :FCHECK:  Â  Â 10.707248409  Â  11.500977675  Â -16.514184574<br>
  Â  Â :FCHECK:  Â  Â 10.289877753  Â  10.742969135  Â -10.026151044<br>
  Â  Â :FCHECK:  Â  Â  7.566779454  Â  Â 6.128096382  Â  -3.725471547<br>
  Â  Â :FCHECK:  Â  Â  3.572302068  Â  Â 1.087140350  Â  Â 6.011501945<br>
  Â  Â :FCHECK:  Â  Â -4.275691705  Â  -7.052038873  Â  31.588782325<br>
  Â  Â :FCHECK:  Â  Â -6.385420320  Â  -8.403344519  Â  33.874614737<br>
  Â  Â :FCHECK:  Â  Â -6.106651082  Â  -8.096761064  Â  35.566398852<br>
  Â  Â :FCHECK:  Â  Â -5.627436387  Â  -7.565248768  Â  32.843232377<br>
  Â  Â :FCHECK:  Â  Â -3.074299958  Â  -4.674255384  Â  27.987425362<br>
<br>
14.2 :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46215529<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46926026<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46462435<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46263442<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46223935<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46214440<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46218336<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46213858<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46236760<br>
  Â  Â :ENE  : *WARNING** TOTAL ENERGY  Â  Â  -222784.46181283<br>
<br>
13.1 :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46212826<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46212374<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46212843<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46207306<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46209495<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46204241<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46199458<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46200681<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46200118<br>
  Â  Â :ENE  : ********** TOTAL ENERGY  Â  Â  -222784.46201120<br>
<br>
(It appears that the warnings come from the forces.  As mentioned<br>
before there is not :WAR tag.)<br>
<br>
14.2 :DIS  :  ( 0.0048283 for atom  Â  7 spin 1)  Â  Â  0.0015628<br>
  Â  Â :DIS  :  ( 0.1175699 for atom  Â  5 spin 1)  Â  Â  0.0137938<br>
  Â  Â :DIS  :  ( 0.0692919 for atom  Â  5 spin 1)  Â  Â  0.0083222<br>
  Â  Â :DIS  :  ( 0.0098508 for atom  Â  7 spin 1)  Â  Â  0.0017036<br>
  Â  Â :DIS  :  ( 0.0167530 for atom  Â  5 spin 1)  Â  Â  0.0022822<br>
  Â  Â :DIS  :  ( 0.0278296 for atom  Â  7 spin 1)  Â  Â  0.0038446<br>
  Â  Â :DIS  :  ( 0.0366785 for atom  Â  5 spin 1)  Â  Â  0.0048138<br>
  Â  Â :DIS  :  ( 0.0371346 for atom  Â  5 spin 1)  Â  Â  0.0047958<br>
  Â  Â :DIS  :  ( 0.0559135 for atom  Â  5 spin 1)  Â  Â  0.0068198<br>
  Â  Â :DIS  :  ( 0.0157960 for atom  Â  5 spin 1)  Â  Â  0.0035943<br>
<br>
13.1 :DIS  :  ( 0.0048362 for atom  Â  7 spin 1)  Â  Â  0.0015593<br>
  Â  Â :DIS  :  ( 0.0022134 for atom  Â 12 spin 2)  Â  Â  0.0012247<br>
  Â  Â :DIS  :  ( 0.0052633 for atom  Â  5 spin 1)  Â  Â  0.0014252<br>
  Â  Â :DIS  :  ( 0.0015732 for atom  Â  5 spin 1)  Â  Â  0.0007705<br>
  Â  Â :DIS  :  ( 0.0046946 for atom  Â  7 spin 1)  Â  Â  0.0010640<br>
  Â  Â :DIS  :  ( 0.0018880 for atom  Â  8 spin 2)  Â  Â  0.0011590<br>
  Â  Â :DIS  :  ( 0.0017283 for atom  Â  5 spin 1)  Â  Â  0.0008169<br>
  Â  Â :DIS  :  ( 0.0009736 for atom  Â  8 spin 2)  Â  Â  0.0006944<br>
  Â  Â :DIS  :  ( 0.0018888 for atom  Â  5 spin 1)  Â  Â  0.0007013<br>
  Â  Â :DIS  :  ( 0.0009991 for atom  Â  7 spin 1)  Â  Â  0.0003925<br>
<br>
14.2 :MIX  :  Â PRATT  REG: 1.00E-06  GREED: 0.200<br>
  Â  Â :MIX  :  Â PRATT  REG: 1.00E-06  GREED: 0.200  Reduce 0.25  0.20<br>
  Â  Â :MIX  :  Â MSE1a  REG: 1.88E-06  GREED: 0.120  Newton 1.00  0.01<br>
  Â  Â :MIX  :  Â MSE1a  REG: 2.17E-06  GREED: 0.083  Newton 1.00  0.01<br>
  Â  Â :MIX  :  Â MSE1a  REG: 3.19E-06  GREED: 0.050  Newton 1.00  0.01<br>
  Â  Â :MIX  :  Â MSE1a  REG: 4.08E-06  GREED: 0.034  Newton 1.00  0.03<br>
  Â  Â :MIX  :  Â MSE1a  REG: 5.17E-06  GREED: 0.022  Newton 1.00  0.00<br>
  Â  Â :MIX  :  Â MSE1a  REG: 5.20E-06  GREED: 0.021  Newton 1.00  0.06<br>
  Â  Â :MIX  :  Â MSE1a  REG: 7.51E-06  GREED: 0.018  Newton 1.00  0.14<br>
  Â  Â :MIX  :  Â MSE1a  REG: 4.52E-06  GREED: 0.017  Newton 1.00  0.10<br>
<br>
13.1 :MIX  :  Â PRATT  REG: 1.00E-06  GREED: 0.025<br>
  Â  Â :MIX  :  Â MSR1a  REG: 1.00E-06  GREED: 0.030  Newton 1.00  0.03<br>
  Â  Â :MIX  :  Â MSE1a  REG: 1.59E-06  GREED: 0.030  Newton 1.00  0.18<br>
  Â  Â :MIX  :  Â MSE1a  REG: 2.61E-06  GREED: 0.042  Newton 1.00  0.29<br>
  Â  Â :MIX  :  Â MSE1a  REG: 3.72E-06  GREED: 0.057  Newton 1.00  0.54<br>
  Â  Â :MIX  :  Â MSE1a  REG: 3.58E-06  GREED: 0.093  Newton 1.00  0.46<br>
  Â  Â :MIX  :  Â PRATT  REG: 1.00E-06  GREED: 0.025<br>
  Â  Â :MIX  :  Â MSR1  Â REG: 1.00E-06  GREED: 0.050  Newton 1.00  0.05<br>
  Â  Â :MIX  :  Â MSE1  Â REG: 1.60E-06  GREED: 0.050  Newton 1.00  0.22<br>
  Â  Â :MIX  :  Â MSE1  Â REG: 2.41E-06  GREED: 0.055  Newton 1.00  0.40<br>
<br>
On 11/03/2015 02:47 PM, Laurence Marks wrote:<br>
> I am not sure exactly what the problem(s) were, but from you case.scf it<br>
> was "right" that there were Warnings<br>
><br>
> a) :FCHECK (bottom of case.scf) was large. This is the sum of all the<br>
> forces, and should be small. Particularly for cells without inversion<br>
> one can get bad, highly asymmetric densities in which case MSR1a can<br>
> have problems.<br>
><br>
> b) The greed is small. Too small a value can be as bad as too large. I<br>
> have struggled with this for years and failed to find a strong ansatz<br>
> for this, although I believe the next release of the mixer will be better.<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
--<br>
Elias Assmann<br>
Institute of Theoretical and Computational Physics<br>
TU Graz  Â  Â  Â  Â  Â  Â  Â  Â  Â âŸ¨<a href="https://itp.tugraz.at/" rel="noreferrer" target="_blank">https://itp.tugraz.at/</a>⟩<br>
<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" target="_blank">MURI4D.numis.northwestern.edu</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>
</div>