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
Moderator: bguster
About nph1l
- Attachments
-
- phonon dispersion curve
-
trf2_5.in- (4.67 KiB) Downloaded 209 times
Re: About nph1l
Dear hmao,
nph1l is the number of q-points in which the phonon frequencies will be computed and plotted, then you have to give the list of q-point 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 q-points 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 LO-TO 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 non-analytical correction is applied. A shortcut tip for this LO-TO 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 q-points in which the phonon frequencies will be computed and plotted, then you have to give the list of q-point 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 q-points 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 LO-TO 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 non-analytical correction is applied. A shortcut tip for this LO-TO 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 q-points, but I want to know if there is a script that generates q-points with abinit.
Best wishes,
Mao
Thanks a lot to reply me.
I am a new, and I am trying to add more q-points, but I want to know if there is a script that generates q-points with abinit.
Best wishes,
Mao
Re: About nph1l
Dear Mao,
You're welcome!
For q-point grid you can either use nqpath and qpath, but if you want to control exactly the number of q-point 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 q-point grid you can either use nqpath and qpath, but if you want to control exactly the number of q-point 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