About nph1l

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
hmao
Posts: 2
Joined: Mon Aug 10, 2020 6:35 am

About nph1l

Post by hmao » Mon Aug 10, 2020 6:50 am

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.
Attachments
phonon dispersion curve
phonon dispersion curve
trf2_5.in
(4.67 KiB) Downloaded 182 times

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: About nph1l

Post by ebousquet » Wed Aug 12, 2020 11:53 am

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

hmao
Posts: 2
Joined: Mon Aug 10, 2020 6:35 am

Re: About nph1l

Post by hmao » Tue Aug 18, 2020 11:20 am

Dear Eric,
Thanks a lot to reply me. :D
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

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: About nph1l

Post by ebousquet » Sat Aug 29, 2020 5:52 pm

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

Locked