Q-Chem-oniom使用求助

柒月小鱼

大家好，小弟这里有几个Q-chem中oniom的使用问题想要求助大家，至于为什么要用Q-Chem，当然是看中了他的CDFT和ALMO版块，我使用的是我们合作者的Q-Chem，所以没有售后我主是遇到了如下几个问题
1.oniom中的电荷设置是每种原子对应与一个电荷值，但是其实同一种原子的电荷差别还是有一些的，为此我对与同一种原子的电荷值取了平均
1.1这种方法可行么？？？
（如amber力场的氮的NB类型原子对应于一个电荷，但是我在自己拟合力场的时候，为了方便有些原子就写成了一种原子，没有很细分）（gaussian的oniom是每个原子对应一个电荷，Q-Chem能和gaussian一样每个原子对应一个电荷么？）
1.2我把MM层其实是完全冻结的，不然的话我不用oniom的方法直接使用电荷嵌入的方法怎么样？？？（因为我的原子数目有点多，重新再拟合参数再调oniom的参数实在太累了，想偷懒了）
2.我现在用oniom时算正常的单点就SCF 不收敛，这算CDFT或者ALMO可怎么办啊，以下是我各种查手册后的关键词，还是不行，大家觉得我的关键词怎么样，
$rem
qm_mm_interface         oniom
force_field             amber99
user_connect            true
molden_format           true
JOBTYPE         sp
EXCHANGE        b3lyp
correlation     none
BASIS           def2-svp
mem_total       40000
mem_static      12000
SCF_PRINT       TRUE
MAX_SCF_CYCLES = 200
scf_algorithm diis
xc_grid 1
SYMMETRY false
SYM_IGNORE true
$end



希望各位路过给出些意见，拜托了

alonewolfyang

1，1.1具体力场参数部分设置(你说的每种原子所对电荷问题)看Tinker中相关力场的设置的说明，别把QChem扯进去，QChem指管总的电荷 多重度；
     1.2我没具体用QM/MM方法计算，你可以测试一下完全冻结和不冻结的区别；
2，仔细看看ONIOM的说明书相关介绍，
Analytic QM/MM gradients are available for QM subsystems described with density functional theory
(DFT) or Hartree-Fock (HF) electronic structure theory, allowing for geometry optimizations and
QM/MM molecular dynamics.
• Single-point QM/MM energy evaluations are available for QM subsystems described with most post-
HF correlated wave functions.
• Single-point QM/MM calculations are available for excited states of the QM subsystem, where the
latter may be described using CIS, TDDFT, or correlated wave function models. Analytic gradients
for excited states are available for QM/MM calculations if the QM subsystem is described using CIS.
• Single-point MM or QM/MM energy evaluations are available using periodic boundary conditions
with Ewald summation.
• Implicit solvation for both Janus QM/MM calculations as well as MM-only calculations is available
using the Polarizable Continuum Models (PCMs) discussed in Section 11.2.2.
• Gaussian blurring of MM external charges is available for Janus QM/MM calculations.
• The user may add new MM atoms types and MM parameters.
• The user may define his/her own force field.
个人建议更改一下格点大小试试，其次，你在QM里面要使用CDFT或者ALMO，我也不敢保证能顺利进行，假如你这里能跑下来，你可以尝试一下

柒月小鱼

Using fractional MO occupation numbers can improve the convergence for small-gap systems. In this approach,
the occupation numbers of MOs around the Fermi level are allowed to assume non-integer values. This
“occupation smearing” allows one to include multiple electron configurations in the same optimization,
which improves the stability of the optimization.  

柒月小鱼

alonewolfyang 发表于 2019-5-17 10:43
1，1.1具体力场参数部分设置(你说的每种原子所对电荷问题)看Tinker中相关力场的设置的说明，别把QChem扯进 ...

1.谢谢您对第一点的回复，我去看看Tinker手册
关键SCF 不收敛的问题可能是我对Pt用的全电子基组def2-svp，没有用赝势（使用赝势之后，会这样报错 Cannot use atomic densities guess with general basis）
不使用oniom的话，加赝势当然是没有问题的，  不行的话 我换台机器试试（可能这台机器的qchem安装不太好，上次就是莫名其妙的死了）
我再多多尝试把
2.对于oniom使用CDFT是没有问题的，我对于小体系亲测可行

最后的最后 谢谢您的详细的回复 比心