应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于计算之后个别分子仍然受力很大请教各位前辈……
51/1返回列表
查看: 797  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

甲斐之虎

银虫 (正式写手)

[求助] 关于计算之后个别分子仍然受力很大请教各位前辈……

问题1：首先上INCAR，EDIFFG设置的是-5E-2，照理说这个标准不高，但是收敛极慢，计算起来都着急，不知道这个INCAR有什么问题没有？
问题2：不知道OUTCAR的这个TOTAL-FORCE怎么看，前三列是原子xyz坐标，后三列是原子在xyz方向受的力吗？
问题3：运行过程的二张截图，一个是最开始的，一个是最后的，可以发现从1F到27F，F的数值一直都是.296E+03，这个应该是原子受到最大的力吧，从头到尾一直都是这个数值，一直到收敛结束，不知道怎么回事……
最后，得到的结果与文献上相比，总是有5%的误差，不知道是不是这里出问题了
关于计算之后个别分子仍然受力很大请教各位前辈……
INCAR.jpg


关于计算之后个别分子仍然受力很大请教各位前辈……-1
OUTCAR-Total-Force.jpg


关于计算之后个别分子仍然受力很大请教各位前辈……-2
运行过程.jpg


关于计算之后个别分子仍然受力很大请教各位前辈……-3
运行过程2.jpg



[ 来自小组 石墨烯 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
1楼2013-10-26 21:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)
引用回帖:
2楼: Originally posted by sungjen at 2013-10-27 15:19:05
这个收敛速度也跟体系有关吧，不知道先优化了没，第二个问题应该你说的是对的，这个total drift 目前还不知道是什么意思，第三个问题我也等待高手解答，包括DAV:1 这一行整个都需要解释，能力有限啊 。

多谢参与啊，大侠，我们一起坐等大神出现吧
赞一下 回复此楼
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
3楼2013-10-27 16:05:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

sungjen

银虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2013-10-28 07:36:50
这个收敛速度也跟体系有关吧，不知道先优化了没，第二个问题应该你说的是对的，这个total drift 目前还不知道是什么意思，第三个问题我也等待高手解答，包括DAV:1 这一行整个都需要解释，能力有限啊 。
赞一下(1人) 回复此楼
希望有一天成为大神
2楼2013-10-27 15:19:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
甲斐之虎: 金币+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
赞一下(2人) 回复此楼
4楼2013-11-04 02:30:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)
引用回帖:
4楼: Originally posted by jpchou at 2013-11-04 02:30:19
问题1：首先上INCAR，EDIFFG设置的是-5E-2，照理说这个标准不高，但是收敛极慢，计算起来都着急，不知道这个INCAR有什么问题没有？

你指的"慢" 是指 "跑了很多ionic loops"
最簡單最直覺的 ...

感謝前輩指教！
就第二個問題，如果?奶友Y面的圖?砜矗瑐�別原子的受力似乎仍然很大，比如a31(矩陣第三行第一列的元素)是-0.126333，這個力仍然很大，如果我需要收斂到每個原子受力小於5E-2（EDIFFG=-5E-2），那麼仍然是沒有達到標準，是這樣的嗎，前輩？
第三個問題，是該用F還是E0計算吸附能呢？這個問題我比較頭疼，現在我的靜態計算都採用ISMEAR=-5，這樣就沒有entropy了，迴避了這個問題
赞一下 回复此楼
别人笑我太疯颠，我笑他人看不穿。不见五陵豪杰墓，无花无酒锄作田。
5楼2013-11-04 12:14:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭