应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 分享与交流CASTEP 计算中截断能和K值的选择经验
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 8939  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

spphins510

新虫 (小有名气)

[交流] 分享与交流CASTEP 计算中截断能和K值的选择经验 已有6人参与

一直在小木虫索取，这次也想分享一下~虽然很基础的内容，但感觉对初学者，尤其是基础比较差的（我自己）还是会有帮助的。这些内容来源于软件说明书、小木虫的讨论等资源。而且如果有不对的地方，还请路过的大家能够指点一下~谢谢啦。

截断能
Energy cutoff：控制平面波基组的大小，平衡准确性和计算速度。在做吸附计算的时候，要分别计算气体分子，洁净表面和吸附后表面的，为了消除最大误差，该值不能够采取默认值，而应该统一。
一般400eV对于计算表面科学足够了。

K点
k-points： 确定布里渊区域采样的准确性；k点的质量对于金属系统尤其重要，在金属系统中，电子结构的快速变化可能发生在跨越费米能级的能带上。
经验：实空间越大，K空间就越小。金属体系需要精细的网格，绝缘体使用很少的k 点通常就可以。小单胞需要精细格点，大单胞很可能不需要。因此：单位晶胞内原子数很多（比如40-60 个）的绝缘体，可能仅需要一个（移动后的）k 点。另一方面，面心立方的铝可能需要上万个k点以获得好的DOS。对于孤立原子或分子的超晶胞，仅需要在Gamma点计算。对于表面（层面）的超晶胞计算，仅需要（垂直于表面）z 方向上有1 个k 点。
取值依据：（1）a×ka ≈ b×kb ≈ c×kc 或者说 1/a : 1/b : 1/c ≈ ka : kb : kc
                 （2）ka~15 Å 为绝缘体
                          ka~20 Å 为半导体
                          ka~25 Å 为simple metals
                          ka~30 Å 为d band metals

截断能和k值的收敛性测试：
（1）方法：首先测试Ecut，把K点设置成足够高的值，以步长50eV逐步增加Ecut，查看总能的变化，当变化小于0.001 eV/atom（dE/dEcut/atom）的时候就已经足够了。然后固定Ecut为测试好的值，逐步增加K值，以同样的收敛标准测试就行了。
（2）对于固体表面计算，和在切面之前对晶体做收敛性分析相比，在切面后做收敛性测试更合理。
（3）但要注意改变k值能量可能不是平稳的收敛。
回复此楼

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

MS学习积累

» 本帖已获得的红花（最新10朵）

kecy

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2019-04-21 21:27:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

网易云之夜

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
收敛性测试那部分解释，我看不太懂
赞一下 回复此楼
2楼2019-05-05 13:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kecy

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
（dE/dEcut/atom）这个具体数值在哪里看呢？比如dE， dEcut，atom分别指什么。atom是体系总原子数吗？
赞一下 回复此楼
never say die
3楼2020-03-12 20:56:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by kecy at 2020-03-12 20:56:34
（dE/dEcut/atom）这个具体数值在哪里看呢？比如dE， dEcut，atom分别指什么。atom是体系总原子数吗？

手动做，调整Ecut看能量变化。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
4楼2020-03-13 13:21:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DR_Zhou

新虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
请教一下，做收敛性测试，固定截断能，改变k值，在从校到大增加k值的过程中，k值该怎样设置呢？比如，现在已经算了k为1x1x1，那么下一步k应该设置成2x2x2还是1x1x2（2x1x1或者1x2x1）？
赞一下(1人) 回复此楼
5楼2020-05-23 14:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fgmsyutong

木虫 (文坛精英)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
截断能，改为custom（定值），多设几个，做出总能量随截断能关系曲线，斜率趋于0即可

发自小木虫Android客户端
赞一下 回复此楼
3d104s1
6楼2020-05-23 17:01:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fgmsyutong

木虫 (文坛精英)
引用回帖:
6楼: Originally posted by fgmsyutong at 2020-05-23 17:01:11
截断能，改为custom（定值），多设几个，做出总能量随截断能关系曲线，斜率趋于0即可

k值，也是custom，找能量最小值对应的k

发自小木虫Android客户端
赞一下(1人) 回复此楼
3d104s1
7楼2020-05-23 17:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fgmsyutong

木虫 (文坛精英)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
还有，研究表面，最后一个数字改成1，比如，如果晶体用k值999得到能量最小，表面就改为991

发自小木虫Android客户端
赞一下(1人) 回复此楼
3d104s1
8楼2020-05-23 17:06:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmaobi

铁虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
The value of the dEtot/d lnEcut gives a good indication of the convergence...
... the value reported in the CASTEP output file divided by the number of atoms...
这里是不是应该dEtot/d lnEcut/Atoms才对啊？在官方手册上看到的
赞一下 回复此楼
长路漫漫，与贤为伴。
9楼2022-02-07 13:49:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 spphins510 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +4 李青山山山 2026-04-07 4/200 2026-04-07 16:59 by 啵啵啵0119
[考研] 本科郑州大学，一志愿华东师范大学282求调剂 +18 熊哥xtk 2026-04-07 21/1050 2026-04-07 15:54 by 乔哒哒哒
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 13/650 2026-04-06 23:48 by 蓝云思雨
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 10/500 2026-04-06 19:47 by jean5056
[考研] 材料专硕322 +11 哈哈哈吼吼吼哈 2026-04-05 11/550 2026-04-06 14:07 by lqwchd
[考研] 专硕304找调剂，一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 考研生物学考A区211，初试322，科目生化和生物综合，求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 求 +5 化工专硕323分 2026-04-04 5/250 2026-04-05 08:02 by 544594351
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊！
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 化学070300-总分378-求调剂 +5 挪椅子的泡泡糖 2026-04-02 5/250 2026-04-02 22:20 by ZXlzxl0425
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
信息提示
关闭
请填处理意见
关闭