| 查看: 1433 | 回复: 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，总能量当然是不一样的。如果算键能（分子总能量减去原子总能量），结果可想而知。 |
回复此楼» 收录本帖的淘帖专辑推荐
 基组的介绍总结 |
» 猜你喜欢
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有146人回复
-
304求调剂
已经有5人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
304求调剂
已经有6人回复
-
304求调剂
已经有2人回复
-
一志愿985中药学306，求调剂
已经有0人回复
-
河南一本高校化学，材料，化工专业招收研究生
已经有0人回复
-
大连海事大学船舶洁净能源研究中心2026年博士研究生招生启事
已经有10人回复
-
2026年湘潭大学接收化学工程与技术相关方向学术型硕士研究生调剂生
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
