VASP-DFPT & phonopy calculation 计算声子谱求助
1. 第一次计算声子谱，想重复vasp手册上面的NaCl例子，但是做出来的图线有点多，
以下是我的band,conf里面的参数
ATOM_NAME = Na Cl
FORCE_CONSTANTS = READ
DIM = 2 2 2
BAND = 0.00 0.00 0.00 0.50 0.00 0.00 0.50 0.50 0.00 0.00 0.00 0.00 0.50 0.50 0.50 波矢路径
BAND_POINTS = 200
INCAR里面的参数如下：
PREC = Accurate
ENCUT = 500
IBRION = 8
EDIFF = 1.0e-08
IALGO = 38
ISMEAR = 0; SIGMA = 0.1
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.
2. 文献上面说如果晶体中有N个原胞，每个原胞中有n个原子，那么一共有3nN种不同频率模式。这3nN个独立振动分为3n支，其中3支是声学波，其余3（n-1）支是光学波，请问波支是什么，是波的条数？还是什么？还有声学波与光学波怎么区分？
TZp.png
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
算例中的NaCl 一个原胞里面有三个原子吧，所以有6个振动模式。但是你这个怎么会这么多振动模式？看输入参数没看出有什么问题？NaCl 分不分原胞跟惯用胞？
声子谱最后画出来的是原胞的色散关系 所以振动模式应该是3n 支（ n表示原胞中的原子数）至于怎么区分声学波 光学波 建议你看看黄昆的固体物理 晶格振动那一章
我的原胞和教程里面的一模一样，就是如下：
Na Cl
1.00000000000000
5.6903014761756712 0.0000000000000000 0.0000000000000000
0.0000000000000000 5.6903014761756712 0.0000000000000000
0.0000000000000000 0.0000000000000000 5.6903014761756712
Na Cl
4 4
Direct
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.5000000000000000 0.5000000000000000
0.5000000000000000 0.0000000000000000 0.5000000000000000
0.5000000000000000 0.5000000000000000 0.0000000000000000
0.5000000000000000 0.5000000000000000 0.5000000000000000
0.5000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.5000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.5000000000000000
还有，您说的NaCl胞有三个原子，那应该画出来是九条线，而不是六条呀
，
这不得了，你用的是惯用胞。8个原子。
第一次计算声子谱，所以问题比较小白：
1. PRIMITIVE这个参数是在band.conf中设置吗，加入PRIMITIVE= 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0，后执行% phonopy --dim="2 2 2" -c POSCAR-unitcell band.conf，老是报错。
2. 还有POSCAR中初胞原子的坐标格式也是上面的POSCAR一致吗？只是相应的原子数减小，坐标个数减少？
3. 最主要一点，第一次求助，如何给你金币
你应该根据phonopy手册上面的算例来输入命令。PRIMITIVE是将惯用胞转换成原胞的矩阵，这个矩阵你是否有输对? 另外我知道为什么处理band.conf的时候还需要用到--dim这个扩胞命令? 其次 我不知道你第二个问题是什么意思 最后 不用给我金币了 好好看看手册，根据手册上面说的一步步来应该没有问题才是啊。