24小时热门版块排行榜    

查看: 1160  |  回复: 6

Gina88

木虫 (正式写手)


[交流] 【求助完毕】CASTEP怎样选择ENCUT参数?

在用CASTEP进行计算的时候要进行参数测试吗?
在进行参数测试的时候怎样看总能随参数的变化啊?

我是这样进行的,不知道对不对啊。
将energy cutoff选为300eV,进行单点能计算,之后看.castep文件中的NB est. 0K energy (E-0.5TS) 并记录下来。
然后将energy cutoff改为350eV,400eV,450eV,500eV,550eV,600eV,分别进行单点能计算,并记录NB est. 0K energy (E-0.5TS) 的值。
如果将energy cutoff按照从小到大的顺序排列,NB est. 0K energy (E-0.5TS) 也按照energy cutoff的顺序排列,去比较相邻NB est. 0K energy (E-0.5TS) 的差值,按理说,当这个差值小于0.001eV的时候,energy cutoff就是计算所需的参数。

可是我计算的响铃NB est. 0K energy (E-0.5TS) 差值并不是逐步减小的,问题出在哪里了?

[ Last edited by cenwanglai on 2011-3-27 at 21:56 ]
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
youzhizhe(金币+2): 谢谢交流。 2011-03-28 11:15:52
看波动幅度。能量收敛相关参数不止一个。测试截断能的时候K点要选大一点。
2楼2011-03-27 21:55:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)


引用回帖:
Originally posted by cenwanglai at 2011-03-27 21:55:52:
看波动幅度。能量收敛相关参数不止一个。测试截断能的时候K点要选大一点。

首先非常感谢区长,本应以金币表示致谢的,现在不能了。
但还是想继续问区长几个问题:

这个能量的波动幅度要小到什么值为止,有没有一个评判的标准呢?

进行截断能测试时,K点大概要取多大呢(下面是我截断能测试时的k点,结构是一维纳米线)?

-------------------------------
                              k-Points For BZ Sampling
                           -------------------------------
                       MP grid size for SCF calculation is  1  1 12
                         Number of kpoints used =           6

我在进行参数测试时,计划的是先测试截断能,后测试K点,最后测试sigma,(如果sigma变化,再测试K点,测试sigma,直到sigma不变为止)。请问:版主平时计算时还测试其它参数吗?都测试些什么参数呢?测试的顺序应该是什么样的呢?

我在进行截断能测试时,其它的参数都没有变化。

[ Last edited by Gina88 on 2011-3-28 at 10:09 ]
3楼2011-03-28 09:05:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+2): 谢谢交流。 2011-03-28 11:16:02
测试顺序这样可以。二维结构没有算过。看计算什么,一般分析总能(差)0.01eV就可以了。
4楼2011-03-28 09:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)


引用回帖:
Originally posted by cenwanglai at 2011-03-28 09:42:25:
测试顺序这样可以。二维结构没有算过。看计算什么,一般分析总能(差)0.01eV就可以了。

谢谢区长:)

那您在计算时,除了能量截断,K点,sigma外,其它参数还测试吗?
您在测试截断能时,k点一个方向取多少值才够呢?一般一个方向是去奇数个K点还是取偶数个k点呢?
5楼2011-03-28 10:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+2): 谢谢交流。 2011-03-28 11:17:13
11埃的长度上最少取3个。去点密度>=0.4/埃。尺寸越小取点越多。奇数偶数都可以,保证Gamma点能够取到。具体参看手册.

取点的方式是:尽量取高对称点,这样可以减少K点密度却能达到较高的计算精度。根据POSCAR的对称性来取K点,可以比较一下3-4种取点方法,测试一下,看看哪种方法好取点最经济。

[ Last edited by cenwanglai on 2011-3-28 at 11:19 ]
6楼2011-03-28 11:02:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awmc2008

至尊木虫 (文坛精英)


有收获,谢谢。
7楼2011-03-28 18:31:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 21/1050 2026-08-31 18:04 by 异氰酸酯
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见