24小时热门版块排行榜    

查看: 2271  |  回复: 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 的主题更新
信息提示
请填处理意见