24小时热门版块排行榜    

查看: 3090  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

pleasant1789

木虫 (初入文坛)

[求助] G09中计算AIM计算过渡金属原子和氢原子间作用的问题

以前从没进行过AIM计算,这次投的文章里审稿人提出用AIM计算一下,应该是想用AIM计算证实过渡金属原子和氢原子间存在相互作用。

从小木虫上现有的与AIM计算相关的文章里看,http://muchong.com/bbs/viewthread.php?tid=1767755
这个帖子里给出了在Gaussian中计算AIM的方法,用AIM=all关键字实现,然而我计算的时候出错,出错信息为“ Error termination in NtrErr:
NtrErr Called from FileIO.”
换用较低的基组,则可以进行到输出“ I. ATTRACTORS”的部分,然后“FATAL ERROR: DUPLICATE CRITICAL POINTS FOUND:   50  41
”609错误结束
计算的体系是IMes配位的Ni与另外一个分子之间的配位,重原子数34。对Ni原子使用了lanl2dz的ECP。

搜解决原因的时候看到各种说法
1. 计算软件不合适?有人提到Gaussian的AIM有问题,故应该用专业的AIM软件计算,然而免费的AIMPac需要编译(这个太复杂,搞不定);AIMall免费的只能算10个以下的原子,AIM2000也是商业软件,故都不可行,平时并不用AIM分析所以也不可能因为这个去购买商业软件。
2.“芳香环片段高斯不可以做,会出错,高斯不可以算环状结构,可以使用aimall做 ”出自 http://muchong.com/bbs/viewthread.php?tid=1781656
是否为真?如果是如此那就真是没办法了,是否有相关的文献说明这点呢。
3.也有提出ECP不能解决,但是http://muchong.com/bbs/viewthread.php?tid=2511503 提到“所有AIM软件都能支持ECP计算产生的wfn文件”,而且G09 的B01版本其实是可以产生可用的wfx文件
4.有人提出高斯中计算AIM不可信?

我的问题总结来说就是
1. 能否用Gaussian09计算这个34个重原子含过渡金属Ni,且使用ECP的体系?
2. 如果不能,是否有可用的免费软件且不太麻烦的可以进行计算?
3.具体如何操作和读取相关信息?

问题非常琐碎和麻烦,望哪位同仁能够不吝回答,非常感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pleasant1789

木虫 (初入文坛)

引用回帖:
2楼: Originally posted by sobereva at 2011-12-16 13:40:15:
无论高斯自带的AIM模块能不能正常分析某个体系,我都不建议用,这个模块年久失修,速度、可靠性、稳定性都比较低。
如果不需要做盆积分,使用Multiwfn做AIM分析是最佳选择,操作与AIM2000一样方便,而速度比AIMA ...

非常感谢!
我只需要一个很简单的功能,能通过AIM说明金属和H原子之间有作用就可以了(其实从键长是可以看出来的,不过审稿人要求做AIM所以不得不用一下)。于是我去看看手册说明怎么用

另外想问一下,现在那个讲Multiwfn的文章还是在ASAP状态,没有页码,那么我引用的时候是否还是引用网址?或者等我定稿的时候卷期号已经出来了,引用文章就可以了吗?
3楼2011-12-16 14:57:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★ ★
gmy1990(金币+3): 2011-12-16 18:52:25
无论高斯自带的AIM模块能不能正常分析某个体系,我都不建议用,这个模块年久失修,速度、可靠性、稳定性都比较低。
如果不需要做盆积分,使用Multiwfn做AIM分析是最佳选择,操作与AIM2000一样方便,而速度比AIMALL还要快得多。
对于过渡金属使用赝势的情况,对于LANL2这样小核赝势,直接用各种AIM软件去分析价层性质都是可以的,例如Journal of Computational Chemistry, Vol. 18, No. 3, 416]429 (1996)文中给出了确切的结论(SC代表small core):
For transition metals, the electron density obtained
from a SC-ECP calculation is good enough
to reproduce not only the correct topology of the
metal]ligand bonds but the values of the local
properties as well

在G09 B01及之后版本中支持wfx格式,会自动将事先拟合好的内核密度写进去,使得分析结果和全电子计算十分一致,见J. Phys. Chem. A 2011, 115, 12879–12894。目前只有Multiwfn和AIMALL支持wfx格式下的拓扑分析。而至于G09自带的AIM模块,我虽然没试过,但是估计它自己都不支持内核密度的处理。
2楼2011-12-16 13:40:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by pleasant1789 at 2011-12-16 14:57:03:
非常感谢!
我只需要一个很简单的功能,能通过AIM说明金属和H原子之间有作用就可以了(其实从键长是可以看出来的,不过审稿人要求做AIM所以不得不用一下)。于是我去看看手册说明怎么用

另外想问一下,现在那 ...

引用doi号不就行了
我就是这样引用的
4楼2011-12-16 15:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pleasant1789

木虫 (初入文坛)

2楼的不是应助回帖我没法评分啊orz
怎么办,要通知版主处理吗
5楼2011-12-16 15:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见