24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

后天一

木虫 (小有名气)

我们的自由！

1ST强帖: 3
应助: 9 (幼儿园)
金币: 2544.6
红花: 4
帖子: 204
在线: 311小时
虫号: 519608
注册: 2008-03-06
专业: 凝聚态物性 II ：电子结构

[交流] 【讨论】结构优化能量最低但存在较大的应力

最近优化结构发现，能量最低对于的结构其应力有可能很大
比如
  FORCE on cell =-STRESS in cart. coord.  units (eV/reduce length):
  Direction X       Y       Z       XY    YZ    ZX
  --------------------------------------------------------------------------------------
  Alpha Z 1791.70 1791.70 1791.70
  Ewald -7349.25  -7349.29  -7349.25    0.00    0.00    0.01
  Hartree 1670.82 1670.80 1670.82    0.00    0.00    0.00
  E(xc) -1457.85  -1457.85  -1457.85    0.00    0.00    0.00
  Local    665.82 665.87 665.81    0.00    0.00    -0.02
  n-local -616.31 -616.31 -616.31    0.00    0.00    0.00
  augment 745.38 745.38 745.38    0.00    0.00    0.00
  Kinetic 4544.65 4544.65 4544.65    0.00    0.00    0.00
  -------------------------------------------------------------------------------------
  Total    -5.05    -5.05    -5.05    0.00    0.00    0.00
  in kB    -22.90 -22.91 -22.90    0.00    0.01    0.00
  external pressure =    -22.90 kB  Pullay stress =       0.00 kB

我想知道如何能使stress变为0?

回复此楼

» 猜你喜欢

26年申博自荐-计算机视觉已经有4人回复
考博已经有4人回复
药化及相关博士的申请已经有3人回复
一篇MDPI论文改变了学习工作和生活已经有4人回复
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇已经有4人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有3人回复
收到国自然专家邀请后几年才会有本子送过来评已经有3人回复

1楼 2010-01-01 10:38:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 3590.1
散金: 2141
红花: 37
帖子: 1432
在线: 681.4小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

★ ★ ★ ★ ★ ★
aylayl08(金币+5,VIP+0):谢谢解答 1-4 18:49
后天一(金币+1,VIP+0):谢谢您的意见，我试试看！ 1-5 13:42

22.90kB(2.29GPa)对于有些问题来说可能并不会影响太大，本身这种量级的压力在理论模拟中还算是比较小的（和理论误差相比）。所以楼主得看具体问题是什么。
还有，variable cell可能是个好办法，但是由于这种方法大大增加了优化的自由度，所以可能最终收敛的结果和我们所期望的会有些差距。在体系结构比较了解的情况下，建议手动调节原胞体积来消除压力，而且这样计算量也会小一些。

赞一下(11人)

回复此楼

先进使役材料计算（https://aidme.nimte.ac.cn）

4楼2010-01-04 14:53:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 13 个回答

xbelinda

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 391.7
红花: 1
帖子: 160
在线: 6.4小时
虫号: 587329
注册: 2008-08-22
性别: MM
专业: 无机合成和制备化学

★ ★ ★ ★ ★ ★ ★ ★ ★
aylayl08(金币+8,VIP+0):谢谢详细解答 1-4 10:48
后天一(金币+1,VIP+0):谢谢解答，不过这个真的是结构优化后的结果！ 1-5 13:41

你应该做的是scf吧！
你如果想把stress变成0.实际就是在做优化。
让结构relax到0GPa。
你需要在INCAR中加入：
NSW=100 #其实我觉得大于50就可以了。
EDIFFG=1E-5 # 通常是EDIFF的值的10倍。
ISIF=3 #改变胞的体积和形状。
IBRION=2 #这个可以看看manual，1，2，3 应该都可以。选择哪个取决于自己。
PSTRESS=0 #单位是0Kbar， 10Kbar=1GPa。这个是优化的最终压力。负的stress=压力。
这些参数应该差不多了。

还有，“最近优化结构发现，能量最低对于的结构其应力有可能很大” 你这句话说的其实有问题。第一，你根本就没做结构优化。第二，如果是哪一个轴的stress大还可能有问题，如果一样的话说明是受的净水压，应该没问题。

赞一下(11人)

回复此楼

2楼2010-01-03 23:06:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖