24小时热门版块排行榜    

查看: 2756  |  回复: 7

enola

捐助贵宾 (正式写手)

[求助] 晶格常数优化

晶格常数优化:
按照候老师的步骤,我对一个cubic的晶体的晶格常数进行优化:
1等间距增加的取10个晶格常量,计算能量;
2计算完成后,提取能量和体积参数;
3用Birch-Murnaghan拟合,得到晶格参数,
4用获得的晶格常数再计算一次能量,查看OUTCAR里的,发现
  external pressure =      -47.30 kB  Pullay stress =        0.00 kB
压力依然比较大。
所以想请教下大家,可可还需要在微调,还是这样的压力是可以接受的,谢谢。。
回复此楼
timeflies..
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


uuv2010(金币+1): 多谢提示 2011-09-12 20:34:53
enola(金币+10): 谢谢。 2011-09-13 08:52:48
如果你是晶体压力是肯定有的,一般情况下只看原子的受力吧,我是这样处理的,对于压力那一块还不懂
嘟啦啦嘟啦啦嘟
2楼2011-09-12 15:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

引用回帖:
2楼: Originally posted by 贺仪 at 2011-09-12 15:52:20:
如果你是晶体压力是肯定有的,一般情况下只看原子的受力吧,我是这样处理的,对于压力那一块还不懂

恩,是会有压力。我的情况是,我用一个原子半径比较大的原子替换了原先半径比较小的原子,所以引入了自身的压力,我想通过优化结构,尽量减少这样的压力。
timeflies..
3楼2011-09-13 08:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)

本帖内容被屏蔽

4楼2011-09-13 09:03:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
zzy870720z(金币+2): 谢谢提示 2011-09-13 15:50:02
引用回帖:
4楼: Originally posted by gougou531 at 2011-09-13 09:03:51:
贺仪版主,哪看原子受力的还,一般什么要求啊?能否截个图示范一下,还有看原子受力的话,那计算的时候参数怎么设置?谢谢

受力查看
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      2.84235      2.44083      5.39857         0.000000      0.000000      0.000000


EDIFFG控制原子受力
嘟啦啦嘟啦啦嘟
5楼2011-09-13 10:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by enola at 2011-09-13 08:51:08:
恩,是会有压力。我的情况是,我用一个原子半径比较大的原子替换了原先半径比较小的原子,所以引入了自身的压力,我想通过优化结构,尽量减少这样的压力。

期待牛人来解答!
嘟啦啦嘟啦啦嘟
6楼2011-09-13 10:34:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)

本帖内容被屏蔽

7楼2011-09-13 14:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

★ ★
zzy870720z(金币+2): 鼓励讨论 2011-09-13 15:50:23
引用回帖:
1楼: Originally posted by enola at 2011-09-12 15:30:09:
晶格常数优化:
按照候老师的步骤,我对一个cubic的晶体的晶格常数进行优化:
1等间距增加的取10个晶格常量,计算能量;
2计算完成后,提取能量和体积参数;
3用Birch-Murnaghan拟合,得到晶格参数,
4用获得 ...

再补充下问题:
我用ISIF=3作结构优化,发现很矛盾的结果:
用Birch-Murnaghan拟合,得到晶格参数是比初始的大(我觉得合理,因为我替换的元素半径比原先的大)
ISIF=3的结果是晶格常数比初始的小。
也许是我的计算精度不够,但是大家有什么其他好的建议吗?
还有就是关于
external pressure =      -47.30 kB  Pullay stress =        0.00 kB
这个量,其实于晶格常量,以及内部原子的位置有关。。。不知道其表达的是是什么物理含义,等我查阅,也欢迎大家讨论阿。
timeflies..
8楼2011-09-13 15:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 enola 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 跳槽后在研项目怎么办? +4 简单化xn 2026-08-22 8/400 2026-08-23 07:20 by ilovexiaobin
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见