| 查看: 1263 | 回复: 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? |
回复此楼» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有14人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ice_rain(金币+1,VIP+0):thanks very much 1-13 17:32
小木虫(金币+0.5):给个红包，谢谢回帖交流
ice_rain(金币+1,VIP+0):thanks very much 1-13 17:32
|
随便抄一段：The Pulay stress arises from the fact that the plane wave basis set is not complete with respect to changes of the volume. Thus, unless absolute convergence with respect to the basis set has been achieved - the diagonal components of the stress tensor are incorrect. This error is often called ``Pulay stress''. The error is almost isotropic (i.e. the same for each diagonal component), and for a finite basis set it tends to decrease volume compared to fully converged calculations (or calculations with a constant energy cutoff). 人算出来所谓的external pressure 是因为Pulay Stress。我没时间去分辨文献对错，只对基本概念有点了解。 [ Last edited by wei.ji on 2010-1-13 at 00:03 ] |
13楼2010-01-12 23:48:47
★ ★ ★ ★ ★ ★ ★ ★ ★
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 (小学生)
- 金币: 3651.6
- 散金: 2321
- 红花: 37
- 帖子: 1432
- 在线: 681.5小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
4楼2010-01-04 14:53:49