24小时热门版块排行榜    

查看: 4415  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +41 医学老男孩 2026-08-13 94/4700 2026-08-18 10:20 by 医学老男孩
[论文投稿] 投稿咨询 +3 wwm09 2026-08-17 4/200 2026-08-18 09:52 by xs74101122
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 5/250 2026-08-18 08:03 by 基诺咪客
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见