| 5 | 1/1 | 返回列表 |
| 查看: 5035 | 回復(fù): 15 | ||||
| 本帖產(chǎn)生 1 個(gè) 1ST強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
vasp計(jì)算帶電粒子吸附計(jì)算的問題 已有3人參與
|
||||
|
最近想計(jì)算帶電粒子在催化劑表面的吸附性質(zhì),但是看了手冊(cè),和壇子里的介紹,發(fā)現(xiàn)對(duì)于slab模型,對(duì)整個(gè)體系加入電荷(NELECT控制),好像計(jì)算是有問題的.... 這個(gè)帖子有過說明http://www.gaoyang168.com/bbs/viewthread.php?tid=8266545&fpage=1&target=blank http://cms.mpi.univie.ac.at/vasp/vasp/FAQ.html最后五段 vasp官方FAQ中也有過解答。 問題: 1、針對(duì)slab模型,調(diào)節(jié)NELECT去結(jié)構(gòu)優(yōu)化,是不是完全沒有意義 對(duì)一個(gè)電中性條件下已優(yōu)化好的結(jié)構(gòu)(CONTCAR_1),其能量為E1,然后對(duì)CONTCAR_1加一個(gè)電子并結(jié)構(gòu)優(yōu)化得到帶電下的穩(wěn)定結(jié)構(gòu)(CONTCAR_2),計(jì)算CONTCAR_2在電中性下的能量E2,比較E1和E2來說明材料帶電對(duì)其結(jié)構(gòu)的影響。這樣做有沒有意義。對(duì)于官方教程,我的理解是帶電后,會(huì)產(chǎn)生一個(gè)靜電場,使得能量不準(zhǔn)。如果我不橫向比較能量(比較帶電前后),有沒有意義 2、如果給體系(slab)帶上一個(gè)電子,然后去吸附電中性的小分子,在這種條件下,所有的能量全是建立在帶一個(gè)電子的基礎(chǔ)上,這樣比較吸附的過程能量變化有沒有意義。 3、最近也看到過 對(duì)體系(slab)加電子模擬,在催化劑表面上超氧負(fù)離子O2-還原過程,這么做有沒有意義 4、仔細(xì)分析了一下啊,產(chǎn)生這些問題的原因都是調(diào)控NELECT不能精準(zhǔn)地把電子加到某個(gè)原子上,有沒有其他的軟件固定原子的所帶的電荷?如果有請(qǐng)告知一下。 PS:曾經(jīng)在壇子里看到有說改POTCAR中電子的數(shù)目,來調(diào)控帶電,仔細(xì)分析,這樣做確實(shí)可以調(diào)節(jié)外來帶電粒子A-(slab結(jié)構(gòu)總不含A元素)的電荷,并且把電荷局域到A上,但是POTCAR太深?yuàn)W,不敢改,也就沒有試過...有沒有哪位試過,或者看到過類似的做法,希望告知。 非計(jì)算出身,自學(xué)的vasp,對(duì)于理論基礎(chǔ)把握不到位,請(qǐng)大家指導(dǎo)。 1.jpg |
弛豫 |


專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |
|
VASP那個(gè)FAQ第五段似乎給出了方案。 1、不見得沒意義,以往無法處理過剩電荷原因是因?yàn)閱伟麅?nèi)帶有電荷會(huì)導(dǎo)致庫倫發(fā)散,F(xiàn)在不少程序允許使用凝膠模型,引入背景電荷與過剩電荷相互作用來克服發(fā)散的問題。如果不能夠使用凝膠模型,那應(yīng)該是無法進(jìn)行計(jì)算的。但即便如此,過剩電荷可能會(huì)引起體系的兩個(gè)slab之間出現(xiàn)靜電作用,比較合適的做法是足夠的真空層(最好也考慮work function的校正)。當(dāng)然,這樣的做法可能多少還是存在一些爭議。 2、你指的是Mol(-)+Slab=Mol-Slab(-)?這樣應(yīng)該還是可行的吧。3也類似,前提是都是要符合1中提到的條件。 4、這個(gè)問題必須搞清楚原理,化合價(jià)只是個(gè)經(jīng)典形式電荷的問題,在量子物理中,電子是離域的,強(qiáng)行使之局域在某個(gè)原子根本無意義(由于方法缺陷導(dǎo)致的該局域而不局域又是另外一碼事)。 不必調(diào)整POTCAR,直接指定NELECT結(jié)果原則上應(yīng)該沒什么不同吧。還望進(jìn)一步討論。 |


專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|