cenwanglai

[交流] 【討論】孤立O2分子的DFT總能量

氧空位的形成能的計(jì)算與O2分子的能量有關(guān)。它的三線態(tài)能量最低。

我用邊長(zhǎng)13的box, PAW-GGA(PW91)，ENCUT= 480 ISPIN=2 MAGMOM=2*1 Gamma點(diǎn) 來(lái)優(yōu)化，鍵長(zhǎng)為1.23.

Gamma單點(diǎn)計(jì)算自洽磁矩為2，能量為-9.81eV.計(jì)算O空位形成能的時(shí)候，就是用這個(gè)能量的一半的。

單個(gè)O原子能量-1.61eV. O2形成能為-6.59eV.

如果要指定單線態(tài)，設(shè)置updown = 0是不是可以實(shí)現(xiàn)呢？要指定三線態(tài)，指定updown=1?

看到文獻(xiàn)有中有單獨(dú)用ENCUT = 1000eV算O2形成能的。問(wèn)問(wèn)這樣做是否有必要？

歡迎評(píng)價(jià)我的計(jì)算細(xì)節(jié)。

[ Last edited by cenwanglai on 2011-3-14 at 22:29 ]

回復(fù)此樓

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

第一原理資料匯編

» 猜你喜歡

334分一志愿武理材料求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有18人回復(fù)
環(huán)境工程 085701，267求調(diào)劑已經(jīng)有6人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有12人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有12人回復(fù)
277跪求調(diào)劑已經(jīng)有10人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有19人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有15人回復(fù)
085600 286分材料求調(diào)劑已經(jīng)有9人回復(fù)
南京大學(xué)化學(xué)調(diào)劑已經(jīng)有9人回復(fù)

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

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

IgG是一種抗體，那anti IgG呢？已經(jīng)有10人回復(fù)
PCR退火溫度態(tài)度會(huì)導(dǎo)致擴(kuò)增不出目的條帶么？已經(jīng)有4人回復(fù)
復(fù)試中遇到的一個(gè)問(wèn)題已經(jīng)有28人回復(fù)
請(qǐng)教藥物設(shè)計(jì)中的帶負(fù)電性基團(tuán)和帶正電性基團(tuán)的概念已經(jīng)有3人回復(fù)
求助，為什么SSR擴(kuò)增產(chǎn)物瓊脂糖跑出的條帶會(huì)比maker最大分子量還要大呢？？已經(jīng)有7人回復(fù)
3.22文獻(xiàn)解讀：三元材料在高能量密度鋰離子電池中的應(yīng)用已經(jīng)有236人回復(fù)
DFT Scan之后使用MP2計(jì)算能量，如何建立批量輸入文件？已經(jīng)有6人回復(fù)
【求助】O2分子的鍵能很高（494kJ/mol）,為什么還很活潑？已經(jīng)有11人回復(fù)
【求助】如何利用分子動(dòng)力學(xué)加速DFT計(jì)算的收斂已經(jīng)有6人回復(fù)
【求助】小分子用DFT能算，大分子用AM1也能算，但換DFT就會(huì)出錯(cuò) 已經(jīng)有12人回復(fù)
【neweroica個(gè)人文集】DFT-D能量校正項(xiàng)計(jì)算程序（EDisp v0.2）已經(jīng)有29人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-14 06:02:20

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

sunyang1988

至尊木蟲 (著名寫手)

★ ★ ★
cenwanglai(金幣+1):謝謝參與
zzy870720z(金幣+2): 鼓勵(lì)討論 2011-03-14 12:57:46

很巧我也在算氧氣的結(jié)合能，我選的NECUT=600，選的15*15*15，E=-9.77eV
沒有設(shè)置初始磁矩，但會(huì)自動(dòng)弛豫到有磁矩為2
按理說(shuō)差距不該這么大啊，看來(lái)我算的有問(wèn)題，你的INCAR怎么設(shè)的？

贊一下(2人)

回復(fù)此樓

2樓2011-03-14 09:21:29

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

cenwanglai

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

geo_opt INCAR:
ENCUT = 480
EDIFF= 1E-5
IBRION = 1
ALGO = F
NFREE = 4
NSW = 30
NELMIN = 5
ISPIN = 2
MAGMON = 2*1
EDIFFG = -0.02

scf INCAR:
ISTART = 0
ICHARG = 2
EDIFF = 1E-6
NELMIN = 5
ISPIN = 2
VOSKOWN = 1

贊一下

回復(fù)此樓

3樓2011-03-14 18:57:06

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

cenwanglai

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

引用回帖:

Originally posted by sunyang1988 at 2011-03-14 09:21:29:
很巧我也在算氧氣的結(jié)合能，我選的NECUT=600，選的15*15*15，E=-9.77eV
沒有設(shè)置初始磁矩，但會(huì)自動(dòng)弛豫到有磁矩為2
按理說(shuō)差距不該這么大啊，看來(lái)我算的有問(wèn)題，你的INCAR怎么設(shè)的？

》》回復(fù)了。

回復(fù)此樓

4樓2011-03-14 18:58:40

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

sunyang1988

至尊木蟲 (著名寫手)

★
zzy870720z(金幣+1): 鼓勵(lì)討論 2011-03-15 08:32:35
cenwanglai(金幣+2): 謝謝回復(fù) 2011-05-16 05:47:09

我用的IBRION=2，POTIM=0.1，是不是有問(wèn)題啊
POTIM再大一點(diǎn)就很容易不收斂，看來(lái)是不能用IBRION=2

贊一下(1人)

回復(fù)此樓

5樓2011-03-14 20:01:08

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

sunyang1988

至尊木蟲 (著名寫手)

★
zzy870720z(金幣+1): 鼓勵(lì)討論 2011-03-15 08:32:44
cenwanglai(金幣+2): 謝謝回復(fù) 2011-05-16 05:47:22

我看到文獻(xiàn)上說(shuō)GGA算分子是很不準(zhǔn)的，對(duì)于O2來(lái)說(shuō)，計(jì)算得出的O-O鍵over-binding，結(jié)合能偏大

贊一下(1人)

回復(fù)此樓

6樓2011-03-14 20:04:46

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

cenwanglai

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

引用回帖:

Originally posted by sunyang1988 at 2011-03-14 20:04:46:
我看到文獻(xiàn)上說(shuō)GGA算分子是很不準(zhǔn)的，對(duì)于O2來(lái)說(shuō)，計(jì)算得出的O-O鍵over-binding，結(jié)合能偏大

是偏大的。實(shí)驗(yàn)值是5.12eV.

贊一下

回復(fù)此樓

7樓2011-03-14 20:22:47

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

vimtex

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 443.2
帖子: 34
在線: 17.5小時(shí)
蟲號(hào): 481120

★
cenwanglai(金幣+1):謝謝參與

引用回帖:

Originally posted by cenwanglai at 2011-03-14 18:58:40:
》》回復(fù)了。

I prefer using kpts=1x1x1 to 15x15x15, because we are trying to mimic the O2 molecular which is not a periodic system

贊一下(1人)

回復(fù)此樓

8樓2011-03-20 05:20:19

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

cenwanglai

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

引用回帖:

Originally posted by vimtex at 2011-03-20 05:20:19:
I prefer using kpts=1x1x1 to 15x15x15, because we are trying to mimic the O2 molecular which is not a periodic system

We used kpts 1x1x1 too.

贊一下

回復(fù)此樓

9樓2011-03-20 09:08:36

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

franch

版主 (著名寫手)

★
cenwanglai(金幣+1):謝謝參與

不一樣的參數(shù)，，一樣的結(jié)果，，用MS的路過(guò)。。

贊一下(1人)

回復(fù)此樓

10樓2011-04-06 11:26:36

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

ellsaking

金蟲 (正式寫手)

★
cenwanglai(金幣+1):謝謝參與
cenwanglai(金幣+1): 謝謝回復(fù)！ 2011-05-16 05:48:26

氧分子我也算過(guò)，dft能量為-9.85eV

贊一下(1人)

回復(fù)此樓

11樓2011-04-06 14:00:36

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

zhoudayan2

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1550.3
帖子: 109
在線: 67小時(shí)
蟲號(hào): 1098641

★
cenwanglai(金幣+1):謝謝參與

引用回帖:

Originally posted by ellsaking at 2011-04-06 14:00:36:
氧分子我也算過(guò)，dft能量為-9.85eV

我用PAW-PBE算的，ENCUT=400，得到的結(jié)果死-9.85ev

贊一下(1人)

回復(fù)此樓

12樓2011-04-25 22:19:31

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

韓D小希

鐵蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 193.3
帖子: 440
在線: 157.7小時(shí)
蟲號(hào): 1249416

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

樓主您好，請(qǐng)問(wèn)是不是開自旋優(yōu)化得磁矩為2，就是三線態(tài)？謝謝

贊一下

回復(fù)此樓

13樓2015-05-14 16:42:30

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

umklapp

銀蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 463.6
帖子: 123
在線: 59.1小時(shí)
蟲號(hào): 1181441

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

Singlet 可以通過(guò)設(shè)定 NUPDOWN＝0 得到，你可以在NUPDOWN＝0 下收斂電子態(tài)密度，然后去掉constrain，再收斂一次；
同時(shí)注意symmetry breaking 對(duì)singlet 能量有影響， higher symmetry, 能量也高

贊一下

回復(fù)此樓

14樓2015-05-15 04:58:40

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

zhouhao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 130.8
帖子: 19
在線: 20小時(shí)
蟲號(hào): 7183686

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

10樓: Originally posted by franch at 2011-04-06 11:26:36
不一樣的參數(shù)，，一樣的結(jié)果，，用MS的路過(guò)。。

請(qǐng)問(wèn)你用MS算出來(lái)了嗎？怎樣算的？

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

15樓2018-01-01 21:42:39

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

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

[考研] 環(huán)境工程 085701，267求調(diào)劑 +6	minht 2026-03-29	6/300	2026-03-29 16:21 by 學(xué)員8dgXkO
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 總分293求調(diào)劑 +6	加一一九 2026-03-25	9/450	2026-03-29 13:40 by 加一一九
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 332求調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 考研化學(xué)308分求調(diào)劑 +10	你好明天你好 2026-03-23	12/600	2026-03-27 14:43 by shangxh
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[考研] 求調(diào)劑 +6	研研，接電話 2026-03-24	7/350	2026-03-24 17:01 by barlinike
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte