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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【原创】做重元素post－HF计算的注意了，Gaussian的一个BUG！
71/1返回列表
查看: 1422  |  回复: 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的回帖

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的回帖

recoli

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

fwenling0506

木虫 (小有名气)
拜读一下！
回复此楼
4楼2011-04-16 08:09:41
已阅   回复此楼   关注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的回帖

yongleli

木虫 (正式写手)
引用回帖:
Originally posted by beefly at 2011-05-11 21:41:17:
最近做了更多的测试，发现在计算镧系、锕系体系时，不会有f轨道部分冻结部分活性的问题，只有算5d、6d金属时存在这个BUG。但是，如果你用的基组已经把Hg的4f电子放进了ECP，就不会有问题了。

谢谢！
回复此楼
7楼2011-05-12 12:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 beefly 的主题更新
71/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料284求调剂，一志愿郑州大学英一数二 +5 想上岸的土拨鼠 2026-02-28 5/250 2026-03-01 11:23 by gaoxiaoniuma
[考研] 0856材料专业298分有科研经历 硕士研究生调剂自荐信 +3 zyf上岸 2026-03-01 3/150 2026-03-01 11:21 by gaoxiaoniuma
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 高分子化学与物理调剂 +5 好好好1233 2026-02-28 9/450 2026-03-01 10:59 by fengyu211
[考研] 材料化工调剂 +5 今夏不夏 2026-03-01 5/250 2026-03-01 10:46 by ms629
[考研] 291分工科求调剂 +7 science饿饿 2026-03-01 8/400 2026-03-01 10:43 by sunny81
[考研] 0856化工专硕求调剂 +5 董boxing 2026-03-01 5/250 2026-03-01 10:41 by sunny81
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 272求调剂 +4 材紫有化 2026-02-28 4/200 2026-03-01 09:20 by L135790
[考研] 材料类求调剂 +7 wana_kiko 2026-02-28 7/350 2026-03-01 07:55 by ms629
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭