[Wien] ‘mixer’ crashes when forces too large

[Laurence Marks]

In principle this is a bug in either mixer or sumpara - we did not expect
FSU to be so large. Unless you are using MSR1a you probably need to change
FOR to TOT at the top of case.in2. If you are using mini there is an option
to switch the forces off for the start of the next run; please check the UG.

If you are using MSR1a you may need to turn it off for some cycles to get a
better density, then do MSR1a again (removing broyd files as needed).

Last, but not least, you might have what I called "ghost forces". These are
due to (probably) something going wrong with the physical model or
linearization energies. Exactly what, I dont know.
On May 24, 2013 10:59 AM, "Elias Assmann" <elias.assmann at gmail.com> wrote:

> A quick follow-up on my last e-mail.
>
> On 05/24/2013 01:12 PM, Elias Assmann wrote:
> > Upon investigation, I found that the cause of the crash was a line
> >
> >     :FSU033:  33.ATOM 120503.840520617-118793.163618910  -4395.235042983
> >   19749.475095455
> >
> > in the ‘scf’ file, which it had tried to read as [scfana.f:365]
> >
> >     READ(MARGN(18:120),*)FHELP(1:4) !READ(MARGN,740) (FHELP(j),j=1,4)
>
> >   * change scfana.f to formatted input,
> >
> >   * change whatever writes this line to list-directed output, or even
> >
> >   * emit an error when the forces are so large (?)
>
> As I said, I changed the ‘scf’ file (essentially the second option
> above), then ‘x mixer’ ran through.  However, the problem reappeared in
> the next iteration.
>
>         Elias
>
> _______________________________________________
> Wien mailing list
> Wien at 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20130524/96c263ef/attachment.htm>