VASP直接计算频率与用phonopy计算声子谱的区别

persei

卡开发发 发表于 2020-2-29 07:32
VASP直接计算的话没办法处理声子的BZ超过Gamma的情况，Phonopy的过程其实就是对超晶格算Gamma点的声子谱 ...

"Phonopy的过程其实就是对超晶格算Gamma点的声子谱再反折叠的过程"
麻烦你再说一说band unfolding，怎样决定supercell的尺寸合适与否？supercell 的布里渊区扩展成两倍，四倍，八倍等等，对应于supercell的能带图中，两个能带合成一个，四个能带合成一个，八个能带合成一个较小晶胞的能带图中的能带，是不是？可以三个supercell 的能带合成两个较小晶胞的能带吗？我想可以，但是还没想清楚怎么操作

卡开发发

persei 发表于 2023-2-21 10:19
"Phonopy的过程其实就是对超晶格算Gamma点的声子谱再反折叠的过程"
麻烦你再说一说band unfolding，怎样 ...

多大严格说要测试。

persei

卡开发发 发表于 2023-2-21 10:40
多大严格说要测试。

没明白你说的意思，你说了，VASP 计算声子谱只能算原点的色散关系，那么得到整个原胞对应的布里渊区的色散关系，比如说原胞含有8个原子，那么应该选多大的supercell 才合适？如果只是band unfolding，那么supercell布里渊区边界上的微扰的作用就自然带进了原胞布里渊区内，而这个本来不应该存在在那里，是不是？

sobereva

persei 发表于 2023-2-21 13:20
没明白你说的意思，你说了，VASP 计算声子谱只能算原点的色散关系，那么得到整个原胞对应的布里渊区的色 ...

做类似这种收敛性测试

persei

sobereva 发表于 2023-2-21 13:27
做类似这种收敛性测试

右下角的legend 表示 1*1*1  2*2*2  3*3*3  4*4*4 5*5*5  6*6*6 supercell， 是不是？那你自然以为6*6*6 的结果就是标准了，符合它自然就是正确的，你看到 3*3*3 和 6*6*6 的结果差不多，所以3*3*3 的supercell 就可以了，是不是这样的？另外，关于那个图，有个问题，声子能带上有40THz 的能量 （160meV），为什么声子态密度上只有 1THz 的信号？

卡开发发

persei 发表于 2023-2-21 13:20
没明白你说的意思，你说了，VASP 计算声子谱只能算原点的色散关系，那么得到整个原胞对应的布里渊区的色 ...

多大你可以不断扩胞去测，楼上给的图大概接近4x4x4才基本上收敛。声子态密度的图你看的不对，横轴是单位能量态的数目而不是能量区间。supercell上的扰动不会引入多余的东西。