<div dir="ltr">You cannot minimize with any of the canned routines available with -so as the forces are wrong.<div><br></div><div>To my knowledge the "NOSE" option in Wien2k is more than 10 (perhaps 20) years old, and I am not sure it works. I doubt that it has been used for decades. If you want to do MD I suspect that other codes such as QE are much better, as they focus on that.</div><div><br></div><div>In principle MSR1a could be modified to do MD, but this would take about 1 person year and I have not heard of any volunteers.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 17, 2015 at 11:50 PM, Gavin Abo <span dir="ltr"><<a href="mailto:gsabo@crimson.ua.edu" target="_blank">gsabo@crimson.ua.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I haven't tried MD with min_lapw.<br>
<br>
The case.finM should be created by the script min_lapw [<br>
<a href="http://zeus.theochem.tuwien.ac.at/pipermail/wien/2006-May/007135.html" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/pipermail/wien/2006-May/007135.html</a> ].<br>
<br>
One of the other users might know better, but it looks like your<br>
"min_lapw -p -so" command does not look right.<br>
<br>
I suggest you run in the terminal the command: min_lapw -h<br>
<br>
I'm not seeing an option for '-so' in the output from that command.<br>
<br>
Since you did not use '-j', it looks like the script uses "run_lapw -I<br>
-fc 1. -i 40" by default (whereas '-sp' should use "runsp_lapw -I -fc 1.<br>
-i 40" by default).<br>
<br>
For SO, you probably need to use the "-j" option:<br>
<br>
min_lapw -p -j "runsp_lapw -I -so -fc 1. -i 40"<br>
<br>
or it looks like there is the '-m' option that will extract force-input<br>
and execute mini (without JOB) and exit.<br>
<br>
On 10/23/2003 8:03 PM, Guo-ping Zhang wrote:<br>
> Dear Wien2k users,<br>
><br>
> I am interested in using Mini to run MD. After reading the manual and<br>
> old mailing list, I could not find a workable example. For instance,<br>
> the manual has an example for NOSE (case.inM) not for MOLD. In the<br>
> structure optimization section,  all the examples are used for<br>
> minimization (which works well for my problem using MSR1a option in<br>
> case.inm though the mixer is very slow).<br>
><br>
> I tried to use min_lapw -p -so (after a prior successful SCF)  but<br>
> failed, since it needs case.finM (but the manual does not say how to<br>
> get case.finM). I already looked at haupt_cad.f, which seems to work.<br>
><br>
> Has anyone succeeded using Wien2k to run MD?<br>
><br>
> If so, is it possible to share with me some details how this is done?<br>
><br>
> I would appreciate any hints.<br>
><br>
><br>
> Thanks a lot in advance!<br>
><br>
> Guoping<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</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>