24小时热门版块排行榜    

查看: 1624  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 6/300 2026-08-11 13:19 by JACS/Angew
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 帮忙看看fileCode +6 wwncly 2026-08-10 12/600 2026-08-11 13:01 by wwncly
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +5 布布和一二 2026-08-10 6/300 2026-08-11 10:15 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +15 Lanmanbaby 2026-08-09 23/1150 2026-08-11 09:51 by demonstreets
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +6 Tide man 2026-08-10 6/300 2026-08-11 09:45 by tangmnt
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见