24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 1079  |  回复: 12
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

后天一

木虫 (小有名气)

我们的自由!

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

最近优化结构发现,能量最低对于的结构其应力有可能很大
比如
  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?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)

★ ★
wuchenwf(金币+1,VIP+0):谢谢 1-9 22:01
后天一(金币+1,VIP+0):谢谢指点 1-10 11:21
我实在看不下去了,你不要说人家概念不清楚,这是pullay stress。。。
引用回帖:
Originally posted by xbelinda at 2010-1-5 19:23:
晕。
murn拟合的话怎么会有stress。
肯定是你的参数设置有问题。
你可以把INCAR和KPOINTS粘一下。
还有,你最需要补一下基本概念。目前最关键的问题,是你概念不清!

http://sim.phys.ruc.edu.cn
10楼2010-01-08 02:17:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

xbelinda

银虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★
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大还可能有问题,如果一样的话说明是受的净水压,应该没问题。
2楼2010-01-03 23:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyjwx

木虫 (小有名气)

★ ★
aylayl08(金币+1,VIP+0):谢谢提示 1-4 10:48
后天一(金币+1,VIP+0):谢谢! 1-5 13:41
楼上讲得很好,更正一点,应该是“静水压”不是“净水压”
夜半追忆往昔
3楼2010-01-04 10:44:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

★ ★ ★ ★ ★ ★
aylayl08(金币+5,VIP+0):谢谢解答 1-4 18:49
后天一(金币+1,VIP+0):谢谢您的意见,我试试看! 1-5 13:42
22.90kB(2.29GPa)对于有些问题来说可能并不会影响太大,本身这种量级的压力在理论模拟中还算是比较小的(和理论误差相比)。所以楼主得看具体问题是什么。
还有,variable cell可能是个好办法,但是由于这种方法大大增加了优化的自由度,所以可能最终收敛的结果和我们所期望的会有些差距。在体系结构比较了解的情况下,建议手动调节原胞体积来消除压力,而且这样计算量也会小一些。
先进使役材料计算(https://aidme.nimte.ac.cn)
4楼2010-01-04 14:53:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见