 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2798 | 回復(fù): 5 | |||||||
| 【獎勵】 本帖被評價5次，作者beefly增加金幣 3.8 個 | |||||||
beefly專家顧問 (職業(yè)作家)
|
[資源]
一個用ADF進(jìn)行能量分解計算的例子
|
||||||
|
ADF的能量分解類似于Morokuma的能量分解。后者應(yīng)用范圍很窄，只適用于閉殼層體系的Hartree-Fock計算，而前者可用于ADF中的各種RDFT和UDFT，甚至還能考慮旋軌耦合。 這個能量分解算例是以前做練習(xí)用的，今天偶然從箱底翻出來，既然老有人問類似的問題，就貼在下面。計算需要人工創(chuàng)建分子片，然后組合成HArF分子。本例要分析H-Ar鍵，所以需要創(chuàng)建的分子片是H和ArF。 #! /bin/sh $ADFBIN/dirac -n1 <$ADFRESOURCES/Dirac/H $ADFBIN/dirac -n1 <$ADFRESOURCES/Dirac/Ar $ADFBIN/dirac -n1 <$ADFRESOURCES/Dirac/F mv TAPE12 t12.rel $ADFBIN/adf -n1 < Create H file=$ADFRESOURCES/ZORA/TZ2P/H XC GGA pbe End Relativistic Scalar zora CorePotentials t12.rel & H 1 End End Input eor mv TAPE21 H.t21 $ADFBIN/adf -n1 < Create Ar file=$ADFRESOURCES/ZORA/TZ2P/Ar XC GGA pbe End Relativistic Scalar zora CorePotentials t12.rel & Ar 2 End End Input eor mv TAPE21 Ar.t21 $ADFBIN/adf -n1 < Create F file=$ADFRESOURCES/ZORA/TZ2P/F XC GGA pbe End Relativistic Scalar zora CorePotentials t12.rel & F 3 End End Input eor mv TAPE21 F.t21 $ADFBIN/adf < title ArF atoms Ar 0.000000 0.000000 0.000000 F 0.000000 0.000000 -1.950000 end fragments Ar Ar.t21 F F.t21 end XC gga pbe End relativistic scalar zora corepotentials t12.rel & Ar 2 F 3 end end input eor mv TAPE21 ArF.t21 $ADFBIN/adf << eor >> harf.out Title test atoms H 0.000000 0.000000 1.330000 Ar 0.000000 0.000000 0.000000 f=ArF F 0.000000 0.000000 -1.950000 f=ArF end restricted CHARGE 0 0 relativistic scalar zora corepotentials t12.rel & H 1 Ar 2 F 3 end integration 6.0 XC gga pbe End fragments H H.t21 ArF ArF.t21 end End input eor 計算完成后，在HArF計算的輸出部分找到如下結(jié)果： hartree eV kcal/mol kJ/mol -------------------- ----------- ---------- ----------- Pauli Repulsion Kinetic (Delta T^0): 1.782470985375970 48.5035 1118.52 4679.88 Delta V^Pauli Coulomb: -1.136710157595204 -30.9315 -713.30 -2984.43 Delta V^Pauli LDA-XC: -0.183123868092539 -4.9831 -114.91 -480.79 Delta V^Pauli GGA-Exchange: 0.012690279451145 0.3453 7.96 33.32 Delta V^Pauli GGA-Correlation: -0.003832183382116 -0.1043 -2.40 -10.06 -------------------- ----------- ---------- ----------- Total Pauli Repulsion: 0.471495055757255 12.8300 295.87 1237.91 (Total Pauli Repulsion = Delta E^Pauli in BB paper) Steric Interaction Pauli Repulsion (Delta E^Pauli): 0.471495055757255 12.8300 295.87 1237.91 Electrostatic Interaction: -0.117045546884956 -3.1850 -73.45 -307.30 (Electrostatic Interaction = Delta V_elstat in the BB paper) -------------------- ----------- ---------- ----------- Total Steric Interaction: 0.354449508872299 9.6451 222.42 930.61 (Total Steric Interaction = Delta E^0 in the BB paper) Orbital Interactions SIGMA: -0.454298601408872 -12.3621 -285.08 -1192.76 PI: -0.006896798127161 -0.1877 -4.33 -18.11 DELTA: 0.000000000000000 0.0000 0.00 0.00 PHI: 0.000000000000000 0.0000 0.00 0.00 -------------------- ----------- ---------- ----------- Total Orbital Interactions: -0.461195399536032 -12.5498 -289.40 -1210.87 Alternative Decomposition Orb.Int. Kinetic: -1.726978833271678 -46.9935 -1083.70 -4534.18 Coulomb: 1.236130927590834 33.6368 775.68 3245.46 XC: 0.029652506144804 0.8069 18.61 77.85 -------------------- ----------- ---------- ----------- Total Orbital Interactions: -0.461195399536040 -12.5498 -289.40 -1210.87 Residu (E=Steric+OrbInt+Res): -0.000000043854616 0.0000 0.00 0.00 Total Bonding Energy: -0.106745934518349 -2.9047 -66.98 -280.26 Summary of Bonding Energy (energy terms are taken from the energy decomposition above) ====================================================================================== Electrostatic Energy: -0.117045546884956 -3.1850 -73.45 -307.30 Kinetic Energy: 0.055492152104291 1.5100 34.82 145.69 Coulomb (Steric+OrbInt) Energy: 0.099420726141014 2.7054 62.39 261.03 XC Energy: -0.144613265878707 -3.9351 -90.75 -379.68 -------------------- ----------- ---------- ----------- Total Bonding Energy: -0.106745934518357 -2.9047 -66.98 -280.26 H-Ar鍵的分析過程非常簡單，對照這篇文章看就行了（因為手頭沒有電子版，不清楚具體頁碼，應(yīng)該是在比較靠后的位置），需要的數(shù)據(jù)都在上面列出。由于計算方法不完全一樣，所以這個計算的結(jié)果可能稍微不同于文獻(xiàn)值。 Matthias Lein, Jan Frunzke and Gernot Frenking, Christian Klixbull Jørgensen and the Nature of the Chemical Bond in HArF http://www.springerlink.com/content/e57dnhpudynhw1yv/ [ Last edited by beefly on 2011-5-4 at 09:27 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化軟件學(xué)習(xí) | 我的收藏 | 水和水的研究 |
» 猜你喜歡
cif求助下載
已經(jīng)有1人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有5人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有157人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有2人回復(fù)
-
化學(xué)化工學(xué)院 招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）
已經(jīng)有1人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學(xué)新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學(xué)方向）
已經(jīng)有6人回復(fù)
-
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有2人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- DFT Scan之后使用MP2計算能量，如何建立批量輸入文件？
已經(jīng)有6人回復(fù)
- vasp計算單個氧原子和分子的能量
已經(jīng)有20人回復(fù)
- 大家?guī)兔匆幌�，ADF計算帶有未成對電子體系的能量分解分析
已經(jīng)有3人回復(fù)
- ADF進(jìn)行能量分解計算
已經(jīng)有3人回復(fù)
- 【求助】有一個286個原子的體系，在ADF上使用什么方法和基組
已經(jīng)有8人回復(fù)
- 【討論】ADF能否z在考慮spin orbit相對論效應(yīng)下計算吸收光譜
已經(jīng)有8人回復(fù)
- 【求助】同一個結(jié)構(gòu)，G09和G03 opt計算的能量差別大嗎？
已經(jīng)有7人回復(fù)
- 【求助】如何計算一個原子的能量
已經(jīng)有10人回復(fù)
- 【求助】adf怎么算總能量，如果不能，用什么可以替代呢
已經(jīng)有3人回復(fù)
- 【求助】Gaussian計算完后能量看哪個值��？
已經(jīng)有23人回復(fù)
- 【sobereva個人文集】電子定域性的圖形分析
已經(jīng)有56人回復(fù)
- 【求助】計算雜質(zhì)形成能時 如何算單個原子的能量
已經(jīng)有11人回復(fù)
2樓2011-05-09 08:07:12
4樓2011-07-25 17:33:57
5樓2011-11-15 19:56:03
簡單回復(fù)
2011-05-15 15:37
回復(fù)
五星好評 
lilachi6樓
2013-03-14 19:57
回復(fù)
五星好評 謝謝分享~
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 beefly 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 一定有學(xué)上- 2026-03-14 | 8/400 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +9 | Copy267 2026-03-10 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 工科，求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[考研] 【考研調(diào)劑求收留】 +3 | Ceciilia 2026-03-11 | 3/150 |
|
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123！ 2026-03-11 | 4/200 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | 柯淮然 2026-03-10 | 8/400 |
|