<div dir="ltr">There is an algorithmic error in MSR1a that effects release 14. Non centro-symmetric cases may converge poorly or not at all, and fixing an atom may not help. This is associated with an addition I made to try and avoid the atoms drifting across the cell which may fail.<div><br></div><div>Two, equivalent patches can be used:</div><div>a) A "soft" patch is to add the word "DRIFT" at the bottom of case.inm</div><div>b) A "hard" patch is to change in GetOpts.f the variable "NoDrift" to false, e.g. in line 32 have</div><div><div>       NoDrift  =.false.</div><div><br></div><div>Caveat. In rare cases it may happen that the contents of the unit cell will drift along one or more axes, because there is nothing to fix the origin along it. In that case I suggest editing by hand case.inM to fix the relevant atom displacement by changing a 1.0 to 0.0. For instance, to fix the z-component of atom 1 (if there is no inversion symmetry along z) the first line of case.inM would be changed from</div><div><br></div><div><div>1.0 1.0 1.0 1.0   #Atom    1 Generated by pairhess</div></div><div><br></div><div>to</div><div><br></div><div><div>1.0 1.0 0.0 1.0   #Atom    1 Generated by pairhess</div></div><div><br></div><div>N.B., in many cases this may not be an issue. If you are concerned that your calculation did not properly converge check the forces at the end of case.scf or case.scfm.</div><div><br></div><div>N.N.B., the technical explanation is that with the trap that was introduced the algorithm can converge to a spurious fixed-point solution cwhere all the forces are the same (along the unfixed direction), and :FCHECK is large.</div><div><br></div><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></div>