| 查看: 792 | 回复: 5 | ||
[求助]
关于计算之后个别分子仍然受力很大请教各位前辈……
|
|
问题1：首先上INCAR，EDIFFG设置的是-5E-2，照理说这个标准不高，但是收敛极慢，计算起来都着急，不知道这个INCAR有什么问题没有？ 问题2：不知道OUTCAR的这个TOTAL-FORCE怎么看，前三列是原子xyz坐标，后三列是原子在xyz方向受的力吗？ 问题3：运行过程的二张截图，一个是最开始的，一个是最后的，可以发现从1F到27F，F的数值一直都是.296E+03，这个应该是原子受到最大的力吧，从头到尾一直都是这个数值，一直到收敛结束，不知道怎么回事…… 最后，得到的结果与文献上相比，总是有5%的误差，不知道是不是这里出问题了  INCAR.jpg  OUTCAR-Total-Force.jpg  运行过程.jpg  运行过程2.jpg [ 来自小组 石墨烯 ] |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有9人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
关于NBANDS、NEDOS、MAXMIX的设定请教各位前辈？
已经有7人回复- 关于美国签证：急！！请教各位前辈：3个人同一张签证费收据能否分开预约面签？
已经有7人回复
- 有关一个henry反应 有问题请教各位前辈
已经有7人回复
- 请教各位前辈，同一个高分子物质如何分离成不同分子量区间的几个部分？
已经有7人回复
- 请教各位前辈，专科工作两年后读研，可以申请出国留学吗？会不会受专科学历影响？
已经有10人回复
- 请教各位前辈，我14年考研，需不需要报辅导班？
已经有16人回复
- 请教各位前辈，药物化学硕士毕业后如果不做实验，都做什么了？
已经有52人回复
- 请教各位前辈关于EMSA（gel shift）的问题！
已经有6人回复
- 请教各位前辈，关于套磁回信的问题请教
已经有12人回复
- 请教各位前辈，用木质素制备离子交换树脂的问题。
已经有4人回复
- 请教各位前辈 三电极体系测试充放电和CV时，铂片上产生大量气泡 ，改如何解决
已经有10人回复
- 请教各位前辈，强酸性溶液如何过滤
已经有12人回复
- 关于工作户口怎么解决的问题，请教各位前辈
已经有12人回复
- 【求助】请教各位前辈，上海交大的药剂怎么样？
已经有15人回复
- 【求助】向各位前辈请教
已经有31人回复
2楼2013-10-27 15:19:05
3楼2013-10-27 16:05:10
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
甲斐之虎: 金币+30, ★★★★★最佳答案, 這位前輩太?�?�?�?�?姶罅� 2013-11-04 12:07:34
fzx2008: 金币+2, 谢谢指导 2013-11-13 09:30:52
甲斐之虎: 金币+30, ★★★★★最佳答案, 這位前輩太?�?�?�?�?姶罅� 2013-11-04 12:07:34
fzx2008: 金币+2, 谢谢指导 2013-11-13 09:30:52
|
问题1：首先上INCAR，EDIFFG设置的是-5E-2，照理说这个标准不高，但是收敛极慢，计算起来都着急，不知道这个INCAR有什么问题没有？ 你指的"慢" 是指 "跑了很多ionic loops" 最簡單最直覺的答案就是 POSCAR "不好" POSCAR 給的結構不好造成 relaxation 過程需要許多功夫 至於 INCAR 的問題~~不少 我有回信在你別的帖 问题2：不知道OUTCAR的这个TOTAL-FORCE怎么看，前三列是原子xyz坐标，后三列是原子在xyz方向受的力吗？ yes 问题3：运行过程的二张截图，一个是最开始的，一个是最后的，可以发现从1F到27F，F的数值一直都是.296E+03，这个应该是原子受到最大的力吧，从头到尾一直都是这个数值，一直到收敛结束，不知道怎么回事…… NO, F 不是受力 F 是能量， 同一行 E0 也是能量 F 是 free energy, E0 是 total energy 兩者之間差一個 entropy |
4楼2013-11-04 02:30:19
5楼2013-11-04 12:14:52
★ ★
fzx2008: 金币+2, 谢谢回帖！ 2013-11-13 10:52:01
fzx2008: 金币+2, 谢谢回帖！ 2013-11-13 10:52:01
|
就第二個問題，如果?奶友Y面的圖?砜矗瑐�別原子的受力似乎仍然很大，比如a31(矩陣第三行第一列的元素)是-0.126333，這個力仍然很大，如果我需要收斂到每個原子受力小於5E-2（EDIFFG=-5E-2），那麼仍然是沒有達到標準，是這樣的嗎，前輩？ 首先 圖太不清楚了 老人家眼睛不好~看得真的很吃力 再者 這個 TOTAL-FORCE data 到底是第一個 ionic loop 還是最後一個 ionic loop ?? 第三個問題，是該用F還是E0計算吸附能呢？這個問題我比較頭疼，現在我的靜態計算都採用ISMEAR=-5，這樣就沒有entropy了，迴避了這個問題 用 F or E0 都可以 只要一致就好 不要 case A 用 F case B 用 E0 然後又把兩個拿去相減 (個人喜歡 E0 ... 沒有任何原因) 如果你沒有計算 MD ~ 那就不會有 entropy 的問題 跟你是不是用 ISMEAR -5 好像是沒有關連性 |
6楼2013-11-04 17:24:14