24小时热门版块排行榜    

查看: 1967  |  回复: 9
【悬赏金币】回答本帖问题,作者nan1026feng将赠送您 16 个金币

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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 你都取的太低了。。  精度问题
遇弱则弱
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点的网格选取有关,如果选取的越高越接近精确
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 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 和导师撕破脸了 +17 圆满123456 2024-11-13 21/1050 2024-11-16 00:30 by ZZZemmm
[论文投稿] 国内期刊审稿人数量 +3 新时代核动力驴 2024-11-13 5/250 2024-11-15 20:37 by 新时代核动力驴
[论文投稿] 找不到审稿人 +6 吉祥如意141913 2024-11-13 6/300 2024-11-15 16:48 by LEyblyt
[基金申请] 博后面上太难了 +14 不安定因素 2024-11-14 17/850 2024-11-15 16:05 by lupbao
[教师之家] 大学老师 +11 考研一路顺风 2024-11-13 14/700 2024-11-15 15:24 by LNP@mRNA
[基金申请] 博后资助名单出来了 +12 Shxjjxjkx 2024-11-14 14/700 2024-11-15 14:40 by fanfenggui
[硕博家园] 实验进度 +9 麻辣烫。 2024-11-11 14/700 2024-11-15 13:14 by hanweiguo
[有机交流] 二甲胺的使用 20+5 太阳谷 2024-11-14 13/650 2024-11-15 13:11 by xiaomei1031
[基金申请] 求助 +4 Enenenene 2024-11-15 4/200 2024-11-15 09:36 by 榨菜拌青椒
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+3 猪小耍 2024-11-13 12/600 2024-11-15 08:41 by 北京莱茵润色
[基金申请] 博后基金分组排名 +7 攻城2025 2024-11-14 7/350 2024-11-14 21:19 by 实验小白ha
[论文投稿] 投稿系统中的通讯作者和文章中的通讯作者不一样,文章目前被录用了? +4 babybabygo 2024-11-12 5/250 2024-11-14 19:13 by 走了002
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[有机交流] 硫谜用双氧水氧化出现的困惑(老板说一定要用双氧水氧化做好)求指教 10+4 luowenwu 2024-11-10 15/750 2024-11-14 10:52 by luolinfeng
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[论文投稿] 论文接受后,编辑部发邮件提醒缺少附件 15+4 上善若水明泪 2024-11-10 4/200 2024-11-12 14:24 by 北京莱茵润色
[论文投稿] 爱思唯尔投稿系统里的通讯作者可以和文章里的通讯作者标注不同吗 +7 Omnissiah 2024-11-10 7/350 2024-11-12 14:07 by holypower
[有机交流] 求助NMR +5 苯巴比妥! 2024-11-09 5/250 2024-11-12 11:20 by 88817753
信息提示
请填处理意见