| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1408 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題,作者plout..將贈(zèng)送您 5 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
plout..新蟲 (初入文壇)
|
[求助]
對(duì)于SrTiO3這個(gè)鈣鈦礦體系,如何算鍵解離能? 已有1人參與
|
||
|
有沒有MS大佬來教我咋算鍵解離能呀 發(fā)自小木蟲手機(jī)客戶端 |
新蟲 (初入文壇)
|
我用DMOl3模塊,先進(jìn)行超胞(2×2×2)然后計(jì)算總能量。再去掉一個(gè)氧原子,同樣用DMol3模塊計(jì)算,我對(duì)這個(gè)缺陷的晶體也加了一個(gè)自旋1,然后單獨(dú)算了一個(gè)氧原子的能量加了自旋3。最后得到差值是氧空位形成能。但是在這里我有個(gè)疑問,就是斷掉去掉一個(gè)氧原子那斷掉的是不是就不止一個(gè)鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機(jī)客戶端 |
木蟲 (正式寫手)
|
用 Quantum ESPRESSO 來計(jì)算 SrTiO₃ 中的鍵解離能(Bond Dissociation Energy, BDE) 是完全可行的。下面是詳細(xì)的 DFT 計(jì)算流程,假設(shè)你想算的是 Ti–O 鍵的 BDE。 1. 使用大一點(diǎn)的超胞,例如:2x2x2 SrTiO₃ 擴(kuò)展晶胞,共40個(gè)原子. 目的是避免空位之間的相互作用影響能量計(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)建一個(gè)新的結(jié)構(gòu),移除一個(gè)Ti附近的O原子,形成一個(gè)氧空位,保存為 SrTiO3-Vo.scf.in。再做一次 scf 計(jì)算。 結(jié)構(gòu)弛豫(使用 relax)是可選的,但建議進(jìn)行: calculation = 'relax' 4. 氧氣分子的總能(O₂)模擬一個(gè)單獨(dú)的 O₂ 分子作為化學(xué)勢(shì)參考,文件名:O2.scf.in 建一個(gè)大真空盒(~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,近似為單個(gè)O原子的能量。 5. 計(jì)算氧空位形成能 ≈ Ti–O 鍵解離能 𝐸form𝑉𝑂=𝐸defect+1/2𝐸O2−𝐸perfect𝐸perfect 這個(gè)值大致就是斷裂一個(gè) Ti–O 鍵所需的能量(會(huì)比真實(shí)的 BDE 高一些,因?yàn)榭紤]了固體環(huán)境和晶格松弛)。 |
新蟲 (初入文壇)
送紅花一朵 |
我用DMOl3模塊,先進(jìn)行超胞(2×2×2)然后計(jì)算總能量。再去掉一個(gè)氧原子,同樣用DMol3模塊計(jì)算,我對(duì)這個(gè)缺陷的晶體也加了一個(gè)自旋1,然后單獨(dú)算了一個(gè)氧原子的能量加了自旋3。最后得到差值是氧空位形成能。這是不是可行的?但是在這里我有個(gè)疑問,就是斷掉去掉一個(gè)氧原子那斷掉的是不是就不止一個(gè)鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機(jī)客戶端 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +7 | 倘若起風(fēng)了呢 2026-03-05 | 9/450 |
|
|---|---|---|---|---|
|
[考博] 售瀆SCI文章,我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-06 | 5/250 |
|
|
[考研] 初試297求調(diào)劑 +3 | 一個(gè)蘿蘿 2026-02-28 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 考研282分求調(diào)劑,接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
[考研] 0856材料與化工,270求調(diào)劑 +17 | YXCT 2026-03-01 | 20/1000 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
[考研] 298求調(diào)劑 +7 | axyz3 2026-02-28 | 8/400 |
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩(shī)0927 2026-03-02 | 5/250 |
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|