单壁碳纳米管的声子谱问题求助

vayne

图1是我用vasp+phonopy花了17天算的单壁碳纳米管的声子谱..
图2是一篇博士论文里的碳纳米管声子谱
怎么差别这么大呀 我的图里声子支都是平缓的..
这是为什么呢
有懂的老师指导一下吗


我用Seekpath网站得到该模型有18个高对称点，如图3所示。 但是他的k-path路径却只含11个点，如图4，我有点不明白这个路径跟高对称点对我计算声子谱有什么影响吗？
所以在通过力常数文件生成band.data文件时
我的band.conf文件该如何编写?图5是我的模型，是MS里能建的最少原子的5，5-SWNT模型（20个C原子）

vayne

卡开发发 发表于 2019-11-23 23:00
C方向不需要真空。应该把A、B扩大。算声子的时候，C方向需要给dim 1 1 X，X最后足够大。

额 好的老师- - 谢谢您

卡开发发

vayne 发表于 2019-11-23 22:34
我感觉不是原胞，毕竟我不了解这方面的知识   我挂在一楼了，图5就是模型。老师，不是原胞应该也能做 ...

C方向不需要真空。应该把A、B扩大。算声子的时候，C方向需要给dim 1 1 X，X最后足够大。

vayne

卡开发发 发表于 2019-11-23 19:38
嗯，所以你还得看下你的模型。

我刚刚看到一篇硕士论文 说这种模型是原胞，不知道靠不靠谱..

vayne

卡开发发 发表于 2019-11-23 19:38
嗯，所以你还得看下你的模型。

我感觉不是原胞，毕竟我不了解这方面的知识   我挂在一楼了，图5就是模型。老师，不是原胞应该也能做计算吧，扩胞什么的无非就是加大计算量，还会有别的影响吗？

卡开发发

vayne 发表于 2019-11-23 17:52
最小重复单元好像是原胞吗，比如石墨烯的原胞就两个原子。我是不知道单壁碳纳米管的最小重复单元是什么样 ...

嗯，所以你还得看下你的模型。

vayne

get-it 发表于 2019-11-23 13:20
一维单原子链的布里渊区是[-\pi/a,\pi/a]，又因为w(k)=w(-k),所以只需要画[0,\pi/a]的色散关系，其实就是 ...

嗷嗷 懂了呢  我知道文献横坐标轴啥意思了 - - 谢谢您

vayne

卡开发发 发表于 2019-11-23 14:55
一维高对称点是很固定的，上面也有人提到了，其实是和一维氢原子链一样，就G-X的路径没必要找程序去生成 ...

最小重复单元好像是原胞吗，比如石墨烯的原胞就两个原子。我是不知道单壁碳纳米管的最小重复单元是什么样的...

卡开发发

vayne 发表于 2019-11-23 08:34
我结构优化就是用的单胞，40个原子的SWNT。所以，老师，高对称点这个东西应该怎么确定呢？网上还有vaspki ...

一维高对称点是很固定的，上面也有人提到了，其实是和一维氢原子链一样，就G-X的路径没必要找程序去生成。我不确定你怎么定义单胞，按道理应该是最小重复单元，我不确定生成出来的结构是否一定是单胞。

get-it

vayne 发表于 2019-11-23 08:38
你好，我看了看一维单原子链那节内容。我太愚钝了，没明白您说的什么意思。还请您指教。

一维单原子链的布里渊区是[-\pi/a,\pi/a]，又因为w(k)=w(-k),所以只需要画[0,\pi/a]的色散关系，其实就是文献中的路径。

vayne

get-it 发表于 2019-11-22 21:31
你看下固体物理中的一维单原子链就知道了。碳纳米管只是原子更多而已。

你好，我看了看一维单原子链那节内容。我太愚钝了，没明白您说的什么意思。还请您指教。

vayne

卡开发发 发表于 2019-11-22 11:36
1、应该是使用超晶格导致色散关系被折叠。所以分析的时候应该指定dim去折叠。比较合理的做法是结构优化也 ...

我结构优化就是用的单胞，40个原子的SWNT。所以，老师，高对称点这个东西应该怎么确定呢？自动生成的结构真空应该不充分，那看起来我得保证每次都加真空层了，谢谢老师。

get-it

你看下固体物理中的一维单原子链就知道了。碳纳米管只是原子更多而已。

卡开发发

vayne 发表于 2019-11-22 09:34
老师，上图中的声子谱，我的band.conf文件里BAND=0.0 0.0 0.0  0.0 0.5 0.0  0.0 0.5 0.125  0.0 0.0 0.1 ...

1、应该是使用超晶格导致色散关系被折叠。所以分析的时候应该指定dim去折叠。比较合理的做法是结构优化也单胞做，声子计算用phonopy指定dim扩成超晶格做，最后分析再指定dim去折叠。
2、seekpath我觉得应该可能针对3维周期来考虑的吧。
3、程序自动生成的结构是有真空的，大小是否充分我不好说。

vayne

卡开发发 发表于 2019-11-22 04:22
声子计算我不算特别了解，但我想大概是这样的问题：
（1）一维结构，高对称点可能只是Gamma-X这样（X沿着 ...

感觉我对真空层这个设置不太重视..是不是只要关于晶体的计算都得设置真空层呢老师