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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于用VASP计算功函数的疑问
91/1返回列表
查看: 4573  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 20:36 by xxniao123
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
关闭
请填处理意见
关闭