[Wien] LDA+U
Tomohiko Tasaka
tasaka at affinity-science.com
Tue Jun 30 14:57:59 CEST 2009
Hello,
Did you execute LDA+U calculation with different
configuration before(e.g. single orbital u)?
Maybe, case.dmat* files are damaged. please try to update them.
For example,
$ runsp -dm -i 1
$ runsp -orb
or
$ rm *.dmat*
$ runsp -orb
With best regards,
Tomo
-----Original Message-----
From: wien-bounces at zeus.theochem.tuwien.ac.at [mailto:wien-bounces at zeus.theochem.tuwien.ac.at] On Behalf Of
Tulika Maitra
Sent: Tuesday, June 30, 2009 8:51 PM
To: A Mailing list for WIEN2k users
Subject: Re: [Wien] LDA+U
Dear Dr. Wu,
No the case.dmatup file is not empty. The density matrices are written in
it. I am not doing a complex calculation.
Thanks for any help.
-best regards
Tulika Maitra
------------------------------------------------------
On Tue, 30 Jun 2009, wu at ph2.uni-koeln.de wrote:
> Is your case.dmatup file void ?
> Did you run a complex calculation which needs
> case.indmc file?
>
> regards -- H. Wu
>
> >
> > Dear Wien-users,
> >
> > I wanted to know if it is possible to apply U in two orbitals of the same
> > atom (for example d and f orbitals of the same atom). From the user guide
> > and mailing list discussions it seems to be possible. But when I tried to
> > do such a calculation I encountered the following problem. I am using
> > Wien2k_08.2 version.
> >
> > The first scf cycle was OK. In the second cycle it stops at 'orb -up'
> > after lapw0 saying 'error in Vorb' and in the dayfile it says 'can't open
> > file 'case.dmatup' whereas in the case.dmatup the density matrices for
> > both d and f orbitals are written. I don't really understand where is the
> > problem.
> >
> > My input files are given below:
> >
> > case.inorb
> >
> > 1 1 0 nmod, natorb, ipr
> > PRATT 1.0 BROYD/PRATT, mixing
> > 1 2 2 3 iatom nlorb, lorb
> > 1 nsic 0..AFM, 1..SIC, 2..HFM
> > 0.37 0.00 U J (Ry) Note: we recommend to use U_eff =
> > 0.52 0.00 U J
> >
> > case.indm
> > -9. Emin cutoff energy
> > 1 number of atoms for which density matrix is
> > calculated
> > 1 2 2 3 index of 1st atom, number of L's, L1
> > 0 0 r-index, (l,s)index
> >
> >
> > Thanks for any help.
> >
> > With regards
> >
> > Tulika Maitra
> >
> > _______________________________________________
> > Wien mailing list
> > Wien at zeus.theochem.tuwien.ac.at
> > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> >
>
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
More information about the Wien
mailing list