应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 会用VASP计算单个原子能量的高手进来！
41/1返回列表
查看: 3006  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小小鑫000

铁虫 (小有名气)

[求助] 会用VASP计算单个原子能量的高手进来！已有2人参与

各位前辈大家好：
过渡金属Ni"(n=1一13)团簇"团簇中每个原子的平均束缚能EB采用如下公式计算:
EB=(NxE atom一Ecluster)/n     ；Eamot和Ecluster分别代表单个Ni原子的自由能和所研究原子团簇的自由能。
这是我从文献上看到的，我想计算一下单个Ni 原子的自由能。但不知道怎么算，于是从侯老师的手册上得知了两种计算原子能量的方法。
他说单个原子的计算有两个目的：1) 检验赝势的好坏；2) 对称性被破坏后自旋极化情况下的原子基态能量，对结合能进行修正。
然后就第一个目的，他INCAR设置如下：
SYSTEM = Al: atom
ENCUT = 250.00 eV   
NELMDL = 5
ISMEAR = 0; SIGMA=0.1
就第二个目的，INCAR设置为：
SYSTEM = Al: atom
ENCUT = 250  
ISYM = 0 ! no symmetry
ISPIN = 2 ! allow for spin polarisation
VOSKOWN = 1    ! this is important, in particular for GGA
ISMEAR = 0 ! Gaussian smearing, otherwise negative occupancies
SIGMA = 0.1 ! intermid. smearing width
AMIX = 0.2 ! mixing set manually
BMIX = 0.0001
NELM = 20 ! 20 electronic steps
ICHARG = 1
我不明白上面的意思，为什要分两种情况，主要是不知道我所要算的Ni原子的自由能到底该按上面哪种情况设置，还是采用别的设置方法。我这方面真是一塌糊涂，还望前辈细心解答，不胜感激！
回复此楼

» 本帖@通知

@jpchou @前田庆次

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
好好研究，天天向上！
1楼2014-05-03 19:46:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
小小鑫000: 金币+5 2014-05-04 09:31:33
若是你不做自旋极化直接自洽计算得到总能即可，若是自旋极化加上自旋项，然后自洽得到总能即可。
赞一下 回复此楼
2楼2014-05-03 23:38:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by mywai520 at 2014-05-03 23:38:18
若是你不做自旋极化直接自洽计算得到总能即可，若是自旋极化加上自旋项，然后自洽得到总能即可。

谢谢前辈，前辈的意思是说无自旋按照第一种设置和有自旋按第二种方法设置INCAR 吗？

[ 发自手机版 http://muchong.com/3g ]
赞一下(1人) 回复此楼
好好研究，天天向上！
3楼2014-05-04 09:35:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
小小鑫000: 金币+5 2014-05-05 10:06:45
提問 問題 跟寫 PAPER ，兩者其�?很接近
如果你能講一下所謂侯老師的手冊在那裡 (reference)
這樣其他沒有看過侯老師手冊的人會更容易了解前因後果


我?木W路上找了一下你講的內容，找到如下網頁
http://wenku.baidu.com/view/a98cfc2c4b73f242336c5ff3
?倪@講義內容?砜� (假設這個講義的確就是侯老師上課用講義)
就你目前想要做的內容?碚f (計算金屬團簇束縛能)
跟你所謂的  "  1) 检验赝势的好坏" 完全沒有關係...... 不是嬤 ??
為什麼會  "主要是不知道我所要算的Ni原子的自由能到底该按上面哪种情况设置" 呢 ?


你的問題就是單純地就是要求束縛能
束縛能的計算就是要計算出兩個能量
團簇總能 & 單一顆原子總能 (這裡用總能似乎是很蠢的事....也不過就一顆...)
你知道如何計算團簇總能
同理
用 "同樣的方式"  去計算單一顆原子總能，即可
(自旋的部份需小心)

或~更直接的講法就是
用 "同樣的 INCAR" 去計算單一顆原子總能，即可
赞一下 回复此楼
4楼2014-05-05 03:13:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小鑫000 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭