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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【讨论】结构优化能量最低但存在较大的应力
北京石油化工学院2026年研究生招生接收调剂公告
131/2返回列表12下一页
查看: 1142  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 一志愿安徽大学085601 +11 zyx上岸！ 2026-04-04 11/550 2026-04-06 22:55 by qlm5820
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 327考研调剂推荐 +6 呜呜呜呜呢 2026-04-06 6/300 2026-04-06 21:39 by 啵啵啵0119
[考研] 0703化学调剂325分 +12 15771691647 2026-04-04 13/650 2026-04-06 12:00 by lijunpoly
[考研] 086000生物与医药298调剂求助 +9 元元青青 2026-03-31 12/600 2026-04-05 21:03 by 学员8dgXkO
[考研] 296求调剂 +3 汪！？！ 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 材料工程085601数二英一335求调剂 +6 双马尾痞老板2 2026-03-31 6/300 2026-04-04 22:29 by hemengdong
[考研] 一志愿南农090401，268，求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 求调剂：085600材料与化工，考材科基，总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 334求调剂 +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by 琢珥丶
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 土木304求调剂 +4 兔突突突， 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突，
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 化学工程专硕324分，一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
信息提示
关闭
请填处理意见
关闭