| 查看: 2084 | 回复: 16 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
wenxuan5788金虫 (正式写手)
|
[交流]
[求助]如何用vasp计算Gd单个原子能量？
|
|||
|
请问如何用vasp计算Gd的单个原子能量？INCAR里面如何设置？ [ Last edited by zeoliters on 2009-6-13 at 17:00 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 化学电源之锂离子电池 |
» 猜你喜欢
286求调剂
已经有9人回复
-
一志愿南昌大学324求调剂
已经有6人回复
-
0703化学/290求调剂/本科经历丰富/工科也可
已经有10人回复
-
291求调剂
已经有5人回复
-
308求调剂
已经有7人回复
-
085404求调剂，总分309，本科经历较为丰富
已经有4人回复
-
压国家一区线，求导师收留，有恩必谢！
已经有3人回复
-
一志愿中南大学化学0703总分337求调剂
已经有4人回复
-
083000学硕274求调剂
已经有6人回复
-
304求调剂
已经有4人回复
wenxuan5788
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 170.1
- 散金: 191
- 红花: 1
- 帖子: 658
- 在线: 43.5小时
- 虫号: 581528
- 注册: 2008-07-18
- 专业: 凝聚态物性I:结构、力学和
4楼2009-01-17 16:09:07
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
老虎大王
木虫 (著名写手)
- 应助: 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
