<DIV>&nbsp;</DIV>
<DIV>&gt;Message:&nbsp;2<BR>&gt;Date:&nbsp;Thu,&nbsp;15&nbsp;Oct&nbsp;2009&nbsp;15:42:34&nbsp;+0200<BR>&gt;From:&nbsp;"Lorenzo&nbsp;Paulatto"&nbsp;&lt;paulatto@sissa.it&gt;<BR>&gt;Subject:&nbsp;Re:&nbsp;[Pw_forum]&nbsp;total&nbsp;force&nbsp;ascend&nbsp;in&nbsp;the&nbsp;course&nbsp;of&nbsp;VC-relax<BR>&gt;To:&nbsp;"PWSCF&nbsp;Forum"&nbsp;&lt;pw_forum@pwscf.org&gt;<BR>&gt;Message-ID:&nbsp;&lt;op.u1ugc8i1a8x26q@paulax&gt;<BR>&gt;Content-Type:&nbsp;text/plain;&nbsp;charset=utf-8;&nbsp;format=flowed;&nbsp;delsp=yes<BR>&gt;<BR>&gt;In&nbsp;data&nbsp;15&nbsp;ottobre&nbsp;2009&nbsp;alle&nbsp;ore&nbsp;15:38:39,&nbsp;Q.J.Wang&nbsp;&lt;wangqj1@126.com&gt;&nbsp;ha&nbsp;&nbsp;<BR>&gt;scritto:<BR>&gt;<BR>&gt;&gt;&nbsp;Dear&nbsp;all<BR>&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;I&nbsp;did&nbsp;geometry&nbsp;optimization&nbsp;with&nbsp;vc-relax&nbsp;in&nbsp;a&nbsp;supercell.&nbsp;When&nbsp;I&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;increase&nbsp;the&nbsp;concentration&nbsp;of&nbsp;doped&nbsp;atoms&nbsp;,the&nbsp;value&nbsp;of&nbsp;total&nbsp;force&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;ascend&nbsp;instead&nbsp;of&nbsp;declining.Before&nbsp;I&nbsp;did&nbsp;not&nbsp;increase&nbsp;the&nbsp;concentration&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;of&nbsp;doped&nbsp;atoms&nbsp;,It&nbsp;converges&nbsp;well&nbsp;.Anyone&nbsp;who&nbsp;encounter&nbsp;the&nbsp;same&nbsp;problem&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;as&nbsp;me&nbsp;?<BR>&gt;<BR>&gt;I&nbsp;sounds&nbsp;normal&nbsp;to&nbsp;me,&nbsp;when&nbsp;you&nbsp;perturb&nbsp;the&nbsp;system&nbsp;you&nbsp;will&nbsp;have&nbsp;to&nbsp;redo&nbsp;&nbsp;<BR>&gt;the&nbsp;relaxation...&nbsp;why&nbsp;do&nbsp;you&nbsp;find&nbsp;it&nbsp;puzzling?</DIV>
<DIV>The problem is when I redo the relaxation ,the&nbsp;value&nbsp;of&nbsp;'Total force'and the 'P' become large in every bfgs cycle .And at last shows error :</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<EM>&nbsp;&nbsp; from bfgs : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<BR></EM><I>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; bfgs history already reset at previous step</I></DIV>
<DIV><I>The input file as following :</I></DIV>
<DIV><I>&nbsp;&nbsp;  &amp;CONTROL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'vc-relax' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/root/work3/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/root/pseudo/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'GaAs' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nstep = 250 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .true. ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true. ,<BR>&nbsp;/<BR>&nbsp;&amp;SYSTEM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;nat = 48,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 4,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 30 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 300 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nosym = .true. ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = 'smearing' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.005 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = 'gaussian' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin = 2 ,<BR>&nbsp;&nbsp; starting_magnetization(1) = 0.5,<BR>&nbsp;&nbsp; starting_magnetization(2) = 0.5,<BR>&nbsp;&nbsp; starting_magnetization(3) = 0.5,<BR>&nbsp;&nbsp; starting_magnetization(4) = 0.5,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lda_plus_u = .false. ,<BR>&nbsp;/<BR>&nbsp;&amp;ELECTRONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.D-7 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'plain' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.3 ,<BR>&nbsp;/<BR>&nbsp;&amp;IONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ion_dynamics = 'bfgs' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; phase_space = 'full' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pot_extrapolation = 'second_order' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; wfc_extrapolation = 'second_order' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; upscale = 100 ,<BR>&nbsp;/<BR>&nbsp;&amp;CELL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cell_dynamics = 'bfgs' ,<BR>&nbsp;/<BR>ATOMIC_SPECIES</I></DIV>
<DIV><EM>Ga&nbsp;&nbsp; ...........</EM></DIV>
<DIV><I>ATOMIC_POSITIONS angstrom&nbsp;</I></DIV>
<DIV><EM>...............</EM></DIV>
<DIV><I>&nbsp; K_POINTS automatic <BR>&nbsp; 3 3 3&nbsp;&nbsp; 1 1 1 <BR></I><I><BR></I></DIV>
<DIV>&gt;<BR>&gt;<BR>&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The&nbsp;second&nbsp;problem&nbsp;is&nbsp;when&nbsp;I&nbsp;increase&nbsp;the&nbsp;relative&nbsp;distance&nbsp;of&nbsp;the&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;two&nbsp;doped&nbsp;atoms&nbsp;in&nbsp;the&nbsp;supercell&nbsp;and&nbsp;did&nbsp;geometry&nbsp;optimization&nbsp;with&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;vc-relax,and&nbsp;the&nbsp;error&nbsp;turn&nbsp;up&nbsp;as&nbsp;like&nbsp;:&nbsp;from&nbsp;c_bands&nbsp;:&nbsp;error&nbsp;#&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;1<BR>&gt;&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;too&nbsp;many&nbsp;bands&nbsp;are&nbsp;not&nbsp;converged<BR>&gt;&gt;&nbsp;So&nbsp;I&nbsp;am&nbsp;puzzled&nbsp;why&nbsp;before&nbsp;I&nbsp;increase&nbsp;relative&nbsp;position&nbsp;of&nbsp;two&nbsp;doped&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;atoms&nbsp;the&nbsp;vc-relax&nbsp;converged&nbsp;and&nbsp;the&nbsp;latter&nbsp;turns&nbsp;up&nbsp;error&nbsp;in&nbsp;the&nbsp;same&nbsp;&nbsp;<BR>&gt;&gt;&nbsp;parameter&nbsp;.<BR>&gt;<BR>&gt;It&nbsp;happens&nbsp;with&nbsp;peculiar&nbsp;configurations,&nbsp;the&nbsp;easiest&nbsp;way&nbsp;to&nbsp;get&nbsp;rid&nbsp;of&nbsp;it&nbsp;&nbsp;<BR>&gt;is&nbsp;normally&nbsp;to&nbsp;increase&nbsp;the&nbsp;number&nbsp;of&nbsp;bands.<BR>&gt;<BR>&gt;cheers<BR><BR>&gt;--&nbsp;</DIV>
<DIV>Best regards</DIV>
<DIV>&nbsp;</DIV>
<DIV>Q.J.Wang</DIV>
<DIV>&nbsp;</DIV>
<DIV>XiangTan University </DIV>
<DIV><BR><BR></DIV><br><br><span title="neteasefooter"/><hr/>
<a href="http://allyes.nie.163.com/main/adfclick?db=afanie&bid=1260,614,23&cid=148,4,1&sid=1357&show=ignore&url=http://tx2.163.com/fab.html">看陆川杨幂新片《琴棋书画》，品网易3D国韵网游《天下贰》</a>
</span>