| 查看: 2609 | 回复: 10 | |||
lulutiantian新虫 (小有名气)
|
[求助]
Au 自由原子计算的结果，大家帮忙看下对不对？
|
||
|
因为要计算一个体系的结合能，so计算一个金原子在自由态的能量，建立一个15*15*15A的空格子，将一个Au原子放在中心。未做截断能测试，选取截断能ecut为360ev,进行结构优化，然后静态计算，k点设为3*3*3，以Gamma点为中心产生，静态计算后结果 E=-0.10956204 而候博士的入门指南中以Al为例说˜„œ单原子能量在1mev-10mev之间。也就是vasp计算的单个原子能量与原子参考组态的能量差。在vasp所计算的总能是扣除了原子参考组态的能量，即POTCAR中EATOM的值。 问题 1.做自由原子能量计算不测试截断能可以吗？ 2.我的这个E=-0.10956204很明显不在1mev-10mev之间。这个值可以作为Au自由原子的能量吗？ 3.我在计算中的发现 Au 自由原子优化结算结果 E= -0.116817 Au自由原子静态计算结果E=-0.10956204 发现这两个值是有差别的。那么应该是以静态计算的结果为结合能。 希望做过这方面计算的虫友赐教！ |
回复此楼» 猜你喜欢
070300化学279求调剂
已经有19人回复
-
材料考研调剂
已经有19人回复
-
272分材料子求调剂
已经有41人回复
-
求调剂，262机械专硕
已经有7人回复
-
266求调剂，一志愿哈工程电子信息，本科获多项国奖和省奖
已经有6人回复
-
0860004 求调剂 309分
已经有8人回复
-
270求调剂
已经有5人回复
-
347求调剂
已经有3人回复
-
280求调剂
已经有9人回复
-
药学305求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一篇PRL文献讨论，电子自旋与电子运动产生磁场的相互作用
已经有48人回复
- SPR 增强是指增强什么，一般有哪些增强方法？？
已经有20人回复
- H取代O，是SnO2本征n型的原因！《文献学习》PRL 101, 055502 (2008)
已经有76人回复
- 求Au原子模守恒赝势
已经有8人回复
- 亲电吸附基团协助下Au面上CO的稳定化（PRL, 2005）
已经有13人回复
- 关于“密度泛函（DFT）”，帮忙翻译一段话！
已经有13人回复
- 测试几种常用化学键分析方法的普适性
已经有8人回复
- 原子转移自由基聚合
已经有8人回复
- 【求助】请教用高斯09算溶剂化单点时，怎样得到溶剂中的Gibbs自由能？
已经有9人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
- 【资源】有用的链接@厦门大学理论化学研究中心
已经有8人回复
- 【求助】Castep中如何计算自由原子的能量？
已经有10人回复
- 【求助】请问用cc-pvtz-pp基组算一个C两个Au原子，怎么写输入文件，谢谢。
已经有12人回复
- 【求助】纳米线拉伸性时，截面的面积大小
已经有20人回复
wanglikun08
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1261.8
- 散金: 25
- 红花: 2
- 帖子: 212
- 在线: 114.1小时
- 虫号: 1411874
- 注册: 2011-09-21
- 专业: 凝聚态物性 II ：电子结构
2楼2013-11-22 21:11:30
lulutiantian
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 328.3
- 散金: 1
- 红花: 1
- 帖子: 116
- 在线: 137.2小时
- 虫号: 2544166
- 注册: 2013-07-13
- 专业: 凝聚态物性I:结构、力学和
3楼2013-11-22 21:20:08
mika
金虫 (正式写手)
纳米宝贝
- 应助: 6 (幼儿园)
- 金币: 992.4
- 散金: 76
- 帖子: 402
- 在线: 121.1小时
- 虫号: 836617
- 注册: 2009-08-29
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
4楼2013-11-23 15:49:39
lulutiantian
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 328.3
- 散金: 1
- 红花: 1
- 帖子: 116
- 在线: 137.2小时
- 虫号: 2544166
- 注册: 2013-07-13
- 专业: 凝聚态物性I:结构、力学和
5楼2013-11-23 16:32:15
wanglikun08
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1261.8
- 散金: 25
- 红花: 2
- 帖子: 212
- 在线: 114.1小时
- 虫号: 1411874
- 注册: 2011-09-21
- 专业: 凝聚态物性 II ：电子结构
6楼2013-11-23 16:34:25
lulutiantian
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 328.3
- 散金: 1
- 红花: 1
- 帖子: 116
- 在线: 137.2小时
- 虫号: 2544166
- 注册: 2013-07-13
- 专业: 凝聚态物性I:结构、力学和
7楼2013-11-23 17:00:53
superdonnie
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 150.1
- 红花: 1
- 帖子: 62
- 在线: 44小时
- 虫号: 2819306
- 注册: 2013-11-22
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★
感谢参与，应助指数 +1
lulutiantian: 金币+1 2013-11-24 15:41:33
感谢参与，应助指数 +1
lulutiantian: 金币+1 2013-11-24 15:41:33
|
本人用的不是vasp，所以针对一般情况随便说说，不对的话楼主指正： 1）求结合能应该分别得出Au原子的能量和Au晶体的能量，然后做差。 2）选定赝势以后，对于Au原子和Au晶体，至少有一个必须测定截断能，截断能必然不能随便乱选的，最好自己测测，或者有前人的经验也行。Au原子和Au晶体可以使用同样大小的截断能。 3）K点跟正格子的大小有关系，对于Au原子来说，需要使用一个大格子模拟，所以一个K点基本是够的，Au晶体格子小，需要密集一些的K点。 4）Au原子直接放在格子什么地方都行，一个格子里一个Au原子，不需要做结构优化，但格子应该取足够大，以减少在周期性条件下各个格子之间的相互作用。Au晶体的话，对于一般的fcc的Au来说，需要求出最低能量的晶格常数，可以做varible-cell relaxation，或者直接看各个晶格常数下Au晶体的能量，然后选取对应最低能量的晶格常数。 |
8楼2013-11-24 14:23:23
lulutiantian
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 328.3
- 散金: 1
- 红花: 1
- 帖子: 116
- 在线: 137.2小时
- 虫号: 2544166
- 注册: 2013-07-13
- 专业: 凝聚态物性I:结构、力学和
9楼2013-11-24 15:36:39
superdonnie
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 150.1
- 红花: 1
- 帖子: 62
- 在线: 44小时
- 虫号: 2819306
- 注册: 2013-11-22
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
10楼2013-11-24 16:00:44
