24小时热门版块排行榜    

查看: 1192  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 03:40 by k06BKNblrGRH
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-17 02:14 by Tql5LQhh5rLK
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 02:04 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 01:52 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-16 21:39 by MHf40yA0Uhyb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 18:37 by rRvxEj3zorCy
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 QUjhNVAcOSff 2026-09-13 9/450 2026-09-16 18:36 by rRvxEj3zorCy
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:27 by LH5NkK5Ud4bh
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-16 02:38 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
请填处理意见