求教Qchem软件里，如何对分子复合体系在CDFT下进行结构优化？

lczs

我想对分子复合体系在CDFT下进行结构优化，用的是Qchem软件。下面是我的关键字，但总是优化失败。我尝试了许多次，改了收敛范围，都不行。我检查了输出文件，发现离收敛标准很远时都error了，根本没有好好做结构优化。
我用google查询，去Qchem的手册找，但没有相关的信息。请问是根本做不了这个结构优化，还是关键字有误呢？

下面是关键字和错误输出信息。



$rem
symmetry false
sym_ignore  true
JOBTYPE opt
METHOD B3LYP
DFT_D D3_BJ
BASIS  6-31+G(d)
SCF_CONVERGENCE 8
THRESH 12
MAX_SCF_CYCLES  200
SCF_PRINT TRUE
XC_GRID  000120000194
UNRESTRICTED TRUE
CDFT TRUE
CDFT_THRESH  7
CDFT_PREDIIS TRUE
CDFT_BECKE_POP TRUE
SOLVENT_METHOD PCM
$end

$pcm
Theory CPCM
Method SWIG
Solver Inversion
HeavyPoints 194
HPoints 194
Radii Bondi
vdwScale 1.2
$end

$solvent
Dielectric 78.39
$end

$cdft
1
1.0  46 47
1
1.0  46 47  s
$end



Q-Chem fatal error occurred in module scfman/scfman.C, line 5612:

SCF failed to converge


Please submit a crash report at q-chem.com/reporter


Error in Q-Chem run part 1
remove work dirs /tmp/qchem-10970/qchem184218.0 -- /tmp/qchem-10970/qchem184218.-1
rm -rf /tmp/qchem-10970/qchem184218
Error in the Q-Chem run

wzkchem5

检查一下最后一步的结构是否合理。

lczs

wzkchem5 发表于 2023-3-30 00:26
检查一下最后一步的结构是否合理。

非常不合理，结构完全乱了，什么都不是。
Step 14        :
        Energy is -1277.0906043898
                            Maximum         Tolerance       Converged?
         Gradient           6.60667392e+00        3.00000000e-04        false
         Displacement       2.33084841e-01        1.20000000e-03        false
         Energy change      2.23539534e-01        1.00000000e-06        false

lczs

有大佬能看看嘛

wzkchem5

lczs 发表于 2023-3-30 02:42
非常不合理，结构完全乱了，什么都不是。
Step 14        :
        Energy is -1277.0906043898

那把结构调整成合理的，继续跑

lczs

wzkchem5 发表于 2023-3-30 16:02
那把结构调整成合理的，继续跑

好滴 我先试试