 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2817 | 回復(fù): 5 | |||||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)5次，作者beefly增加金幣 3.8 個(gè) | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||
beefly專家顧問 (職業(yè)作家)
|
[資源]
一個(gè)用ADF進(jìn)行能量分解計(jì)算的例子
|
||||||
|
ADF的能量分解類似于Morokuma的能量分解。后者應(yīng)用范圍很窄，只適用于閉殼層體系的Hartree-Fock計(jì)算，而前者可用于ADF中的各種RDFT和UDFT，甚至還能考慮旋軌耦合。 這個(gè)能量分解算例是以前做練習(xí)用的，今天偶然從箱底翻出來(lái)，既然老有人問類似的問題，就貼在下面。計(jì)算需要人工創(chuàng)建分子片，然后組合成HArF分子。本例要分析H-Ar鍵，所以需要?jiǎng)?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鍵的分析過(guò)程非常簡(jiǎn)單，對(duì)照這篇文章看就行了（因?yàn)槭诸^沒有電子版，不清楚具體頁(yè)碼，應(yīng)該是在比較靠后的位置），需要的數(shù)據(jù)都在上面列出。由于計(jì)算方法不完全一樣，所以這個(gè)計(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í) | 我的收藏 | 水和水的研究 |
» 猜你喜歡
一志愿南大理論與計(jì)算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計(jì)算機(jī)背景的考研咨詢
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有290人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團(tuán)隊(duì)招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎(jiǎng)學(xué)金
已經(jīng)有0人回復(fù)
-
【211博士招生】環(huán)境化學(xué)、地學(xué)、毒理方向，擅長(zhǎng)污染物遷移轉(zhuǎn)化降解、計(jì)算模擬者優(yōu)先
已經(jīng)有17人回復(fù)
-
求助如何利用vasp計(jì)算磁光克爾效應(yīng)
已經(jīng)有1人回復(fù)
-
070300總分298求調(diào)劑
已經(jīng)有2人回復(fù)
-
一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
已經(jīng)有11人回復(fù)
-
理論與計(jì)算化學(xué)方向調(diào)劑——初試388分（數(shù)學(xué)144）
已經(jīng)有0人回復(fù)
-
調(diào)劑招生
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- DFT Scan之后使用MP2計(jì)算能量，如何建立批量輸入文件？
已經(jīng)有6人回復(fù)
- vasp計(jì)算單個(gè)氧原子和分子的能量
已經(jīng)有20人回復(fù)
- 大家?guī)兔匆幌�，ADF計(jì)算帶有未成對(duì)電子體系的能量分解分析
已經(jīng)有3人回復(fù)
- ADF進(jìn)行能量分解計(jì)算
已經(jīng)有3人回復(fù)
- 【求助】有一個(gè)286個(gè)原子的體系，在ADF上使用什么方法和基組
已經(jīng)有8人回復(fù)
- 【討論】ADF能否z在考慮spin orbit相對(duì)論效應(yīng)下計(jì)算吸收光譜
已經(jīng)有8人回復(fù)
- 【求助】同一個(gè)結(jié)構(gòu)，G09和G03 opt計(jì)算的能量差別大嗎？
已經(jīng)有7人回復(fù)
- 【求助】如何計(jì)算一個(gè)原子的能量
已經(jīng)有10人回復(fù)
- 【求助】adf怎么算總能量，如果不能，用什么可以替代呢
已經(jīng)有3人回復(fù)
- 【求助】Gaussian計(jì)算完后能量看哪個(gè)值��？
已經(jīng)有23人回復(fù)
- 【sobereva個(gè)人文集】電子定域性的圖形分析
已經(jīng)有56人回復(fù)
- 【求助】計(jì)算雜質(zhì)形成能時(shí) 如何算單個(gè)原子的能量
已經(jīng)有11人回復(fù)
5樓2011-11-15 19:56:03
2樓2011-05-09 08:07:12
4樓2011-07-25 17:33:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 林深溫瀾 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 265求調(diào)劑 +11 | yelck 2026-04-01 | 12/600 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +15 | 柳溪邊 2026-03-26 | 16/800 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +10 | hanamiko 2026-04-01 | 10/500 |
|
|
|
[考研] 283求調(diào)劑 +9 | A child 2026-03-28 | 9/450 |
|
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
|
[碩博家園] 博一被送出聯(lián)培感覺不適應(yīng)怎么辦 +3 | 全村的狗 2026-03-31 | 3/150 |
|
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +25 | lllllos 2026-03-30 | 26/1300 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +12 | Koxui 2026-03-28 | 12/600 |
|
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑 +9 | s1145 2026-03-31 | 9/450 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-30 | 6/300 |
|
|
|
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|