24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 3765  |  回复: 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 的主题更新
信息提示
请填处理意见