[Wien] calculation quality

Peter Blaha pblaha at theochem.tuwien.ac.at
Tue Jan 20 08:32:31 CET 2009


No, the atomic number is not really a good parameter.
The needed value of RKmax depends much more on the type of
valence electrons and an "old" rule said:
for s,p wavefunctions you need RKmax 6-7
     d                                7-8
     f                                8-9

So a "heavier" element, which has only s,p electrons (eg. S, Se, Te) needs
a much samller RKmax than a "light" Fe atom.

To some extend we try to "correct" this automatically in a compound and the
setup make a sphere of sp-atoms smaller than the sphere of 3d or even 4f elements.

Of course, within a row in the periodic table, your arguments would be correct,
i.e. a O or F atom needs a larger RKmax than Be or B. Thus for a CO molecule
I'd make the C sphere smaller than the O-sphere (this is NOT done yet by setrmt).

H, which has only s-electrons, is even more special and a much smaller RKmax is enough.

k-points: sure, the size of the unit cell matters and for one and the same material
the k-mesh scales with volume (eg. when creating a supercell).
But other properties may be even more important. For certain magnetic metallic
materials you may need 100-10000 times more k-points than for an insulator,
although the volume might be identical.

There is NO other choice than: start with a realistic (but small) k-mesh and test
the result of the desired properties with respect to k-mesh.

What "realistic, but small" means, needs a bit experience, but there were many
discussions about that in the mailing-list or in WIEN2k-workshops.

Pablo de la Mora schrieb:
> Dear Wien-users,
> There are two parameters that meassure the quality of a calculation; RxK
> and the number of k-points.
> 
> With respect to RxK; for hydrogen R is quite small but it also has one 
> nuclear charge, so a small RxK parameter is acceptable.
> So my question is: would a more complex parameter be more meaningful?
> Such as
>         RxK/n
> where n is the atomic number?
> The idea behind is that the nuclear potential increases with n, so a
> larger R is needed and R/n would be more or less constant.
> 
> With respecto to the number of k-points (Nk) with the size of the cell 
> the number of k-points should be reduced to maintain the same quality, 
> so this parameter could be modified to VxNk.
> 
>        Pablo de la Mora
> 
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien

-- 

                                       P.Blaha
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-15671             FAX: +43-1-58801-15698
Email: blaha at theochem.tuwien.ac.at    WWW: http://info.tuwien.ac.at/theochem/
--------------------------------------------------------------------------


More information about the Wien mailing list