[Wien] orbital pot (DFT+U)
tran at theochem.tuwien.ac.at
tran at theochem.tuwien.ac.at
Thu Nov 26 17:53:41 CET 2015
Did you adapt natorb, iatom, etc. to your particular case?
Also, 0.07 Ry is a very small U which is one order of
magnitude smaller than usual values.
F. Tran
On Thu, 26 Nov 2015, Jihoon Park wrote:
> Dear users,
>
> This is my case.inorb file.
>
> 1 2 0 nmod, natorb, ipr
> PRATT 1.0 BROYD/PRATT, mixing
> 1 1 2 iatom nlorb, lorb
> 2 1 2 iatom nlorb, lorb
> 1 nsic 0..AMF, 1..SIC, 2..HFM
> 0.07 0.00 U J (Ry) Note: we recommend to use U_eff = U-J and J=0
> 0.07 0.00 U J
>
> If there is any problem, please give me some advice.
> And if my case.inorb is correct, I think the step-by-step procedure for this
> calculation is wrong.
> Please give me any possible reasons for the same calculated total energies
> for different U values.
>
>
> All my best,
> Jihoon Park
>
> On Thu, Nov 26, 2015 at 9:37 AM, Jihoon Park
> <maximumenergyproduct at gmail.com> wrote:
> Yes, I used it.
>
> On Wed, Nov 25, 2015 at 1:15 PM, <tran at theochem.tuwien.ac.at>
> wrote:
> Did you use run calculation with the -orb flag:
> runsp_lapw -orb
>
> On Wed, 25 Nov 2015, Jihoon Park wrote:
>
> Dear users,
>
> I have tried to add GGA+U with different
> U values, but always get the same total
> energies.
> The steps that I did are as follows:
> 1. Do spin-polarized calculation.
> 2. setup spin-orbit coupling and orbital
> pot (DFT+U) together.
> 3. Run
>
> This procedure is very simple and
> straightforward, but I don't know why
> there is no difference for different U
> values.
> Would you please let me know what
> mistakes I have been making?
>
>
> All my best,
> Jihoon Park
>
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> SEARCH the MAILING-LIST at:
> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
>
>
>
>
>
More information about the Wien
mailing list