Hi Eyvaz.,<br><br>I've tried this method (doing relaxation while holding some atoms fixed) long ago when I started learning using PWSCF. But the problem was that the number of iterations needed in this case is much much more than if I set all atoms free!! I have this problem in my mind since then. Do any one know why?<br><br>Thank you!<br><br>Lawrence<br><br><b><i>Eyvaz Isaev <eyvaz_isaev@yahoo.com></i></b> 說：<blockquote class="replbq" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> <style type="text/css"><!-- DIV {margin:0px;} --></style><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">Hi Prithwish,<br><br>You can do it easily. From INPUT_PW:<br>    X  x y z   {if_pos(1) if_pos(2) if_pos(3)}<br>  where :<br>    X        
 Character: label of the atom as specified in ATOMIC_SPECIES<br>    x, y, z   Real: atomic positions<br>    if_pos:   Integer, optional ( default = 1 ): component i of the force for<br>              this atom is multiplied by if_pos(i), which must be either 0 or 1.<br>              Used to keep selected atoms and/or selected components  fixed in<br>              meta-dynamics, neb, smd, MD dynamics or structural optimization<br>              run.<br>I.e. like this where Al atom is supposed to be relaxed, but As atom not.<br><br>Al 0.5 0.5 0.5 1 1 1<br>...........................<br>As  0 0 .5   0 0 0<br><br>Bests,<br>Eyvaz.<br>P.S. By the way there is no
 need to duplicate your message.<br><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">----- Original Message ----<br>From: Prithwish Nandi <pnspeaks@sify.com><br>To: pw_forum@pwscf.org<br>Sent: Thursday, April 5, 2007 1:38:55 PM<br>Subject: [Pw_forum] How to calculate the relaxation calculation by setting one of the atomic coordinate fixed???<br><br><pre>  Hi, <br>     I want to do a relaxation calculation of a 32<br>  atom<br>cell by keeping one atom at a fixed position and other 31 atoms to be<br>relaxed.Is it possible to do the problem by setting X-iforce,Y-iforce<br>and Z-iforce equals to 1 for 31 atoms and by not setting for the atom I<br>want to be fixed ? Or, there is some other procedures for doing the job?<br>Please help.<br><br><br>  Prithwish Nandi<br>  Research Scholar<br>  IGCAR,<br>Kalpakkam, India<br>   </pre>  <br><br><a rel="nofollow" target="_blank"
 href="http://ads.sify.com/RealMedia/ads/click_nx.ads/mail.sify.com/sentmail@Bottom"> <img border="0" height="40" width="300"></a> <br> <div><font color="#ff3300" face="arial, helvetica, sans-serif"><strong>Now book your  Railway Tickets by cash at Sify Iway</strong></font>. For more details contact  our <a rel="nofollow" target="_blank" href="http://customercare.sify.com/">Customer Care</a> <br><br> Complete Coverage of the ICC World Cup '07! Log on to <a rel="nofollow" target="_blank" href="http://www.sify.com/khel">www.sify.com/khel</a> for latest updates, expert columns, schedule, desktop scorecard, photo galleries and more! <br><br> For the Expert view of the ICC World Cup log on to <a rel="nofollow" target="_blank" href="http://www.sify.com/khel">www.sify.com/khel</a>. Read exclusive interviews with Sachin, Ganguly, Yuvraj, Sreesanth, Expert Columns by Gavaskar, Web chat with Dhoni and more! </div> </div><br></div></div><br>  <hr size="1"><a
 href="http://us.rd.yahoo.com/evt=49981/*http://advision.webevents.yahoo.com/mailbeta/features_spam.html">Sucker-punch spam</a> with award-winning protection.<br> Try the <a href="http://us.rd.yahoo.com/evt=49981/*http://advision.webevents.yahoo.com/mailbeta/features_spam.html">free Yahoo! Mail Beta.</a></blockquote><br><p> _______________________________________<br> YM - 離線訊息<br> 就算你沒有上網，你的朋友仍可以留下訊息給你，當你上網時就能立即看到，任何說話都冇走失。<br> http://messenger.yahoo.com.hk