应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于用VASP计算功函数的疑问
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 4257  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小小鑫000

铁虫 (小有名气)

[求助] 关于用VASP计算功函数的疑问 已有1人参与

各位前辈：
大家好！我首先计算了clean Ag(111)表面的功函数，也用同样的算法计算了单个K原子吸附在Ag（111）表面后，Ag表面的功函，发现二者的值基本相等，也就是说K的加入对表面功函几乎没影响，但是看文献，发现文献上clean Ag(111)表面的功函数的数值与我的基本相同（说明我的计算方法没问题），但是文献上加了K以后Ag（111）表面的功函数值明显减少了很多。因此，我猜想是不是在金属表面加了吸附物以后，某些计算方法发生了变化，还是别的什么原因？有做过类似反应的前辈，希望解惑！不胜感激！
回复此楼

» 猜你喜欢
好好研究，天天向上！
1楼 2018-08-23 11:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
小小鑫000: 金币+5 2018-08-24 08:39:23
小红豆: 金币+5, 回答帖子 2018-08-26 12:57:33
吸附K的Ag(111)表面积是多大，也就是说你考虑的K吸附度是多少？吸附K之后，功函数应该减少。肯定是你的K/Ag(111)计算输入文件某个地方设置错了。贴出POSCAR和INCAR以及OUTCAR。
赞一下 回复此楼
2楼2018-08-23 21:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by valenhou001 at 2018-08-23 21:47:59
吸附K的Ag(111)表面积是多大，也就是说你考虑的K吸附度是多少？吸附K之后，功函数应该减少。肯定是你的K/Ag(111)计算输入文件某个地方设置错了。贴出POSCAR和INCAR以及OUTCAR。

多谢前辈的指导，我的是在3x3x3的Ag（111）表面上吸附一个K原子。OOSCAR,INCAR和KPOINTS贴在下面了，OUTCAR太大了，在上传文件里。我最后计算出来的值是4.4，文献里是2.5到3之间。
我的POSCAR为：
Ag-K-1
   1.00000000000000
     8.8670997620000005    0.0000000000000000    0.0000000000000000
    -4.4335498810000002    7.6791336518000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000   17.2399997710999990
   Ag   K
    36     1
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000   F   F   F
-0.0005481086812443  0.0002945031603394  0.4175492366596158   T   T   T
  0.2224693444537190  0.1114479022138608  0.2790869698130677   T   T   T
  0.1111100019999967  0.2222200040000004  0.1399800029999980   F   F   F
  0.3333300050000005  0.0000000000000000  0.0000000000000000   F   F   F
  0.3341327414374464  0.0001280312862744  0.4175463669582398   T   T   T
  0.5553160205813739  0.1112440984261727  0.2780573198840947   T   T   T
  0.4444400069999972  0.2222200040000004  0.1399800029999980   F   F   F
  0.6666700239999983  0.0000000000000000  0.0000000000000000   F   F   F
  0.6673032516525320  0.0013079748918270  0.4152925386186909   T   T   T
  0.8893146928716145  0.1112679000560649  0.2780453264638525   T   T   T
  0.7777799959999996  0.2222200040000004  0.1399800029999980   F   F   F
  0.0000000000000000  0.3333300050000005  0.0000000000000000   F   F   F
  0.0003925343029520  0.3340696476839030  0.4164943124703328   T   T   T
  0.2221918001795668  0.4442549078529263  0.2781175516406809   T   T   T
  0.1111100019999967  0.5555599930000028  0.1399800029999980   F   F   F
  0.3333300050000005  0.3333300050000005  0.0000000000000000   F   F   F
  0.3340775428172547  0.3347727721386630  0.4174089661883459   T   T   T
  0.5554867366757332  0.4443452076217984  0.2781780990472212   T   T   T
  0.4444400069999972  0.5555599930000028  0.1399800029999980   F   F   F
  0.6666700239999983  0.3333300050000005  0.0000000000000000   F   F   F
  0.6669443419471347  0.3340591778941079  0.4164930981790806   T   T   T
  0.8891001515026636  0.4449081052628588  0.2783154894142121   T   T   T
  0.7777799959999996  0.5555599930000028  0.1399800029999980   F   F   F
  0.0000000000000000  0.6666700239999983  0.0000000000000000   F   F   F
  0.0001272820991698  0.6670793815056536  0.4169189999347203   T   T   T
  0.2223020531412813  0.7784377144467273  0.2783328997846240   T   T   T
  0.1111100019999967  0.8888900279999987  0.1399800029999980   F   F   F
  0.3333300050000005  0.6666700239999983  0.0000000000000000   F   F   F
  0.3326354675994421  0.6664957561090066  0.4154936342111973   T   T   T
  0.5555505813843156  0.7776899831079608  0.2772388255386020   T   T   T
  0.4444400069999972  0.8888900279999987  0.1399800029999980   F   F   F
  0.6666700239999983  0.6666700239999983  0.0000000000000000   F   F   F
  0.6670675078105478  0.6667129740449578  0.4156738581628971   T   T   T
  0.8894851509113517  0.7785105442632916  0.2784068013330614   T   T   T
  0.7777799959999996  0.8888900279999987  0.1399800029999980   F   F   F
  0.5547645634826891  0.7803550920131926  0.5854364333345364   T   T   T

INCAR为：
SYSTEM = metal                  
ENCUT = 400                     
ISTART = 0;ICHARG = 2            
ISMEAR = 0; SIGMA = 0.2         
IBRION = -1                     
ALGO = fast                     
ISPIN = 1                        
LREAL = Auto                     
LVHAR = .TRUE.                  
SYMPREC = 0.001                  
IDIPOL = 3
LDIPOL = .TURE.  

KPOINTS为：
Automatic generation
0
Monhkorst-Pack
3 3 1
0.0 0.0 0.0
赞一下(1人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : OUTCAR
    • 2018-08-24 08:48:46, 1.07 M
好好研究，天天向上！
3楼2018-08-24 08:49:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
小小鑫000: 金币+5 2018-08-25 15:00:46
引用回帖:
3楼: Originally posted by 小小鑫000 at 2018-08-24 09:49:00
多谢前辈的指导，我的是在3x3x3的Ag（111）表面上吸附一个K原子。OOSCAR,INCAR和KPOINTS贴在下面了，OUTCAR太大了，在上传文件里。我最后计算出来的值是4.4，文献里是2.5到3之间。
我的POSCAR为：
Ag-K-1
   1. ...

存在几个问题：
1.真空层厚度只有10 Angstrom左右，略小点。适当增加到12到15angstrom。
2. k-mesh有点过小，只有在第一布里渊区只有5个k点，对金属体系，这个有点过少。适度地增加，至少5x5x1以上。
3. ISMEAR推荐为-5。

另外你在对LOCPOT处理时，看真空层区域的静电势是不是趋于平坦。

从你的OUTCAR来看Fermi能级是E-fermi :  -1.5835。
赞一下(1人) 回复此楼
4楼2018-08-24 19:04:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1437570483

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by valenhou001 at 2018-08-24 19:04:58
存在几个问题：
1.真空层厚度只有10 Angstrom左右，略小点。适当增加到12到15angstrom。
2. k-mesh有点过小，只有在第一布里渊区只有5个k点，对金属体系，这个有点过少。适度地增加，至少5x5x1以上。
3. ISMEAR ...

大神可以加个qq吗，我最近也在做表面吸附的第一性原理计算，想请教请教
赞一下 回复此楼
5楼2018-11-14 10:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

surelv

新虫 (小有名气)
可能需要加力吧
回复此楼
6楼2019-03-05 10:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Wq一路向北

禁虫 (正式写手)
本帖内容被屏蔽
7楼2019-05-30 18:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

希望能上up

新虫 (小有名气)
引用回帖:
3楼: Originally posted by 小小鑫000 at 2018-08-24 08:49:00
多谢前辈的指导，我的是在3x3x3的Ag（111）表面上吸附一个K原子。OOSCAR,INCAR和KPOINTS贴在下面了，OUTCAR太大了，在上传文件里。我最后计算出来的值是4.4，文献里是2.5到3之间。
我的POSCAR为：
Ag-K-1
   1. ...

能请教一下吗？我也遇到这种情况了，clean Pt表面和文献相差不大，但是吸附CO和乙醇之后发现并没有任何变化，不知道是什么原因？您可以帮助一下吗？
赞一下 回复此楼
8楼2019-08-01 15:19:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

弗伦德里希

新虫 (初入文坛)
楼主最后调整了哪些参数，计算结果与文献吻合？有做K的覆盖度 功函计算吗？
赞一下 回复此楼
9楼2022-09-27 11:12:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小鑫000 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +10 辰沐5211314 2026-03-26 10/500 2026-03-29 21:01 by 唐沐儿
[考研] 一志愿武汉理工，总分321，英一数二，求老师收留。 +11 nnnnnnn5 2026-03-25 11/550 2026-03-29 20:42 by 无际的草原
[考研] 085600 材料与化工 329分求调剂 +14 Mr. Z 2026-03-25 14/700 2026-03-29 17:27 by 学员8dgXkO
[考研] 349求调剂 +6 李木子啊哈哈 2026-03-25 6/300 2026-03-29 12:47 by 无际的草原
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 071000生物学求调剂，初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 308求调剂 +7 墨墨漠 2026-03-27 7/350 2026-03-28 07:43 by 热情沙漠
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 086000调剂 +3 7901117076 2026-03-26 3/150 2026-03-27 21:34 by Jianing_Mi
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 化学调剂 +4 爱吃番茄的旭 2026-03-24 5/250 2026-03-27 17:50 by kiokin
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 12/600 2026-03-27 14:43 by shangxh
[考研] 一志愿华东理工大学081700，初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 一志愿陕师大生物学071000，298分，求调剂 +5 SYA！ 2026-03-23 5/250 2026-03-27 09:29 by 不吃魚的貓
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
信息提示
关闭
请填处理意见
关闭