计算化学公社

标题: 求教Qchem软件里，如何对分子复合体系在CDFT下进行结构优化？ [打印本页]

作者
Author: lczs 时间: 2023-3-29 21:56
标题: 求教Qchem软件里，如何对分子复合体系在CDFT下进行结构优化？
我想对分子复合体系在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

作者
Author: wzkchem5 时间: 2023-3-30 00:26
检查一下最后一步的结构是否合理。

作者
Author: lczs 时间: 2023-3-30 09:42

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

作者
Author: lczs 时间: 2023-3-30 11:18
有大佬能看看嘛

作者
Author: wzkchem5 时间: 2023-3-30 16:02

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

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

作者
Author: lczs 时间: 2023-3-30 16:16

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

好滴我先试试

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