| 查看: 1316 | 回复: 6 | |||
| 【奖励】 本帖被评价1次，作者beefly增加金币 0.8 个 | |||
| 本帖产生 1 个 QC强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[资源]
【原创】做重元素post－HF计算的注意了，Gaussian的一个BUG！已有2人参与
|
|||
|
在Gaussian所有的版本中，默认的冻芯设置一旦遇到f轨道占据的情况，很可能发生问题。例如汞化合物的计算，Gaussian都倾向于把4f(-3)和4f(+3)分量的轨道冻结，只对4f(-2,-1,0,+1,+2)轨道考虑相关能。如果用来优化结构，这么做也许还行。可是千万别做原子计算，因为4f的7简并个轨道分量会把Gaussian逼疯！同样的输入文件，腊月三十是一个结果，正月初一是另外一个结果。所以做post－HF计算不要依赖Gaussian，最好自己设置冻芯，例如：CCSD(T,window=n)，其中的n值为芯轨道数+1。 最后，亲身体会一下疯狂的感觉吧： #p mp2/gen 5d 7f Int=dkh2 DKH2-MP2 for Hg 0 1 Hg Hg 0 s 6 1.0 1778058.5043390000 0.8962883239D-01 790248.2241510000 -0.2283432867D-01 351221.4329560000 0.1802289753D+00 156098.4146470000 0.7672072391D-01 69377.0731760000 0.3713806894D+00 30834.2547450000 0.4601707434D+00 s 1 1.0 13704.1132200000 1.0 s 1 1.0 6090.7169870000 1.0 s 1 1.0 2706.9853270000 1.0 s 1 1.0 1203.1045900000 1.0 s 1 1.0 534.7131510000 1.0 s 1 1.0 237.6502890000 1.0 s 1 1.0 105.6223510000 1.0 s 1 1.0 46.9432670000 1.0 s 1 1.0 20.8636740000 1.0 s 1 1.0 9.2727440000 1.0 s 1 1.0 4.1212200000 1.0 s 1 1.0 1.8316530000 1.0 s 1 1.0 0.8140680000 1.0 s 1 1.0 0.3618080000 1.0 s 1 1.0 0.1608040000 1.0 s 1 1.0 0.0714680000 1.0 p 5 1.0 24956.6090260000 0.1712344575D-01 9982.6436100000 0.1579005785D-01 3993.0574440000 0.9011193001D-01 1597.2229780000 0.2528001579D+00 638.8891910000 0.7338856627D+00 p 1 1.0 255.5556760000 1.0 p 1 1.0 102.2222710000 1.0 p 1 1.0 40.8889080000 1.0 p 1 1.0 16.3555630000 1.0 p 1 1.0 6.5422250000 1.0 p 1 1.0 2.6168900000 1.0 p 1 1.0 1.0467560000 1.0 p 1 1.0 0.4187020000 1.0 p 1 1.0 0.1674810000 1.0 p 1 1.0 0.0669920000 1.0 d 4 1.0 1928.0943100000 0.8055355574D-02 701.1252040000 0.4429306498D-01 254.9546200000 0.2455594297D+00 92.7107710000 0.8056186267D+00 d 1 1.0 33.7130080000 1.0 d 1 1.0 12.2592750000 1.0 d 1 1.0 4.4579180000 1.0 d 1 1.0 1.6210610000 1.0 d 1 1.0 0.5894770000 1.0 d 1 1.0 0.2143550000 1.0 d 1 1.0 0.0779470000 1.0 f 4 1.0 96.6910160000 0.5840546881D-01 32.2303390000 0.2860558584D+00 10.7434460000 0.5719063200D+00 3.5811490000 0.3785524869D+00 f 1 1.0 1.1937160000 1.0 f 1 1.0 0.3979050000 1.0 **** 然后察看冻芯轨道数（搜NFC=）。有时候是27，有时候是28，总能量当然是不一样的。如果算键能（分子总能量减去原子总能量），结果可想而知。 |
回复此楼» 收录本帖的淘帖专辑推荐
 基组的介绍总结 |
» 猜你喜欢
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有177人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
Gaussian计算NMR个别原子误差大该怎么调整
已经有8人回复- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 在Gaussian 09 计算中怎么做核冻结处理？
已经有6人回复
- gaussian 03Revison-E01版在Window 下不能并行计算啊!
已经有9人回复
- 【求助】请问Gaussian计算带电荷有机物时电荷和自旋多重度如何选择
已经有18人回复
- 【求助】关于Gaussian计算中的能量问题
已经有14人回复
- 【求助】用Gaussian计算分子的熵时，不同分子的熵的计算基准是一样的吗？谢谢
已经有11人回复
- 【求助】利用Gaussian一次计算多个压力温度下的热力学常数
已经有3人回复
- 【求助】Gaussian计算HOMO-LUMO问题
已经有17人回复
- 【求助】gaussian频率计算问题
已经有9人回复
- 【求助】Gaussian 03 计算错误
已经有4人回复