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

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


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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習經(jīng)歷 +8 | 材料學257求調(diào)劑 2026-03-23 | 9/450 |
|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
[考研] 生物學071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|
|
[考研] 生物學調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 化學求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 材料學求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 本科鄭州大學物理學院,一志愿華科070200學碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|