24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3411  |  回复: 7

scpchaos

新虫 (初入文坛)

[求助] VASP引入空穴以及电子的一些问题

关于VASP中如何引入空穴和电子论坛上有过讨论:
http://muchong.com/html/201106/3322747.html
但我还是有一些问题:
如果我要研究Au/TiO2 的光催化,那么我需要引入一个空穴和一个电子。如果我只关心电子,用NELECT来引入电子的话,根据
http://cms.mpi.univie.ac.at/vasp ... LDIPOL_EPSILON.html
Total energies from charged slab calculations are hence useless, and can not be used to determine relative energies.
如果是这样的话,原文为何还能计算能量呢?

如果我希望引入的电子停留在Au上的话,有什么办法呢?
谢谢!
回复此楼

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

Photochemistry

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

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
scpchaos: 金币+10, ★★★很有帮助 2014-01-17 04:06:11
sunyang1988: 金币+2, 谢谢交流 2014-01-17 13:54:41
scpchaos: 金币+20, ★★★★★最佳答案 2014-01-18 07:50:03
只是说vasp里面没有方法修正能量charged slab得能量 有些其他方法可以的 但是都比较难弄 需要比较好的理论和编程基础 楼主可以看看这两片文章:
PRL 102 016402
PRL 110 095505
尤其第二篇讲了表面的电荷修正

还有一种方法就是用cluster model 没有周期性边界条件的 所以可以随便加charge 但是vasp做不了 像FHI-aims bigdft之类的可以 用cluster model的问题就是边界上的原子的性质和表面上的是不一样的 也会造成很多问题
2楼2014-01-17 00:54:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scpchaos

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by gump_813 at 2014-01-17 00:54:42
只是说vasp里面没有方法修正能量charged slab得能量 有些其他方法可以的 但是都比较难弄 需要比较好的理论和编程基础 楼主可以看看这两片文章:
PRL 102 016402
PRL 110 095505
尤其第二篇讲了表面的电荷修正

...

非常感谢回复~~
不知道是否可以这么认为:可以通过charged slab 来引入hole,但是这样计算出来的能量是不对的,需要correction?
如果是这样的话,那之前提到的那篇算co2- 在tio2上,并没有correction,是不是认为这样计算得到的结果不准确呢?
J. Phys. Chem. C 2010, 114, 21474–21481

在VASP 的FAQ里有这样一段:
Q:I adsorb, an ionic species e.g. O$ ^-$ on an insulating surface. To select a specific charge state, I have increased the number of electrons by one compared to the neutral system. Now, I have no clue how to evalute the total energy properly (i.e. are there convergence corrections).
A:Actually, you MUST NOT set the number of electrons manually for a slab calculation. I.e., when you calculate the slab-O$ ^-$ system you are not allowed to select a specific charge state for the oxygen ion, by increasing the number of electrons manually. Specific charge state calculations make sense only in 3D systems and for cluster calculations.

If you conduct the calculations properly, i.e. if your slab is large enough and the lateral dimension (x,y) of your surface is large enough the energy should converge to the proper value, i.e. the O should acquire the correct charge state automatically.

Reason: If you set the number of electrons in the INCAR file for a slab calculation you end up with a charged slab. The electrostatic energy of such a slab is however only conditionally convergent and worse, in practice, even infinite (BASIC, BASIC ELECTROSTATICS). Therefore, no method whatsoever exists to correct the error in the electrostatic energy. E.g. the energy converges towards infinity, when the vacuum width is increased. You can try to validate this, by simply increasing the vacuum width in VASP for a charged slab. You will find that the energy increases or decreases linearly with the vacuum width.

Well, there is maybe one method that can surmount the aforementioned problem. You can charge the slab and increase systematically the distance between the O- species (by increasing the lateral dimensions of your supercell) at a fixed vacuum width, and finally extrapolate the energies towards infinite lateral distances. The energy should converge towards the correct value as 1/$ d$, where $ d$ is the distance between the adsorbed species. This might yield a converged value. The point is that, as I mentioned above, the electrostatic energy is only conditionally convergent for the case of a charged slab/system, and results depend on how you evaluate the limit towards infinity. However, to the best of my knowledge, this has not been done or attempted hereto (and therefore we can not assist you on that issue).

这样的话,如果我只是要计算Au/TiO2 中的charge transfer应该可以,但是如果在光催化下,还有什么方法可以加入electron 和hole 呢?
谢谢!
3楼2014-01-17 04:05:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813

金虫 (小有名气)

★ ★
sunyang1988: 金币+2, 谢谢交流 2014-01-17 13:54:53
引用回帖:
3楼: Originally posted by scpchaos at 2014-01-17 04:05:35
非常感谢回复~~
不知道是否可以这么认为:可以通过charged slab 来引入hole,但是这样计算出来的能量是不对的,需要correction?
如果是这样的话,那之前提到的那篇算co2- 在tio2上,并没有correction,是不是认 ...

不是那么容易做的 比如说可以引入一个碱金属 因为电负性很低 所以electron很容易出来 总体系还是保持中性的 (hole的话可以引入一个Cl?)不过这样做也会引入其他的问题 而且电子的transfer到什么地方怎么就能保证和实际体系一样呢? 不过倒是可以试一试 这个总的来说是很难做的 哈哈
可能也又其他方法吧 再查查文献吧! 我看到有人就只做一个au cluster 加电子

你说的那篇jpcc 我没看 但是如果就加个电子就说自己是有CO2-的话 而且是周期性便捷条件 应该是不对的 没看过 不知道到底是什么
4楼2014-01-17 05:54:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813

金虫 (小有名气)


sunyang1988: 金币+1, 谢谢交流 2014-01-17 13:55:03
引用回帖:
3楼: Originally posted by scpchaos at 2014-01-17 04:05:35
非常感谢回复~~
不知道是否可以这么认为:可以通过charged slab 来引入hole,但是这样计算出来的能量是不对的,需要correction?
如果是这样的话,那之前提到的那篇算co2- 在tio2上,并没有correction,是不是认 ...

刚看了一眼 他们那个是用cluster model做的阿 不是周期的 这个加电荷是没问题的 至于是不是CO2-就要电荷分析了
5楼2014-01-17 05:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scpchaos

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by gump_813 at 2014-01-17 05:54:58
不是那么容易做的 比如说可以引入一个碱金属 因为电负性很低 所以electron很容易出来 总体系还是保持中性的 (hole的话可以引入一个Cl?)不过这样做也会引入其他的问题 而且电子的transfer到什么地方怎么就能保证 ...

嗯 看了好多文献感觉这类问题现阶段的package只能qualitative的算些东西~~
准备用cluster model算下,看看结果和cluster size是不是很关系很大
6楼2014-01-18 07:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

送红花一朵
本帖仅楼主可见
7楼2018-10-08 08:35:21
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)

本帖仅楼主可见
8楼2018-10-08 08:35:56
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 scpchaos 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂求收留 +33 果然有我 2026-04-10 34/1700 2026-04-15 00:37 by zzzggc
[考研] 297工科调剂? +13 河南农业大学-能 2026-04-13 13/650 2026-04-14 16:46 by Art1977
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[考研] 求调剂 +20 MAX怅惘 2026-04-09 22/1100 2026-04-14 14:57 by 独醉梦孤城
[考研] 机械工程313分找工科调剂 +4 双一流本科机械 2026-04-08 4/200 2026-04-14 07:32 by Abskk
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 考研二轮调剂 +10 故人?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 280求调剂 +7 兮兮夜夜 2026-04-09 10/500 2026-04-12 00:33 by 蓝云思雨
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 282,电气工程专业,求调剂,不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
信息提示
请填处理意见