| 4 | 1/1 | 返回列表 |
| 查看: 1421 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題,作者plout..將贈(zèng)送您 5 個(gè)金幣 | |||
plout..新蟲 (初入文壇)
|
[求助]
對(duì)于SrTiO3這個(gè)鈣鈦礦體系,如何算鍵解離能? 已有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個(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ī)客戶端 |
新蟲 (初入文壇)
送紅花一朵 |
我用DMOl3模塊,先進(jìn)行超胞(2×2×2)然后計(jì)算總能量。再去掉一個(gè)氧原子,同樣用DMol3模塊計(jì)算,我對(duì)這個(gè)缺陷的晶體也加了一個(gè)自旋1,然后單獨(dú)算了一個(gè)氧原子的能量加了自旋3。最后得到差值是氧空位形成能。這是不是可行的?但是在這里我有個(gè)疑問,就是斷掉去掉一個(gè)氧原子那斷掉的是不是就不止一個(gè)鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機(jī)客戶端 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856化工原理 +5 | z2839474511 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +4 | 步川酷紫123 2026-03-11 | 4/200 |
|
|
[考博] 讀博申請(qǐng) +5 | 感dd 2026-03-10 | 7/350 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +8 | mingmingry 2026-03-09 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-08 | 7/350 |
|
|
[考研] 收調(diào)劑 +7 | 調(diào)劑的考研學(xué)生 2026-03-10 | 7/350 |
|
|
[考研] 一志愿:武漢理工,材料工程,英二數(shù)二 總分314 +3 | 2202020125 2026-03-10 | 4/200 |
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼?~+123456 2026-03-08 | 8/400 |
|
|
[考研] 考研一志愿長安大學(xué)材料與化工309分請(qǐng)求調(diào)劑 +6 | dtdxzxx 2026-03-06 | 8/400 |
|
|
[考研] 293求調(diào)劑 +4 | 上班不著吉 2026-03-09 | 4/200 |
|
|
[考研] 考研調(diào)劑,一志愿山東大學(xué)材料與化工,328分,政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 8/400 |
|
|
[教師之家] 交大前校長王樹國:現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室,而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
[考研] 288求調(diào)劑(0703)一志愿東北大學(xué) +5 | 好好- 2026-03-07 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 誠求調(diào)劑,323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研]
|
程晴之 2026-03-06 | 6/300 |
|