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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » CASTEP几何优化问题
101/1返回列表
查看: 2990  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

nan1026feng

木虫 (小有名气)

[求助] CASTEP几何优化问题

请问下，我在做bulk Fe 几何优化时，始终不能得到与实验值相差不大的估算值。我改变了多次参数，还是不行啊，而且当晶格参数与实验值较为接近时，体弹性模量却相差很大，截断能为260ev，k点网格为2,2,1时，晶格参数为2.858327,（实验值为2.86埃），但体弹性模量估算为500GPa,实验值却为174GPa，神奇的是，无论我怎么改变k点网格和截断能，变的是晶格参数，体弹性模量始终没变啊，咋个回事？另外顺带问下，优化后的体弹性模量，是不是看*.castep中前面的Geometry Optimization Parameters下的estimated bulk modulus                         :  500.0       GPa这一项啊？我刚刚开始学castep，很多都不懂，请教哈大家啊，希望大家能帮帮我啊
JTW8@NOORVF)CV7`T73CDF.jpg[/filename]




[ Last edited by nan1026feng on 2011-12-26 at 21:56 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-12-26 21:49:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
WDD880227(金币+1): 感谢交流~~~有帮助~~ 2011-12-27 09:22:54
nan1026feng(金币+2): ★有帮助 谢谢你对这个问题的关注，感谢啊 2011-12-28 11:16:11
金属？

Cutoff， K-point 你都取的太低了。。  精度问题
赞一下(1人) 回复此楼
遇弱则弱
2楼2011-12-26 21:59:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444239852

银虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
WDD880227: 2011-12-27 09:23:05
nan1026feng(金币+2): ★有帮助 谢谢，我试了，不过感觉截断能和k-point设不了多大样啊，我设置大了后，与实验值差的更大了啊 2011-12-28 11:17:51
感觉也跟K点的网格选取有关，如果选取的越高越接近精确
赞一下(1人) 回复此楼
3楼2011-12-27 09:11:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nan1026feng

木虫 (小有名气)
这个低了啊？要多高才高啊
赞一下 回复此楼
4楼2011-12-27 09:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nan1026feng

木虫 (小有名气)
引用回帖:
楼: Originally posted by franch at 2011-12-26 21:59:54:
金属？

Cutoff， K-point 你都取的太低了。。  精度问题

到底是越大越高吗？还是越小越精细啊？我刚学，不太明白啊
赞一下 回复此楼
5楼2011-12-27 09:38:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nan1026feng

木虫 (小有名气)
引用回帖:
楼: Originally posted by 444239852 at 2011-12-27 09:11:17:
感觉也跟K点的网格选取有关，如果选取的越高越接近精确

哦，谢谢你，要的，那我再改改试试看
赞一下 回复此楼
6楼2011-12-27 09:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nan1026feng

木虫 (小有名气)
引用回帖:
2楼: Originally posted by franch at 2011-12-26 21:59:54:
金属？

Cutoff， K-point 你都取的太低了。。  精度问题

我改了cutoff，k-point啊，我现在把cutoff取为450，k-point取为15*15*15，但是算出来，晶格参数2.747（实验值2.86），体弹性模量也不对，算出来还是500GPa（实验值174GPa），到底咋个回事啊，望解救啊，谢谢你啦
赞一下 回复此楼
7楼2011-12-27 10:07:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nan1026feng

木虫 (小有名气)
引用回帖:
2楼: Originally posted by franch at 2011-12-26 21:59:54:
金属？

Cutoff， K-point 你都取的太低了。。  精度问题

大侠，算出来有很多体弹性模量，到底看哪个呢，是下面这个吗？
BFGS: Final Enthalpy     = -8.64900225E+002 eV
BFGS: Final bulk modulus =    1785.41620 GPa，我这算出来咋个那么大啊，实验值才174GPa，这差距也太大了，咋个回事哦
赞一下 回复此楼
8楼2011-12-27 10:25:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaq2800

铁杆木虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
cutoff小了点，体弹性模量要另外单算的把，不是这个选项么？
赞一下 回复此楼
9楼2011-12-27 14:37:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nan1026feng

木虫 (小有名气)
引用回帖:
9楼: Originally posted by aaq2800 at 2011-12-27 14:37:58:
cutoff小了点，体弹性模量要另外单算的把，不是这个选项么？

是吗？体弹性模量要单独重新算啊？cutoff还小了吗？那到底如何判断我所设置的参数正确呢？难道就只是看优化后的晶格参数和实验值接不接近吗？我在截断能取260，k-point取2*2*1的时候算出来，晶格参数为2.858327啊，跟实验值差不多，那难道就能说明我设置的参数是正确的吗？
赞一下 回复此楼
10楼2011-12-27 18:39:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nan1026feng 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-31 02:45 by hZiFeudZyoGR
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-31 02:35 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-31 02:13 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-31 01:38 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:37 by hZiFeudZyoGR
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:07 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 12/600 2026-08-31 00:34 by hZiFeudZyoGR
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:21 by hZiFeudZyoGR
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 AVl0JAsg9hyo 2026-08-29 5/250 2026-08-30 19:57 by hZiFeudZyoGR
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:36 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 14:03 by GkCVyODlx4mR
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭