24小时热门版块排行榜    

查看: 1249  |  回复: 13

Brilliantly

捐助贵宾 (正式写手)

[求助] 求Zr的热电子激发对自由能的贡献(Thermal electron contributions)100金币奉上

求高手帮我计算一下Zr的热电子激发对自由能的贡献,态密度已经用VASP计算出来了,在Excel文件中,Zr的费米能级为8.9633,所需要的公式在下面的word文件中,实在是由于水平有限,都不会去积分啊,求助帮忙算一下,顺便说一下所需要的积分软件和计算的结果,谢谢了!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Zr的态密度_-_副本.xls
  • 2015-03-15 16:00:21, 244 K
  • 附件 2 : 热电子计算公式.doc
  • 2015-03-15 16:00:46, 27 K
  • 附件 3 : 计算公式.jpg
  • 2015-03-15 16:02:45, 32.29 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

加油!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Brilliantly

捐助贵宾 (正式写手)

自己先顶一下!
求Zr的热电子激发对自由能的贡献(Thermal electron contributions)100金币奉上
计算公式.jpg

加油!
2楼2015-03-15 16:04:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Brilliantly

捐助贵宾 (正式写手)

给自己顶一下!
加油!
3楼2015-03-15 19:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
fzx2008: 金币+2, 谢谢指导 2015-03-17 09:16:13
Brilliantly: 金币+100, ★★★很有帮助, 请问到底用什么软件进行编程啊?还有怎么去学啊?请指导,本人对此一窍不通啊!我明白之后再给你加100金币,谢谢了! 2015-03-17 13:53:34
既然你自己知道公式了,那么就可以编程计算了。不过如果楼主比较懒的话,那我就给你一个懒人的方法
(请一定给我足够的金币哦):
在VASP INCAR里面,设置:
ISMEAR = -1
然后离散地调节SIGMA,从0调到k_B*T,得到不同地自由能(OUTCAR or OSZICAR里面有),然后用简单地函数去拟合出一个连续地变化关系。一般对于Zr,1000 K内地电子激发导致的自由能下降在10 meV/atom左右。
先进使役材料计算(https://aidme.nimte.ac.cn)
4楼2015-03-16 22:44:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

引用回帖:
4楼: Originally posted by lfhuang at 2015-03-16 22:44:51
既然你自己知道公式了,那么就可以编程计算了。不过如果楼主比较懒的话,那我就给你一个懒人的方法
(请一定给我足够的金币哦):
在VASP INCAR里面,设置:
ISMEAR = -1
然后离散地调节SIGMA,从0调到k_B*T, ...

另外,很多人都拿这种计算来验证通过公式计算得到的结果。
先进使役材料计算(https://aidme.nimte.ac.cn)
5楼2015-03-16 22:46:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Brilliantly

捐助贵宾 (正式写手)

引用回帖:
4楼: Originally posted by lfhuang at 2015-03-16 22:44:51
既然你自己知道公式了,那么就可以编程计算了。不过如果楼主比较懒的话,那我就给你一个懒人的方法
(请一定给我足够的金币哦):
在VASP INCAR里面,设置:
ISMEAR = -1
然后离散地调节SIGMA,从0调到k_B*T, ...

还有一个问题,你所说的k_B*T中的k_B是不是指玻尔兹曼常数呢?调节SIGMA,从0调到k_B*T时一般间隔是多少呢?多谢了
加油!
6楼2015-03-17 15:41:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Brilliantly

捐助贵宾 (正式写手)

引用回帖:
5楼: Originally posted by lfhuang at 2015-03-16 22:46:20
另外,很多人都拿这种计算来验证通过公式计算得到的结果。...

我用你的方法取ISMEAR=-1,SIGMA=1.3806488E-20(这样应该算是T=1000K吧?)计算了Zr的自由能,发现自由能变化dE=0,这是怎么回事呢?下面是我的INCAR和OSZICAR不知道哪里出错了呢?
求Zr的热电子激发对自由能的贡献(Thermal electron contributions)100金币奉上-1
2.png


求Zr的热电子激发对自由能的贡献(Thermal electron contributions)100金币奉上-2
1.png

加油!
7楼2015-03-17 19:01:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

【答案】应助回帖

(1)自己编程用什么编程语言都行的,比如Fortran, C, C++, python ...
(2)计算SIGMA的时候,请注意单位。300 K对应的是SIGMA=0.025875 eV。
先进使役材料计算(https://aidme.nimte.ac.cn)
8楼2015-03-18 03:14:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

哦,对了,看自由能的变化不是看dE,那是判定收敛的能量差别。你要分别计算不同SIGMA下的能量,看OSZICAR里F=后面那个数值的变化。
PS:你不是专门搞理论的吧?
先进使役材料计算(https://aidme.nimte.ac.cn)
9楼2015-03-18 03:20:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Brilliantly

捐助贵宾 (正式写手)

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

我不是搞理论的,我想计算有限温度下的弹性常数,审稿人说自由能要加上这一项才算正确。我按照您的方法,算了当sigma=0时,Zr的F=16.955219;当sigma=0.08625 (应该是1000K对应的sigma吧)时,F=17.066694,这样的话就相差了111.475meV,居然差那么大,不知道哪里搞错了
加油!
10楼2015-03-18 16:38:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Brilliantly 的主题更新
信息提示
请填处理意见