<div dir="ltr"><div>Dear Ali</div>
<div>First check your input file with xcrysden&nbsp;. then check ecut and ecutrho(8-10 times of ecut is good).then use</div>
<div>&amp;ions</div>
<div>pot_extrapolation=&#39;second_order&#39;,</div>
<div>wfc_extrapolation=&#39;second_order&#39;,</div>
<div>good luck</div>
<div>Mansoureh <br><br></div>
<div class="gmail_quote">On Mon, Sep 8, 2008 at 11:41 AM, ali kazempoor <span dir="ltr">&lt;<a href="mailto:kazempoor2000@yahoo.com">kazempoor2000@yahoo.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Dear all<br>&nbsp;For test ,I constructed a simple Ni nanowire along [001] and I want to calculate the magnetic anisotropy energy. I get the full relativistic pseudopotential and lattice constant from example 26 of espresso.But for relaxation with nanocoliin and lspinorbit switch it does&#39;nt converge . where is the error or which paprameter must be modified?<br>
thanks<br>Ali Kazempour,Isfahan univ .Of .Tech,Iran<br>&nbsp; &amp;control<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;title = ni ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; calculation = &#39;relax&#39; ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;restart_mode = &#39;from_scratch&#39; ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; tstress = .true.<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;tprnfor = .true.<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;outdir = &#39;/scratch19/user14/&#39; ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;pseudo_dir = &#39;./&#39; ,<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;prefix = ni111 ,<br>&nbsp;/<br>&nbsp;&amp;system<br>&nbsp; &nbsp;ibrav= &nbsp;6, celldm(1)=19.95, celldm(3)=0.3333333,nat=5, ntyp=1,<br>
&nbsp; &nbsp;ecutwfc = 40 , &nbsp;ecutrho =400,<br>&nbsp; &nbsp;lspinorb=.true.,<br>&nbsp; &nbsp;noncolin=.true.,<br>&nbsp; &nbsp;starting_magnetization(1)=0.8,<br>&nbsp; &nbsp;degauss=0.001,<br>&nbsp; &nbsp;smearing=&#39;mp&#39;,<br>&nbsp; &nbsp;angle1(1)=60<br>&nbsp; &nbsp;angle2(1)=90<br>&nbsp; &nbsp;occupations=&#39;smearing&#39;,<br>
&nbsp; &nbsp;ecutwfc =40<br>&nbsp; &nbsp;ecutrho =500.0<br><br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp; &nbsp; &nbsp; &nbsp; conv_thr = &nbsp;1.0d-4<br>&nbsp; &nbsp; &nbsp; &nbsp; mixing_beta = 0.8<br>/<br>&amp;IONS<br>ion_dynamics=&#39;bfgs&#39;<br>/<br>ATOMIC_SPECIES<br>&nbsp;Ni &nbsp;28.0 &nbsp;NirelPBE.RRKJ3.UPF<br>
ATOMIC_POSITIONS crystal<br>&nbsp; Ni &nbsp; &nbsp; &nbsp;0.5 &nbsp; &nbsp; &nbsp; &nbsp; 0.5 &nbsp; &nbsp; &nbsp; &nbsp; 0.0<br>&nbsp; Ni &nbsp; &nbsp; &nbsp;0.3333333 &nbsp; 0.5 &nbsp; &nbsp; &nbsp; &nbsp; 0.5<br>&nbsp; Ni &nbsp; &nbsp; &nbsp;0.5 &nbsp; &nbsp; &nbsp; &nbsp; 0.66666666 &nbsp;0.5<br>&nbsp; Ni &nbsp; &nbsp; &nbsp;0.6666666 &nbsp; 0.5 &nbsp; &nbsp; &nbsp; &nbsp; 0.5<br>&nbsp; Ni &nbsp; &nbsp; &nbsp;0.5 &nbsp; &nbsp; &nbsp; &nbsp; 0.33333333 &nbsp;0.5<br>
K_POINTS automatic<br>&nbsp;1 1 50 &nbsp; 0 0 1<br><br><br><br><br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br></div>