24小时热门版块排行榜    

查看: 4266  |  回复: 15
本帖产生 1 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

无所谓109

银虫 (正式写手)

[求助] vasp计算带电粒子吸附计算的问题已有3人参与

最近想计算带电粒子在催化剂表面的吸附性质,但是看了手册,和坛子里的介绍,发现对于slab模型,对整个体系加入电荷(NELECT控制),好像计算是有问题的....
这个帖子有过说明http://muchong.com/bbs/viewthread.php?tid=8266545&fpage=1&target=blank
http://cms.mpi.univie.ac.at/vasp/vasp/FAQ.html最后五段
vasp官方FAQ中也有过解答。

问题:
1、针对slab模型,调节NELECT去结构优化,是不是完全没有意义
      对一个电中性条件下已优化好的结构(CONTCAR_1),其能量为E1,然后对CONTCAR_1加一个电子并结构优化得到带电下的稳定结构(CONTCAR_2),计算CONTCAR_2在电中性下的能量E2,比较E1和E2来说明材料带电对其结构的影响。这样做有没有意义。对于官方教程,我的理解是带电后,会产生一个静电场,使得能量不准。如果我不横向比较能量(比较带电前后),有没有意义  

2、如果给体系(slab)带上一个电子,然后去吸附电中性的小分子,在这种条件下,所有的能量全是建立在带一个电子的基础上,这样比较吸附的过程能量变化有没有意义。

3、最近也看到过 对体系(slab)加电子模拟,在催化剂表面上超氧负离子O2-还原过程,这么做有没有意义

4、仔细分析了一下啊,产生这些问题的原因都是调控NELECT不能精准地把电子加到某个原子上,有没有其他的软件固定原子的所带的电荷?如果有请告知一下。
      PS:曾经在坛子里看到有说改POTCAR中电子的数目,来调控带电,仔细分析,这样做确实可以调节外来带电粒子A-(slab结构总不含A元素)的电荷,并且把电荷局域到A上,但是POTCAR太深奥,不敢改,也就没有试过...有没有哪位试过,或者看到过类似的做法,希望告知。

非计算出身,自学的vasp,对于理论基础把握不到位,请大家指导。

vasp计算带电粒子吸附计算的问题
1.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

弛豫

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

向金字塔顶攀登
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无所谓109

银虫 (正式写手)

引用回帖:
4楼: Originally posted by 卡开发发 at 2015-03-30 18:29:28
研究一下IDIPOL相关的关键字,偶极校正的信息论坛能找出来不少。...

好的,谢谢您!
向金字塔顶攀登
5楼2015-03-30 19:00:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
franch: 金币+2, 谢谢回帖交流。。。。 2015-03-31 10:18:55
franch: 1ST强帖+1, 应助指数+1, 鼓励交流。。。。 2015-03-31 10:22:58
无所谓109: 金币+20 2015-07-25 20:55:15
VASP那个FAQ第五段似乎给出了方案。
1、不见得没意义,以往无法处理过剩电荷原因是因为单胞内带有电荷会导致库伦发散。现在不少程序允许使用凝胶模型,引入背景电荷与过剩电荷相互作用来克服发散的问题。如果不能够使用凝胶模型,那应该是无法进行计算的。但即便如此,过剩电荷可能会引起体系的两个slab之间出现静电作用,比较合适的做法是足够的真空层(最好也考虑work function的校正)。当然,这样的做法可能多少还是存在一些争议。

2、你指的是Mol(-)+Slab=Mol-Slab(-)?这样应该还是可行的吧。3也类似,前提是都是要符合1中提到的条件。

4、这个问题必须搞清楚原理,化合价只是个经典形式电荷的问题,在量子物理中,电子是离域的,强行使之局域在某个原子根本无意义(由于方法缺陷导致的该局域而不局域又是另外一码事)。

不必调整POTCAR,直接指定NELECT结果原则上应该没什么不同吧。还望进一步讨论。

» 本帖已获得的红花(最新10朵)

不一定挂在论坛,计算问题问题欢迎留言。
2楼2015-03-29 23:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无所谓109

银虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by 卡开发发 at 2015-03-29 23:45:56
VASP那个FAQ第五段似乎给出了方案。
1、不见得没意义,以往无法处理过剩电荷原因是因为单胞内带有电荷会导致库伦发散。现在不少程序允许使用凝胶模型,引入背景电荷与过剩电荷相互作用来克服发散的问题。如果不能够 ...

谢谢您的解答,现在对这个问题弄明白了些。早上就分析了一下我的结果,把表面静电势导出来看了看。发现加了电子后,我体系两个slab之间确实存在比较明显的静电场的影响。

您提到的work function校正,能不能简要说说校正方法,之前没处理过work function校正,在坛子里也没有找到相关的计算方法。
vasp计算带电粒子吸附计算的问题-1
1.jpg


vasp计算带电粒子吸附计算的问题-2
12.jpg

向金字塔顶攀登
3楼2015-03-30 16:06:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by 无所谓109 at 2015-03-30 16:06:34
谢谢您的解答,现在对这个问题弄明白了些。早上就分析了一下我的结果,把表面静电势导出来看了看。发现加了电子后,我体系两个slab之间确实存在比较明显的静电场的影响。

您提到的work function校正,能不能简要 ...

研究一下IDIPOL相关的关键字,偶极校正的信息论坛能找出来不少。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2015-03-30 18:29:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 求助!请问海外博后依托内地单位申请青基的优劣? 100+3 2280999712 2024-11-15 9/450 2024-11-15 23:08 by dxcharlary
[教师之家] 这种人痛苦吗 +12 2671 2024-11-12 14/700 2024-11-15 22:15 by 仇恨毁灭财富
[论文投稿] 国内期刊审稿人数量 +3 新时代核动力驴 2024-11-13 5/250 2024-11-15 20:37 by 新时代核动力驴
[考博] 随缘读博 一篇中科院一区Top 一篇中科院二区Top 两篇专利 +6 Ricoch4t 2024-11-13 15/750 2024-11-15 19:24 by 青青之洲
[基金申请] 博后资助名单出来了 +12 Shxjjxjkx 2024-11-14 14/700 2024-11-15 14:40 by fanfenggui
[硕博家园] 实验进度 +9 麻辣烫。 2024-11-11 14/700 2024-11-15 13:14 by hanweiguo
[有机交流] 二甲胺的使用 20+5 太阳谷 2024-11-14 13/650 2024-11-15 13:11 by xiaomei1031
[硕博家园] reject后小感 +6 sioc-sunj 2024-11-14 8/400 2024-11-15 11:50 by 畅21
[硕博家园] 课题和专业一点都不相关怎么办 +7 imt4n9 2024-11-11 10/500 2024-11-15 01:55 by 成长虫
[硕博家园] 研究生的生活该是什么样 +4 lqy0719 2024-11-14 4/200 2024-11-14 16:45 by 阿荣喝酒
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[硕博家园] 大龄已婚想读博如何 +15 米娅阳 2024-11-11 18/900 2024-11-14 14:18 by xiaomi0401
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[基金申请] 浙江省青基怎么看A还是B? +8 ToallGoof 2024-11-11 14/700 2024-11-13 09:11 by Praybest
[论文投稿] 二审审一年的佛系编辑 +10 thefoxrain 2024-11-09 15/750 2024-11-12 19:27 by lide966
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见