<br>Dear Gabriele Sclauzero,<br><br>Thanks,<br><br>Jie<br><br><div class="gmail_quote">On Wed, Jun 1, 2011 at 3:05 AM, Gabriele Sclauzero <span dir="ltr"><<a href="mailto:sclauzer@sissa.it">sclauzer@sissa.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap: break-word;">Dear Jie Jiang,<div><br></div><div> Are you sure that your atomic structure is correct? We cannot judge since you did not provide the atomic positions. You can check by yourself using xcrysden, for instance. Do the standard PBE calculations work well and give meaningful results? For PBE0, you can maybe try to reduce the mixing_beta, or use another mixing scheme. Finally, are you sure about your k-point sampling? Your cell has celldm(2)=b/a < 1, then b<a, so why do you use more k-points along the first direction? I would rather use 3 6 1, instead of 6 3 1.</div>
<div><br></div><div>Regards,</div><div><br></div><div><br></div><div>in the behalf of Mr. PWscf,</div><div><br></div><div> GS</div><div><br></div><div><br></div><div><div><div><div>Il giorno 31/mag/2011, alle ore 15.10, Jie Jiang ha scritto:</div>
<br><blockquote type="cite"><div><div></div><div class="h5">Dear PWscf,<br>
<br>
I get a problem in PBE0 convergence. Below is part of the input file,<br>
<br>
&CONTROL<br>
calculation = "scf",<br>
/<br>
&SYSTEM<br>
ibrav = 12,<br>
celldm(1) = 15.73981,<br>
celldm(2) = 0.62519,<br>
celldm(3) = 4.39673,<br>
celldm(4) = 0.62519,<br>
nat = 58,<br>
ntyp = 5,<br>
ecutwfc = 80.0,<br>
input_dft = 'pbe0',<br>
nqx1 = 2, nqx2 = 1, nqx3 = 1,<br>
exxdiv_treatment= 'gygi-baldereschi'<br>
ecutvcut = 0.7<br>
x_gamma_extrapolation= .TRUE.<br>
occupations = "smearing",<br>
smearing = "methfessel-paxton",<br>
degauss = 0.01,<br>
/<br>
&ELECTRONS<br>
conv_thr = 1.D-6,<br>
mixing_beta = 0.5D0,<br>
/<br>
ATOMIC_SPECIES<br>
C 12.0107 <br>
S 32.066 <br>
H 1.00794 <br>
O 15.9994 <br>
Zn 65.409 <br>
/<br>
K_POINTS {automatic}<br>
6 3 1 0 0 0<br>
<br>
The output 'dexx' is<br>
<br>
est. exchange err (dexx) = 0.00059720 Ry<br>
est. exchange err (dexx) = -0.02431459 Ry<br>
est. exchange err (dexx) = -0.02596414 Ry<br>
est. exchange err (dexx) = -0.02720608 Ry<br>
est. exchange err (dexx) = -0.02821286 Ry<br>
est. exchange err (dexx) = -0.02946368 Ry<br>
est. exchange err (dexx) = -0.03148666 Ry<br>
est. exchange err (dexx) = -0.03481403 Ry<br>
est. exchange err (dexx) = -0.03964928 Ry<br>
est. exchange err (dexx) = -0.04531623 Ry<br>
<br>
and Fock energy is<br>
<br>
-177.9185854736362<br>
-177.9185854736362 -177.8610717829301 -177.8047524875259<br>
-177.8047524875259 -177.8356769333793 -177.8179722072562<br>
-177.8179722072562 -177.8457748551563 -177.8216492155755<br>
-177.8216492155755 -177.8490273372818 -177.8219932983469<br>
-177.8219932983469 -177.8500300290520 -177.8216410368603<br>
-177.8216410368603 -177.8505498629616 -177.8205313239625<br>
-177.8205313239625 -177.8511288881755 -177.8187531280260<br>
-177.8187531280260 -177.8519180436694 -177.8154548907789<br>
-177.8154548907789 -177.8530355735678 -177.8113176980133<br>
-177.8113176980133 -177.8542927054368 -177.8066352532881<br>
<br>
Can you please help me in this issue.<br>
<br>
Thanks,<br>
<br>
Jie Jiang<br>
<br>
Applied Physics<br>
Yale University <br></div></div>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org" target="_blank">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><font color="#888888"><br><div>
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
<span style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic;"><br>§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font color="#7e7e7e"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div>
</span>
</div>
<br></font></div></div></div><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>
<br></blockquote></div><br>