| 查看: 2820 | 回复: 5 | |||||
| 【奖励】 本帖被评价5次，作者beefly增加金币 3.8 个 | |||||
[资源]
一个用ADF进行能量分解计算的例子
|
|||||
|
ADF的能量分解类似于Morokuma的能量分解。后者应用范围很窄，只适用于闭壳层体系的Hartree-Fock计算，而前者可用于ADF中的各种RDFT和UDFT，甚至还能考虑旋轨耦合。 这个能量分解算例是以前做练习用的，今天偶然从箱底翻出来，既然老有人问类似的问题，就贴在下面。计算需要人工创建分子片，然后组合成HArF分子。本例要分析H-Ar键，所以需要创建的分子片是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计算的输出部分找到如下结果： 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键的分析过程非常简单，对照这篇文章看就行了（因为手头没有电子版，不清楚具体页码，应该是在比较靠后的位置），需要的数据都在上面列出。由于计算方法不完全一样，所以这个计算的结果可能稍微不同于文献值。 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 量化软件学习 | 我的收藏 | 水和水的研究 |
» 猜你喜欢
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有17人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有278人回复
-
070300总分298求调剂
已经有2人回复
-
一志愿厦门大学化学学硕307求调剂
已经有13人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有0人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
-
硕士研究生调剂
已经有0人回复
-
中国电机工程学报的投稿系统
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- DFT Scan之后使用MP2计算能量，如何建立批量输入文件？
已经有6人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- 大家帮忙看一下，ADF计算带有未成对电子体系的能量分解分析
已经有3人回复
- ADF进行能量分解计算
已经有3人回复
- 【求助】有一个286个原子的体系，在ADF上使用什么方法和基组
已经有8人回复
- 【讨论】ADF能否z在考虑spin orbit相对论效应下计算吸收光谱
已经有8人回复
- 【求助】同一个结构，G09和G03 opt计算的能量差别大吗？
已经有7人回复
- 【求助】如何计算一个原子的能量
已经有10人回复
- 【求助】adf怎么算总能量，如果不能，用什么可以替代呢
已经有3人回复
- 【求助】Gaussian计算完后能量看哪个值啊？
已经有23人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
- 【求助】计算杂质形成能时 如何算单个原子的能量
已经有11人回复
2楼2011-05-09 08:07:12
4楼2011-07-25 17:33:57
5楼2011-11-15 19:56:03
简单回复
2011-05-15 15:37
回复
五星好评 
lilachi6楼
2013-03-14 19:57
回复
五星好评 谢谢分享~
