<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><P>Dear Laurence.</P>
<P>&nbsp;</P>
<P>I want to optimizar my structure with min_lapw. I read the instructions in the userguide, but is not clear because I'm new in Wien2k. I run first </P>
<P>&nbsp;</P>
<P>run_lapw -fc 1.0 &amp;&nbsp; and the results are OK.</P>
<P>&nbsp;</P>
<P>I make the case.inM </P>
<P>&nbsp;</P>
<P>PORT 1.0 0.25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #(PORT, NEWT, BFGS, MOLD, NOSE;&nbsp; tolf (a4,f5.2))<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom1 (NOSE, MOLD:Masse, delta t, T, nose-frequency)<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom2 (NEWT: 1,2,3:delta, 4:eta(1=MOLD))<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom3 (PORT,BFGS, NEWT: [ 1 || 2 || 3 ] = 0 constraint)<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom4 (PORT,BFGS: delta1,2,3=preconditioning (eta not used)<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom5<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom6<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom7<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom8<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom9<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom10<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom11<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom12<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom13<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom14<BR>0.0 0.0 0.0
 0.0&nbsp;&nbsp;&nbsp; # Atom15<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom16<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom17<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom18<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom19<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom20<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom21<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom22<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom23<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom24<BR>0..0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom25<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom26<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom27<BR>0.0 0.0 0.0 0..0&nbsp;&nbsp;&nbsp; # Atom28<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom29<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom30<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom31<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom32<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom33<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom34<BR>1.0 1.0 1.0 1.0&nbsp;&nbsp;&nbsp; # Atom35<BR>0.0 0.0 0.0
 0.0&nbsp;&nbsp;&nbsp; # Atom36<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom37<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom38<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom39<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom40<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom41<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom42<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom43<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom44<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom45<BR>0..0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom46<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom47<BR>0.0 0.0 0.0 0.0&nbsp;&nbsp;&nbsp; # Atom48<BR></P>
<P>I wanto optimizas some positions.</P>
<P>&nbsp;</P>
<P>before I make min_lapw &amp;</P>
<P>&nbsp;</P>
<P>after sometime the job stop. In the case..dayfile I have:</P>
<P>&nbsp;</P>
<P>&nbsp;</P>
<P>&nbsp;&nbsp;&nbsp; start&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (Mon Apr 14 14:36:39 AMT 2008) with mixer (40/20 to go)<BR>&gt;&nbsp;&nbsp; mixer&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (14:36:39) 51.370u 2.740s 0:54.20 99.8% 0+0k 0+0io 203pf+0w<BR>:ENERGY convergence:&nbsp; 0 0 0<BR>:CHARGE convergence:&nbsp; 0 0 .0000009</P>
<P>&nbsp;&nbsp;&nbsp; cycle 2&nbsp;&nbsp;&nbsp;&nbsp; (Mon Apr 14 14:37:36 AMT 2008)</P>
<P>&nbsp;&nbsp;&nbsp; cycle 2&nbsp;&nbsp;&nbsp;&nbsp; (39/19 to go)<BR>&gt;&nbsp;&nbsp; lapw0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (14:37:36) 1148.980u 1.680s 19:11.29 99.9%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0+0k 0+0io 263pf+0w<BR>:FORCE convergence:&nbsp; 1 1. 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 <BR>XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO 0 XCO<BR>&gt;&nbsp;&nbsp; lapw1&nbsp; -c&nbsp;&nbsp; (15:00:04) 5333.340u 12.350s 1:29:15.11 99.8%&nbsp;&nbsp; 0+0k 0+0io 311pf+0w<BR>&gt;&nbsp;&nbsp; lapw2 -c&nbsp;&nbsp;&nbsp; (16:29:19) 906.960u 11.740s 15:20.64 99.7%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0+0k 0+0io 308pf+0w<BR>&gt;&nbsp;&nbsp; lcore&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (16:44:40) 11.560u 0.260s 0:11.88 99.4% 0+0k 0+0io 169pf+0w<BR>&gt;&nbsp;&nbsp;
 mixer&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (16:44:54) 50.440u 1.720s 0:53.61 97.2% 0+0k 0+0io 203pf+0w<BR>:ENERGY convergence:&nbsp; 0 0 .0000120000000000<BR>:CHARGE convergence:&nbsp; 0 0 .0000018</P>
<P>&nbsp;&nbsp;&nbsp; cycle 3&nbsp;&nbsp;&nbsp;&nbsp; (Mon Apr 14 16:45:50 AMT 2008)</P>
<P>&nbsp;&nbsp;&nbsp; cycle 3&nbsp;&nbsp;&nbsp;&nbsp; (0/18 to go)</P>
<P>&gt;&nbsp;&nbsp; stop<BR>&gt;&nbsp;&nbsp; mini&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (16:45:58) 5.390u 0.030s 0:05.48 98.9%&nbsp; 0+0k 0+0io 176pf+0w</P>
<P>&gt;&nbsp;&nbsp; stop forces small, minimization stopped</P>
<P>&gt;&nbsp;&nbsp; relaxation steps summarized in bfesi2_vFeI_si.scf_mini</P>
<P>&gt;&nbsp;&nbsp; relaxed structure and results in bfesi2_vFeI_si.struct, scf, clmsum</P>
<P>&gt;&nbsp;&nbsp; please use&nbsp;&nbsp; save_lapw save_file_name&nbsp;&nbsp; to save the results<BR></P>
<P>&nbsp;</P>
<P>&nbsp;</P>
<P>I use&nbsp;&nbsp; save_lapw save_file_name&nbsp;&nbsp; to save the results,&nbsp;OK.</P>
<P>&nbsp;</P>
<P>But the old case struct is the&nbsp;same new case.struct, and still, in the case.outputnn after make x nn the near-neighbor is the same for old and new case.struct.</P>
<P>&nbsp;</P>
<P>I need your help.</P>
<P>&nbsp;</P>
<P>tanks</P>
<P>&nbsp;</P>
<P>Luiz Claudio</P>
<P>&nbsp;</P>
<P>&nbsp;</P>
<P>&nbsp;</P>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif"><BR><BR>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">----- Mensagem original ----<BR>De: Laurence Marks &lt;L-marks@northwestern.edu&gt;<BR>Para: A Mailing list for WIEN2k users &lt;wien@zeus.theochem.tuwien.ac.at&gt;<BR>Enviadas: Terça-feira, 8 de Abril de 2008 11:24:58<BR>Assunto: Re: [Wien] positions optimizations with min_lapw and Newt.<BR><BR>N.B., I just noticed that the optimization notes don't mention<BR>pairhess. Please look at the user guide for this code. (I will put<BR>updating the optimization notes onto my "to do" list.)<BR><BR>On Tue, Apr 8, 2008 at 9:22 AM, Laurence Marks &lt;<A href="mailto:L-marks@northwestern.edu" ymailto="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</A>&gt; wrote:<BR>&gt; You are approaching this wrong. Please run through the example in the<BR>&gt;&nbsp; user guide for optimizing the positions for TiO2 and also read the<BR>&gt;&nbsp; optimization FAQ at<BR>&gt;&nbsp; <A
 href="http://www.wien2k.at/reg_user/textbooks/Optimization-Notes.pdf" target=_blank>http://www.wien2k.at/reg_user/textbooks/Optimization-Notes.pdf</A><BR>&gt;<BR>&gt;<BR>&gt;<BR>&gt;&nbsp; On Tue, Apr 8, 2008 at 9:05 AM, Luiz Cláudio Carvalho<BR>&gt;&nbsp; &lt;<A href="mailto:lccfisica@yahoo.com.br" ymailto="mailto:lccfisica@yahoo.com.br">lccfisica@yahoo.com.br</A>&gt; wrote:<BR>&gt;&nbsp; &gt; Dear all,<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt; I´m new in Wien2k. I want to optimizar the positions in a unit cell with 48<BR>&gt;&nbsp; &gt; atoms. I want to using the option NEWT. I run_lapw -fc 1.0 &amp;. The result is<BR>&gt;&nbsp; &gt; Ok. Then I create the file case.inM and change the position of the atoms<BR>&gt;&nbsp; &gt; that I want to optimizar. I Have 3 questions that the userguide is not<BR>&gt;&nbsp; &gt; clear.<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt; 1) How I create the file case.finM?<BR>&gt;&nbsp; &gt; 2) After creating the file case.finM it is
 only to run min_lapw?<BR>&gt;&nbsp; &gt; 3) How to known if is it varying the wanted positions and converging<BR>&gt;&nbsp; &gt; correctly?<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt; Thanks for help<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt; Luiz Cláudio<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp; &gt;&nbsp; ________________________________<BR>&gt;&nbsp; &gt; Abra sua conta no Yahoo! Mail, o único sem limite de espaço para<BR>&gt;&nbsp; &gt; armazenamento!<BR>&gt;&nbsp; &gt; _______________________________________________<BR>&gt;&nbsp; &gt;&nbsp; Wien mailing list<BR>&gt;&nbsp; &gt;&nbsp; <A href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac..at</A><BR>&gt;&nbsp; &gt;&nbsp; <A href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target=_blank>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</A><BR>&gt;&nbsp; &gt;<BR>&gt;&nbsp;
 &gt;<BR>&gt;<BR>&gt;<BR>&gt;<BR>&gt;&nbsp; --<BR>&gt;&nbsp; Laurence Marks<BR>&gt;&nbsp; Department of Materials Science and Engineering<BR>&gt;&nbsp; MSE Rm 2036 Cook Hall<BR>&gt;&nbsp; 2220 N Campus Drive<BR>&gt;&nbsp; Northwestern University<BR>&gt;&nbsp; Evanston, IL 60208, USA<BR>&gt;&nbsp; Tel: (847) 491-3996 Fax: (847) 491-7820<BR>&gt;&nbsp; email: L-marks at northwestern dot edu<BR>&gt;&nbsp; Web: www.numis.northwestern.edu<BR>&gt;&nbsp; Commission on Electron Diffraction of IUCR<BR>&gt;&nbsp; www.numis.northwestern.edu/IUCR_CED<BR>&gt;<BR><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: www.numis.northwestern.edu<BR>Commission on Electron Diffraction of
 IUCR<BR>www.numis.northwestern.edu/IUCR_CED<BR>_______________________________________________<BR>Wien mailing list<BR><A href="mailto:Wien@zeus.theochem.tuwien.ac.at" ymailto="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" target=_blank>http://zeus.theochem.tuwien.ac..at/mailman/listinfo/wien</A><BR></DIV><BR></DIV></div><br>


between 0000-00-00 and 9999-99-99 
      <hr size=1>Abra sua conta no <a href="http://br.rd.yahoo.com/mail/taglines/mail/*http://br.mail.yahoo.com/">Yahoo! Mail</a>, o único sem limite de espaço para armazenamento! 
</body></html>