cenwanglai

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

[交流] 【求助完畢】硫酸根離子的結(jié)構(gòu)優(yōu)化，偶極修正及勢能修正。已有2人參與

SO4放在1nm邊長的立方盒子中間,中性分子簇價電子數(shù)為30.

對于離子的計算，按照結(jié)構(gòu)優(yōu)化和自洽計算兩個階段是否帶電可分為一下四種情況

帶電（NELECT=32）
不帶電(NELECT=30)

      結(jié)構(gòu)優(yōu)化       靜態(tài)自洽計算
A)    帶電                帶電             E1=-35.994eV
B）不帶電          帶電             E2=-35.765eV
C）帶電             不帶電          --
D）  不帶電          不帶電          E4=-26.539eV

以上結(jié)果說明，硫酸根是可以穩(wěn)定單獨(dú)存在的。而且，結(jié)構(gòu)優(yōu)化是否帶電對體系基態(tài)是有顯著影響的，二者能量差|E1-E2|=0.23eV.

但是這里沒有加偶極修正。

在結(jié)構(gòu)優(yōu)化和自洽計算都帶電且考慮偶極修正（NELECT=32, IDIPOL=4)后，在靜態(tài)自洽文件夾中：
grep "TOTEN" OUTCAR

free energy TOTEN =  -34.6451eV
*********** adding dipol energy to TOTEN NOW *************
free energy TOTEN = -28.091eV

看起來考慮偶極修正后，差異還是很大的。

問題：
（1）對于slab結(jié)構(gòu)吸附帶電離子，結(jié)構(gòu)優(yōu)化和靜態(tài)計算的時候需要加電荷。在結(jié)構(gòu)優(yōu)化的同時就需要考慮偶極修正（slab IDIPOL =1~3,帶電小分子 IDIPOL=4）？這樣計算好像要慢很多，各位一般是怎么做的？

（2）還有一項(xiàng)勢能修正（LDIPOL=.TRUE.），這個手冊上說對于帶電體系只對立方盒子有效。因?yàn)閹щ婓w系要做四極矩修正，這個只能在立方盒子超胞下才能做的。但是大部分超胞吸附模型都不是立方盒子，比如我的是六方的。吸附結(jié)構(gòu)引入高介電常數(shù)的底材，勢能修正項(xiàng)是否可以忽略呢？各位在計算的時候一般是怎么考慮勢能修正這一項(xiàng)的呢？

（3）偶極修正（dipol correction）好理解，勢能修正(potential correction)是什么意思呢？

（4）手冊上講：For charged cells or for calculations of molecules and surfaces with a large dipole moment, the energy converges very slowly with respect to the size L of the supercell.

...for charged systems with excess electrons more reliable results can be obtained if the energy after correction of the linear error (1/L) is plotted against 1/L^3 to extrapolates results manually for L->infinite.
這個多級矩修正與超胞尺寸L的關(guān)系怎么理解呢？

完全蒙。請過路高手跟帖解惑或者提供參考資料。

[ Last edited by cenwanglai on 2010-10-22 at 19:45 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

缺陷穩(wěn)定性分析

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

進(jìn)口色譜純級的甲醇、乙腈有沒有可能帶有很多的雜離子，大家交流交流經(jīng)驗(yàn) 已經(jīng)有15人回復(fù)
MS里怎么用原胞的坐標(biāo)建結(jié)構(gòu) 已經(jīng)有6人回復(fù)
“聲子液體”優(yōu)化熱電性能已經(jīng)有9人回復(fù)
gaussian 是如何進(jìn)行部分優(yōu)化的？已經(jīng)有5人回復(fù)
請教一個問題，優(yōu)化結(jié)構(gòu)鍵長過長怎么解決已經(jīng)有5人回復(fù)
結(jié)構(gòu)-熱分析時如何施加溫度梯度載荷？已經(jīng)有6人回復(fù)
非水體系銀/銀離子電極的制作方法已經(jīng)有7人回復(fù)
GaAs等異質(zhì)結(jié)構(gòu)太赫茲輻射材料已經(jīng)有6人回復(fù)
用Dmol優(yōu)化收斂的結(jié)構(gòu)，用vasp不收斂？已經(jīng)有9人回復(fù)
先進(jìn)鋰離子電池正極材料LiFePO4的研究進(jìn)展已經(jīng)有237人回復(fù)
過渡元素離子都能形成配合物嗎？已經(jīng)有1人回復(fù)
化學(xué)中可以通過什么作用來結(jié)合重金屬離子？已經(jīng)有11人回復(fù)
請教用molpro打印氫分子離子的重疊積矩陣、單電子積分矩陣、偶極矩陣已經(jīng)有3人回復(fù)
【求助】平坦勢能面下的結(jié)構(gòu)優(yōu)化已經(jīng)有4人回復(fù)

1樓 2010-09-20 18:19:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hgye

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1406
帖子: 180
在線: 46.1小時
蟲號: 461024
注冊: 2007-11-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
mazuju028(金幣+1):謝謝交流 2010-09-20 21:36:57

是的，必須的

回復(fù)此樓

2樓2010-09-20 19:04:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

問題更新了，希望大家給給予幫助！

贊一下

回復(fù)此樓

3樓2010-09-21 09:12:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

求助。

回復(fù)此樓

4樓2010-09-25 17:34:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunyang1988

至尊木蟲 (著名寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 24 (小學(xué)生)
貴賓: 0.984
金幣: 9397.7
散金: 1593
紅花: 23
沙發(fā): 1
帖子: 1744
在線: 1255.8小時
蟲號: 965316
注冊: 2010-03-09
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

結(jié)構(gòu)優(yōu)化與自洽計算的區(qū)別是什么呢，結(jié)構(gòu)優(yōu)化不需要自洽計算么？

贊一下

回復(fù)此樓

心之所向無懼無悔

5樓2010-09-25 17:59:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

引用回帖:

Originally posted by sunyang1988 at 2010-09-25 17:59:11:
結(jié)構(gòu)優(yōu)化與自洽計算的區(qū)別是什么呢，結(jié)構(gòu)優(yōu)化不需要自洽計算么？

呵呵，我的說法過于簡化了。

結(jié)構(gòu)優(yōu)化需要自洽。我在一樓說的自洽計算是指結(jié)構(gòu)優(yōu)化后的靜態(tài)自洽計算，k點(diǎn)加密了的。對于分子或者小的離子，我還是在Gamma點(diǎn)結(jié)構(gòu)優(yōu)化的基礎(chǔ)上加密K點(diǎn)做靜態(tài)自洽計算來算能量的。

[ Last edited by cenwanglai on 2010-9-25 at 19:00 ]

贊一下

回復(fù)此樓

6樓2010-09-25 18:58:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunyang1988

至尊木蟲 (著名寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 24 (小學(xué)生)
貴賓: 0.984
金幣: 9397.7
散金: 1593
紅花: 23
沙發(fā): 1
帖子: 1744
在線: 1255.8小時
蟲號: 965316
注冊: 2010-03-09
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

cenwanglai(金幣+2):謝謝回復(fù)。呵呵。 2010-09-25 21:08:23

原來是這樣，這樣算出來精度差了點(diǎn)吧
我算過LFPO（001）表面帶電荷的加dipole correction后，根本就不收斂

贊一下

回復(fù)此樓

心之所向無懼無悔

7樓2010-09-25 19:38:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

我去查查文獻(xiàn)，研究研究先。回復(fù)的人太少，不知道是什么原因。

贊一下

回復(fù)此樓

8樓2010-09-25 21:11:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fanchen021

鐵桿木蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 7994.9
散金: 1844
紅花: 4
帖子: 776
在線: 569小時
蟲號: 773188
注冊: 2009-05-18
專業(yè): 理論和計算化學(xué)

VASP 可以計算帶電離子的表面吸附么？是5.2 的新功能么
INCAR要如何設(shè)置？

贊一下

回復(fù)此樓

9樓2010-09-26 13:58:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

不是新功能。手冊6.57節(jié)有。vasp4.6是可以上算帶電離子吸附的。我和侯老師郵件交流的結(jié)果也是這個意思。

[ Last edited by cenwanglai on 2010-9-26 at 14:08 ]

贊一下

回復(fù)此樓

10樓2010-09-26 14:01:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

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

[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	3/150	2026-03-18 00:13 by PY-Hzb
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 070300化學(xué)319求調(diào)劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 274求調(diào)劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 308 085701 四六級已過求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考博] 2026年博士申請 +3	QwQwQW10 2026-03-11	3/150	2026-03-12 17:58 by gxch43