[Wien] How to include d and f-states in calculation ??
tran at theochem.tuwien.ac.at
tran at theochem.tuwien.ac.at
Thu Jul 21 13:39:28 CEST 2011
The best is to give an example:
case.inorb for the case LDA+U is applied to the electrons d and f on
atoms number 2, 4 and 5:
1 3 0 nmod, natorb, ipr
PRATT,1.0 mixmod, amix
2 2 2 3 iatom nlorb, lorb
4 2 2 3 iatom nlorb, lorb
5 2 2 3 iatom nlorb, lorb
0.400 0.95 U J U and J for d electrons on atom 2
0.500 0.95 U J U and J for f electrons on atom 2
0.400 0.95 U J U and J for d electrons on atom 4
0.500 0.95 U J U and J for f electrons on atom 4
0.400 0.95 U J U and J for d electrons on atom 5
0.500 0.95 U J U and J for f electrons on atom 5
I chose the values of U and J randomly. To know which value
The indices iatom are in case.struct.
The two main ways of choosing U are
1) to search for values in the literature
2) to do calculations with several values of U and to see which value
gives the best agreement with experiment.
If you use LDA+U you need to create case.inorb and case.indm.
For EECE/hybrid you need only case.ineece (case.inorb and case.indm will
be created automatically)
On Wed, 20 Jul 2011, yedu kondalu wrote:
> Dear Wien2k experts and users,
>
> I am running wien version 11. Iam new to wien2k, I want to treat both d
> and f electrons of X atom. I have some questions in this regard.
>
> *1*. How can I modify the case.inorb file to include d anf f - states ??
> *2*. If the cell contains more than one formula unit, then how to
> treat(include d and f - states for the same type of atom) again ??
>
> As an example, for two atoms of same type to include d and f-states
>
> *case.inorb*
> 1 2 0 nmod, natorb, ipr
> PRATT,1.0 mixmod, amix
> 1 2 2 3 iatom nlorb, lorb
> 2 2 2 3 iatom nlorb, lorb
> 0.5000
> 1
> 0.500 0.0 U J
> 0.500 0.0 U J does this file works good or ??
>
> *3*. How to fix the U parameters ??
> *4*. Is it necessary to change the files called* case.ineece* and *case.indm
> * ??
> *5*. If needed, What are all the changes should be ??
>
> Can u please some body explain me in detail because I just started doing
> calculations of strongly correlated systems.
> Thank you very much in advance
>
> KONDAL
>
More information about the Wien
mailing list