给Znpc 分子加赝势优化，分子散架

好多于

在优化 Znpc 的过程中，在B3LYP方法上，对C，H，N原子选用6-31G*，给 Zn 加赝势LANL2DZ，结果总是出现距离Zn最近的四个 N原子 向中心的 Zn 靠拢，造成分子散架从而难以达到收敛限。
之后，用相同的输入文件，去掉赝势基组，优化结果收敛。

请问各位老师，这种情况可能由什么原因导致？

zjxitcc

一般是由于自定义基组格式写错了不信可以上传输入文件（指上传gjf文件，不是指将其复制粘贴在帖子里）

好多于

zjxitcc 发表于 2022-12-22 20:54
一般是由于自定义基组格式写错了不信可以上传输入文件（指上传gjf文件，不是指将其复制粘贴在帖子里）

附上了

chands

好多于 发表于 2022-12-22 21:00
附上了

将gen改成genecp并加上赝势。
去掉pop=full
建议开对称性。

LittlePupil

除了把“gen”改成“genecp”外，也可以保留“gen”但加上“Pseudo=Read”；
设置赝势时不止要指定赝势基组，还要指定赝势，参看：
《详解Gaussian中混合基组、自定义基组和赝势基组的输入》
http://sobereva.com/60；
如果有特殊需求可以保留“pop=full”，但依然不建议优化时用这个，会使输出文件巨大。
最后，给的内存未免太小了，虽然Gaussian的普通DFT计算对内存需求不大。

好多于

chands 发表于 2022-12-22 21:23
将gen改成genecp并加上赝势。
去掉pop=full
建议开对称性。

好的谢谢。不过我这个后边要进行扫描，所以还是留着nosymm把

好多于

LittlePupil 发表于 2022-12-22 21:32
除了把“gen”改成“genecp”外，也可以保留“gen”但加上“Pseudo=Read”；
设置赝势时不止要指定赝势基 ...

好的谢谢，已经在看帖的路上了

sobereva

好多于 发表于 2022-12-22 21:38
好的谢谢。不过我这个后边要进行扫描，所以还是留着nosymm把

这和做不做扫描没有任何联系
什么时候才需要用nosymm，看
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297（http://bbs.keinsci.com/thread-1355-1-1.html）
你这个体系明显有对称性，让Gaussian利用对称性算得又快又容易收敛

好多于

sobereva 发表于 2022-12-23 13:31
这和做不做扫描没有任何联系
什么时候才需要用nosymm，看
谈谈Gaussian中的对称性与nosymm关键词的使用 ...

好的谢谢老师回复。我需要将分子固定在xy平面，后期需要做z（垂直分子平面）方向的拉曼强度，还要做分子的拉曼图像。所以为了后期方便起见，还是加了nosymm，

sobereva

好多于 发表于 2022-12-24 19:14
好的谢谢老师回复。我需要将分子固定在xy平面，后期需要做z（垂直分子平面）方向的拉曼强度，还要做分子 ...

对于纯平面体系，大多数情况标准朝向本来就是XY平面，所以加nosymm依然是多余的。即便最后发现实际被置于了诸如YZ平面上，也可以之后再用下文提到的功能令其精确摆在XY平面上

Multiwfn中非常实用的几何操作和坐标变换功能介绍
http://sobereva.com/610（http://bbs.keinsci.com/thread-24674-1-1.html）

好多于

sobereva 发表于 2022-12-25 23:58
对于纯平面体系，大多数情况标准朝向本来就是XY平面，所以加nosymm依然是多余的。即便最后发现实际被置于 ...

收到，谢谢老师的回帖