24小时热门版块排行榜    

查看: 2995  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天大材料与化工(085600)总分338 +3 蔡大美女 2026-03-13 3/150 2026-03-18 15:53 by jhhcooi
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +7 rare12345 2026-03-18 7/350 2026-03-18 14:31 by laoshidan
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +9 yangfz 2026-03-17 9/450 2026-03-18 12:38 by 尽舜尧1
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 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%、
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见