[Wien] Problem in LDA+U calculation

[Robert Laskowski]

Are you sure you used wien2k orb. You can copy orb binary from wien2k to 
wienncm directory.

Robert


On Wednesday 27 June 2012 09:22:32 tran at theochem.tuwien.ac.at wrote:
> Dear Hena,
> 
> Actually, it should work with WIEN2k (I don't know with Wienncm).
> I have tried myself on NiO with
> l=1 and 2 on nickel and l=1 on oxygen and it's running properly.
> The first thing you should do is to try with WIEN2k and
> see if it works or not. If it works with WIEN2k, but not with
> Wienncm, it could be that also a c-shell script of Wienncm has to be
> modified.
> 
> But just in case, you can also try this:
> 
> replace
> 
>   1 2 3 2
>   2 2 3 2
> 
> by
> 
>   1 2 2 3
>   2 2 2 3
> 
> in case.indmc and case.inorb
> 
> F. Tran
> 
> On 25.06.2012 17:45, Hena Das wrote:
> >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

-- 
Dr Robert Laskowski
Vienna University of Technology, Institute of Materials Chemistry, 
Getreidemarkt 9/165-TC, A-1060 Vienna, Austria
tel. +43 1 58801 165303               Fax  +43 1 58801 15698