[QE-users] phonon frequencies at Gamma using qplot = .true. along K-Gamma-L
Pietro Davide Delugas
pdelugas at sissa.it
Mon Jun 24 11:11:36 CEST 2019
Thanks for reporting
this is a bug, I started an issue on the GitLab repository. If you are
interested you may follow the discussion here.
<https://gitlab.com/QEF/q-e/issues/129>
Thanks again
Pietro
On 6/23/19 8:55 PM, Hari Paudyal wrote:
> Dear experts,
>
> To calculate the phonon along a specific path of BZ (K-Gamma-L) for
> Si, I used qplot = .true. and q_in_band_form = .true. (taking 2
> points in between). Surprisingly, I got the following frequencies at
> Gamma,
> q = ( 0.000000000 0.000000000 0.000000000 )
>
> **************************************************************************
> freq ( 1) = 11.535533 [THz] = 384.783979 [cm-1]
> freq ( 2) = 11.535533 [THz] = 384.783979 [cm-1]
> freq ( 3) = 11.535533 [THz] = 384.783979 [cm-1]
> freq ( 4) = 18.554595 [THz] = 618.914662 [cm-1]
> freq ( 5) = 18.554595 [THz] = 618.914662 [cm-1]
> freq ( 6) = 18.554595 [THz] = 618.914662 [cm-1]
> **************************************************************************
>
> On the other hand, If I calculate only at Gamma point, I get
> q = ( 0.000000000 0.000000000 0.000000000 )
>
> **************************************************************************
> freq ( 1) = -0.117104 [THz] = -3.906180 [cm-1]
> freq ( 2) = -0.117104 [THz] = -3.906180 [cm-1]
> freq ( 3) = -0.117104 [THz] = -3.906180 [cm-1]
> freq ( 4) = 15.035771 [THz] = 501.539319 [cm-1]
> freq ( 5) = 15.035771 [THz] = 501.539319 [cm-1]
> freq ( 6) = 15.035771 [THz] = 501.539319 [cm-1]
> **************************************************************************
>
> Why am I getting the large frequencies at Gamma if I calculate along
> K-Gamma-L? Do I need to specify Gamma always in the first?Why?
>
> My input files for the scf and phonon for Si are as the following.
>
> &control
> calculation = 'scf'
> prefix = 'Si'
> restart_mode = 'from_scratch'
> wf_collect = .true.
> pseudo_dir = './'
> outdir = './'
> /
> &system
> ibrav = 2
> celldm(1) = 10.261217
> nat = 2
> ntyp = 1
> ecutwfc = 45.0
> /
> &electrons
> diagonalization = 'david'
> mixing_beta = 0.3
> conv_thr = 1.0d-10
> /
> ATOMIC_SPECIES
> Si 28.0855 Si.pz-vbc.UPF
> ATOMIC_POSITIONS alat
> Si 0.00 0.00 0.00
> Si 0.25 0.25 0.25
> K_POINTS automatic
> 8 8 8 0 0 0
>
> Si
> &inputph
> outdir = './',
> prefix = 'Si',
> fildyn = 'Si.dyn',
> fildvscf = 'dvscf',
> ldisp = .true.,
> tr2_ph = 1.0d-14
> qplot = .true.
> q_in_band_form = .true.
> /
> 3
> -0.7500000 0.7500000 0.0000000 2 ! K
> 0.0000000 0.0000000 0.0000000 2 ! G
> -0.5000000 0.5000000 0.5000000 1 ! L
>
> Thank You,
> Hari Paudyal
> SUNY Binghamton
>
> _______________________________________________
> Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso)
> users mailing list users at lists.quantum-espresso.org
> https://lists.quantum-espresso.org/mailman/listinfo/users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20190624/8d6dbf32/attachment.html>
More information about the users
mailing list