please vary the % of change or volume -10,-5.0,5,10 like this.If your systeem is optimised one i think lattice constants will not change too much..<br><br>
<div class="gmail_quote">On Wed, Jul 9, 2008 at 6:25 PM, zhs064 &lt;<a href="mailto:zhs064@mail.usask.ca">zhs064@mail.usask.ca</a>&gt; wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Dear Wien2k users,<br><br>&nbsp; &nbsp; I&#39;m working on Ca_cmca (space group is 64) at high pressure. I got<br>all cell constants<br>
and atomic positions from the optimized values generated by the other<br>program. The initial structure file looks like:<br>**************************************************************<br>Title Ca &nbsp; &nbsp; &nbsp; &nbsp; cmca<br>CXY &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1<br>
MODE OF CALC=RELA<br>&nbsp;8.435170 &nbsp;8.861303 &nbsp;8.446130 90.000000 90.000000 90.000000<br>ATOM &nbsp;-1: X=0.00000000 Y=0.34493000 Z=0.19252000<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;MULT= 4 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ISPLIT=<br>&nbsp; &nbsp; &nbsp;-1: X=0.00000000 Y=0.15507000 Z=0.69252000<br>
&nbsp; &nbsp; &nbsp;-1: X=0.00000000 Y=0.84493000 Z=0.30748000<br>&nbsp; &nbsp; &nbsp;-1: X=0.00000000 Y=0.65507000 Z=0.80748000<br>Calcium &nbsp; &nbsp;NPT= &nbsp;781 &nbsp;R0=.000500000 RMT= &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Z: 20.0<br>LOCAL ROT MATRIX: &nbsp; &nbsp;1.0000000 0.0000000 0.0000000<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 0.0000000 1.0000000 0.0000000<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 0.0000000 0.0000000 1.0000000<br>&nbsp; 0<br>***************************************************************<br>I wanna do the optimization on both cell constants and atomic positions<br>with constant<br>volumes. The question is I found the optimized cell constants don&#39;t<br>
change and the atomic positions changed<br>very small (&lt; 0.001), I think these results doesn&#39;t make sense. Anybody<br>can tell me some ideas about them?<br>What I did is the following commands:<br>-- instgen_lapw<br>
-- init_lapw &nbsp;(I used PBE -6.0 Ry )<br>-- run_lapw (scf)<br>-- x optimize &nbsp;(actually in this step, I chose &quot;[3] &nbsp;VARY C/A RATIO with<br>CONSTANT VOLUME and B/A (orthorh lattice)&quot;<br>&nbsp; at -4,-2,0,2,4, then at each C/A variation I run &quot;x optimize&quot; again<br>
to get 5 B/A variations and renamed them, i.e. chose<br>&nbsp; &quot;[4] &nbsp;VARY B/A RATIO with CONSTANT VOLUME and C/A (orthorh lattice)&quot;,<br>finally I got 25 structure files)<br>-- I changed the script &quot; optimize.job &quot;, see below.<br>
****************************************************<br>#!/bin/csh -f<br>&nbsp;# &nbsp; Modify this script according to your needs:<br>&nbsp;# &nbsp; &nbsp; &nbsp;Uncomment one of the lines ...<br>&nbsp;# &nbsp; &nbsp; &nbsp;Change run_lapw to runsp_lapw or use different convergence criterium<br>
&nbsp;# &nbsp; &nbsp; &nbsp;Change save_lapw -d XXX<br>foreach i ( \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___0.0_0 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___0.0_-4 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-4.0_-4 &nbsp;\<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-2.0_-4 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___2.0_-4 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___4.0_-4 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-4.0_-2 &nbsp;\<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-2.0_-2 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___0.0_-2 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___2.0_-2 &nbsp;\<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___4.0_-2 &nbsp;\<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-4.0_0 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-2.0_0 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___2.0_0 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___4.0_0 &nbsp; \<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-4.0_2 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-2.0_2 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___0.0_2 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___2.0_2 &nbsp; \<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___4.0_2 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-4.0_4 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa__-2.0_4 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___0.0_4 &nbsp; \<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___2.0_4 &nbsp; \<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Ca_cmca_1_5_boa___4.0_4 &nbsp; \<br>&nbsp;)<br>&nbsp; &nbsp; rm Ca_cmca_1_5.struct &nbsp; &nbsp; &nbsp; &nbsp; # NFS-bug<br>&nbsp; &nbsp; cp &nbsp;$i.struct Ca_cmca_1_5.struct<br>
<br>&nbsp;# Please uncomment and adapt any of the lines below according to your needs<br>&nbsp;# &nbsp; &nbsp;cp &nbsp;$i.clmsum Ca_cmca_1_5.clmsum<br>&nbsp; &nbsp; x dstart<br>&nbsp;# &nbsp; &nbsp;x dstart -c<br>&nbsp;# &nbsp; &nbsp;run_lapw -ec 0.0001 -in1new 3 -in1orig -renorm<br>&nbsp;# &nbsp; &nbsp;runsp_lapw -ec 0.0001<br>
&nbsp; &nbsp; min -I -j &quot;run_lapw -p -fc 0.1 -i 200 -NI &quot;<br>&nbsp;# &nbsp; &nbsp;run_lapw -ec 0.0001<br><br>&nbsp; &nbsp; set stat = $status<br>&nbsp; &nbsp; if ($stat) then<br>&nbsp; &nbsp; &nbsp; &nbsp;echo &quot;ERROR status in&quot; $i<br>&nbsp; &nbsp; &nbsp; &nbsp;exit 1<br>&nbsp; &nbsp; endif<br>&nbsp;# &nbsp; &nbsp;save_lapw &nbsp;$i<br>
&nbsp; &nbsp; save_lapw &nbsp;-f -d /home/victor/wien2k_07/Ca_cmca_1_5/final $i<br>&nbsp;end<br>**********************************************************<br><br>-- I also changed *.inM file, changed the force tolerance form 2 to 0.1<br>mRy/a.u., I found the default value is<br>
&nbsp; too high.<br><br>-- finally I run the script &quot;optimize.job&quot;<br><br>After the calculations, I got the results as I stated before.<br><br>&nbsp; &nbsp;Thank you very much!<br><br>Yours<br><br>Victor<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" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>SUSANTA K MOHANTA<br>RESEARCH SCHOLAR<br>DEPT.OF PHYSICS<br>IIT KHARAGPUR<br>*****************************