[Wien] Error in spaghetti (irrep?)

Peter Blaha pblaha at theochem.tuwien.ac.at
Fri Jun 10 14:51:44 CEST 2005


TOLDG of 0.5E-3 is meaningless.

Degenerate eigenvalues should be degenerate to about 2.E-8, thus maybe my
recommendation to increase TOLDG slightly.

A few more points: ONLY at Gamma and A (maybe I missed some others) there
are 2-dimensional representations (degenerate eigenvalues). At your knum=7
all eigenvalues are non-degenerate and do not need to be identical.
At Gamma you have eg:
                   
     -0.0384385   -0.0275091   -0.0275090   -0.0252905   -0.0242101
     -0.0242100   -0.0212862   -0.0212862   -0.0190025   -0.0186596

2+3, 5+6 and 7+8 should most likely be degenerate, but only 7+8 are.
TOLDG should be set such that also 2+3 are "degenerate".  

In order to test spaghetti, I need case.struct, case.output1up, 
case.outputirrup,case.irrepup and case.insp

Please send some "non-parallel" test.

> Increasing TOLDG up to 0.5E-3, the number of errors in spaghetti are reduced
> systematically. However, even for TOLDG =0.5E-3, all errors do not disapperar.
> Is this a reasonable value for TOLDG?
> 
> For each k-point at which an error occurs in the new plot, there are two
> eigenvalues close together which 'float' alone and are not 'joined' to other
> eigenvalues in the spaghetti plot.
> In outputir, these points are skipped (they correspond to the incorrect points
> in spaghetti):
> Spin up:   knum=7, bnd 17,   knum=14, bnd 17,   in outputirup_1
> Spin dn:   knum=6, bnd 17,   knum=9, bnd 17,    in  outputirdn_1  and
> knum=10, bnd17   in outputirdn_4
> Other bands have =?? (eg SpinUP band 44 of knum 7 and band 48 of knum 14) most
> likely due to the earlier point's problem.
> 
> For knum=7 (eigenvalues below) band 17 (-0.9725160) differs by 0.0004 from its
> previous eigenvalue, while band 44 (-0.0969005) differs by -0.0000012 which
> seems strange with TOLDG=0.5E-3?
> All the other differences are of similar size and smaller than 0.5E-3. For the
> spin down point knum=10, bnd 17 in outputirdn_4, the difference is 0.4E-6.
> 
>     K=   0.15000   0.00000   0.00000             (knum=7)
>     -3.8395122   -2.1050534   -2.1050530   -2.1023162   -1.0487342
>     -1.0139718   -1.0130222   -1.0122362   -1.0075969   -0.9951106
>     -0.9878561   -0.9783544   -0.9770046   -0.9769033   -0.9748384
>     -0.9729166  *-0.9725160*  -0.9616254   -0.9609163   -0.9583208
>     -0.1549498   -0.1527483   -0.1464105   -0.1450850   -0.1439287
>     -0.1379788   -0.1367145   -0.1312673   -0.1306861   -0.1283666
>     -0.1260235   -0.1256493   -0.1245272   -0.1151976   -0.1142609
>     -0.1137208   -0.1116372   -0.1089792   -0.1078344   -0.1063183
> 
>     -0.1057823   -0.1011095   -0.1005828  *-0.0969005*  -0.0968993
>     -0.0921372   -0.0910378   -0.0881868   -0.0881706   -0.0838515
> 
> I am sending the new spaghetti files (for TOLDG=0.5E-3) and output files
> directly to your e-mail.


                                      P.Blaha
--------------------------------------------------------------------------
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-15671             FAX: +43-1-58801-15698
Email: blaha at theochem.tuwien.ac.at    WWW: http://info.tuwien.ac.at/theochem/
--------------------------------------------------------------------------



More information about the Wien mailing list