应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » castep晶格优化问题
161/2返回列表12下一页
查看: 3898  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

菲顶顶

铜虫 (正式写手)

[求助] castep晶格优化问题

我想问一下用castep做晶格优化时我选了task里的geometry optimization，还勾选了Optimize cell，band structure，density of states和里面的calculate PDOS，然后点run。这个步骤应该是正确的吧，
问题一：我用的是ICSD里查出的数据，必须要进行结构优化才可以进行以后的分析吗？
问题二：查看优化以后的结果是在哪里看啊，是.castep文件吗？里面有好多东西，晶格常数那些的出现了好几次，是看最后那一部分的吗？
问题三：我看到有人说要多次优化，两次优化的什么能量差<0.001收敛才算有话成功，这个是什么意思啊 ，我要改变参数并重复上面的步骤两次吗？
问题四：我还看到有人说优化结构时改变截断能，在不同的截断能下优化得到不同的结构参数，选择合适的什么的，是什么意思啊？
望请各位高手多多帮助讲解，先在这里谢过了。
回复此楼

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

实验技术 ms中castep模块疑难杂症

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-05 16:42:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

书万里

木虫 (著名写手)

历史控+原著控

【答案】应助回帖

★ ★ ★ ★
liliangfang: 金币+2, 谢谢指教 2012-09-10 08:26:31
菲顶顶: 金币+2 2012-09-10 17:10:38
引用回帖:
5楼: Originally posted by 菲顶顶 at 2012-09-08 10:55:31
对于第三个问题的回答我还有以下问题：
1，是在哪里设定能量差啊，是软件自动设定的吗？找不到；
2，还要计算力和位移吗，那个也需要设定吗？在哪里？有什么标准吗？ ；
3，看下面那个图，左边的paramenter都指 ...

你在task中选择optimization后，旁边的more按钮会成为活动的，点击后里面就有设置四个收敛标准的。你图中的parameter四个选项就是四个收敛标准(分别是能量差、力、位移和压力），tolerance就是你设定的值，value是当前计算的值，只有当value小于tolerance时才会出现yes（你的图中出现是个yes，表明你的计算已收敛）。不同体系采用的收敛标准肯定不同，具体要自己测试或查阅文献。
    收敛测试应在计算之前，他是为了保证你的计算准确。
赞一下(4人) 回复此楼

» 本帖已获得的红花（最新10朵）

菲顶顶
喂人民服雾。
6楼2012-09-09 02:07:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

书万里

木虫 (著名写手)

历史控+原著控

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xueht987: 金币+3, 鼓励交流！ 2012-09-05 20:29:49
菲顶顶: 金币+5 2012-09-08 09:51:43
liliangfang: 金币+2, 谢谢指教 2012-09-08 10:31:54
我也是新手，来尝试回答楼主的问题。
1.还是必须优化下，结构优化时一切性质计算的基础。
2.是的.castep文件中在optimization completely后面的数据就是优化最终的数据。
3.能量差是收敛标准，比如你设定为0.001，如果优化迭代过程中两次得到的能量相差小于0.001，就表示能量差标准达到，就会产生一个YES（见.castep文件）。一般计算还有另外两个标准力和位移。
4.他们说的应该是收敛测试。平面波截断能和k点是两个重要的参数，关系到你计算的精度和时间。适合你体系的截断能和k点就需要你自己先测试。一般这样，固定一个大的截断能，改变k点，观察总能的变化；再固定一个大的k点，改变截断能，观察总能变化。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

菲顶顶
喂人民服雾。
2楼2012-09-05 18:28:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
送鲜花一朵
引用回帖:
6楼: Originally posted by 书万里 at 2012-09-09 02:07:06
你在task中选择optimization后，旁边的more按钮会成为活动的，点击后里面就有设置四个收敛标准的。你图中的parameter四个选项就是四个收敛标准(分别是能量差、力、位移和压力），tolerance就是你设定的值，value是 ...

非常感谢 讲解的很清楚
赞一下(1人) 回复此楼
7楼2012-09-09 09:29:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

alienlj

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-09-06 11:52:29
菲顶顶: 金币+3 2012-09-08 09:52:02
刚刚接触ms 尝试回答下
1.文献里面都是优化的，所以咱们也优化下吧 反正也不麻烦
2.对 final figuration里面的就是最终结果了
3.我看过是改变截断能 然后对比不同截断能下的总能 差距小于一定值就可以认为收敛了
4.这个问题我也在纠结 一般都是多做几组 然后对比找出适合的吧
赞一下 回复此楼
3楼2012-09-06 09:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
菲顶顶: 回帖置顶 2012-09-09 09:30:00
第一个和第二个问题明白了 我再问一下  我上面那个优化的步骤是对的吧   为什么优化后我得到的晶格常数和文献里的相差很大呢？ 而且结构也变了  在文献里查到的是a=b=6，不等于c=4点多，且角度为90，优化后变成a=b=c=4点多了，角度变成104了，这是怎么回事啊
赞一下 回复此楼
4楼2012-09-08 10:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by 书万里 at 2012-09-05 18:28:46
我也是新手，来尝试回答楼主的问题。
1.还是必须优化下，结构优化时一切性质计算的基础。
2.是的.castep文件中在optimization completely后面的数据就是优化最终的数据。
3.能量差是收敛标准，比如你设定为0.001 ...

对于第三个问题的回答我还有以下问题：
1，是在哪里设定能量差啊，是软件自动设定的吗？找不到；
2，还要计算力和位移吗，那个也需要设定吗？在哪里？有什么标准吗？ ；
3，看下面那个图，左边的paramenter都指的是什么参数啊，看不太懂。
对于第四个问题：
收敛测试是在什么时候进行啊，是在优化以前还是优化完结构以后啊。

CZ9[Z7AR)BVQJ880CQG77)9.jpg
赞一下 回复此楼
5楼2012-09-08 10:55:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
4楼: Originally posted by 菲顶顶 at 2012-09-08 10:28:36
第一个和第二个问题明白了 我再问一下  我上面那个优化的步骤是对的吧   为什么优化后我得到的晶格常数和文献里的相差很大呢？ 而且结构也变了  在文献里查到的是a=b=6，不等于c=4点多，且角度为90，优化后变成a=b= ...

大家帮忙解决一下这个问题啊  没有人回答啊
赞一下 回复此楼
8楼2012-09-09 20:18:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
送鲜花一朵
引用回帖:
2楼: Originally posted by 书万里 at 2012-09-05 18:28:46
我也是新手，来尝试回答楼主的问题。
1.还是必须优化下，结构优化时一切性质计算的基础。
2.是的.castep文件中在optimization completely后面的数据就是优化最终的数据。
3.能量差是收敛标准，比如你设定为0.001 ...

我有个问题哈 就是第四个 先取个大的截断能 然后改变k点，这样会得到一个我们认为合适的k点，，然后再选个大的k点 让偶改变截断能，这样我们得到一个我们认为合适的截断能。然后计算的时候用这两个合适的截断能和k点吗？，总觉得这样不太合适
赞一下 回复此楼
9楼2012-09-13 15:57:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

书万里

木虫 (著名写手)

历史控+原著控
引用回帖:
9楼: Originally posted by 菲顶顶 at 2012-09-13 15:57:08
我有个问题哈 就是第四个 先取个大的截断能 然后改变k点，这样会得到一个我们认为合适的k点，，然后再选个大的k点 让偶改变截断能，这样我们得到一个我们认为合适的截断能。然后计算的时候用这两个合适的截断能和k ...

是的。什么地方不合适了？你的体系如果有人做过，你也可以借用别人的。
赞一下 回复此楼
喂人民服雾。
10楼2012-09-13 20:33:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 菲顶顶 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后，还能找到好工作吗？ +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭