24小时热门版块排行榜    

查看: 3537  |  回复: 13

yqkwww

金虫 (正式写手)

[求助] 求助:用castep 做掺杂计算时,Geometry Optimization的作用?

在用castep进行计算时,遇到一个问题,困扰我很久,想请各位高手帮我解答一下。
计算过程一般我时这样的:
第一步进行结构优化,task项选择Geometry Optimization, Properties项没有选择,第二步计算能量时,task项选择energy,Properties项选择band structure and density of States。请问第一步有什么作用?优化出来时的是cutoff值和K-point 这两个参数么?第一步与第二步之间有什么联系?
回复此楼
change myself if I am stocked!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liu108

木虫 (正式写手)

【答案】应助回帖

上面这两个并不是优化参数,这是为了保证优化的准确性,人为设置的cutoff和k点要准确,这两个参数设置的不合理的话,优化出来的构型也不是最稳定构型。所谓优化参数就是把结构优化的最优构型,即自身能量最小。将优化出的构型参数,例如晶格参数,原子位置,键长键角等和实验值进行比对。
会当凌绝顶,一览众山小。
12楼2012-05-15 19:02:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

liu108

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-15 11:55:15
无论是参杂之前还是参杂后,对体系的结构优化是必须的。优化后使得结构能量收敛的最小值,得到能量稳定的晶型,cutoff和k点是用来判断收敛精度的,这两个值取得大,收敛精度也就越大,得到的结构也就越近于稳定结构。
会当凌绝顶,一览众山小。
2楼2012-05-15 09:06:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qx0813

木虫 (正式写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
fzx2008: 金币+1, 鼓励交流 2012-05-15 11:55:27
结构优化是为了获得稳定结构,结构优化的时候也可以计算性质的。
遇到最好的自己
3楼2012-05-15 09:22:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yqkwww

金虫 (正式写手)

引用回帖:
2楼: Originally posted by liu108 at 2012-05-15 09:06:54:
无论是参杂之前还是参杂后,对体系的结构优化是必须的。优化后使得结构能量收敛的最小值,得到能量稳定的晶型,cutoff和k点是用来判断收敛精度的,这两个值取得大,收敛精度也就越大,得到的结构也就越近于稳定结 ...

那请问下面这两张图的优化出来的cutoff值和k点是不是做结构优化出来的??





change myself if I am stocked!
4楼2012-05-15 09:26:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)

【答案】应助回帖

★ ★
fzx2008: 金币+2, 谢谢指教 2012-05-15 11:56:05
结构优化时k点没有必要取得很大,4*4*4就行,虽然说cutoff取得越大越好,但是取得太大收敛的速度就会慢下来,对月能量的截断一般根据你的体系中的元素来定,一般取能量截断值最大的1.3倍即可,从图上看你的能量截断取曲线中平滑的部分的拐点处即可。
会当凌绝顶,一览众山小。
5楼2012-05-15 10:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yqkwww

金虫 (正式写手)

引用回帖:
5楼: Originally posted by liu108 at 2012-05-15 10:38:22:
结构优化时k点没有必要取得很大,4*4*4就行,虽然说cutoff取得越大越好,但是取得太大收敛的速度就会慢下来,对月能量的截断一般根据你的体系中的元素来定,一般取能量截断值最大的1.3倍即可,从图上看你的能量截 ...

这个图的total energy和k-point、cutoff的值的关系图是怎么得到的?尤其是total energy的值一直不懂怎么得到!
change myself if I am stocked!
6楼2012-05-15 14:10:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2012-05-15 14:56:29
这个问题简单啊,你设置不同的k点和cutoff,在输出文件中就会有不同的total energy与之对已,将这一系列的值用origin作图就可以了。
会当凌绝顶,一览众山小。
7楼2012-05-15 14:37:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yqkwww

金虫 (正式写手)

引用回帖:
7楼: Originally posted by liu108 at 2012-05-15 14:37:59:
这个问题简单啊,你设置不同的k点和cutoff,在输出文件中就会有不同的total energy与之对已,将这一系列的值用origin作图就可以了。

我就是没有找到啊····所以才有此问,能发个截图过来看看在哪找的么?
change myself if I am stocked!
8楼2012-05-15 14:39:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2012-05-15 14:56:40
k点和cutoff是要你自己设置的,total energy 是castep输出文件给的。


会当凌绝顶,一览众山小。
9楼2012-05-15 14:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yqkwww

金虫 (正式写手)

引用回帖:
9楼: Originally posted by liu108 at 2012-05-15 14:49:55:
k点和cutoff是要你自己设置的,total energy 是castep输出文件给的。
af/9b/1210273_1337064591_360.jpg

很是感谢!我知道K和cutoff值是自己设的,但是在不知道的情况和文献中也没有的情况下,结构优化时是不是需要把这两个参数像上面一样,优化出来呢?结构优化我一直不明白优化参数到底是优化什么参数??十分感谢····
change myself if I am stocked!
10楼2012-05-15 16:09:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yqkwww 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:34 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:23 by Xc2FxHsHkBNc
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:06 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见