| 4 | 1/1 | 返回列表 |
| 查看: 1401 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題,作者plout..將贈送您 5 個金幣 | ||
plout..新蟲 (初入文壇)
|
[求助]
對于SrTiO3這個鈣鈦礦體系,如何算鍵解離能? 已有1人參與
|
|
|
有沒有MS大佬來教我咋算鍵解離能呀 發(fā)自小木蟲手機(jī)客戶端 |
木蟲 (正式寫手)
|
用 Quantum ESPRESSO 來計(jì)算 SrTiO₃ 中的鍵解離能(Bond Dissociation Energy, BDE) 是完全可行的。下面是詳細(xì)的 DFT 計(jì)算流程,假設(shè)你想算的是 Ti–O 鍵的 BDE。 1. 使用大一點(diǎn)的超胞,例如:2x2x2 SrTiO₃ 擴(kuò)展晶胞,共40個原子. 目的是避免空位之間的相互作用影響能量計(jì)算。 2. 完美晶體的總能計(jì)算. 用 pw.x 做 self-consistent 計(jì)算 &control calculation = 'scf', prefix = 'sto', outdir = './tmp', pseudo_dir = './pseudo', / &system ibrav = 0, nat = 40, ntyp = 3, ecutwfc = 60, ecutrho = 600, occupations = 'smearing', smearing = 'mp', degauss = 0.01, / &electrons conv_thr = 1.0d-8, / ATOMIC_SPECIES Sr 87.62 Sr.pbe-nsp-van.UPF Ti 47.867 Ti.pbe-sp-van.UPF O 15.999 O.pbe-n-van.UPF CELL_PARAMETERS angstrom ... # lattice vectors here ATOMIC_POSITIONS angstrom ... # 40 atoms here K_POINTS automatic 4 4 4 0 0 0 3. 創(chuàng)建一個新的結(jié)構(gòu),移除一個Ti附近的O原子,形成一個氧空位,保存為 SrTiO3-Vo.scf.in。再做一次 scf 計(jì)算。 結(jié)構(gòu)弛豫(使用 relax)是可選的,但建議進(jìn)行: calculation = 'relax' 4. 氧氣分子的總能(O₂)模擬一個單獨(dú)的 O₂ 分子作為化學(xué)勢參考,文件名:O2.scf.in 建一個大真空盒(~15 ų) 自旋極化(O₂是三重態(tài)) &control calculation = 'scf', prefix = 'o2', outdir = './tmp', pseudo_dir = './pseudo', / &system ibrav = 1, celldm(1) = 30.0, ! cubic box ~15 Å nat = 2, ntyp = 1, ecutwfc = 60, ecutrho = 600, nspin = 2, starting_magnetization(1) = 1.0, / &electrons conv_thr = 1.0d-8, / ATOMIC_SPECIES O 15.999 O.pbe-n-van.UPF ATOMIC_POSITIONS angstrom O 0.0 0.0 0.0 O 1.2 0.0 0.0 K_POINTS gamma 計(jì)算得到的 O₂ 總能 ,然后除以2,近似為單個O原子的能量。 5. 計(jì)算氧空位形成能 ≈ Ti–O 鍵解離能 𝐸form𝑉𝑂=𝐸defect+1/2𝐸O2−𝐸perfect𝐸perfect 這個值大致就是斷裂一個 Ti–O 鍵所需的能量(會比真實(shí)的 BDE 高一些,因?yàn)榭紤]了固體環(huán)境和晶格松弛)。 |
新蟲 (初入文壇)
|
我用DMOl3模塊,先進(jìn)行超胞(2×2×2)然后計(jì)算總能量。再去掉一個氧原子,同樣用DMol3模塊計(jì)算,我對這個缺陷的晶體也加了一個自旋1,然后單獨(dú)算了一個氧原子的能量加了自旋3。最后得到差值是氧空位形成能。但是在這里我有個疑問,就是斷掉去掉一個氧原子那斷掉的是不是就不止一個鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機(jī)客戶端 |
新蟲 (初入文壇)
送紅花一朵 |
我用DMOl3模塊,先進(jìn)行超胞(2×2×2)然后計(jì)算總能量。再去掉一個氧原子,同樣用DMol3模塊計(jì)算,我對這個缺陷的晶體也加了一個自旋1,然后單獨(dú)算了一個氧原子的能量加了自旋3。最后得到差值是氧空位形成能。這是不是可行的?但是在這里我有個疑問,就是斷掉去掉一個氧原子那斷掉的是不是就不止一個鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機(jī)客戶端 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +4 | 啊擺啊擺 2026-03-05 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +6 | 帆船哥 2026-03-04 | 8/400 |
|
|
[考研] 282求調(diào)劑 +6 | 夕~日 2026-03-05 | 7/350 |
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
[考研] 考研282分求調(diào)劑,接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 6/300 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +8 | wengyujian 2026-03-04 | 8/400 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 0855機(jī)械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
[考研] 26考研報考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩0927 2026-03-02 | 5/250 |
|