 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2812 | 回復(fù): 5 | |||||||
| 【獎勵】 本帖被評價5次，作者beefly增加金幣 3.8 個 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||
beefly專家顧問 (職業(yè)作家)
|
[資源]
一個用ADF進(jìn)行能量分解計(jì)算的例子
|
||||||
|
ADF的能量分解類似于Morokuma的能量分解。后者應(yīng)用范圍很窄，只適用于閉殼層體系的Hartree-Fock計(jì)算，而前者可用于ADF中的各種RDFT和UDFT，甚至還能考慮旋軌耦合。 這個能量分解算例是以前做練習(xí)用的，今天偶然從箱底翻出來，既然老有人問類似的問題，就貼在下面。計(jì)算需要人工創(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 計(jì)算完成后，在HArF計(jì)算的輸出部分找到如下結(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àn)槭诸^沒有電子版，不清楚具體頁碼，應(yīng)該是在比較靠后的位置），需要的數(shù)據(jù)都在上面列出。由于計(jì)算方法不完全一樣，所以這個計(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í) | 我的收藏 | 水和水的研究 |
» 猜你喜歡
CSC訪學(xué)博后項(xiàng)目獲批，外方學(xué)校暫停合作該怎么辦？
已經(jīng)有51人回復(fù)
-
[電子教材]天津大學(xué)物理化學(xué)（第七版）上、下冊
已經(jīng)有4人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有72人回復(fù)
-
陰離子交換膜電解二氧化碳還原
已經(jīng)有0人回復(fù)
-
分子動力學(xué)模擬合作需求
已經(jīng)有2人回復(fù)
-
九江學(xué)院2026年最新高層次人才招聘公告
已經(jīng)有0人回復(fù)
-
深圳信息職業(yè)技術(shù)學(xué)院2026年教師招聘公告（最新）
已經(jīng)有0人回復(fù)
-
研究生招生
已經(jīng)有0人回復(fù)
-
山西大學(xué)資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有2人回復(fù)
-
一志愿南師大0703化學(xué) 275求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- DFT Scan之后使用MP2計(jì)算能量，如何建立批量輸入文件？
已經(jīng)有6人回復(fù)
- vasp計(jì)算單個氧原子和分子的能量
已經(jīng)有20人回復(fù)
- 大家?guī)兔匆幌�，ADF計(jì)算帶有未成對電子體系的能量分解分析
已經(jīng)有3人回復(fù)
- ADF進(jìn)行能量分解計(jì)算
已經(jīng)有3人回復(fù)
- 【求助】有一個286個原子的體系，在ADF上使用什么方法和基組
已經(jīng)有8人回復(fù)
- 【討論】ADF能否z在考慮spin orbit相對論效應(yīng)下計(jì)算吸收光譜
已經(jīng)有8人回復(fù)
- 【求助】同一個結(jié)構(gòu)，G09和G03 opt計(jì)算的能量差別大嗎？
已經(jīng)有7人回復(fù)
- 【求助】如何計(jì)算一個原子的能量
已經(jīng)有10人回復(fù)
- 【求助】adf怎么算總能量，如果不能，用什么可以替代呢
已經(jīng)有3人回復(fù)
- 【求助】Gaussian計(jì)算完后能量看哪個值�。�
已經(jīng)有23人回復(fù)
- 【sobereva個人文集】電子定域性的圖形分析
已經(jīng)有56人回復(fù)
- 【求助】計(jì)算雜質(zhì)形成能時 如何算單個原子的能量
已經(jīng)有11人回復(fù)
4樓2011-07-25 17:33:57
2樓2011-05-09 08:07:12
5樓2011-11-15 19:56:03
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 286分 材料求調(diào)劑 +4 | 麻辣魷魚 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +7 | 新時代材料 2026-03-27 | 7/350 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 292求調(diào)劑 +13 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 14/700 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-27 | 7/350 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +8 | zhubinhao 2026-03-27 | 8/400 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 0856求調(diào)劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-24 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|