计算化学公社

标题: Dalton计算旋轨耦合SCF无法收敛(DIIS aborted because max DIIS iteration) [打印本页]

作者
Author: shaowenhao28 时间: 2022-10-24 01:30
标题: Dalton计算旋轨耦合SCF无法收敛(DIIS aborted because max DIIS iteration)
本帖最后由 shaowenhao28 于 2022-10-24 03:25 编辑

各位老师求助！

用Dalton计算镧系金属（Gd）配合物旋轨耦合常数时出错，SCF无法收敛。报错内容如下：

WARNING !!! DIIS aborted because max DIIS iterations reached !

  --- SEVERE ERROR, PROGRAM WILL BE ABORTED ---

@ MPI MASTER, node no.:    0
@ Reason: QC-SCF not implemented/tested for DFT

分子在gaussian中用PBE1PBE优化过，参照sobreva老师的帖子，对Gd使用了MWB53,对其他元素使用6-31G**，优化结果没有问题，波函数稳定无虚频。

分子结构用avogadro提取出来后输入到dalton中，在dalton内对Gd使用了stuttgart_rsc_1997_ecp基组（Basis=stuttgart_rsc_1997_ecp ECP=stuttgart_rsc_1997_ecp）

dalton的输入文件均附在帖子中。希望各位老师能给予帮助！

**DALTON INPUT
.RUN RESPONS
**INTEGRALS
.MNF-SO
**WAVE FUNCTIONS
.DFT
PBE1PBE
**RESPONS
.TRPFLG
*LINEAR
.SINGLE RESIDUE
.PROPRT
X1MNF-SO
.PROPRT
Y1MNF-SO
.PROPRT
Z1MNF-SO
.PRINT
3
.ROOTS
10

**END OF DALTON INPUT

ATOMBASIS
Using the mixed basis
title
Atomtypes=11 Angstrom NoSymmetry
Charge=17.0 Atoms=1 Basis=6-31G**
Cl       0.00094    -0.00018       3.29157
Charge=64.0 Atoms=1 Basis=stuttgart_rsc_1997_ecp ECP=stuttgart_rsc_1997_ecp
Gd       0.00004    -0.00000       0.68588
Charge=7.0 Atoms=1 Basis=6-31G**
N       -2.04840       0.41345    -0.40123
Charge=6.0 Atoms=3 Basis=6-31G**
C       -2.60723       1.65168    -0.55783
C       -1.89875       2.85894    -0.58683
C       -0.51119       3.04364    -0.55788
Charge=7.0 Atoms=1 Basis=6-31G**
N       0.41334       2.04824    -0.40149
Charge=6.0 Atoms=5 Basis=6-31G**
C       1.65156       2.60713    -0.55806
C       1.50753       4.02532    -0.77769
C       0.17168       4.29494    -0.77752
C       2.85885       1.89865    -0.58716
C       3.04355       0.51110    -0.55817
Charge=7.0 Atoms=1 Basis=6-31G**
N       2.04817    -0.41341    -0.40168
Charge=6.0 Atoms=5 Basis=6-31G**
C       4.29483    -0.17178    -0.77790
C       4.02518    -1.50762    -0.77804
C       2.60704    -1.65165    -0.55817
C       1.89852    -2.85894    -0.58712
C       0.51096    -3.04361    -0.55803
Charge=7.0 Atoms=1 Basis=6-31G**
N       -0.41352    -2.04821    -0.40149
Charge=6.0 Atoms=7 Basis=6-31G**
C       -1.65178    -2.60709    -0.55790
C       -1.50778    -4.02526    -0.77763
C       -0.17194    -4.29490    -0.77767
C       -2.85909    -1.89859    -0.58674
C       -3.04377    -0.51106    -0.55767
C       -4.29506       0.17184    -0.77747
C       -4.02543       1.50767    -0.77754
Charge=1.0 Atoms=12 Basis=6-31G**
H       -2.49505       3.75691    -0.72143
H       2.32445       4.71647    -0.94089
H       -0.31314       5.24885    -0.94070
H       3.75679       2.49495    -0.72193
H       5.24874       0.31301    -0.94114
H       4.71631    -2.32454    -0.94133
H       2.49481    -3.75689    -0.72185
H       -2.32473    -4.71639    -0.94076
H       0.31285    -5.24881    -0.94093
H       -3.75705    -2.49491    -0.72131
H       -5.24897    -0.31296    -0.94071
H       -4.71654       2.32461    -0.94076

作者
Author: zjxitcc 时间: 2022-10-24 02:19
本帖最后由 zjxitcc 于 2022-10-24 02:22 编辑

体系的电荷、自旋多重度是？你看看输出文件中Dalton是否正确识别为你想要的电荷和自旋多重度，这个信息千万不能错。另外，有4个C原子采用了cc-pVDZ，而非全部非金属都用了6-31G(d,p)，是写错了还是真的有必要？

作者
Author: shaowenhao28 时间: 2022-10-24 03:24

zjxitcc 发表于 2022-10-24 02:19
体系的电荷、自旋多重度是？你看看输出文件中Dalton是否正确识别为你想要的电荷和自旋多重度，这个信息千万 ...

多谢您的回复。上一版dalton输入文件有些错误，四个C采用了cc-pvdz。已在原贴中修正为6-31G**。不过这并不影响输入结果。
同时我也尝试过增加max DIIS iterations，不过SCF依旧无法收敛。

我想dalton的输出文件加在了帖子中。按照gaussian的几何优化结果，分子不含电荷，自旋为单重态。

dalton的output中显示
Total charge of the molecule 0
Spin multiplicity and 2 M_S 2 1
请问这是表明spin multiplicity 已经被识别为1了么？

作者
Author: zjxitcc 时间: 2022-10-24 10:08
本帖最后由 zjxitcc 于 2022-10-24 10:27 编辑

shaowenhao28 发表于 2022-10-24 03:24
多谢您的回复。上一版dalton输入文件有些错误，四个C采用了cc-pvdz。已在原贴中修正为6-31G**。不过这并 ...

这表示中性、二重态。这个体系中性的时候是电子数是奇数，不可能是单重态。SCF不收敛是次要问题，应先检查结构合理性、电荷和自旋多重度。

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