<div>Thanks very much for your reply!</div>
<div>&nbsp;</div>
<div>The structure file I&#39;m using is exactly the same as that one in the userguide:</div>
<div>&nbsp;</div>
<div>TiO2:</div>
<div>Lattice :</div>
<div>Spacegroup :&nbsp; 136_P42/mnm_</div>
<div>&nbsp;</div>
<div>Splitting of equivalent positions not available.</div>
<div>To split you must select a lattice type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&lt;--why? I thought the equivalent ones have already been splitted.</div>
<div>&nbsp;</div>
<div>Lattice constant:</div>
<div>a&nbsp; 8.682 bohr, b&nbsp; 8.682 bohr, c 8.682 bohr</div>
<div>alpha&nbsp;= beta = gamma = 90 degree</div>
<div>&nbsp;</div>
<div>Atom one : Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Z=22&nbsp;&nbsp;&nbsp; RMT = 2.00</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; Pos 1:&nbsp; x= 0&nbsp;&nbsp; y= 0&nbsp;&nbsp; z=0</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; Pos&nbsp;2:&nbsp; x=0.5 y=0.5&nbsp; z=0.5</div>
<div>&nbsp;</div>
<div>Atom two : O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Z=8.0&nbsp;&nbsp; RMT=1.6</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; Pos&nbsp; 1:&nbsp; x=0.3&nbsp;&nbsp; y=0.3&nbsp; z=0</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; Pos&nbsp; 2:&nbsp; x=0.7&nbsp;&nbsp;&nbsp;y=0.7&nbsp; z=0</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; Pos&nbsp; 3:&nbsp; x=0.2&nbsp;&nbsp; y=0.8&nbsp; z=0.5</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; Pos&nbsp;&nbsp;4:&nbsp; x=0.8&nbsp;&nbsp; y=0.2&nbsp; z=0.5</div>
<div>&nbsp;</div>
<div>After I saved the structure, I didn&#39;t press</div>
<div>&quot;set automatically RMT and continue editing&quot; because I want to follow the RMT value in the userguide. Then I cleaned up and initialize the calculation and so on... that&#39;s what I have done.</div>
<div>&nbsp;</div>
<div>I followed your suggestions and reduced the RMT :&nbsp; Ti : from 2.0--&gt;1.8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O:&nbsp; from 1.6--&gt; 1.5, however, when I initialize the calculation and right after the &quot;x lstart&quot; step, there is a warning saying:
</div>
<div>&nbsp;</div>
<div>WARNING !!!! For good atomic total energies you <br>Check try.outputst for which atom/states the core-leakage occurs and rerun lstart with lower core-seperation energy (or increase RMT)</div>
<div>&nbsp;</div>
<div>Should I increase the RMT again? Thanks a lot again!!!</div>
<div>&nbsp;</div>
<div>Best Wishes!</div>
<div>Sincerely yours,</div>
<div>&nbsp;</div>
<div>Bo Qiu</div>
<div><br>&nbsp;</div>
<div><span class="gmail_quote">On 10/14/07, <b class="gmail_sendername">Laurence Marks</b> &lt;<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Too large RMT (not too small, a typo). Either give a larger reduction<br>when using &quot;setrmt&quot; in the initialization or copy a file with smaller
<br>RMT&#39;s to TiO2.struct_new, use x clminter, then copy TiO2.struct_new to<br>TiO2.struct (and similarly for TiO2.clmsum_new to TiO2.clmsum) are<br>rerun.<br><br>On 10/14/07, <a href="mailto:jiang@fhi-berlin.mpg.de">jiang@fhi-berlin.mpg.de
</a> &lt;<a href="mailto:jiang@fhi-berlin.mpg.de">jiang@fhi-berlin.mpg.de</a>&gt; wrote:<br>&gt; Hi Bo,<br>&gt; I think the problem in your calculations is probably caused by too small RMT<br>&gt; values in the struct file. The muffin-tin spheres should never overlap
<br>&gt; when you optimize the structure. Try reducing RMT value a little bit to<br>&gt; see what happens.<br>&gt;<br>&gt; Good luck!<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;Hong<br>&gt;<br>&gt; &gt; Dear Wien2k users and developers,<br>&gt; &gt;
<br>&gt; &gt; I&#39;m a new user and I&#39;ve been struggling with this problem and find no way<br>&gt; &gt; out for several days and I&#39;m really in despair. I sincerely request for<br>&gt; &gt; your<br>&gt; &gt; help as my last hope, thanks a lot!
<br>&gt; &gt;<br>&gt; &gt; The problem is that:<br>&gt; &gt; When I&#39;m trying to do the optimization with the example TiO2 and try to<br>&gt; &gt; optimize its atomic positions using min_lapw. Firstly I tried to follow
<br>&gt; &gt; the<br>&gt; &gt; userguide in ch 10.3 exactly using w2web, and run scf first. As pointed in<br>&gt; &gt; the userguide, the last force was a jump. Then I proceeded to mini<br>&gt; &gt; positions, however, first of all, before every min_lapw round there is a
<br>&gt; &gt; warning saying:<br>&gt; &gt;<br>&gt; &gt; WARNING: Step size reduced due to overlapping spheres -- check RMT<br>&gt; &gt;<br>&gt; &gt; after several (around 10) rounds, PORT crashed:<br>&gt; &gt;<br>&gt; &gt;&gt;&gt;&nbsp;&nbsp;(mini) PORT stopped -&gt; exit
<br>&gt; &gt;&nbsp;&nbsp;:WARNING you might have inconsistent Forces/Energies<br>&gt; &gt;<br>&gt; &gt; The last few lines on FGL002 inside *mini are (FGL001 is alway zero):<br>&gt; &gt;<br>&gt; &gt; :FGL002:&nbsp;&nbsp; 2.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -
45.621114095&nbsp;&nbsp; -45.621114095<br>&gt; &gt; 0.000000000 t<br>&gt; &gt; :FGL002:&nbsp;&nbsp; 2.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -45.586465863&nbsp;&nbsp; -45.586465863<br>&gt; &gt; 0.000000000 t<br>&gt; &gt; :FGL002:&nbsp;&nbsp; 2.ATOM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -45.575859262&nbsp;&nbsp; -45.575859262
<br>&gt; &gt; 0.000000000 t<br>&gt; &gt;<br>&gt; &gt; Then I removed case.tmpM and modified case.inM as<br>&gt; &gt;<br>&gt; &gt; NEW1 2.0 0.35&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;# PORT/NEWT;&nbsp;&nbsp;tolf, Initial Trust Radius<br>&gt; &gt; 0.0 0.0 0.0 1.0&nbsp;&nbsp; #Atom&nbsp;&nbsp;1 Generated by pairhess
<br>&gt; &gt; 1.0 1.0 0.0 1.0&nbsp;&nbsp; #Atom&nbsp;&nbsp;2 Generated by pairhess<br>&gt; &gt;<br>&gt; &gt; Then rerun min_lapw using w2web<br>&gt; &gt;<br>&gt; &gt; however, this time, in the first round, after 4 cycles I get it<br>&gt; &gt; crashed with NN -Error:
<br>&gt; &gt;<br>&gt; &gt; in cycle 4&nbsp;&nbsp;&nbsp;&nbsp;ETEST: 0&nbsp;&nbsp; CTEST: 0<br>&gt; &gt;&nbsp;&nbsp;LAPW0 END<br>&gt; &gt;&nbsp;&nbsp;LAPW1 END<br>&gt; &gt;&nbsp;&nbsp;LAPW2 END<br>&gt; &gt;&nbsp;&nbsp;CORE&nbsp;&nbsp;END<br>&gt; &gt;&nbsp;&nbsp;MIXER END<br>&gt; &gt; ec cc and fc_conv 1 1 1<br>
&gt; &gt;<br>&gt; &gt;&gt;&nbsp;&nbsp; stop<br>&gt; &gt; NN - Error<br>&gt; &gt;<br>&gt; &gt; The version used: Wien2k_07, compiler ifort, mkl 9.1, in series mode<br>&gt; &gt;<br>&gt; &gt; Could you please help me? Thank you very much!!! If any more information
<br>&gt; &gt; needed, please let me know, thanks!<br>&gt; &gt;<br>&gt; &gt; --<br>&gt; &gt; Thanks and Best Wishes.<br>&gt; &gt;<br>&gt; &gt; Sincerely yours,<br>&gt; &gt; Bo Qiu<br>&gt; &gt;<br>&gt; &gt; <a href="http://web.ics.purdue.edu/~qiub/">
http://web.ics.purdue.edu/~qiub/</a> &lt;<a href="http://web.ics.purdue.edu/%7Eqiub/">http://web.ics.purdue.edu/%7Eqiub/</a>&gt;<br>&gt; &gt; _______________________________________________<br>&gt; &gt; Wien mailing list<br>
&gt; &gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt; &gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
</a><br>&gt; &gt;<br>&gt;<br>&gt; _______________________________________________<br>&gt; Wien mailing list<br>&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>&gt;<br><br><br>--<br>Laurence Marks<br>Department of Materials Science and Engineering<br>MSE Rm 2036 Cook Hall<br>2220 N Campus Drive<br>Northwestern University
<br>Evanston, IL 60208, USA<br>Tel: (847) 491-3996 Fax: (847) 491-7820<br>email: L-marks at northwestern dot edu<br>Web: <a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a><br>Commission on Electron Diffraction of IUCR
<br><a href="http://www.numis.northwestern.edu/IUCR_CED">www.numis.northwestern.edu/IUCR_CED</a><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">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br></blockquote></div><br><br clear="all">