I'm sorry to ask this simple question, but it really bothered me for a long time。。。
I want to plot the phonon dispersion curve, but I don’t understand how to set ‘’nph1l‘’ in the tutorial，It has some explanations about nph1l, but I don’t understand how to generate it. When I used nph1l=71 in the tutorial, I got such a dispersion curve. I think there are some problems in it. Who can you help me? I will be very grateful.
About nph1l
About nph1l
 Attachments

 trf2_5.in
 (4.67 KiB) Downloaded 42 times
Re: About nph1l
Dear hmao,
nph1l is the number of qpoints in which the phonon frequencies will be computed and plotted, then you have to give the list of qpoint directions by qph1l with a user defined path in the reciprocal space, as you have in your input file. I would advise to put more qpoints to have smoother curves.
In your case it seems that you have unstable phonon modes (appearing as negative numbers but they are actually imaginary). If you have an insulator then you'll have to correct the LOTO splitting by, either editing the file where the eigenvalues are written for Gamma point and replace by the corresponding frequencies obtained in the second list of wavectors where the nonanalytical correction is applied. A shortcut tip for this LOTO is to put a qpt value very close to Gamma instead of exactly Gamma in the qph1l, for example:
0.0001 0.0001 0.0001 1.0
Best wishes,
Eric
nph1l is the number of qpoints in which the phonon frequencies will be computed and plotted, then you have to give the list of qpoint directions by qph1l with a user defined path in the reciprocal space, as you have in your input file. I would advise to put more qpoints to have smoother curves.
In your case it seems that you have unstable phonon modes (appearing as negative numbers but they are actually imaginary). If you have an insulator then you'll have to correct the LOTO splitting by, either editing the file where the eigenvalues are written for Gamma point and replace by the corresponding frequencies obtained in the second list of wavectors where the nonanalytical correction is applied. A shortcut tip for this LOTO is to put a qpt value very close to Gamma instead of exactly Gamma in the qph1l, for example:
0.0001 0.0001 0.0001 1.0
Best wishes,
Eric
Re: About nph1l
Dear Eric,
Thanks a lot to reply me.
I am a new, and I am trying to add more qpoints, but I want to know if there is a script that generates qpoints with abinit.
Best wishes,
Mao
Thanks a lot to reply me.
I am a new, and I am trying to add more qpoints, but I want to know if there is a script that generates qpoints with abinit.
Best wishes,
Mao
Re: About nph1l
Dear Mao,
You're welcome!
For qpoint grid you can either use nqpath and qpath, but if you want to control exactly the number of qpoint and their direction, you can do it with a small script or with a spread shit, etc.
Abipy might do that for you too, but I'm not expert with it...
Best wishes,
Eric
You're welcome!
For qpoint grid you can either use nqpath and qpath, but if you want to control exactly the number of qpoint and their direction, you can do it with a small script or with a spread shit, etc.
Abipy might do that for you too, but I'm not expert with it...
Best wishes,
Eric