| 查看: 1200 | 回复: 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? |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
小木虫看见有人已经知道结果了
已经有22人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有11人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有21人回复
-
2026年国自然面上资助率
已经有19人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
★ ★ ★ ★ ★ ★ ★ ★ ★
aylayl08(金币+8,VIP+0):谢谢详细解答 1-4 10:48
后天一(金币+1,VIP+0):谢谢解答，不过这个真的是结构优化后的结果！ 1-5 13:41
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
happyjwx
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3099.5
- 散金: 70
- 帖子: 178
- 在线: 176.9小时
- 虫号: 934744
- 注册: 2009-12-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2010-01-04 10:44:53
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 3640.6
- 散金: 2141
- 红花: 37
- 帖子: 1432
- 在线: 681.5小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
4楼2010-01-04 14:53:49
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 3640.6
- 散金: 2141
- 红花: 37
- 帖子: 1432
- 在线: 681.5小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
5楼2010-01-04 14:58:53
后天一
木虫 (小有名气)
我们的自由！
- 1ST强帖: 3
- 应助: 9 (幼儿园)
- 金币: 2544.6
- 红花: 4
- 帖子: 204
- 在线: 311小时
- 虫号: 519608
- 注册: 2008-03-06
- 专业: 凝聚态物性 II ：电子结构
6楼2010-01-05 13:44:43
7楼2010-01-05 19:23:34
后天一
木虫 (小有名气)
我们的自由！
- 1ST强帖: 3
- 应助: 9 (幼儿园)
- 金币: 2544.6
- 红花: 4
- 帖子: 204
- 在线: 311小时
- 虫号: 519608
- 注册: 2008-03-06
- 专业: 凝聚态物性 II ：电子结构
8楼2010-01-05 23:30:16
★ ★ ★ ★ ★ ★ ★ ★
aylayl08(金币+8,VIP+0):感谢再次指点 1-7 15:36
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拟合还是很准的。 |
9楼2010-01-07 12:37:06
10楼2010-01-08 02:17:13
15