[Pw_forum] specify correct q-point grid!

陳 少華 chen_shao_hua197 at yahoo.com.tw
Mon Dec 8 20:04:44 CET 2008


Dear forum's members,

I do an el-ph calculation.
The error message is
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         0
     from matdyn : error #         1
     specify correct q-point grid!
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

So I examine the  espresso-4.0.3/PH/matdyn.f90
...
IF (dos) THEN
        IF (nk1 < 1 .OR. nk2 < 1 .OR. nk3 < 1) &
             CALL errore  ('matdyn','specify correct q-point grid!',1)
        ntetra = 6 * nk1 * nk2 * nk3
...

And I set

K_POINTS {automatic}
12 12 8 0 0 0           ( in scf.in )

  nq1=3, nq2=3, nq3=2    ( in elph.in )

I don't know why that error message show up and how to solve it ?

Thanks in advance
Max
Candidate of PhD
Physics department, National Taiwan University, Taipei, Taiwan


P.S. cat *.dyn0

   3   3   2
   6
   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00
   0.000000000000000E+00   0.000000000000000E+00  -0.337153960894342E+00
   0.000000000000000E+00   0.338170830884862E+00   0.000000000000000E+00
   0.000000000000000E+00   0.338170830884862E+00  -0.337153960894342E+00
   0.338170830884862E+00   0.338170830884862E+00   0.000000000000000E+00
   0.338170830884862E+00   0.338170830884862E+00  -0.337153960894342E+00





      ______________________________________________________________________________________________________
付費才容量無上限?Yahoo!奇摩電子信箱2.0免費給你,信件永遠不必刪! http://tw.mg0.mail.yahoo.com/dc/landing



More information about the users mailing list