| 4 | 1/1 | 返回列表 |
| 查看: 1402 | 回復: 3 | ||
| 【懸賞金幣】回答本帖問題,作者plout..將贈送您 5 個金幣 | ||
plout..新蟲 (初入文壇)
|
[求助]
對于SrTiO3這個鈣鈦礦體系,如何算鍵解離能? 已有1人參與
|
|
|
有沒有MS大佬來教我咋算鍵解離能呀 發(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ā)自小木蟲手機客戶端 |
新蟲 (初入文壇)
送紅花一朵 |
我用DMOl3模塊,先進行超胞(2×2×2)然后計算總能量。再去掉一個氧原子,同樣用DMol3模塊計算,我對這個缺陷的晶體也加了一個自旋1,然后單獨算了一個氧原子的能量加了自旋3。最后得到差值是氧空位形成能。這是不是可行的?但是在這里我有個疑問,就是斷掉去掉一個氧原子那斷掉的是不是就不止一個鍵了,是不是還要除以鍵數(shù)呢? 發(fā)自小木蟲手機客戶端 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 080500材料科學與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | forgman95 2026-03-05 | 5/250 |
|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|
[考研]
材料學學碩308分/本科東北大學/一志愿西北工業(yè)大學/
5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
[考研] 290求調(diào)劑 +5 | Bananaiy 2026-03-04 | 5/250 |
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
[考研] 化學 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲 2026-03-03 | 11/550 |
|
|
[考研] 成績276,專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
[考研] 能動297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級已過,英語能力強 +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
[考研] 一志愿中科大能動297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
[基金申請]
剛錄用,沒有期刊號,但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|