[Wien] case.indm or case.indmc

Ahmad Ranjbar ranjbar at imr.edu
Tue Jul 27 04:22:52 CEST 2010


Dear WIEN2K users,

I have some small questions.I am wondering if anyone can reply me.


1- My system is graphene sheet. The following is the case.in1 file.


----------------- top of file: case.in1------------------------------------------------------------------------------
WFFIL        (WFPRI, SUPWF)
  7.00       10    4 (R-MT*K-MAX; MAX L IN WF, V-NMT
  0.30    3  0      (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)
 0   -0.78      0.010 CONT 1
 0    0.30      0.000 CONT 1
 1    0.30      0.000 CONT 1
K-VECTORS FROM UNIT:4   -9.0       2.0    7   red   emin/emax/nband
------------------bottom of file--------------------------------------------------------------------------------------

I put some part of the last iteration of case.scf file here:

------------------------------------------------------------------------------------------------------------
          ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM  C 1
:e__0001: OVERALL ENERGY PARAMETER IS    0.3000
          OVERALL BASIS SET ON ATOM IS LAPW
:E0_0001: E( 0)=   -0.7800   E(BOTTOM)=   -1.760   E(TOP)= -200.000
             APW+lo
:E0_0001: E( 0)=    0.3000
             LOCAL ORBITAL
:E1_0001: E( 1)=    0.3000
             APW+lo

       K=   0.00000   0.00000   0.00000            1
:RKM  : MATRIX SIZE  221LOs:  10  RKM= 6.81  WEIGHT= 1.00  PGR:
       EIGENVALUES ARE:
:EIG00001:      -0.7960186    0.4854188    0.4854189    0.4943940    1.3050301
:EIG00006:       1.3050303    1.4322420    1.5433034    1.7005535
------------------------------------------------------------------------------------------------------------

As you see, E(TOP) is smaller than E(BOTTOM)!!!!!! and it remain at the initial value set to -200.
I think that during scf cycle and searching up and down in energy starting with El= -0.78,  both Etop and Ebottom could not to be found.

I am just wondering is it right or not? How can i fix this problem. Is WIEN2K results sensitive  to parameters located in case.in1?




2) When i use w2web to initialize a LAPWDM SCF calculation, i need to edit an input file for LAPWDM as follow:

------------------ top of file: case.indm --------------------
-9. Emin cutoff energy
1 number of atoms for which density matrix is calculated
1 1 2 index of 1st atom, number of L’s, L1
0 0 r-index, (l,s)-index
------------------- bottom of file ------------------------


But when i use terminal to submit my job, there is an error which says there is no case.indmc file.

What is the problem? Is there any difference between case.indm  and  case.indmc?



Thank you very much for your cooperation?



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20100727/6b8d663a/attachment.htm>


More information about the Wien mailing list