计算化学公社

标题: 加入calcall之后几何优化依旧报错 [打印本页]

作者
Author: wwp 时间: 2020-3-15 22:06
标题: 加入calcall之后几何优化依旧报错
各位老师我想请问一下:在opt=calcall freq b3lyp/genecp scrf=(iefpcm,solvent=acetonitrile) geom=connectivity这样条件下几何优化之后依旧报错,优化的结果我看走后一直走平了,是怎么回事呢?
(, 下载次数 Times of downloads: 30)

谢谢各位老师了.

作者
Author: liyuanhe211 时间: 2020-3-15 22:09
Displacement已经是0了，优化过程中结构不会变的。稍微调整一下结构试试

作者
Author: wwp 时间: 2020-3-15 23:06

liyuanhe211 发表于 2020-3-15 22:09
Displacement已经是0了，优化过程中结构不会变的。稍微调整一下结构试试

好的,谢谢老师

作者
Author: sobereva 时间: 2020-3-16 11:05
如果换了初猜结构后还是这个问题，改用其它优化算法，比如opt=gdiis、opt=rfo再试试，或者改用其它坐标优化，比如opt=cartesian

作者
Author: wwp 时间: 2020-3-16 11:57

sobereva 发表于 2020-3-16 11:05
如果换了初猜结构后还是这个问题，改用其它优化算法，比如opt=gdiis、opt=rfo再试试，或者改用其它坐标优化 ...

好,老师我想请问一下在冻结坐标坐标的情况下还可以用opt=cartesian这样优化吗

作者
Author: snljty 时间: 2020-3-16 12:37
其实我挺好奇是什么结构位移会比受力收敛早这么多的，方便您给一下结构文件用来学习么？或者说一下大概特征可以么？谢谢！

作者
Author: wwp 时间: 2020-3-16 15:30

snljty 发表于 2020-3-16 12:37
其实我挺好奇是什么结构位移会比受力收敛早这么多的，方便您给一下结构文件用来学习么？或者说一下大概特征 ...

好的,没有问题.

作者
Author: snljty 时间: 2020-3-16 16:09

wwp 发表于 2020-3-16 15:30
好的,没有问题.

谢谢！顺便，你这个基组有些问题，Ga的LanL2DZ太小了还不带极化函数，搭配其他原子的6-311G**不太合理。建议至少提升到LanL08。可以看看卢老师的《谈谈赝势基组的选用》http://sobereva.com/373。
这个团簇感觉有点儿小，而且边缘是不是要考虑饱和一下。我不是太了解团簇计算，建议咨询一下别人，或者看看卢老师的《使用量子化学程序基于簇模型计算金属表面吸附问题》http://sobereva.com/540。

作者
Author: wwp 时间: 2020-3-16 16:33

snljty 发表于 2020-3-16 16:09
谢谢！顺便，你这个基组有些问题，Ga的LanL2DZ太小了还不带极化函数，搭配其他原子的6-311G**不太合理。 ...

是的，算到后面是感觉团簇小了一点，我把其他原子的基组换到6-31+G**了，这样搭配可以吗？

作者
Author: snljty 时间: 2020-3-16 17:02

wwp 发表于 2020-3-16 16:33
是的，算到后面是感觉团簇小了一点，我把其他原子的基组换到6-31+G**了，这样搭配可以吗？

主要是LanL2DZ太小了，都不带极化函数，实在无法保证可靠性。

作者
Author: zjxitcc 时间: 2020-3-16 17:15
本帖最后由 zjxitcc 于 2020-3-16 17:16 编辑

取的团簇实在是太小了，连一小半层感觉都没到。这算下去意义不大啊

Ga可以用LANL2DZdp，其他原子用6-31G(d)就行了，如果有负离子考虑给负离子加弥散。

作者
Author: wwp 时间: 2020-3-16 17:19

snljty 发表于 2020-3-16 17:02
主要是LanL2DZ太小了，都不带极化函数，实在无法保证可靠性。

好，我换成LanL08

作者
Author: wwp 时间: 2020-3-16 17:19

zjxitcc 发表于 2020-3-16 17:15
取的团簇实在是太小了，连一小半层感觉都没到。这算下去意义不大啊
Ga可以用LANL2DZdp，其他原子用6- ...

嗯，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)