24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 3030  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 288环境专硕,求调材料方向 +35 lllllos 2026-04-04 39/1950 2026-04-07 23:24 by 一只好果子?
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 土木水利专硕276分求调剂 +4 我想上学!!6 2026-04-05 7/350 2026-04-07 22:10 by Delta2012
[考研] 材料专硕(0856) 339分求调剂 +14 哈哈哈鹅哈哈哈 2026-04-05 14/700 2026-04-07 22:02 by lijunpoly
[考研] 085602调剂 初试总分335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
[考研] 280求调剂 +7 李rien 2026-04-04 7/350 2026-04-07 17:33 by 蓝云思雨
[考研] 22408 调剂材料 +7 我叫ez 2026-04-06 8/400 2026-04-07 17:12 by 蓝云思雨
[考研] 一志愿国科大信工所,英二数二408总分293分求调剂 +3 ilcyuan 2026-04-02 4/200 2026-04-06 16:35 by likeihood
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 调剂 +3 好好读书。 2026-04-02 3/150 2026-04-05 13:02 by arrow8852
[考研] 085601,一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 085600,320分求调剂 +6 大馋小子 2026-04-02 6/300 2026-04-02 21:54 by dongzh2009
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 食品学硕362求调剂 +3 xuanxianxian 2026-04-01 3/150 2026-04-01 21:05 by 啊李999
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
信息提示
请填处理意见