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