24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 3035  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

rustboy

铜虫 (初入文坛)

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

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

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

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

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

化学电源之锂离子电池

» 猜你喜欢

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

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

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

铜虫 (初入文坛)

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326分,一志愿沪9,求生物学调剂 +3 刘墨墨 2026-04-05 3/150 2026-04-08 00:57 by 尼尼宁宁
[考研] 一志愿吉大化学327求调剂 +10 王王白石 2026-04-06 11/550 2026-04-07 23:54 by JourneyLucky
[考研] 材料专硕调剂 +9 哈哈哈吼吼吼哈 2026-04-07 10/500 2026-04-07 22:48 by 来看流星雨10
[考研] 土木水利专硕276分求调剂 +4 我想上学!!6 2026-04-05 7/350 2026-04-07 22:10 by Delta2012
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 材料调剂 +10 一样YWY 2026-04-06 10/500 2026-04-06 21:05 by lbsjt
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] 346分的生物与医药08600求调剂 +6 常雨阳上岸 2026-04-05 7/350 2026-04-06 12:36 by lys0704
[考研] 0855求调剂材料 +11 红桃灼灼 2026-04-04 12/600 2026-04-06 10:26 by 蓝云思雨
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +15 路痴小琪 2026-04-05 15/750 2026-04-06 09:27 by cql1109
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 13/650 2026-04-06 07:24 by houyaoxu
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 材料与化工306分找调剂 +23 沧海轻舟e 2026-04-02 27/1350 2026-04-04 21:52 by laoshidan
[考研] 求调剂 +6 朔朔话 2026-04-02 7/350 2026-04-04 19:16 by 蓝云思雨
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦! 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
信息提示
请填处理意见