24小时热门版块排行榜    

查看: 4001  |  回复: 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 界面的问题,根据实验值建立结构是否能成为直接计算能量的模型,不去优化,可以吗?如有解答能且附上文献,重谢。
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 4/200 2026-03-16 23:52 by zbzihdhd
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
信息提示
请填处理意见