the steps to remove the strange peak in phonon bands.

(1) Calculate the dielectric constant with the following INCAR (it's an 
example for PBE functional):

general:dielectric ele
 System =
 KSPACING = 0.1
 KGAMMA = .TRUE.
 ISTART = 0 ; ICHARG = 2
 ENCUT? =??? 600
 ISMEAR = 0; SIGMA = 0.01
 IBRION = 2
 EDIFF = 1E-8
 EDIFFG = -0.01
 LEPSILON = .TRUE.
 LWAVE = .F.
 LCHARG = .F.
 #LPEAD=.TRUE.

(2) After finishing the calculation;

run the script born.sh as

bash born.sh (atoms number) (filename);

then you can get (filename).born

(3) copy this (filename).born to the alamode folder.

Then

add NONANALYTIC = 2; BORNINFO = (filename).born in phband.in

run anphon phband.in again




the scripts to generate the alamode input files with a relaxed structure file.

Put POSCAR(a relaxed structure file), POTCAR,INCAR and attached two 
scripts (revised from Prof. Tadano's scripts) in your calculation folder 
and then follow the steps:

1. python to_stdcell.py POSCAR? > pri.vasp

2. vi makedisp.py;

to modify ALAMODE_root to your install pathway, prefix, and 
scaling_matrix (supercell size you use).

3. python makedisp.py
Then you can get structures with displacement for force calculations in 
VASP.