应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【原创】做重元素post－HF计算的注意了，Gaussian的一个BUG！
51/1返回列表
查看: 1466  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价1次，作者beefly增加金币 0.8 个
本帖产生 1 个 QC强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

beefly

专家顾问 (职业作家)

[资源] 【原创】做重元素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，总能量当然是不一样的。如果算键能（分子总能量减去原子总能量），结果可想而知。
回复此楼

» 收录本帖的淘帖专辑推荐

基组的介绍总结

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-04-13 23:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)
LZ很犀利，感谢分享
回复此楼
3楼2011-04-14 05:07:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

fooo

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by beefly at 2011-04-13 23:05:00:
在Gaussian所有的版本中，默认的冻芯设置一旦遇到f轨道占据的情况，很可能发生问题。例如汞化合物的计算，Gaussian都倾向于把4f(-3)和4f(+3)分量的轨道冻结，只对4f(-2,-1,0,+1,+2)轨道考虑相关能。如果用来优化结 ...

学习了。最近算的体系中就遇到过问题，不知道是不是LZ说的这种情况，反正结合能很离谱，这几天一直在考虑这事，正找原因呢。
赞一下(1人) 回复此楼
2楼2011-04-14 02:06:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

引用回帖:
Originally posted by beefly at 2011-04-13 23:05:00:
在Gaussian所有的版本中，默认的冻芯设置一旦遇到f轨道占据的情况，很可能发生问题。例如汞化合物的计算，Gaussian都倾向于把4f(-3)和4f(+3)分量的轨道冻结，只对4f(-2,-1,0,+1,+2)轨道考虑相关能。如果用来优化结 ...

哦！

最近做的体系恰好含有Hg。
请问，用pseudo potential会出现这个问题吗？
赞一下 回复此楼
5楼2011-05-11 21:22:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)
引用回帖:
Originally posted by yongleli at 2011-05-11 21:22:41:
哦！

最近做的体系恰好含有Hg。
请问，用pseudo potential会出现这个问题吗？

最近做了更多的测试，发现在计算镧系、锕系体系时，不会有f轨道部分冻结部分活性的问题，只有算5d、6d金属时存在这个BUG。但是，如果你用的基组已经把Hg的4f电子放进了ECP，就不会有问题了。
赞一下 回复此楼
6楼2011-05-11 21:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学求调剂 +12 喽哈加油 2026-04-14 14/700 2026-04-15 23:49 by 陈皮皮
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 291分调剂 +11 上岸小莹加油 2026-04-09 12/600 2026-04-15 21:55 by noqvsozv
[考研] 22408 312求调剂 +20 门路摸摸 2026-04-14 21/1050 2026-04-15 20:39 by cuisz
[考研] 307中医考研调剂 +4 于以采蘩 2026-04-14 4/200 2026-04-15 19:14 by AN流800
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[基金申请] 2026 WR青拔 +3 冬日阳光CAS 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 085410 273求调剂 +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by 逆水乘风
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 本人女孩 +7 吼吼， 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
信息提示
关闭
请填处理意见
关闭