| 查看: 3033 | 回复: 6 | |||
[求助]
vasp 计算单分子/原子 时晶胞大小的选取
|
|
我觉得，原则上讲应该是使相邻晶胞内分子/原子距离足够大，以消除计算过程中互相的影响。满足这一条件后，在相同的精度下，所得到的单分子/原子的能量大致应该是相当的。 但是本人在计算单个金属Ni时发现，使用10*10*10的晶胞以及15*15*15晶胞，其它参数一样。计算所得单原子的能量相差较大。 不知道大家平时计算类似的体系时，是如何设置晶胞大小的？ |
回复此楼» 收录本帖的淘帖专辑推荐
 化学电源之锂离子电池 |
» 猜你喜欢
机械工程264学硕求调剂
已经有3人回复
-
264求调剂
已经有9人回复
-
求调剂
已经有6人回复
-
22408 266求调剂
已经有7人回复
-
求调剂
已经有14人回复
-
273求调剂
已经有43人回复
-
307求调剂
已经有13人回复
-
327求调剂
已经有7人回复
-
338求调剂
已经有6人回复
-
326分，一志愿沪9，求生物学调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用VASP计算H原子的振动频率
已经有16人回复
- 态密度计算 vasp 4.6 计算Zn原子的PDOS 关于Zn（3d104s2）的4p对价带的贡献
已经有4人回复
- 怎么用vasp的计算结果分析分子中各原子带电多少情况
已经有6人回复
- 画VASP计算的原子结构图的一个方法（自己瞎想的，如和他人方法雷同请见谅）
已经有4人回复
- vasp计算界面结构时原子弛豫的问题
已经有12人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- VASP 四方晶系（TiO2 anatase ）单胞建模 原子数的问题 麻烦大家指导
已经有19人回复
- 在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？
已经有8人回复
- vasp可以固定单个原子的磁矩吗？
已经有5人回复
- 如何用VASP计算体系中每一类原子的DOS？
已经有10人回复
- vasp中如何构建超晶胞
已经有3人回复
- 求助：vasp 计算氢原子能量是-1.16ev 而不是 -13.6ev，why ?
已经有18人回复
- 【求助】如何显示晶胞结构中的原子编号与空间坐标？
已经有10人回复
- 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？
已经有5人回复
- 【求助】vasp为什么采用周期性边界条件（或超原胞模型）来处理原子分子表面体系等，有
已经有8人回复
- 【求助】能否发一份用vasp计算超晶胞某个表面的输入文件给我
已经有7人回复
- 【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
- 【求助】关于锂原子vasp计算后的1s2s2p分波态密度处理
已经有12人回复
emilyoyang
木虫 (正式写手)
- 应助: 232 (大学生)
- 金币: 4340.4
- 红花: 42
- 帖子: 667
- 在线: 189.9小时
- 虫号: 1630291
- 注册: 2012-02-20
- 性别: MM
- 专业: 理论和计算化学
2楼2012-12-30 20:19:41
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼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
5楼2013-01-01 10:51:27
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
sunyang1988: 金币+1, 谢谢交流 2013-01-03 12:10:18
sunyang1988: 金币+1, 谢谢交流 2013-01-03 12:10:18
|
我自己算了一下，结果如下： 15*15*15 -1.840572eV 10*10*10 -1.844656eV INCAR： general: SYSTEM = general ENCUT = 400 ISTART = 0 ; ICHARG = 2 ISMEAR = 0 ; SIGMA = 0.05 IALGO = 38 NELM = 200 #PREC = Accurate EDIFF = 1E-5 ISPIN=2 MAGMOM=3.0 Write flags LWAVE = .FALSE. LCHARG = .FALSE. Cluster NPAR = 8 ; NSIM = 4 LPLANE = .TRUE. LSCALU = .FALSE. |
6楼2013-01-03 10:24:41
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2013-01-03 10:36:48
