| 查看: 1078 | 回复: 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人回复
★ ★ ★ ★ ★ ★ ★ ★
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
★ ★ ★ ★ ★ ★ ★ ★ ★
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