24小时热门版块排行榜    

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

倔强的小蚊虫

铜虫 (正式写手)

[求助] 如何用VASP计算单个He原子的能量? 已有4人参与

本人系初学者,希望得到前辈的指点。经过了无数次的尝试,都没有得到合适的值,想一想还是关键词设置的问题。现在把输入文件内容贴如下:
INCAR:
SYSTEM = He

LWAVE=.FALSE.
LCHARG =.FALSE.

ISTART = 0
ICHARG = 2
IALGO = 48

LPLANE = .TRUE.
NPAR   = 16
LSCALU = .FALSE.
NSIM   = 4

PREC = ACCURATE
NELMDL=5

#ISPIN = 1
LREAL = Auto

ISMEAR = 0
SIGMA = 0.02
ENCUT = 650


POSCAR:
He
1.0
        10         0.0000000000         0.0000000000
        0          10                   0.0000000000
        0          0                    10
    He
    1
Direct
     0.000000000         0.000000000         0.000000000


KPOINTS:
AUTO GRID
0
M
21 21 21
0 0 0


POTCAR(贴了部分进去):
PAW_PBE He 05Jan2001
2.00000000000000000
parameters from PSCTR are:
   VRHFIN =He: 1s
   LEXCH  = PE
   EATOM  =    78.7251 eV,    5.7861 Ry

   TITEL  = PAW_PBE He 05Jan2001
   LULTRA =        F    use ultrasoft PP ?
   IUNSCR =        0    unscreen: 0-lin 1-nonlin 2-no
   RPACOR =     .000    partial core radius
   POMASS =    4.000; ZVAL   =    2.000    mass and valenz
   RCORE  =    1.100    outmost cutoff radius
   RWIGS  =    1.200; RWIGS  =     .635    wigner-seitz radius (au A)
   ENMAX  =  478.896; ENMIN  =  359.172 eV
   RCLOC  =     .703    cutoff for local pot
   LCOR   =        T    correct aug charges
   LPAW   =        T    paw PP
   EAUG   = 1180.948
   RMAX   =    2.107    core radius for proj-oper
   RAUG   =    1.200    factor for augmentation sphere
   RDEP   =    1.110    radius for radial grids
   QCUT   =   -5.933; QGAM   =   11.866    optimization parameters

……………………………………………………

以及一个脚本文件VASP.PBS。
回复此楼

» 猜你喜欢

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

自制力,创造力,行动力,三力超群,才能卓越。
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lqcata

铜虫 (正式写手)

【答案】应助回帖

As to the He element, ENMAX  =  478.896 , you do not need so high cut off energy for ECUT=650,
By the way, SIGMA = 0.01 is good for moleculars
用知识与勤劳成功自己
3楼2014-03-11 20:29:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

lqcata

铜虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
For single atom, gamma point is enough!
For the box, you can use 7*8*9 to break the symmetry.
For more detail information, you can refer to the manual of vasp.
https://www.vasp.at/vasp-workshop/slides/handsonI.pdf

Be carefull: Do not learn some bad habits when you start to learn how to do calculations!
用知识与勤劳成功自己
2楼2014-03-11 20:02:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
ISPIN=2,需要自旋计算。
4楼2014-03-11 23:12:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by jmlv at 2014-03-11 23:12:16
ISPIN=2,需要自旋计算。

我看错了,你计算的是He(氦)原子,无需自选计算。
我再次看了一下你的INCAR文件,貌似没问题。你说的没得到合适的值是什么意思?
5楼2014-03-11 23:23:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见