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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于用VASP计算功函数的疑问
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 4256  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +3 种圣赐 2026-03-29 3/150 2026-03-29 12:06 by longlotian
[考研] 2026年华南师范大学欢迎化学，化工，生物，生医工等专业优秀学子加入！ +3 llss0711 2026-03-28 6/300 2026-03-29 10:26 by llss0711
[考研] 321求调剂 +7 璞玉~~ 2026-03-25 8/400 2026-03-29 06:41 by 544594351
[考研] 一志愿华理，数一英一285求A区调剂 +8 AZMK 2026-03-25 12/600 2026-03-28 18:15 by AZMK
[考研] 347求调剂 +3 山顶见α 2026-03-25 3/150 2026-03-28 14:13 by 唐沐儿
[考研] 291求调剂 +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[考研] 0856，材料与化工321分求调剂 +12 大馋小子 2026-03-27 13/650 2026-03-28 10:56 by self2008
[考研] 339求调剂，想调回江苏 +6 烤麦芽 2026-03-27 8/400 2026-03-28 10:40 by 烤麦芽
[考研] 291求调剂 +7 孅華 2026-03-22 7/350 2026-03-28 04:02 by fmesaito
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 333求调剂 +3 question挽风 2026-03-23 3/150 2026-03-27 11:29 by 不吃魚的貓
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
[考研] 一志愿郑大085600，310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 292求调剂 +4 求求了收下我吧� 2026-03-26 4/200 2026-03-27 10:37 by zhshch
[考研] 276求调剂。有半年电池和半年高分子实习经历 +10 材料学257求调剂 2026-03-23 11/550 2026-03-27 10:13 by YCIT- LHL
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 0703化学求调剂 +3 丹青奶盖 2026-03-26 5/250 2026-03-26 20:11 by macy2011
[考研] 081700 调剂 267分 +11 迷人的哈哈 2026-03-23 11/550 2026-03-26 15:41 by zzll406
[考研] 考研一志愿苏州大学初始315（英一）求调剂 +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
信息提示
关闭
请填处理意见
关闭