cenwanglai

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

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

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

Gamma單點計算自洽磁矩為2，能量為-9.81eV.計算O空位形成能的時候，就是用這個能量的一半的。

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

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

看到文獻有中有單獨用ENCUT = 1000eV算O2形成能的。問問這樣做是否有必要？

歡迎評價我的計算細節(jié)。

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

回復此樓

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

第一原理資料匯編

» 猜你喜歡

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

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

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

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

查看全部散金貼

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

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

sunyang1988

至尊木蟲 (著名寫手)

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

很巧我也在算氧氣的結合能，我選的NECUT=600，選的15*15*15，E=-9.77eV
沒有設置初始磁矩，但會自動弛豫到有磁矩為2
按理說差距不該這么大啊，看來我算的有問題，你的INCAR怎么設的？

贊一下(2人)

回復此樓

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

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

cenwanglai

榮譽版主 (知名作家)

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

贊一下

回復此樓

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

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

cenwanglai

榮譽版主 (知名作家)

引用回帖:

Originally posted by sunyang1988 at 2011-03-14 09:21:29:
很巧我也在算氧氣的結合能，我選的NECUT=600，選的15*15*15，E=-9.77eV
沒有設置初始磁矩，但會自動弛豫到有磁矩為2
按理說差距不該這么大啊，看來我算的有問題，你的INCAR怎么設的？

》》回復了。

回復此樓

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

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

sunyang1988

至尊木蟲 (著名寫手)

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

我用的IBRION=2，POTIM=0.1，是不是有問題啊
POTIM再大一點就很容易不收斂，看來是不能用IBRION=2

贊一下(1人)

回復此樓

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

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

sunyang1988

至尊木蟲 (著名寫手)

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

我看到文獻上說GGA算分子是很不準的，對于O2來說，計算得出的O-O鍵over-binding，結合能偏大

贊一下(1人)

回復此樓

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

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

cenwanglai

榮譽版主 (知名作家)

引用回帖:

Originally posted by sunyang1988 at 2011-03-14 20:04:46:
我看到文獻上說GGA算分子是很不準的，對于O2來說，計算得出的O-O鍵over-binding，結合能偏大

是偏大的。實驗值是5.12eV.

贊一下

回復此樓

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

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

vimtex

銀蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 443.2
帖子: 34
在線: 17.5小時
蟲號: 481120

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

引用回帖:

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

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

贊一下(1人)

回復此樓

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

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

cenwanglai

榮譽版主 (知名作家)

引用回帖:

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.

贊一下

回復此樓

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

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

franch

版主 (著名寫手)

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

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

贊一下(1人)

回復此樓

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

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

ellsaking

金蟲 (正式寫手)

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

氧分子我也算過，dft能量為-9.85eV

贊一下(1人)

回復此樓

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

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

zhoudayan2

木蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 1550.3
帖子: 109
在線: 67小時
蟲號: 1098641

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

引用回帖:

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

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

贊一下(1人)

回復此樓

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

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

韓D小希

鐵蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 193.3
帖子: 440
在線: 157.7小時
蟲號: 1249416

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

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

贊一下

回復此樓

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

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

umklapp

銀蟲 (小有名氣)

應助: 9 (幼兒園)
金幣: 463.6
帖子: 123
在線: 59.1小時
蟲號: 1181441

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

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

贊一下

回復此樓

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

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

zhouhao

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 130.8
帖子: 19
在線: 20小時
蟲號: 7183686

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

引用回帖:

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

請問你用MS算出來了嗎？怎樣算的？

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

贊一下

回復此樓

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

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

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

[考研] 300求調劑，材料科學英一數(shù)二 +10	leaflight 2026-03-24	10/500	2026-03-29 20:58 by goldfish51
[考研] 一志愿華東師范大學有機化學專業(yè)，初試351分，復試被刷求調劑! +5	真名有冰 2026-03-29	6/300	2026-03-29 20:53 by 唐沐兒
[考研] 296求調劑 +6	彼岸t 2026-03-29	6/300	2026-03-29 20:22 by 永字號
[考研] 考研調劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學員8dgXkO
[考研] 環(huán)境工程 085701，267求調劑 +6	minht 2026-03-29	6/300	2026-03-29 16:21 by 學員8dgXkO
[考研] 070300化學354求調劑 +6	101次希望 2026-03-28	6/300	2026-03-29 12:57 by 無際的草原
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調劑 +8	AZMK 2026-03-25	12/600	2026-03-28 18:15 by AZMK
[考研] 085602 化工專碩 338分求調劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 346求調劑一志愿070303有機化學 +3	蘿卜燉青菜 2026-03-28	3/150	2026-03-28 14:11 by 唐沐兒
[考研] 材料與化工考研調劑 +17	孅華 2026-03-22	17/850	2026-03-28 08:35 by WYUMater
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 324求調劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[考研] 081200-11408-276學碩求調劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 284求調劑 +11	junqihahaha 2026-03-26	12/600	2026-03-27 04:37 by wxiongid
[考研] 321求調劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 一志愿河工大 081700 276求調劑 +4	地球繞著太陽轉 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 一志愿南京郵電大學 288分材料考研求調劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 309求調劑 +4	gajsj 2026-03-25	5/250	2026-03-26 00:27 by Dyhoer
[考研] 290分調劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 求調劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325