[交流] 【求助】请教K点及K点路径的设置

» 收录本帖的淘帖专辑推荐

计算化学

» 猜你喜欢

博士读完未来一定会好吗已经有6人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

1楼2009-08-25 12:04:01

flyfly20069888

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 440.2
散金: 5
帖子: 93
在线: 178.8小时
虫号: 740525
注册: 2009-04-05
性别: MM
专业: 微电子学与固体电子学

在用CASTEP做第一性原理计算时，看到文献中写到：
1、The Brillouin zone integration is performed within the Monkhorst–Pack scheme using 4×4×2 k points.
请问这里的4×4×2 k points是几何优化时的K点，还是算能带结构、态密度等性质时的K点？还是说几何优化、能带结构、态密度时的K点都是4×4×2  ?

2、The optimal atomic positions are determined until the magnitude of the forces acting on all atoms is less than 0.001 eV /A
我理解为这是在几何优化时Max. force应小于0.001 eV /A，但另外三个收敛参数Energy ，  Max. stress 和Max. displacement文献中没说明，该如何设置呢？

3、我计算的是一个面，在上面加真空层，在计算能带结构时，k-point path 的默认为：
Sym From             to
G 0.000 0.000 0.000 0.000  0.500 0.000
F 0.000 0.500 0.000 0.000  0.500 0.500
Q 0.000 0.500 0.500 0.000  0.000 0.500
Z 0.000 0.000 0.500 0.000  0.000 0.000
这样有四条路径，能带结构的横坐标上会有5个特殊点，
我认为虽然加了真空层，但研究的还是个面，k-point path应设为下面三条路径：
G 0.000 0.000 0.000 0.000  0.500 0.000
F 0.000 0.500 0.000 0.000  0.500 0.500
Q 0.000 0.500 0.500 0.000  0.000 0.000
不知道这样想对不对?

请教高人指点，感激不尽!

2楼2009-08-25 12:04:41

flyfly20069888

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 440.2
散金: 5
帖子: 93
在线: 178.8小时
虫号: 740525
注册: 2009-04-05
性别: MM
专业: 微电子学与固体电子学

问题请看1楼中的内容

3楼2009-08-25 12:06:43

guifan

木虫 (著名写手)

应助: 9 (幼儿园)
金币: 2279.1
帖子: 1083
在线: 73小时
虫号: 717940
注册: 2009-03-08
专业: 金属结构材料

★ ★ ★ ★ ★ ★
flyfly20069888(金币+3,VIP+0):谢谢你的解答，对我很有帮助！ 8-25 14:40
wuchenwf(金币+3,VIP+0):感谢你的热心帮助 8-27 20:05

1，是几何优化时的k点。能带结构和态密度的k点可以在优化之后进行Analysis的时候再设，这是因为优化计算之后哈密顿量中已经包含了能带结构和态密度的信息，Analysis时只是根据所设参数往上填充。
2，没错，是几何优化时Max. force，0.001 eV /A已是比较高的设置。若要保证计算精度的话，其他三项也应设得高一点，用默认的Ultrafine应该就可以了。

赞一下(14人)

4楼2009-08-25 13:24:15

xirainbow

木虫 (正式写手)

1ST强帖: 2
应助: 0 (幼儿园)
贵宾: 0.03
金币: 2132.6
红花: 2
帖子: 775
在线: 6.1小时
虫号: 431211
注册: 2007-08-11
专业: 光谱学

★ ★ ★ ★ ★ ★
flyfly20069888(金币+3,VIP+0):谢谢回复！ 8-25 15:07
wuchenwf(金币+3,VIP+0):感谢你的热心帮助 8-27 20:05

你做的真空层的法线是沿着x轴的吧？
你给出的两组k点的设置都是在yz平面内的，所以第一组将给出更多的能带的信息。当然，究竟你选取哪一组，取决于你在分析的时候使用哪个方向而已。所以在计算的时候是越多越好（当然多了更耗时）

赞一下(13人)

5楼2009-08-25 14:03:20

flyfly20069888

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 440.2
散金: 5
帖子: 93
在线: 178.8小时
虫号: 740525
注册: 2009-04-05
性别: MM
专业: 微电子学与固体电子学

引用回帖:

Originally posted by guifan at 2009-8-25 13:24:
1，是几何优化时的k点。能带结构和态密度的k点可以在优化之后进行Analysis的时候再设，这是因为优化计算之后哈密顿量中已经包含了能带结构和态密度的信息，Analysis时只是根据所设参数往上填充。
2，没错，是几何 ...

谢谢你的回复，我还想问一下：
1、在选了Properties里的Density of states 后，会有k-point set 这项，要不要也设成4x4x2, 还是可以重新设置？
2、我是先优化，然后进行性质的计算的，优化时和计算性质时的截止能等是否应保持一致？

期待你的回复

6楼2009-08-25 14:51:48

flyfly20069888

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 440.2
散金: 5
帖子: 93
在线: 178.8小时
虫号: 740525
注册: 2009-04-05
性别: MM
专业: 微电子学与固体电子学

引用回帖:

Originally posted by xirainbow at 2009-8-25 14:03:
你做的真空层的法线是沿着x轴的吧？
你给出的两组k点的设置都是在yz平面内的，所以第一组将给出更多的能带的信息。当然，究竟你选取哪一组，取决于你在分析的时候使用哪个方向而已。所以在计算的时候是越多越好（ ...

我做的真空层的法线是沿着Z轴的,我也不知道默认的为什么在yz平面内，按说是不是应该在xy面才对？

7楼2009-08-25 15:10:21

guifan

木虫 (著名写手)

应助: 9 (幼儿园)
金币: 2279.1
帖子: 1083
在线: 73小时
虫号: 717940
注册: 2009-03-08
专业: 金属结构材料

★ ★ ★ ★ ★ ★
flyfly20069888(金币+3,VIP+0):再次感谢你的热心帮助 8-25 16:11
wuchenwf(金币+3,VIP+0):感谢你的热心帮助 8-27 20:05

引用回帖:

Originally posted by flyfly20069888 at 2009-8-25 14:51:

谢谢你的回复，我还想问一下：
1、在选了Properties里的Density of states 后，会有k-point set 这项，要不要也设成4x4x2, 还是可以重新设置？
2、我是先优化，然后进行性质的计算的，优化时和计算性质时的 ...

1,可以重新设置，根据你的需要来
2，截断能可以不一样，不过推荐设成一样，这样写文章时做一处说明就可以了

赞一下(13人)

8楼2009-08-25 15:16:32

positron

木虫 (职业作家)

应助: 0 (幼儿园)
贵宾: 0.025
金币: 705.6
散金: 91
红花: 1
帖子: 4541
在线: 310.1小时
虫号: 274184
注册: 2006-08-26
专业: 核物理

引用回帖:

Originally posted by flyfly20069888 at 2009-8-25 12:04:
在用CASTEP做第一性原理计算时，看到文献中写到：
1、The Brillouin zone integration is performed within the Monkhorst–Pack scheme using 4×4×2 k points.
请问这里的4×4×2 k points是几何优化时的K点 ...

俺也初学，对你的第一点的软件中的设置都不明白

那个FFT grid set 中的设置和M-P grid中的grid的grid parametersd 的参数设置中的格点对布里渊区的积分—— “The Brillouin zone integration”，用的是哪个格点设置呢？
或者是对应的哪个格点的积分？
（因为我看了计算出的电子密度的文件，感觉电子密度数据对应的是FFT grid，但是它division的却是 a, b, c（应该是正格矢吧））

9楼2009-08-26 06:36:53

flyfly20069888

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 440.2
散金: 5
帖子: 93
在线: 178.8小时
虫号: 740525
注册: 2009-04-05
性别: MM
专业: 微电子学与固体电子学

引用回帖:

Originally posted by guifan at 2009-8-25 15:16:

1、在选了Properties里的Density of states 后，会有k-point set 这项，要不要也设成4x4x2, 还是可以重新设置？

1,可以重新设置，根据你的需要来
2，截断能可以不一样，不过推荐设成一样，这样写文章时做一处说明就可以了

还想请教你一下：
1、文献中没有说明Density of states 的 k-point set ，如果按我自己的设置算的话，做出的图形会不会与文献中的不同？
2、算能带结构时的k-point path不同，是否也会影响做出的图形与文献中的图形进行比较？

10楼2009-08-26 12:49:29