| 查看: 1705 | 回复: 7 | |||
didi5158金虫 (小有名气)
|
[交流]
【求助】结构优化后，stress仍然很大，如何解决？已有5人参与
|
|
对结构进行优化，relax在11个bfgs后，total force 以及 scf 均达到收敛阈值，但是结构所受的stress 很大，达到p= -3443。 采用新的初始原子位置值也是老样子。改用vc-relax 计算后晶格长度变化很大后收敛。请问大家造成这个结果的原因是什么？如何调整！ 先谢谢了 |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有23人回复
-
导师想让我从独立一作变成了共一第一
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
新建碳纳米管，进行结构优化，为什么五天才算完？
已经有6人回复- 用MS的castep做完几何优化后如何查看优化后的结构
已经有8人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 【求助】金属羰基化合物结构优化后CO之间为什么变成了三键？(已解决)
已经有4人回复
- 【求助】一个结构优化的问题
已经有11人回复
- 【求助】castep结构优化与性质计算
已经有5人回复
- 【求助】结构优化时设置ISIF=4时不收敛如何解决
已经有8人回复
- 【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
- 【求助】castep 结构优化中如何根据能量算出体模量B和B‘？
已经有13人回复
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-09-20 16:16:06
guolianshun
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 71.4
- 散金: 193
- 红花: 61
- 帖子: 1403
- 在线: 231小时
- 虫号: 979838
- 注册: 2010-03-23
- 性别: GG
- 专业: 无机材料化学
3楼2010-09-20 16:35:58
didi5158
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 482.4
- 帖子: 250
- 在线: 89.9小时
- 虫号: 273449
- 注册: 2006-08-20
- 性别: GG
- 专业: 无机材料化学
4楼2010-09-20 18:09:25
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2010-09-21 10:04:35
111hzs
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 132.8
- 帖子: 63
- 在线: 44.8小时
- 虫号: 1951933
- 注册: 2012-08-22
- 专业: 凝聚态物性 II ：电子结构
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
在Control部分里，有一个变量forc_conv_thr，缺省值为10^-3，这个值是不是意味着最后计算得到的原子受力要小于0.001？我在做vc-relax计算的时候让forc_conv_thr作为缺省，然后计算，为什么得到的原子受力大于0.001？ atom 1 type 1 force = -0.00112511 -0.00057496 -0.00392619 atom 2 type 1 force = 0.00021704 -0.00010523 0.00140817 atom 3 type 2 force = 0.00348517 -0.00103276 0.00343759 atom 4 type 2 force = -0.00054035 0.00199416 -0.00437423 atom 5 type 2 force = 0.00091431 0.00032076 0.00536982 atom 6 type 2 force = -0.00295107 -0.00060197 -0.00191515 Total force = 0.010488 Total SCF correction = 0.000005 这是怎么回事？？ |
6楼2013-05-20 20:17:46
111hzs
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 132.8
- 帖子: 63
- 在线: 44.8小时
- 虫号: 1951933
- 注册: 2012-08-22
- 专业: 凝聚态物性 II ：电子结构
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
另外一个问题，赝势的名字有什么含义，一般计算如何选取赝势？比如说氧的赝势有那么多，怎么选取呀？： O.blyp-mt.UPF O.pbe-n-kjpaw_psl.0.1.UPF O.pbe-van_bm.UPF O.pz-n-kjpaw_psl.0.1.UPF O.blyp-van_ak.UPF O.pz-n-rrkjus_psl.0.1.UPF O.rel-pbesol-n-rrkjus_psl.0.1.UPF O.bp-van_ak.UPF O.pbe-rrkjus.UPF O.pw91-van_ak.UPF O.pz-rrkjus.UPF O.rel-pz-n-kjpaw_psl.0.1.UPF O.pbe-kjpaw.UPF O.pbesol-n-kjpaw_psl.0.1.UPF O.pz-kjpaw.UPF O.pz-van_ak.UPF O.rel-pz-n-rrkjus_psl.0.1.UPF O.pbe-mt_fhi.UPF O.pbesol-n-rrkjus_psl.0.1.UPF O.pz-mt_fhi.UPF O.rel-pbe-n-kjpaw_psl.0.1.UPF O.rel-pz-rrkjus.UPF |
7楼2013-05-20 20:21:58
8楼2013-06-19 21:27:15