计算化学公社

标题: oniom结合混合基组计算不能运行 [打印本页]

作者
Author: Yoghurt 时间: 2019-4-13 14:24
标题: oniom结合混合基组计算不能运行
1. 用oniom模型计算MOF材料的激发态，选择一个分子为high层，其余部分为Low层，并且考虑混合基组，Eu原子用SDD，其余原子通6-31g*,不知道什么原因总是读完结构就停止，请问如何解决？（输入和输出文件见附件）2. 另外，在这种选择下有化学键被斩断，输入文件中有原子显式L H 两种分层，是否因为这个原因导致的问题呢?

作者
Author: sobereva 时间: 2019-4-13 20:22
对于杂化泛函用/Auto无意义，密度拟合在Gaussian里只能对纯泛函使用
写了Embedcharge，但是既没自己定义low层原子电荷，也没用诸如qeq自动计算，等于根本没考虑高层和低层之间的静电相互作用

作者
Author: Yoghurt 时间: 2019-4-13 20:32

sobereva 发表于 2019-4-13 20:22
对于杂化泛函用/Auto无意义，密度拟合在Gaussian里只能对纯泛函使用
写了Embedcharge，但是既没自己定义lo ...

首先谢谢老师的指点。
按照您的意见，我把命令行改成了# ONIOM(b3lyp/Genecp TD=(NStates=10):UFF=QEq)=Embedcharge geom=connectivity
但是运行仍然很快就结束了，输出文件末尾显式
LnStep:  EOld=  0.000000000D+00 ENew= -1.846020180D+00 DE= -1.846020180D+00
      DOld= -2.467708985D+00 DNew= -1.179085768D+00
LnStep:  A = 1.192587525D-03  4.524560691D-02  6.437153149D-01 -1.857331582D+00  0.000000000D+00
LnStep:  Values at new, mid, old =  -7.620067244D-01  0.000000000D+00  1.084013455D+00
      derivs at new, mid, old =  -1.179085768D+00 -1.857331582D+00 -2.467708985D+00
                              = 1.426745213D+00  1.287430630D+00  1.155271572D+00
LnStep:  XStart =  -5.000000000D-01  0.000000000D+00  5.000000000D-01
LnStep:  X,  F at min = 1.266128129D+00 -1.224792294D+00
      DF, DDF at min = 2.827599266D-16  1.654092808D+00
      FVal =  -2.308805749D+00
ISolve=1 DoInv=T D.G= -2.47D+00 1/2 D.H.D=  5.80D-01 DEPred= -1.89D+00 XLin=  7.66D-01
请问还有可能是哪方面的问题呢？谢谢~

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