24小时热门版块排行榜    

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

rustboy

铜虫 (初入文坛)

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

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

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

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

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

化学电源之锂离子电池

» 猜你喜欢

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

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

已阅   回复此楼   关注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的回帖
查看全部 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的回帖

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 324分 +5 llllkkkhh 2026-03-18 5/250 2026-03-18 15:24 by 星空星月
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 281求调剂(0805) +5 烟汐忆海 2026-03-16 13/650 2026-03-18 14:30 by stone_128
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 085600材料与化工 +4 安全上岸! 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
信息提示
请填处理意见