24小时热门版块排行榜    

查看: 3239  |  回复: 15
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 flyfly20069888 的 8 个金币

flyfly20069888

银虫 (小有名气)

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

请高手留步帮我看下这几个问题,谢了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyfly20069888

银虫 (小有名气)

在用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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyfly20069888

银虫 (小有名气)

问题请看1楼中的内容
3楼2009-08-25 12:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guifan

木虫 (著名写手)

★ ★ ★ ★ ★ ★
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应该就可以了。
4楼2009-08-25 13:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★ ★ ★ ★ ★
flyfly20069888(金币+3,VIP+0):谢谢回复! 8-25 15:07
wuchenwf(金币+3,VIP+0):感谢你的热心帮助 8-27 20:05
你做的真空层的法线是沿着x轴的吧?
你给出的两组k点的设置都是在yz平面内的,所以第一组将给出更多的能带的信息。当然,究竟你选取哪一组,取决于你在分析的时候使用哪个方向而已。所以在计算的时候是越多越好(当然多了更耗时)
5楼2009-08-25 14:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyfly20069888

银虫 (小有名气)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyfly20069888

银虫 (小有名气)

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

我做的真空层的法线是沿着Z轴的,我也不知道默认的为什么在yz平面内,按说是不是应该在xy面才对?
7楼2009-08-25 15:10:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guifan

木虫 (著名写手)

★ ★ ★ ★ ★ ★
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,截断能可以不一样,不过推荐设成一样,这样写文章时做一处说明就可以了
8楼2009-08-25 15:16:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyfly20069888

银虫 (小有名气)

引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 flyfly20069888 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:34 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:28 by Xc2FxHsHkBNc
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:06 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
请填处理意见