24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2085  |  回复: 16
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wenxuan5788

金虫 (正式写手)

[交流] [求助]如何用vasp计算Gd单个原子能量?

请问如何用vasp计算Gd的单个原子能量?INCAR里面如何设置?

[ Last edited by zeoliters on 2009-6-13 at 17:00 ]
回复此楼

» 收录本帖的淘帖专辑推荐

化学电源之锂离子电池

» 猜你喜欢

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

老虎大王

木虫 (著名写手)

★ ★ ★ ★ ★ ★
wuchenwf(金币+3):很好,性情中人啊 1-17 21:09
wenxuan5788(金币+3):非常感谢 1-19 08:15
你气死我了。我这很多天都在进行Gd的运算。用了好几个软件,都算不出来。可能是势函数不好。你用VASP居然一下子算出来了。我早就感到VASP可能有戏,但是我没有这个软件啊。

我来告诉你我的看法:
1.你光建立一个大单胞是不行的。你要建几个,比如15*15*15,20*20*20,25*25*25等,分别进行计算,看看能量是不是收敛到一个值,也就是说,看看元胞大到什么程度的时候,不影响计算结果。
2.元胞尺寸定下来以后,再考察K-point,以及能量截断的影响。办法就是多换几个参数去试,直到能量收敛(一般两次计算能量之差不超过1mhar),也就是这些参数不影响能量值为止。因为只算一个原子,应该很快的。
3.最后你应该得到一个你确信的值。

你问Ga、La单原子的能量是多少,这是没法回答你的。因为大家如果用不同的势或者不同的算法去算,结果绝对值可能有很大不同。你只能以如下方式检验你的结果:

A。计算Gd单个原子能量。
B。对Gd晶体进行结构优化(只优化元胞尺寸),并计算Gd晶体中单个原子能量。
C。求得Gd的结合能(cohesive energy)。
如果你的结合能算出来是对的,那么你的单原子能量就是可以使用的。

Gd的结合能Ecoh,文献值一般在4.14~4.18eV,
Gd的结构如下:
空间群:P6 3/mmC (194#),HCP,Mg型,a=3.615ang,c=5.777ang。

以前也有人计算过如下结果:a=6.873Bohr,c=10.897Bohr, Ecoh=4.344eV,给出的实验值则是a=10.902Bohr,c=6.858Bohr。
这篇文章是PRB, 50(1994),1363。
还有一篇文章在PRB,44(1991),7451。你自己查吧。这两篇文章都是用LSDA算的。

对La和其它原子也是同样的处理和检验。对O,可以通过计算O2分子的原子化能来考察计算的有效性。

我以上说的是比较保险的办法,步步检验。当然这样搞的话比较慢一些,而且需要人很勤奋。呵呵

不过根据我的经验,稀土元素在DFT中是很难处理的。所以以上只供你参考。

如果你算出了好的结果,尤其是如果用GGA算出了好的结果,别忘了跟我联系一下,或者发个帖子讲一下你用的势和Gd的计算结果。谢谢。

[ Last edited by 老虎大王 on 2009-1-17 at 20:18 ]
5楼2009-01-17 20:15:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

wenxuan5788

金虫 (正式写手)

要计算该能量主要为了计算体系的结合能,请高手详细告诉一下,我对这方面不太懂。谢谢
2楼2009-01-16 13:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomagh

铁杆木虫 (正式写手)

★ ★ ★ ★ ★ ★
wuchenwf(金币+4,VIP+0):xiexie讨论
wenxuan5788(金币+2):谢谢 1-19 08:16
引用回帖:
Originally posted by wenxuan5788 at 2009-1-16 12:29:
请问如何用vasp计算Gd的单个原子能量?INCAR里面如何设置?

楼主说的不是很清楚,可以这么回答你,你要是想获得Gd组成某种结构时,单个原子的能量的话,那么你就将Gd按照某种结构进行建模,比如HCP,BCC, FCC等,然后计算总能量除以原子数目即可。
如果你想计算一个孤立原子的能量的话,那么你就要建立一个很大的包,只放一个原子于其中,比如10angstrom的大包。然后计算能量,建议看一下VASP的user guide 上面有。
3楼2009-01-16 17:33:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenxuan5788

金虫 (正式写手)

我计算的体系是La7Gd1Mo8O36。我想计算一下体系的结合能,我看资料上说需要计算单个原子的能量。我计算Gd单个原子能量的步骤是:先建立一个大的单包:15X15X15。只放一个Gd原子,然后计算能量。我计算结果是-0.210eV,我觉得不太对,请问Gd,La的单个原子能量一般是多少?
4楼2009-01-17 16:09:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料求调剂一志愿哈工大324 +6 闫旭东 2026-03-28 8/400 2026-03-28 08:11 by 闫旭东
[考研] 291求调剂 +6 HanBeiNingZC 2026-03-24 6/300 2026-03-28 07:55 by baoball
[考研] 320分,材料与化工专业,求调剂 +7 一定上岸aaa 2026-03-27 9/450 2026-03-28 07:07 by wangy0907
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +6 辰沐5211314 2026-03-26 6/300 2026-03-28 03:31 by fmesaito
[考研] 0703化学 +6 妮妮ninicgb 2026-03-27 6/300 2026-03-27 17:18 by shangxh
[考研] 一志愿北化085600材料专硕275|有文章专利|求调剂 +3 Micky11223 2026-03-25 3/150 2026-03-27 14:52 by 醉在风里
[考研] 348求调剂 +4 小懒虫不懒了 2026-03-27 5/250 2026-03-27 12:47 by 果果妈咪
[考研] 316求调剂 +5 江辞666 2026-03-26 5/250 2026-03-27 08:08 by hypershenger
[考研] 一志愿武汉理工,总分321,英一数二,求老师收留。 +5 nnnnnnn5 2026-03-25 5/250 2026-03-27 04:42 by wxiongid
[考研] 求调剂 +5 芦lty 2026-03-25 6/300 2026-03-26 20:49 by 不吃魚的貓
[考研] 324求调剂 +4 wysyjs25 2026-03-21 4/200 2026-03-26 20:38 by fmesaito
[考研] 一志愿211 初试270分 求调剂 +6 谷雨上岸 2026-03-23 7/350 2026-03-26 18:55 by 不吃魚的貓
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 085602 289分求调剂 +8 WWW西西弗斯 2026-03-24 8/400 2026-03-26 16:33 by 不吃魚的貓
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 打过很多竞赛,085406控制工程300分,求调剂 +3 askeladz 2026-03-26 3/150 2026-03-26 09:08 by 给你你注意休息
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +7 niko- 2026-03-22 7/350 2026-03-25 20:14 by qingfeng258
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
信息提示
请填处理意见