24小时热门版块排行榜    

查看: 4438  |  回复: 7

mazuju028

木虫 (正式写手)

[求助] 单轴施加应变结构优化的方法有哪些?

方法1,使用vasp. 修改constr_cell_relax.F,重新编译vasp,在计算中,我使用基态构型,使用ISIF=3,固定 c 轴,结果a b轴在优化后偏离实验值5.8%, 与不固定c轴,同样的优化,a,b,c偏离实验值4.0%相比,这样的结果合理吗?
方法2,使用abinit,设置参数optcell=9,结果发现 c 轴和 a ,b 轴都在变,怀疑是不是保持对称性的原因,有没有人知道怎么在abinit不保持对称性,有类似vasp中ISYM=0的参数吗?
方法3,使用vasp, 手动调节a,b,c,然后固定abc, 进行内坐标优化(ISIF=2),这种方法太繁琐了吧?
PS: 还有没有人正使用其它更有效的方法呢?

[ Last edited by mazuju028 on 2011-12-8 at 11:34 ]
回复此楼

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

第一性原理相关文档 first principle VASP

» 猜你喜欢

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

锐气藏于胸,和气浮于脸,才气见于事,义气施与人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

mazuju028(金币+2): 谢谢交流 2011-12-09 08:03:49
感觉用第3中方法做就可以了吧,虽然麻烦一点,但是感觉比让程序自己跑更放心一点
2楼2011-12-08 20:01:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
2楼: Originally posted by uuv2010 at 2011-12-08 20:01:04:
感觉用第3中方法做就可以了吧,虽然麻烦一点,但是感觉比让程序自己跑更放心一点

如果是立方或是四方晶系还好调整,如果是正交或是单斜,三边都不等,那调整就有点麻烦了吧。
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
3楼2011-12-09 08:12:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

引用回帖:
3楼: Originally posted by mazuju028 at 2011-12-09 08:12:10:
如果是立方或是四方晶系还好调整,如果是正交或是单斜,三边都不等,那调整就有点麻烦了吧。

好像有一种说法是:计算的晶格常数和实验的误差应该在3%以内。在这个范围内扫描的话计算量会减少很多的
4楼2011-12-09 09:49:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★
mazuju028(金币+3): 谢谢指点。我明白了。 2011-12-09 21:05:50
sunyang1988(金币+3): 谢谢侯老师指点 2011-12-09 21:49:57
直接与实验值进行比较的这种方式,不能说明你固定c轴优化的结果不好。
LDA/GGA在对结构优化时,得到的结构参数一般都与实验值存在一定的偏差。

有意义的比较方式:固定c轴,手动的改变a,b的值,做出E(a,b)的一个三维图,找出这函数的最低点处的a和b值,再把它们与你修改constr_cell_relax.F之后的优化结果相比较。
5楼2011-12-09 18:30:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)

楼主,请问您的这个问题:”单轴施加应变结构优化的方法有哪些?”最后采用了哪种方法?

另外,ABINIT 是如何设置固定一轴而优化其它两轴?能否麻烦您把ABINIT的输入文件都发给我看看,我也遇到了同样的问题。我的邮箱是:hfdong@semi.ac.cn   谢谢您!
好好学习
6楼2012-08-21 09:59:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WXPhoebe

新虫 (初入文坛)

请问您最后采用的是哪种方法?
7楼2019-12-03 21:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)

本帖仅楼主可见
8楼2019-12-04 11:36:41
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 mazuju028 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-31 00:42 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 12/600 2026-08-31 00:34 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 15/750 2026-08-31 00:10 by hZiFeudZyoGR
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:40 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:21 by hZiFeudZyoGR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 jCd0dEvKHShX 2026-08-29 7/350 2026-08-30 19:45 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:23 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 17:47 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:36 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见