| 查看: 2107 | 回复: 16 | ||
| 当前主题已经存档。 | ||
wenxuan5788金虫 (正式写手)
|
[交流]
[求助]如何用vasp计算Gd单个原子能量？
|
|
|
请问如何用vasp计算Gd的单个原子能量？INCAR里面如何设置？ [ Last edited by zeoliters on 2009-6-13 at 17:00 ] |
回复此楼
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★
wuchenwf(金币+3):很好，性情中人啊 1-17 21:09
wenxuan5788(金币+3):非常感谢 1-19 08:15
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
xiaomagh
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 6113.1
- 红花: 5
- 帖子: 652
- 在线: 150.8小时
- 虫号: 224710
- 注册: 2006-03-22
- 性别: GG
- 专业: 金属材料的合金相、相变及
8楼2009-01-19 23:24:10
wenxuan5788
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 170.1
- 散金: 191
- 红花: 1
- 帖子: 658
- 在线: 43.5小时
- 虫号: 581528
- 注册: 2008-07-18
- 专业: 凝聚态物性I:结构、力学和
2楼2009-01-16 13:50:33
xiaomagh
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 6113.1
- 红花: 5
- 帖子: 652
- 在线: 150.8小时
- 虫号: 224710
- 注册: 2006-03-22
- 性别: GG
- 专业: 金属材料的合金相、相变及
3楼2009-01-16 17:33:33
wenxuan5788
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 170.1
- 散金: 191
- 红花: 1
- 帖子: 658
- 在线: 43.5小时
- 虫号: 581528
- 注册: 2008-07-18
- 专业: 凝聚态物性I:结构、力学和
4楼2009-01-17 16:09:07
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
6楼2009-01-17 20:24:56
wenxuan5788
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 170.1
- 散金: 191
- 红花: 1
- 帖子: 658
- 在线: 43.5小时
- 虫号: 581528
- 注册: 2008-07-18
- 专业: 凝聚态物性I:结构、力学和
7楼2009-01-19 08:25:29
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★
fegg7502(金币+2):鼓励交流！ 1-20 16:09
fegg7502(金币+2):鼓励交流！ 1-20 16:09
|
楼上回答的是对的。 结合能＝晶体中每个原子的能量－孤立的单个原子能量。 我说的晶体中的单个原子能量，就是指晶体能量，但由于元胞设置问题，故而强调是平均到每个原子的能量。 你列出的公式也是对的。即晶体总能－孤立的单个原子能量×晶体中原子个数。显然这两个公式结果相差n倍，n＝晶体中原子个数。 但一般所说的结合能的量纲是 能量/每原子。 DFT计算，有意义的不是能量绝对值而是相对的能量值。你算出来的晶体能量和孤立原子能量的值都不重要，关键是二者之差是否符合实验值。 [ Last edited by 老虎大王 on 2009-1-20 at 12:02 ] |
9楼2009-01-20 11:59:43
wenxuan5788
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 170.1
- 散金: 191
- 红花: 1
- 帖子: 658
- 在线: 43.5小时
- 虫号: 581528
- 注册: 2008-07-18
- 专业: 凝聚态物性I:结构、力学和
10楼2009-02-11 09:29:41
