| 查看: 1072 | 回复: 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? |
回复此楼» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
★ ★
小木虫(金币+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 (幼儿园)
- 金币: 3094.5
- 散金: 70
- 帖子: 168
- 在线: 174.5小时
- 虫号: 934744
- 注册: 2009-12-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2010-01-04 10:44:53
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 4306.7
- 散金: 1541
- 红花: 36
- 帖子: 1403
- 在线: 676.6小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
4楼2010-01-04 14:53:49