24小时热门版块排行榜

>论坛更新日志 (4134)
>文献求助 (440)
>虫友互识 (359)
>导师招生 (345)
>考博 (187)
>博后之家 (170)
>硕博家园 (160)
>休闲灌水 (147)
>论文投稿 (127)
>基金申请 (91)
>考研 (60)
>招聘信息布告栏 (59)
>教师之家 (54)
>公派出国 (43)
>绿色求助(高悬赏) (40)
>找工作 (38)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

小小鑫000

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 6.6
帖子: 267
在线: 111.8小时
虫号: 2648268
注册: 2013-09-11
性别: MM
专业: 理论和计算化学

[求助] 会用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原子的自由能到底该按上面哪种情况设置，还是采用别的设置方法。我这方面真是一塌糊涂，还望前辈细心解答，不胜感激！

回复此楼

好好研究，天天向上！

1楼2014-05-03 19:46:35

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

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
小小鑫000: 金币+5 2014-05-05 10:06:45

提問問題跟寫 PAPER ，兩者其�?很接近
如果你能講一下所謂侯老師的手冊在那裡 (reference)
這樣其他沒有看過侯老師手冊的人會更容易了解前因後果

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

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

或~更直接的講法就是
用 "同樣的 INCAR" 去計算單一顆原子總能，即可