| 5 | 1/1 | 返回列表 |
| 查看: 1411 | 回復: 3 | |||
| 【懸賞金幣】回答本帖問題,作者plout..將贈送您 5 個金幣 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
plout..新蟲 (初入文壇)
|
[求助]
對于SrTiO3這個鈣鈦礦體系,如何算鍵解離能? 已有1人參與
|
||
|
有沒有MS大佬來教我咋算鍵解離能呀 發(fā)自小木蟲手機客戶端 |
新蟲 (初入文壇)
送紅花一朵 |
我用DMOl3模塊,先進行超胞(2×2×2)然后計算總能量。再去掉一個氧原子,同樣用DMol3模塊計算,我對這個缺陷的晶體也加了一個自旋1,然后單獨算了一個氧原子的能量加了自旋3。最后得到差值是氧空位形成能。這是不是可行的?但是在這里我有個疑問,就是斷掉去掉一個氧原子那斷掉的是不是就不止一個鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機客戶端 |
木蟲 (正式寫手)
|
用 Quantum ESPRESSO 來計算 SrTiO₃ 中的鍵解離能(Bond Dissociation Energy, BDE) 是完全可行的。下面是詳細的 DFT 計算流程,假設你想算的是 Ti–O 鍵的 BDE。 1. 使用大一點的超胞,例如:2x2x2 SrTiO₃ 擴展晶胞,共40個原子. 目的是避免空位之間的相互作用影響能量計算。 2. 完美晶體的總能計算. 用 pw.x 做 self-consistent 計算 &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)建一個新的結構,移除一個Ti附近的O原子,形成一個氧空位,保存為 SrTiO3-Vo.scf.in。再做一次 scf 計算。 結構弛豫(使用 relax)是可選的,但建議進行: calculation = 'relax' 4. 氧氣分子的總能(O₂)模擬一個單獨的 O₂ 分子作為化學勢參考,文件名: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 計算得到的 O₂ 總能 ,然后除以2,近似為單個O原子的能量。 5. 計算氧空位形成能 ≈ Ti–O 鍵解離能 𝐸form𝑉𝑂=𝐸defect+1/2𝐸O2−𝐸perfect𝐸perfect 這個值大致就是斷裂一個 Ti–O 鍵所需的能量(會比真實的 BDE 高一些,因為考慮了固體環(huán)境和晶格松弛)。 |
新蟲 (初入文壇)
|
我用DMOl3模塊,先進行超胞(2×2×2)然后計算總能量。再去掉一個氧原子,同樣用DMol3模塊計算,我對這個缺陷的晶體也加了一個自旋1,然后單獨算了一個氧原子的能量加了自旋3。最后得到差值是氧空位形成能。但是在這里我有個疑問,就是斷掉去掉一個氧原子那斷掉的是不是就不止一個鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組,低分調劑,招收機械電子信息通信等交叉方向 +52 | lisimayy 2026-03-04 | 65/3250 |
|
|---|---|---|---|---|
|
[考研] 337求調劑 +3 | 睡醒,。 2026-03-09 | 3/150 |
|
|
[考研] 考研調劑,一志愿山東大學材料與化工,328分,政治51 +4 | 關你西紅柿929 2026-03-08 | 5/250 |
|
|
[考研] 求0856調劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 0817一志愿蘇州大學280復試調劑 +3 | kk扛 2026-03-07 | 3/150 |
|
|
[考研] 083000,總分284,求調劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
[考研]
|
Sixuan wang 2026-03-06 | 7/350 |
|
|
[考研] 0307化學求調劑 +6 | 0ok0k 2026-03-07 | 6/300 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學0856材料與化工,前三科206,總分283,求調劑 +7 | 26考研求調劑 2026-03-06 | 7/350 |
|
|
[考研] 085600材料與化工 298 調劑 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|
[考研] 復試調劑 +7 | 呼呼?~+123456 2026-03-05 | 10/500 |
|
|
[考研] 材料與化工304求B區(qū)調劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研]
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
[考博] 2026申博自薦 六級440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
[考研] 304求調劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 一志愿天津大學085600 319分 材料與化工 金屬方向 求調劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
[考研] 材料學碩080500復試調劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
[考研] 322,求調劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
[考研] 0856材料工程,初試313調劑 +7 | 賣個關子吧 2026-03-03 | 7/350 |
|
|
[考研] 085602化學工程350,調劑,有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|