计算化学公社

标题: 请教genecp混合基组计算错误的原因 [打印本页]

作者
Author: heiheihaha 时间: 2020-5-4 22:51
标题: 请教genecp混合基组计算错误的原因
本帖最后由 heiheihaha 于 2020-5-4 22:50 编辑

各位老师好，本人接触gaussian不久，在学习了Sob老师的帖子之后尝试用gaussian09w对铂金属配合物进行几何优化。在初步优化时，尝试用HF函数，混合基组（Pt原子, 赝势基组LanL2DZ，其余原子用3-21G）,但在优化结构出现了明显错误（见附图1）。手动停止后，去掉了3-21G基组，尝试全部原子都用LanL2DZ（虽然Sob老师帖子不建议第四周期之前元素用赝势），竟然优化得到了正确的结构（虽不是最优结果）。请问各位老师，为什么在用混合基组时会出现结构优化错误，反而只用赝势基组时得到了优化结果？

输入文件：
%mem=1GB
%chk=D:\Desktop\1.chk
#p opt hf/genecp scrf=(solvent=water) scf=novaracc

Title Card Required

1 1
Pt                1.51394252 4.62879320 -4.88473076
C                -2.47367289 3.09099251 -4.80237232
C                -1.06287228 3.15498294 -4.81155310
C                -1.18042777 5.47776382 -4.84313055
C                -2.55843057 5.47701599 -4.84359694
C                -3.21628669 4.23972194 -4.82064369
C                -0.48536803 6.67349319 -4.83262394
C                0.88221720 6.58167873 -4.81370883
C                1.66224839 7.74020812 -4.76291735
C                1.00126771 8.99151402 -4.75138024
C                -0.36979699 9.06181305 -4.79447746
C                -1.13524336 7.88535670 -4.83265011
C                -2.54291094 7.89762777 -4.86678071
C                -3.24458267 6.71445979 -4.86410379
H                -4.31380552 6.73220168 -4.87592399
H                -3.06586717 8.83157930 -4.89283641
H                -2.96188823 2.13905859 -4.78202960
H                -4.28559651 4.19470177 -4.81577701
H                1.57457763 9.89340080 -4.70952736
H                -0.85709649 10.01222293 -4.79709595
H                2.72857144 7.68507155 -4.73565595
H                -0.48045248 2.25628237 -4.79521972
C                5.52113463 6.10834297 -4.81633329
C                4.10895714 6.06573310 -4.85280598
C                4.18530121 3.75578721 -4.88270624
C                5.56081793 3.73100893 -4.84446102
C                6.24099331 4.95048750 -4.80097545
C                3.45911318 2.59443054 -4.94291449
C                1.34465466 1.66668622 -5.08511789
C                1.91788510 0.37273220 -5.11599676
C                3.27659129 0.21900163 -5.04598905
C                4.08574305 1.36488057 -4.96629643
C                5.50094117 1.32152859 -4.91061421
C                6.22440717 2.48968077 -4.85178838
H                7.29360291 2.45934487 -4.81304186
H                6.01294846 0.38096895 -4.91571258
H                6.02733676 7.04994808 -4.79904108
H                7.30835432 4.97475928 -4.75763067
H                1.28621021 -0.48837052 -5.19287601
H                3.71335998 -0.75892271 -5.05358222
H                0.28235370 1.79017546 -5.14008689
H                3.54121603 6.97145800 -4.86729421
Cl                1.50817633 4.54607209 -2.60691371
Cl                1.50999419 4.71780728 -7.16157130
N                2.14500883 2.73552580 -4.98544989
N                3.49125162 4.87888208 -4.86743418
N                -0.46194119 4.34950141 -4.84141567

-C -N -H -Cl
3-21G
****
-Pt
LanL2DZ
****

作者
Author: paramecium86 时间: 2020-5-5 05:55
混合基组那个输入文件，只定义了基组，但是没有给Pt 定义赝势。结果就是原子会聚到Pt原子上。

第二种写法虽然不太好（不过其实对于这几个轻原子LanL2DZ 也是全电子的，也不比3-21G差的），但是Pt的赝势至少默认给加上了。

作者
Author: sobereva 时间: 2020-5-5 07:26
没有叫做“用HF函数”的东西，HF那叫理论方法。

用HF优化的结果没有丝毫实际价值，仔细看此文了解关于选择理论方法最基本的知识
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

除非你就是想算着玩，否则拿3-21G做优化毫无意义，对当前体系哪怕是拿这个做预优化也没有意义
仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

用混合基组结果离谱显然是你根本没给Pt设赝势，怎么正确使用赝势基组在此文里已经写得没法更详细了
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

另：强烈建议先把下文末尾列出的必读博文仔细一个一个字看了再做计算，否则你做的计算纯粹是在白浪费计算资源，结果没有任何意义
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

作者
Author: heiheihaha 时间: 2020-5-5 08:59

paramecium86 发表于 2020-5-5 05:55
混合基组那个输入文件，只定义了基组，但是没有给Pt 定义赝势。结果就是原子会聚到Pt原子上。

第二种写 ...

谢谢老师，我明白了~

作者
Author: heiheihaha 时间: 2020-5-5 09:01

sobereva 发表于 2020-5-5 07:26
没有叫做“用HF函数”的东西，HF那叫理论方法。

用HF优化的结果没有丝毫实际价值，仔细看此文了解关于选 ...

谢谢sob老师的建议~

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