首页 > Web开发 > 详细

phonon dispersion spectra VASP 5.3.3 + phonopy 1.10

时间:2015-12-02 20:48:03      阅读:410      评论:0      收藏:0      [点我收藏+]

1. Installation of phonopy

  ubuntu12.04 + python3.5.0

  (1) python 3.5.0

    extract the pack, and then,

    

cd Python-3.5.0/
sudo ./configure --prefix=/opt/python3.5.0
sudo make
sudo make install

 

 

  (2) some pack

      

sudo apt-get install python2.7-dev python-numpy python-matplotlib python-tk python-lxml python-yam

 

     maybe some of them are useless, uhmmmmmm.....even I don‘t know what they mean.

  (3) phonopy 

      

sudo python setup.py install

 

2. VASP & phonopy, DFPT method

  (1) relax the coordinate with a high accuracy by VASP. We will use the POSCAR, maybe the CONTCAR.

  (2) copy and rename the coordinate file to POSCAR-unitcell

  (3) some command

      

phonopy -d --dim="2 2 2" -c POSCAR-unitcell

mv SPOSCAR POSCAR

 

  (4) now we will calculate force constants of the supercell by VASP

    IMPORTANT: IBRION = 8, NSW = 1, LREAL = .FALSE

    NOTICE:    NO NPAR!

    the INCAR may be such like

          

PREC = Accurate
ENCUT = 500
61
IBRION = 8
EDIFF = 1.0e-08
IALGO = 38
ISMEAR = 0; SIGMA = 0.1
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.

 

  (5)after a longtime caculation, we will find the vasprun.xml.

        

phonopy --fc vasprun.xml

 

    we get the FORCE-CONSTANT

    then, we should make a file named xxx.conf(maybe band.conf), which contains some parameters.

        it may like this

ATOM_NAME = Na Cl
DIM = 2 2 2
PRIMITIVE_AXIS = 0.0 0.5 0.5  0.5 0.0 0.5  0.5 0.5 0.0
BAND = 0.0 0.0 0.0  0.5 0.0 0.0  0.5 0.5 0.0  0.0 0.0 0.0  0.5 0.5 0.5
FORCE_CONSTANTS = READ

 

    then, we use this commander

        

phonopy --dim="2 2 2" -c POSCAR-unitcell band.conf

 

    we get the band.yaml, then

bandplot --gnuplot band.yaml > phono.dat

 

  (6)finally, put the phono.dat into origin, matlab, MS excel and etc. 

 

Reference: phonon manual.

phonon dispersion spectra VASP 5.3.3 + phonopy 1.10

原文:http://www.cnblogs.com/HuaNeedsPills/p/5014013.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!