24小时热门版块排行榜    

查看: 729  |  回复: 10
【奖励】 本帖被评价5次,作者beefly增加金币 3.5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] gaussian03的几个已知bug

1. He原子使用ecp的BUG。

He原子和H原子一样,他们的“ECP”基组实际上是全电子基组,这是gaussian为了指定方便而设计的。但使用会发现,如果为He指定了lanl2dz或cep-31g基组,He没有分配任何基函数。导致的结果是:如果体系只含He元素(如He,He+),程序意外终止但没有任何出错信息;如果体系还含有其他元素,就会出一条不起眼的警告信息,继续计算,用户很难发现结果是错的。sdd、sddall、cep-121g则没有这个问题。

2. 容易混淆的full

manual里说,full用来取消post-HF计算的冻芯设置,对全部电子进行相关能计算。以mp2为例,full的用法为mp2=full,或mp2(full)。但是,如果不慎漏掉了等号和括号,写成了mp2 full,程序并不报错,而是进行更昂贵的MP4和BCCD(T)计算。这个用法在manual里并未提到。

3. BD算布居和原子电荷的BUG
BD(包括BCCD,BCCD(T)和BCCD(TQ))是比一般的耦合簇更精确的方法。但是,BD轨道做Mulliken或NBO会有严重问题:分子全部处于高度电离状态!

[ Last edited by beefly on 2008-6-24 at 13:38 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)


引用回帖:
Originally posted by wxjbuilder at 2008-6-24 20:26:
真的吗?孤陋寡闻了,呵呵,那样就没有16G的限制了,如果那样的话MP2+Freq也能很快就算出啦

如果在MP2 opt的时候,需要硬盘突破了16G,直接使用freq是计算不出来的,基本上是硬盘空间上限问题,需要使用freq=numer来计算才能通过
6楼2008-06-25 18:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

wxjbuilder

木虫 (知名作家)


感觉最大的bug应该是rwf文件大小上限问题,如果把这个解决了,很多大的体系就可以算了
2楼2008-06-24 19:21:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)


★★★ 三星级,支持鼓励

牛人啊,学习了
to wxjbuilder : 好像高版本gaussian (如 D01以后的)的rwf文件并不需要分割
3楼2008-06-24 19:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)


真的吗?孤陋寡闻了,呵呵,那样就没有16G的限制了,如果那样的话MP2+Freq也能很快就算出啦
4楼2008-06-24 20:26:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见