24小时热门版块排行榜    

查看: 3997  |  回复: 12

tc1788

铁虫 (小有名气)

[求助] Castep 设置的问题 学习良久,仍有不解,求解!!

有帖子(http://muchong.com/bbs/viewthread.php?tid=2103451&page=2)里说 “选随K增加,总能变化幅度较小的值.有可能K点很少反而能量低, k点多了能量高.K的测试,要看随K加密,总能波动<1meV/atom就可以了.”
“设置了不同的kpoints,但是都满足总能波动<1meV/atom, 这样的话,对于不同的k点,设置一下SCF tolerance 使其收敛”
以下是我测试K点时的结果
1 我设置的计算k点小时候能量较低,k点多了反而能量增加了,能量总在-1460.24 --1460.246间波动,请问这个是总能波动<1meV/atom了吗?若满足那么SCFtolerance 增加越大越好 是吗?
2 针对我的结构 关于k点设置选 separation=0.04时 k点为4 4 2,那么增加k点时我是按照每个数加1的方式成为 553 664 775 等等一直加,这样对吗? separation选0.05 好还是0.04好还是更高更低,这个有点困惑
3 我认为测试k点和截断能的步骤为
31 建立结构
3.2 几何优化,在该过程中先选定一个截断能,之后选小的k ,进行优化结构,再计算能量。之后逐渐增加k点,直到能量几乎不变,就确定k点了。之后增加截断能,再优化结构及计算能量。而这个SCF tolerance 以及以前看到过的searing的设置看到过但不知道什么时候用这两个,求助高人帮忙解答,也可以短消息给我,也会用金币感谢您!!!
4 界面的问题,根据实验值建立结构是否能成为直接计算能量的模型,不去优化,可以吗?如有解答能且附上文献,重谢。
回复此楼

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

第一原理资料汇编 计算MS等 科研软件资料 ms中castep模块疑难杂症
第一性原理 castep问题集合

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hwceng0816

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢交流 2012-08-27 15:05:54
tc1788: 金币+5, 有帮助, 谢谢鼓励! 1 2 3 问题中,您的见解能详细些吗 2012-08-27 16:19:06
你说的大体上都没有错,separation选0.05 好还是0.04好还是更高更低?这个当然是0.04更高,数值越小,说明格子划分的越细,精度越高。界面的问题,根据实验值建立结构是否能成为直接计算能量的模型,不去优化,可以吗?这一定要优化的,在第一性原理里面就算你严格按照实验值要建模也要优化,优化是所有计算的第一步,你会发现不经过优化的计算结果往往都无法接收。
最后你做科研很严谨,好多人根本不做收敛测试。
2楼2012-08-27 14:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

cosfrist

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
xueht987: 金币+1, 鼓励交流! 2012-08-30 09:22:27
1 我设置的计算k点小时候能量较低,k点多了反而能量增加了,能量总在-1460.24 --1460.246间波动,请问这个是总能波动<1meV/atom了吗?若满足那么SCFtolerance 增加越大越好 是吗?

你看看你体系的原子数,你的总能波动现在是6me ,如果体系的原子数大于6个,那么满足了

第二点上面已经回答了

SCF torlence是用于自洽收敛的,我计算的时候用的默认的,如果你对于收敛有特殊要求,一般可以设的更低些,但是过低的话可能导致收敛步数增加。
力学出身搞过生物医学工程的研究计算材料学的医学影像系青椒一枚
3楼2012-08-27 17:45:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tc1788

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by cosfrist at 2012-08-27 17:45:44
1 我设置的计算k点小时候能量较低,k点多了反而能量增加了,能量总在-1460.24 --1460.246间波动,请问这个是总能波动<1meV/atom了吗?若满足那么SCFtolerance 增加越大越好 是吗?

你看看你体系的原子数,你的 ...

ev 和mev换算应该是10^6那么要小数点后六位保持不变的意思吗?可为什么说现在是6me呢 ???困惑中???
4楼2012-08-27 18:11:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tc1788

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by cosfrist at 2012-08-27 17:45:44
1 我设置的计算k点小时候能量较低,k点多了反而能量增加了,能量总在-1460.24 --1460.246间波动,请问这个是总能波动<1meV/atom了吗?若满足那么SCFtolerance 增加越大越好 是吗?

你看看你体系的原子数,你的 ...

可是小数点后6位置不变似乎太精确了吧??困惑???
5楼2012-08-27 18:13:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tc1788

铁虫 (小有名气)

有帖子(http://muchong.com/bbs/viewthread.php?tid=2103451&page=2)里说 “选随K增加,总能变化幅度较小的值.有可能K点很少反而能量低, k点多了能量高.K的测试,要看随K加密,总能波动<1meV/atom就可以了.”
“设置了不同的kpoints,但是都满足总能波动<1meV/atom, 这样的话,对于不同的k点,设置一下SCF tolerance 使其收敛”
以下是我测试K点时的结果
1 我设置的计算k点小时候能量较低,k点多了反而能量增加了,能量总在-1460.24 ev--1460.246ev间波动,请问这个是总能波动<1meV/atom了吗?若满足那么SCFtolerance 增加越大越好 是吗?
2 针对我的结构 关于k点设置选 separation=0.04时 k点为4 4 2,那么增加k点时我是按照每个数加1的方式成为 553 664 775 等等一直加,这样对吗? separation选0.05 好还是0.04好还是更高更低,这个有点困惑
3 我认为测试k点和截断能的步骤为
31 建立结构
3.2 几何优化,在该过程中先选定一个截断能,之后选小的k ,进行优化结构,再计算能量。之后逐渐增加k点,直到能量几乎不变,就确定k点了。之后增加截断能,再优化结构及计算能量。而这个SCF tolerance 以及以前看到过的searing的设置看到过但不知道什么时候用这两个,求助高人帮忙解答,也可以短消息给我,也会用金币感谢您!!!
4 界面的问题,根据实验值建立结构是否能成为直接计算能量的模型,不去优化,可以吗?如有解答能且附上文献,重谢。
6楼2012-08-27 19:39:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cosfrist

木虫 (小有名气)

【答案】应助回帖


xueht987: 金币+1, 鼓励交流! 2012-08-30 09:22:41
xueht987: 应助指数+1 2012-08-30 09:22:51
引用回帖:
5楼: Originally posted by tc1788 at 2012-08-27 18:13:04
可是小数点后6位置不变似乎太精确了吧??困惑???...

这个不存在精确的,要看你使用的什么单位,比如说0.001 kv 感觉很精确了吧,可是其实就是1v 而已,还是1000ev。。所以不能单纯的看小数点后几位的问题
力学出身搞过生物医学工程的研究计算材料学的医学影像系青椒一枚
7楼2012-08-29 11:27:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tc1788

铁虫 (小有名气)

引用回帖:
7楼: Originally posted by cosfrist at 2012-08-29 11:27:16
这个不存在精确的,要看你使用的什么单位,比如说0.001 kv 感觉很精确了吧,可是其实就是1v 而已,还是1000ev。。所以不能单纯的看小数点后几位的问题...

我的单位是ev,改变量为0.006ev 平均每个原子也就是0.001ev。而似乎提到的波动量在1mev换算成ev是1000000ev。那这个能量波动太大了不是吗?我糊涂了?
8楼2012-08-30 09:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunboyccc

铁杆木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-09-06 20:08:14
引用回帖:
8楼: Originally posted by tc1788 at 2012-08-30 09:15:25
我的单位是ev,改变量为0.006ev 平均每个原子也就是0.001ev。而似乎提到的波动量在1mev换算成ev是1000000ev。那这个能量波动太大了不是吗?我糊涂了?...

1meV,这里m是读“兆”吗?还是读“毫”?意思差很远。
9楼2012-08-30 09:27:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cosfrist

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
tc1788: 金币+5, ★★★很有帮助, 非常感谢 2012-08-30 12:49:56
sunyang1988: 金币+1, 谢谢交流 2012-09-06 20:08:21
引用回帖:
8楼: Originally posted by tc1788 at 2012-08-30 09:15:25
我的单位是ev,改变量为0.006ev 平均每个原子也就是0.001ev。而似乎提到的波动量在1mev换算成ev是1000000ev。那这个能量波动太大了不是吗?我糊涂了?...

小写的1mev 应该是指 1*10-3 ev ,所以0.001ev = 1mev
如果大写的1M ev 才是 1*106 ev
力学出身搞过生物医学工程的研究计算材料学的医学影像系青椒一枚
10楼2012-08-30 10:48:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tc1788 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +5 爱擦汗的可乐冰 2026-02-28 6/300 2026-03-02 09:36 by 破风α
[考研] 材料工程269求调剂 +3 白刺玫 2026-03-02 3/150 2026-03-02 09:25 by 一休哥FU
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 085600 英一数二272求调剂 5+3 vida_a 2026-03-01 7/350 2026-03-02 07:51 by ms629
[考研] 求调剂 +5 yunziaaaaa 2026-03-01 6/300 2026-03-01 23:57 by ccp273206157
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 295求调剂 +7 19171856320 2026-02-28 7/350 2026-03-01 18:54 by 18137688336
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
信息提示
请填处理意见