In this case, Mr. Robert Laskowski should know how old is the orb package in Wienncm. I think that we changed also something in one of the lapw scripts.
On Mon, 25 Jun 2012, Hena Das wrote: > Yes, I am using Wienncm code. > > ________________________________________ > From: wien-bounces at zeus.theochem.tuwien.ac.at [wien-bounces at > zeus.theochem.tuwien.ac.at] on behalf of tran at theochem.tuwien.ac.at [tran > at theochem.tuwien.ac.at] > Sent: Monday, June 25, 2012 4:00 PM > To: A Mailing list for WIEN2k users > Subject: Re: [Wien] Problem in LDA+U calculation > > I vaguely remember that around 2006, we changed something such that it > should work, but I am not sure. If this is true, it could be that these > changes were not included in the orb package of the NCM code (you are > mentioning non-collinear spin configuration, right?). > > F. Tran > > On Mon, 25 Jun 2012, Peter Blaha wrote: > > > I know that Pavel Novak programmed orb so that it can read more l-values per > > atom. However, I've never tested this myself because I believe there is no > > physics in it. In other words: I do not know an atom where it makes sense to > > put a U for d and f states. For lanthanides the 5d density will be quite > > delocalized and anyway only some fraction is inside the atomic sphere. > > > > So unless you can convince me that you have a physically sensible system, I > > don't know if I should debug this ... > > > > Am 25.06.2012 17:45, schrieb Hena Das: > > > Dear All, > > > > > > I am doing a scf calculation for a particular non-collinear spin > > > configuration. A am using LDA+U method. When I use the following > > > case.inorb and case.indmc files: > > > case.inorb > > > 1 3 0 nmod, natorb, ipr > > > PRATT 1.0 BROYD/PRATT, mixing > > > 1 1 3 > > > 2 1 3 > > > 3 1 2 > > > 1 nsic 0..AMF, 1..SIC, 2..HFM > > > 0.65 0.07 U J (Ry) > > > 0.65 0.07 > > > 0.45 0.07 > > > case.indmc > > > -9. Emin cutoff energy > > > 3 number of atoms for which density matrix is > > > calculated > > > 1 1 3 index of 1st atom, number of L's, L1 > > > 2 1 3 dtto for 2nd atom, repeat NATOM times > > > 3 1 2 index of 1st atom, number of L's, L1 > > > 0 0 r-index, (l,s)index > > > > > > the program executes properly without any error. However when I use the > > > other set of case.inorb and case.indmc files: > > > case.inorb > > > 1 3 0 nmod, natorb, ipr > > > PRATT 1.0 BROYD/PRATT, mixing > > > 1 2 3 2 > > > 2 2 3 2 > > > 3 1 2 > > > 1 nsic 0..AMF, 1..SIC, 2..HFM > > > 0.65 0.07 > > > 0.45 0.07 U J (Ry) > > > 0.65 0.07 > > > 0.45 0.07 > > > 0.45 0.07 > > > case.indmc > > > -9. Emin cutoff energy > > > 3 number of atoms for which density matrix is > > > calculated > > > 1 2 3 2 index of 1st atom, number of L's, L1 > > > 2 2 3 2 dtto for 2nd atom, repeat NATOM times > > > 3 1 2 index of 1st atom, number of L's, L1 > > > 0 0 r-index, (l,s)index > > > the program stops by giving the error: error in Vorb. It is not writing > > > case.vorbup/dn files. In the next set I just use U at the d states as > > > well as at the f states for atom 1 and 2. Am I giving any wrong input? > > > Help me to solve this problem. > > > > > > Waiting for suggestions. > > > > > > Best, > > > Hena > > > > > > > > > > > > > > > _______________________________________________ > > > Wien mailing list > > > Wien at zeus.theochem.tuwien.ac.at > > > http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien > > > > > > > -- > > Peter Blaha > > Inst.Materials Chemistry > > TU Vienna > > Getreidemarkt 9 > > A-1060 Vienna > > Austria > > +43-1-5880115671 > > > > > > _______________________________________________ > > 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 >