[Pw_forum] specify correct q-point grid!
Eyvaz Isaev
eyvaz_isaev at yahoo.com
Mon Dec 8 22:16:37 CET 2008
Dear Max,
> IF (dos) THEN
> IF (nk1 < 1 .OR. nk2 < 1 .OR. nk3 < 1)
> & CALL errore ('matdyn','specify correct q-point grid!',1)
You correctly noticed that the error is due these lines.
Because of that integration of DOS is done using the tetrahedra method, so, there will be at least 1 tetrahedron with 4 points (vertexes).
> ntetra = 6 * nk1 * nk2 * nk3
This line (total number of tetrahedra) is due to the improved tetrahedra method by Bl\"ochl. You can find more in the original paper.
> nq1=3, nq2=3, nq3=2 ( in elph.in )
This q-points are not connected to DOS calculations in matdyn.x, they are for dynamical matrices D_{a,b}^{i,j}(q) calculations which then used to perform FFT to get real space F_{a,b}^{i,j}(R) force constants matrix.
You should use nk1,nk2,nk3 keywords in matdyn.in file for DOS calculations. By the way, most likely you did not specify these keywords which are zeros by default. Namely this one is stated in the abovementioned error.
Hope this helps,
Bests,
Eyvaz.
-------------------------------------------------------------------
Prof. Eyvaz Isaev,
Theoretical Physics Department, Moscow State Institute of Steel & Alloys, Russia,
Department of Physics, Chemistry, and Biology (IFM), Linkoping University, Sweden
Condensed Matter Theory Group, Uppsala University, Sweden
Eyvaz.Isaev at fysik.uu.se, isaev at ifm.liu.se, eyvaz_isaev at yahoo.com
--- On Mon, 12/8/08, 陳 少華 <chen_shao_hua197 at yahoo.com.tw> wrote:
> From: 陳 少華 <chen_shao_hua197 at yahoo.com.tw>
> Subject: [Pw_forum] specify correct q-point grid!
> To: "Pwscf" <pw_forum at pwscf.org>
> Date: Monday, December 8, 2008, 10:04 PM
> 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
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
More information about the users
mailing list