24小时热门版块排行榜    

查看: 2992  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

rustboy

铜虫 (初入文坛)

[求助] vasp 计算单分子/原子 时晶胞大小的选取

我觉得,原则上讲应该是使相邻晶胞内分子/原子距离足够大,以消除计算过程中互相的影响。满足这一条件后,在相同的精度下,所得到的单分子/原子的能量大致应该是相当的。

但是本人在计算单个金属Ni时发现,使用10*10*10的晶胞以及15*15*15晶胞,其它参数一样。计算所得单原子的能量相差较大。

不知道大家平时计算类似的体系时,是如何设置晶胞大小的?
回复此楼

» 收录本帖的淘帖专辑推荐

化学电源之锂离子电池

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

rustboy

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by hn19870519 at 2012-12-31 22:17:09
你说的能量相差很大是多大?最好把你计算的参数贴出来,这样大家也好帮你。

如下是我的INCAR参数:
  NWRITE =      2;   LPETIM=F    write-flag & timer
   ISTART =      0   job   : 0-new  1-cont  2-samecut
   ICHARG =      2    charge: 1-file 2-atom 10-const  
   INIWAV =      1    electr: 0-lowe 1-rand  2-diag

Electronic Relaxation 1
   PREC = Accurate ! standard precision
   ENCUT =       400.000
   NELM   =      500;   NELMIN= 2; NELMDL= 0     # of ELM steps
   ISPIN  =      2    spin polarized calculation?
   EDIFF  =      0.1E-04  

Ionic Relaxation
   EDIFFG =    -0.03  
   NSW    =      500  number of steps for IOM
   NBLOCK =      1;   
   IBRION =      2   
   ISIF   =      2   
   ISYM   =      2    0-nonsym 1-usesym
   LCORR  =      T    Harris-correction to forces
   POTIM  =     0.5    time-step for ion-motion

Electronic Relaxation 2
   LPLANE =   .TRUE.
   ALGO   =   Fast
   LDIAG  =   T     sub-space diagonalisation
   LWAVE   =  F     write WAVECAR
   LCHARG  = F      write CHGCAR and CHG
   LVDW =   .TRUE.
   LREAL  =   Auto     real-space projection
   NSIM = 4
昨天又算了一个Co原子,情况也是类似,10*10*10的能量是1.480eV,15的是1.93,20的是1.476。
谢前辈帮忙!
4楼2013-01-01 09:51:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

emilyoyang

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2013-01-03 12:09:43
这个能量差别很多  可能是你的cutoff 取得过小。。。不知道你取多少

你试着计算一下9.5*10*10.5 和14.5*15*15.5这两个超胞 中单个Ni原子的能量差 是不是也很大? 也有可能是对称性造成的
2楼2012-12-30 20:19:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢提示 2013-01-03 12:09:52
你说的能量相差很大是多大?最好把你计算的参数贴出来,这样大家也好帮你。
3楼2012-12-31 22:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rustboy

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by emilyoyang at 2012-12-30 20:19:41
这个能量差别很多  可能是你的cutoff 取得过小。。。不知道你取多少

你试着计算一下9.5*10*10.5 和14.5*15*15.5这两个超胞 中单个Ni原子的能量差 是不是也很大? 也有可能是对称性造成的

cutoff我取的是400,具体INCAR参数请见三楼
我新算了一个Co原子。
结果是10*10*10的超胞能量1.480eV,9.5*10*10.5的1.785eV,15*15*15的1.931, 14.5*15*15.5的1.944。貌似没有规律性。
这些能量差距,请问前辈是何原因?

另外,本人菜鸟,前辈提到的 “对称性造成的” 是什么意思?
5楼2013-01-01 10:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见