[Pw_forum] PBE0 convergence
Jie Jiang
jiejiangnc at gmail.com
Wed Jun 1 16:23:23 CEST 2011
Dear Gabriele Sclauzero,
Thanks,
Jie
On Wed, Jun 1, 2011 at 3:05 AM, Gabriele Sclauzero <sclauzer at sissa.it>wrote:
> Dear Jie Jiang,
>
> 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.
>
> Regards,
>
>
> in the behalf of Mr. PWscf,
>
> GS
>
>
> Il giorno 31/mag/2011, alle ore 15.10, Jie Jiang ha scritto:
>
> Dear PWscf,
>
> I get a problem in PBE0 convergence. Below is part of the input file,
>
> &CONTROL
> calculation = "scf",
> /
> &SYSTEM
> ibrav = 12,
> celldm(1) = 15.73981,
> celldm(2) = 0.62519,
> celldm(3) = 4.39673,
> celldm(4) = 0.62519,
> nat = 58,
> ntyp = 5,
> ecutwfc = 80.0,
> input_dft = 'pbe0',
> nqx1 = 2, nqx2 = 1, nqx3 = 1,
> exxdiv_treatment= 'gygi-baldereschi'
> ecutvcut = 0.7
> x_gamma_extrapolation= .TRUE.
> occupations = "smearing",
> smearing = "methfessel-paxton",
> degauss = 0.01,
> /
> &ELECTRONS
> conv_thr = 1.D-6,
> mixing_beta = 0.5D0,
> /
> ATOMIC_SPECIES
> C 12.0107
> S 32.066
> H 1.00794
> O 15.9994
> Zn 65.409
> /
> K_POINTS {automatic}
> 6 3 1 0 0 0
>
> The output 'dexx' is
>
> est. exchange err (dexx) = 0.00059720 Ry
> est. exchange err (dexx) = -0.02431459 Ry
> est. exchange err (dexx) = -0.02596414 Ry
> est. exchange err (dexx) = -0.02720608 Ry
> est. exchange err (dexx) = -0.02821286 Ry
> est. exchange err (dexx) = -0.02946368 Ry
> est. exchange err (dexx) = -0.03148666 Ry
> est. exchange err (dexx) = -0.03481403 Ry
> est. exchange err (dexx) = -0.03964928 Ry
> est. exchange err (dexx) = -0.04531623 Ry
>
> and Fock energy is
>
> -177.9185854736362
> -177.9185854736362 -177.8610717829301 -177.8047524875259
> -177.8047524875259 -177.8356769333793 -177.8179722072562
> -177.8179722072562 -177.8457748551563 -177.8216492155755
> -177.8216492155755 -177.8490273372818 -177.8219932983469
> -177.8219932983469 -177.8500300290520 -177.8216410368603
> -177.8216410368603 -177.8505498629616 -177.8205313239625
> -177.8205313239625 -177.8511288881755 -177.8187531280260
> -177.8187531280260 -177.8519180436694 -177.8154548907789
> -177.8154548907789 -177.8530355735678 -177.8113176980133
> -177.8113176980133 -177.8542927054368 -177.8066352532881
>
> Can you please help me in this issue.
>
> Thanks,
>
> Jie Jiang
>
> Applied Physics
> Yale University
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
>
> § Gabriele Sclauzero, EPFL SB ITP CSEA
> * PH H2 462, Station 3, CH-1015 Lausanne*
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20110601/b48e62d7/attachment.html>
More information about the users
mailing list