Phonon calculations with IBRION=8

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
deepashri
Newbie
Newbie
Posts: 12
Joined: Wed Jan 30, 2013 6:36 am

Phonon calculations with IBRION=8

#1 Post by deepashri » Wed Mar 08, 2017 10:30 am

Dear All,

I am trying to calculate phonon frequencies for a 2-D sheet with the help of DFPT. But the results show that the acoustic sum rule is violated i. e. there are finite values of phonon frequencies at the gamma point in the negative region.

Following are some of the input parameters in the INCAR file,

########################################################
ISTART = 0
INIWAV = 1
NWRITE = 2
LWAVE = .TRUE.
PREC = accurate
EDIFF = 1.0E-6
EDIFFG = -0.001

ISMEAR = 1
SIGMA = 0.001

ICHARG = 2
LCHARG = .TRUE.

ISYM = 1

IBRION = 8
ISIF = 3
NFREE = 2
NSW = 1
POTIM = 0.02

LPETIM = .FALSE.
LREAL = .FALSE.
LORBIT = 1
LPLANE = .TRUE.
LRPA = .FALSE.

###############################################################

The KPOINTS file is as follows,

###############################################################
Automatic mesh
0
Gamma
8 8 1
0. 0. 0.

##############################################################

To perform the phonon calculations, the structure is fully relaxed with ISIF=3 but the z-axis was kept fixed with the help of "selective dynamics" flag. The vacuum in the z-direction was set to 20 Angstrom. Further, the phonon calculations were performed on 2x2x1 supercell of the relaxed structure.

What could be the reason behind the incorrect results that I am getting?

Thanks in advance.

Best,
Deepashri.

User avatar
SKM
Jr. Member
Jr. Member
Posts: 72
Joined: Wed Oct 30, 2019 5:39 am
License Nr.: 5-516

Re: Phonon calculations with IBRION=8

#2 Post by SKM » Tue May 11, 2021 4:25 am

hi
this is very old thread.
but curious to know, how did you resolve the issue, if you did.
Regards

Post Reply