| 查看: 5097 | 回复: 15 | ||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
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，对于理论基础把握不到位，请大家指导。  1.jpg |
回复此楼» 收录本帖的淘帖专辑推荐
 弛豫 |
» 猜你喜欢
085602 307分 求调剂
已经有3人回复
-
340求调剂
已经有4人回复
-
085602 化工专硕 338分 求调剂
已经有8人回复
-
328求调剂
已经有6人回复
-
求调剂
已经有6人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有8人回复
-
0703化学
已经有4人回复
-
289求调剂
已经有6人回复
-
291求调剂
已经有17人回复
-
材料与化工085600，总分304，本科有两篇sci参与，求调剂
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 采用VASP计算得到的禁带宽度为什么比实验值低
已经有11人回复
- 吸附计算如何设计？
已经有4人回复
- 感觉VASP计算缺陷存在很大问题
已经有17人回复
- 关于Vasp TST Tools的NEB计算中的SS NEB计算问题
已经有7人回复
- VASP中断计算之后，续算的准确性问题？
已经有14人回复
- VASP 中离子的能量如何计算
已经有12人回复
- 在用vasp进行表面吸附时，遇到的问题，求高手解惑！
已经有12人回复
- VASP计算H原子 吸附在Pt(111)表面时，最后CONTCAR 里少了H 原子
已经有3人回复
- 关于VASP计算收敛的判断
已经有25人回复
- vasp表面和吸附分子的结构分别优化后，如何将它们放到一起
已经有11人回复
- VASP 中NELECT 的设置
已经有8人回复
- VASP计算表面吸附 但是能量总是不收敛/。。。
已经有9人回复
- vasp 计算 partial charge 时k点选取的问题
已经有19人回复
- vasp 计算 dimer的问题
已经有4人回复
- vasp计算频率后，如何查看OUTCAR中的结果？
已经有8人回复
- vasp计算过渡态问题
已经有13人回复
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- VASP计算能带时是否读取WAVECAR
已经有16人回复
- VASP计算：LORBIT设置出现问题
已经有5人回复
- vasp计算Ag在graphene上的吸附能
已经有7人回复
- vasp计算磁矩的困惑
已经有10人回复
- VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！
已经有23人回复
- 【求助】如何使用VASP计算charged defect的形成焓
已经有12人回复
- 【求助完毕】集群上用vasp计算ELF有ELFCAR而没有DOSCAR是啥问题？
已经有15人回复
kjsun
木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 5557.5
- 散金: 100
- 红花: 5
- 帖子: 434
- 在线: 469.8小时
- 虫号: 285070
- 注册: 2006-10-14
- 性别: GG
- 专业: 催化化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 辛苦 2015-04-03 07:47:59
无所谓109: 金币+10 2015-07-25 20:55:30
感谢参与，应助指数 +1
liliangfang: 金币+2, 辛苦 2015-04-03 07:47:59
无所谓109: 金币+10 2015-07-25 20:55:30
|
我最近也在研究带电粒子的计算问题。 我发现VASP加电荷不是在体系的最低空轨道加入一个电子，而是在整个空间加入一个背景电荷。（说明书上有相信加入公式） VASP这么干目的肯定是为了解决一些问题，但是具体什么，我完全不清楚。 但是我测试的结果表明要想用VASP计算带电系统，必须非常非常慎重。 我测试了很多结构，最简单的如 CO2- (就是CO2+1额外电子） CO2-的结构是非直线型结构，但是VASP优化的结果是直线型。这说明VASP对带电体系的结构优化是不对的。 另外，我的测试结果表明：VASP对带电体系的能量的计算结果也是不对的，至少我测试的那些体系，偏差都极大，完全超出了我的理解范畴。 相应的Gaussian的测试结果则是相对可靠。 |
6楼2015-04-01 13:50:45
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3514.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
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朵）
无所谓1092楼2015-03-29 23:45:56
3楼2015-03-30 16:06:34
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3514.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2015-03-30 18:29:28
