| 查看: 1308 | 回复: 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，总能量当然是不一样的。如果算键能（分子总能量减去原子总能量），结果可想而知。 |
回复此楼» 收录本帖的淘帖专辑推荐
 基组的介绍总结 |
» 猜你喜欢
七嗪类物质合成求助
已经有0人回复
-
软包
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有119人回复
-
软包电池的N/P是多少合适
已经有0人回复
-
招贤纳博（已结束）
已经有28人回复
-
表征测试机构后台私聊我
已经有0人回复
-
Coordination Chemistry Reviews
已经有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人回复
2楼2011-04-14 02:06:31
3楼2011-04-14 05:07:06
4楼2011-04-16 08:09:41
5楼2011-05-11 21:22:41
6楼2011-05-11 21:41:17
7楼2011-05-12 12:19:03