计算化学公社

标题: 单壁碳纳米管的声子谱问题求助 [打印本页]

作者
Author: vayne 时间: 2019-11-21 15:16
标题: 单壁碳纳米管的声子谱问题求助
本帖最后由 vayne 于 2019-11-23 22:32 编辑

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

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

作者
Author: vayne 时间: 2019-11-21 15:48
我感觉是band.conf的问题

作者
Author: vayne 时间: 2019-11-21 15:50
本帖最后由 vayne 于 2019-11-21 15:57 编辑

BAND参数该怎么设置？有人说是DIM*BAND=高对称点坐标

作者
Author: 卡开发发 时间: 2019-11-22 04:22
声子计算我不算特别了解，但我想大概是这样的问题：
（1）一维结构，高对称点可能只是Gamma-X这样（X沿着一维延伸的方向），沿着其他的方向得到的应该都是平的。文献的k∈[0,0.5]的取法就是这样的原因。
（2）这个模型也应该是纳米管周围留有足够的真空，而计算声子谱的时候一维周期的方向应该要给足够的大的超晶格，然后计算的时候要指定dim给定超晶格的大小来将色散关系反折叠到单胞。

作者
Author: vayne 时间: 2019-11-22 09:34

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

老师，上图中的声子谱，我的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.125  0.0 0.0 0.0  -0.5 0.0 0.125  -0.5 0.5 -0.125  0.0 0.5 0.0  -0.5 0.5 0.0  -0.5 0.0 0.0  0.0 0.0 0.0

后来听别人的改成了BAND=0.0 0.0 0.0 0.0 0.0 0.5  ，就跟文献中的图曲线走势差不多，但是k的范围却是0~0.25，我也不是很懂其中的原理..
您的意思是一维结构的高对称点只是（x,y）吗  那这样说来我用seekpath网站算得的高对称点及路径都是错的对不？
而且我计算的时候没有加真空层...

作者
Author: vayne 时间: 2019-11-22 09:36

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

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

作者
Author: 卡开发发 时间: 2019-11-22 11:36

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、程序自动生成的结构是有真空的，大小是否充分我不好说。

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

作者
Author: vayne 时间: 2019-11-23 08:34

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

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

作者
Author: vayne 时间: 2019-11-23 08:38

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

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

作者
Author: get-it 时间: 2019-11-23 13:20

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

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

作者
Author: 卡开发发 时间: 2019-11-23 14:55

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

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

作者
Author: vayne 时间: 2019-11-23 17:52

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

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

作者
Author: vayne 时间: 2019-11-23 17:57

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

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

作者
Author: 卡开发发 时间: 2019-11-23 19:38

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

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

作者
Author: vayne 时间: 2019-11-23 22:34

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

我感觉不是原胞，毕竟我不了解这方面的知识

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

作者
Author: vayne 时间: 2019-11-23 22:52

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

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

作者
Author: 卡开发发 时间: 2019-11-23 23:00

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

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

作者
Author: vayne 时间: 2019-11-24 09:40

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

额好的老师- - 谢谢您

欢迎光临计算化学公社 (http://ccc.keinsci.com/)