应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【讨论】结构优化能量最低但存在较大的应力
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
131/2返回列表12下一页
查看: 1073  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

后天一

木虫 (小有名气)

我们的自由！

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

最近优化结构发现，能量最低对于的结构其应力有可能很大
比如
  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?
回复此楼

» 猜你喜欢
1楼2010-01-01 10:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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大还可能有问题，如果一样的话说明是受的净水压，应该没问题。
赞一下(11人) 回复此楼
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
楼上讲得很好，更正一点，应该是“静水压”不是“净水压”
赞一下(11人) 回复此楼
夜半追忆往昔
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可能是个好办法，但是由于这种方法大大增加了优化的自由度，所以可能最终收敛的结果和我们所期望的会有些差距。在体系结构比较了解的情况下，建议手动调节原胞体积来消除压力，而且这样计算量也会小一些。
赞一下(11人) 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
4楼2010-01-04 14:53:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)
★ ★
aylayl08(金币+2,VIP+0):谢谢解答 1-4 18:49
附加一点：根据xbelinda的分析，楼主这个例子应该是静水压情况。所以个人认为手动调节体系原胞体积会比较可取一些。
赞一下(10人) 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
5楼2010-01-04 14:58:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

后天一

木虫 (小有名气)

我们的自由！
谢谢各位大佬，貌似我说的不是很明确。
这样的无论我采用ISIF=3或者是murn拟合，最后的结果stress都不为零！是不是计算精度的问题？
赞一下 回复此楼
6楼2010-01-05 13:44:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)
★ ★ ★
后天一(金币+1,VIP+0):thanks! 1-5 23:30
ice_rain(金币+1,VIP+0):谢谢交流 1-6 11:53
aylayl08(金币+1,VIP+0):谢谢交流讨论 1-7 15:36
晕。
murn拟合的话怎么会有stress。
肯定是你的参数设置有问题。
你可以把INCAR和KPOINTS粘一下。
还有，你最需要补一下基本概念。目前最关键的问题，是你概念不清！
赞一下(16人) 回复此楼
7楼2010-01-05 19:23:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

后天一

木虫 (小有名气)

我们的自由！
★
ice_rain(金币+1,VIP+0):讨论奖励 1-6 11:54
引用回帖:
Originally posted by xbelinda at 2010-1-5 19:23:
晕。
murn拟合的话怎么会有stress。
肯定是你的参数设置有问题。
你可以把INCAR和KPOINTS粘一下。
还有，你最需要补一下基本概念。目前最关键的问题，是你概念不清！

呵呵，是这样的，我用ISIF=3 relax后的CONTCAR再静态计算！
或者手动改变晶格常数，然后murn拟合后的晶格参数，再静态计算！

然后参考静态计算的stress。

不知道这样行不?
赞一下(6人) 回复此楼
8楼2010-01-05 23:30:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★
aylayl08(金币+8,VIP+0):感谢再次指点 1-7 15:36
用isif=3 relax后的CONTCAR 静态计算后，你要看stress是否相等，如果相等后就可以用murn去拟合。因为实际拟合需要的是 E和V。
之所以要用isif=3  relax, 是因为如果不是立方结构的话， 晶格abc的比例在改变体积后可能会变化。 如果是立方体系，则不用relax，那样很慢。
你用murn拟合后，得到的0压晶格常数（即平衡体积），再做静态计算，有可能三轴stress不为0，这就是为什么要用murn去拟合，而不直接用relax的0压晶格常数。
这个问题主要是因为stress算的不准。毕竟能量如果说是2阶错误的话，那么力就是1阶错误。  所以要想得到精确的受力还需要特别精确的能力。
目前dft的一个问题，就是压力算的不是很准。但随着压力提高， 有所好转。   
所以对于一般的结构，murn拟合还是很准的。
赞一下(10人) 回复此楼
9楼2010-01-07 12:37:06
 已阅   回复此楼   关注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粘一下。
还有，你最需要补一下基本概念。目前最关键的问题，是你概念不清！

赞一下(10人) 回复此楼
http://sim.phys.ruc.edu.cn
10楼2010-01-08 02:17:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 后天一 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭