24小时热门版块排行榜

返回列表

Brilliantly

捐助贵宾 (正式写手)

应助: 7 (幼儿园)
金币: 420
散金: 2788
红花: 8
帖子: 768
在线: 312.2小时
虫号: 1360339
注册: 2011-08-03
性别: GG
专业: 功能陶瓷

引用回帖:

9楼: Originally posted by lfhuang at 2015-03-18 03:20:58
哦，对了，看自由能的变化不是看dE，那是判定收敛的能量差别。你要分别计算不同SIGMA下的能量，看OSZICAR里F=后面那个数值的变化。
PS:你不是专门搞理论的吧？

对了，我的结构只含有2个Zr原子！

赞一下

回复此楼

加油！

11楼2015-03-18 16:39:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 4318.2
散金: 1541
红花: 36
帖子: 1403
在线: 676.7小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
Brilliantly: 金币+100, ★★★★★最佳答案, 我确实用了两个原子的模型，感谢您的回答，总算算起来差不多了，再送您100金币作为回报！ 2015-03-19 09:15:42

引用回帖:

10楼: Originally posted by Brilliantly at 2015-03-18 16:38:20
我不是搞理论的，我想计算有限温度下的弹性常数，审稿人说自由能要加上这一项才算正确。我按照您的方法，算了当sigma=0时，Zr的F=16.955219；当sigma=0.08625 （应该是1000K对应的sigma吧）时，F=17.066694，这样的 ...

SIGMA取0.0的时候，由于Fermi-Dirac分布的公式，会出现一定的数值误差，建议用一个非常小的温度作为参考，比如10 K。你这个结果（数值误差110 meV，两个原子的吧？）挺有意思的，我在Ti（和Zr类似）上也计算了一下，确实也发现了45 meV per atom的误差，应该和你的情况类似。实际上，1000K的温度会造成-13 meV per atom的自由能下降，和我前面说的类似。看了附件中的图，你就明白了。