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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于用VASP计算功函数的疑问
51/1返回列表
查看: 4056  |  回复: 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

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
小小鑫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的回帖
查看全部 9 个回答

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

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的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭