Gaussian/gview计算单个原子结合能时报错

lomon

我在计算团簇中单个原子的结合能是高斯报错，单原子out文件已上传附件，原子和团簇的关键字如图，目前我做了如下的尝试：
报错l1.exe时，提示我从第二行的gemo=connectivity开始有语法错误，我删除了第二行的geom=connectivity，并空一行再运行，报错l101.exe，我删除了第五行的geom=connectivity，并在第二行加上了的geom=connectivity，删除空行后运行，仍报错l1.exe，之后我修改第二行的内容为# HF/3-21G BOND(Phase(1,2)),仍报错l1.exe,提示我从3-21开始有语法错误。
请大佬指点这里究竟为什么出问题

zjxitcc

做计算前多思考：结构优化的本质是啥？单个原子没有键长、键角、二面角，无法做结构优化。如果想得到热力学量，直接freq即可。如果连热力学量都不需要，那就是简单的算个单点。

lomon

zjxitcc 发表于 2022-5-20 11:24
做计算前多思考：结构优化的本质是啥？单个原子没有键长、键角、二面角，无法做结构优化。如果想得到热力学 ...

谢谢您的回复。单点能我已经算完了，我现在需要计算单个原子的能量来计算平均结合能。热力学数据我现在不需要

snljty2

lomon 发表于 2022-5-20 15:54
谢谢您的回复。单点能我已经算完了，我现在需要计算单个原子的能量来计算平均结合能。热力学数据我现在不 ...

什么叫“单点能我已经算完了，我现在需要计算单个原子的能量”然后又“热力学数据现在不需要”，你这句话本来就是矛盾的，不带热力学数据的单个原子的能量就属于单点能。另外不要用CEP-121G这么非主流的基组。换成def2-SVP更合理。如果体系就这么小，直接def2-TZVP优化即可。HF/3-21G完全不用考虑。优化通常不要带Population=(Full,NBO)，不然只会多一堆你基本不看的信息。对过渡金属做NBO分析本来结果也往往不明。M06L这里也未必是最佳选择。看卢老师博文
常见的多余的和被滥用的Gaussian关键词
简谈量子化学计算中DFT泛函的选择
谈谈赝势基组的选用

lomon

snljty2 发表于 2022-5-20 16:38
什么叫“单点能我已经算完了，我现在需要计算单个原子的能量”然后又“热力学数据现在不需要”，你这句话 ...

谢谢您的回复。我把单点能和分子总结合能的概念弄混了，在网上查到的资料没看明白。基组这个因为是老师跑的时候用的就是这个基组，师兄说算单点能不要修改关键字所以我就这么用了。我按照您说的把pop这一个删去再算算看。

lomon

snljty2 发表于 2022-5-20 16:38
什么叫“单点能我已经算完了，我现在需要计算单个原子的能量”然后又“热力学数据现在不需要”，你这句话 ...

另外我想请教一个问题：像我现在需要算单个原子的能量，那么基组我可以选择跟计算团簇时不同的基组吗

sobereva

lomon 发表于 2022-5-22 13:15
另外我想请教一个问题：像我现在需要算单个原子的能量，那么基组我可以选择跟计算团簇时不同的基组吗

若要和团簇的能量求差，当然基组必须相同

lomon

sobereva 发表于 2022-5-30 06:26
若要和团簇的能量求差，当然基组必须相同

谢谢sob老师的回复，这个地方我弄明白了。