计算化学公社

标题: 请问碳纳米管的KPOINTS该如何取值。 [打印本页]

作者
Author: vayne 时间: 2019-10-11 09:47
标题: 请问碳纳米管的KPOINTS该如何取值。
本帖最后由 vayne 于 2019-10-11 09:49 编辑

KPOINTS的建议取值如下图，即每个晶格矢量的长度乘以这个方向上的k点数。

但碳纳米管根据导电性能的不同，可分为金属型和半导体型。
请问各位老师金属型碳纳米管和半导体型碳纳米管该如何取KPOINTS？

作者
Author: 卡开发发 时间: 2019-10-11 10:29
拿半导体的猜，然后做收敛性测试。

作者
Author: vayne 时间: 2019-10-11 11:04

卡开发发发表于 2019-10-11 10:29
拿半导体的猜，然后做收敛性测试。

好的 - - ，谢谢您。再次请教您，我算了a,b,c三个的晶胞长度分别为16.907A，16.907A，2.4595A；假设为金属型，那么K点是否是3 3 20呢？

作者
Author: 卡开发发 时间: 2019-10-11 11:16

vayne 发表于 2019-10-11 11:04
好的 - - ，谢谢您。再次请教您，我算了a,b,c三个的晶胞长度分别为16.907A，16.907A，2.4595A；假设为金 ...

那要看是否考虑纳米管之间的相互作用，如果考虑的话，a、b两个轴向的k也应当纳入测试；若不考虑，则a、b方向应该取一个比较大的值避免管子之间相互作用，此时a、b轴向的k取1即可。如果只看尺寸，a、b轴向的k取1大致是够了，测试一下c的就行。

作者
Author: vayne 时间: 2019-10-11 14:34

卡开发发发表于 2019-10-11 11:16
那要看是否考虑纳米管之间的相互作用，如果考虑的话，a、b两个轴向的k也应当纳入测试；若不考虑，则a、b ...

好深奥啊老师。管内碳碳原子间的相互作用肯定要考虑的，但是如图所示的取K值方法，不是根据尺寸的比例来吗？

作者
Author: 卡开发发 时间: 2019-10-11 14:37

vayne 发表于 2019-10-11 14:34
好深奥啊老师。管内碳碳原子间的相互作用肯定要考虑的，但是如图所示的取K值方法，不是根据尺寸的比例来 ...

有些材料并不是各向同性的。这个值只是参考值，严格来说还是应该做收敛性测试，这些值作为k收敛性测试的初始参考比较合适。

作者
Author: vayne 时间: 2019-10-11 14:44
好的老师 - - 谢谢您的答疑

作者
Author: vayne 时间: 2019-10-11 14:44

卡开发发发表于 2019-10-11 14:37
有些材料并不是各向同性的。这个值只是参考值，严格来说还是应该做收敛性测试，这些值作为k收敛性测试的 ...

好的老师谢谢您的答疑

作者
Author: vayne 时间: 2019-10-18 15:37

卡开发发发表于 2019-10-11 11:16
那要看是否考虑纳米管之间的相互作用，如果考虑的话，a、b两个轴向的k也应当纳入测试；若不考虑，则a、b ...

老师，K点设置1X1X14会不会太少了。我现在准备用VASP结合Phonony计算碳纳米管的声子谱，结构优化的时候看别人的案例K点都几百个，弄的我有点慌。。

作者
Author: 卡开发发 时间: 2019-10-18 17:34

vayne 发表于 2019-10-18 15:37
老师，K点设置1X1X14会不会太少了。我现在准备用VASP结合Phonony计算碳纳米管的声子谱，结构优化的时候看 ...

收敛性测试做好就没有必要慌。包括两个“非周期性”的轴向的晶格尺寸，“周期性”轴向的k点，动能截断，尽量看一下体系的最大力也能随着这些量收敛到足够高的精度，如果满足这些条件，声子谱按道理不会有虚频。

注意，计算声子谱的时候“周期性”的那个方向最好延伸足够大。

作者
Author: vayne 时间: 2019-10-18 21:12

卡开发发发表于 2019-10-18 17:34
收敛性测试做好就没有必要慌。包括两个“非周期性”的轴向的晶格尺寸，“周期性”轴向的k点，动能截断， ...

好的老师谢谢您

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