24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4032  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖

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的回帖

前田庆次

木虫 (著名写手)

【答案】应助回帖

NPAR   = 16 难道你有256个核心???!!!太强悍了!羡慕嫉妒恨一下
个人有几个建议:
1. SIGMA最好取-5
2. 加上ISPIN=2
3. 我一般都是用Gamma的,不知道您的爱好如何
个人拙见
长枪一横花飘零,松风追月伴我行。无双人间世如梦,倾奇万世永留名。
17楼2014-03-15 10:20:24
已阅   关注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的回帖

倔强的小蚊虫

铜虫 (正式写手)

引用回帖:
5楼: Originally posted by jmlv at 2014-03-11 23:23:43
我看错了,你计算的是He(氦)原子,无需自选计算。
我再次看了一下你的INCAR文件,貌似没问题。你说的没得到合适的值是什么意思?...

我目前在计算钨基体中He的缺陷形成能。根据相关公式,氦的哪一项不应该太小(每次计算都是-0.000XXXX),这样的话,对整个缺陷形成能来说,氦的影响就不大,甚至可以忽略,所以这样的结果是不合适的。我想的结果应该是小数点后一位或者两位都可以接受。
自制力,创造力,行动力,三力超群,才能卓越。
6楼2014-03-13 13:21:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

倔强的小蚊虫

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by lqcata at 2014-03-11 20:02:05
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/ ...

谢谢你的应助。
自制力,创造力,行动力,三力超群,才能卓越。
7楼2014-03-13 13:22:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)

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

请问大金虫,我如果算一个H原子,加入自旋的话,只需要写入ISPIN=2吗,MAGMOM和NUPDOWN,难道不用设置吗。。。
慢慢的思考。。。
8楼2014-03-13 16:36:17
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

【答案】应助回帖

引用回帖:
8楼: Originally posted by 刘小伟110 at 2014-03-13 16:36:17
请问大金虫,我如果算一个H原子,加入自旋的话,只需要写入ISPIN=2吗,MAGMOM和NUPDOWN,难道不用设置吗。。。...

计算单个H原子,只需将ISPIN设为2即可,VASP会给出MAGMOM。
9楼2014-03-14 00:06:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

引用回帖:
6楼: Originally posted by 倔强的小蚊虫 at 2014-03-13 13:21:37
我目前在计算钨基体中He的缺陷形成能。根据相关公式,氦的哪一项不应该太小(每次计算都是-0.000XXXX),这样的话,对整个缺陷形成能来说,氦的影响就不大,甚至可以忽略,所以这样的结果是不合适的。我想的结果应该 ...

我用VASP.5.2.2计算He原子的总能量是-0.0003806 eV。这个值应该是对的,因为He原子是电子填满壳层的原子且无极化,VASP就该给出接近于参考值的总能量。
10楼2014-03-14 00:45:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 倔强的小蚊虫 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +11 新时代材料 2026-03-27 11/550 2026-03-28 18:09 by YK。。
[考研] 调剂310 +6 温柔的晚安 2026-03-25 7/350 2026-03-28 16:04 by qq2000
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +6 芋泥小铃铛 2026-03-28 6/300 2026-03-28 14:15 by Ainin_
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +8 辰沐5211314 2026-03-26 8/400 2026-03-28 11:25 by xuxiang
[考研] 085602 307分 求调剂 +7 不知道叫什么! 2026-03-26 7/350 2026-03-28 09:57 by 神马都不懂
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂 +3 wywy.. 2026-03-27 4/200 2026-03-28 03:32 by fmesaito
[考研] 一志愿上海理工能源动力(085800)310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 279 分 求调剂 +4 睡个好觉_16 2026-03-24 4/200 2026-03-27 15:05 by 醉在风里
[考研] 085600,材料与化工321分,求调剂 +9 大馋小子 2026-03-27 9/450 2026-03-27 14:30 by mmm just
[考研] 334求调剂 +3 雨清天晴 2026-03-21 3/150 2026-03-27 11:30 by 不吃魚的貓
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 341求调剂 +7 青柠檬1 2026-03-26 7/350 2026-03-27 00:19 by wxiongid
[考研] 321求调剂 +6 Ymlll 2026-03-24 6/300 2026-03-26 20:50 by 不吃魚的貓
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
信息提示
请填处理意见