q-chem中CDFT-CI输入文件$cdft部分写法问题

锤子

我使用q-chem程序中的CDFFT-CI方法来计算两个透热态之间的电子耦合，其中两个透热态分别是三重态局域在第一个片段（1-23）上，和三重态局域在第二个片段（24-51）上。命令部分和限制部分是这么写的：

$rem
symmetry off
sym_ignore = true
method = b3lyp
basis = 6-31g*
unrestricted = true
scf_convergence = 8
max_scf_cycles = 200
cdftci = true
cdftci_print = 2
cdft_thresh = 7
$end

$cdft
2.0
1.0   1 23 s
0.0
1.0   24 51 s

-------------
0.0
1.0   1  23 s
2.0
1.0   24 51 s
$end


但得到的两个透热态的becke布局分析显示都是三重态集中在第一个片段上。得到如下形式的overlap矩阵：
               1               2
    1    1.0000000001    1.0000000001
    2    1.0000000001    1.0000000001

很明显两个态是一样的，都是三重态局域在第一个片段上。
我想请教一下我的输入部分哪里存在问题。十分感谢各位前辈的指导。

喵星大佬

我觉得这种东西并不应该用CDFT算

锤子

喵星大佬 发表于 2021-3-24 19:55
我觉得这种东西并不应该用CDFT算

我使用的是D-A这种体系，看一下分子内的能量转移，现在是先限制自旋局域在不同的片段上，学习一下这种方法的使用。下一步要限制分子内某个原子上的电荷，重点是后面限制电荷后的两个态的耦合，使用CDFT-CI也是为了后面能够限制电荷。但是我第一步都没能成功。
原理上这种方法是应该能算出来两个态之间的耦合的，但我现在输入文件存在问题，导致结果明显异常。希望能够有前辈帮我找找原因在哪，我试了好多种限制的方式都没能成功。

alonewolfyang

Qchem这玩意有问题就写信问开发者吧

锤子

alonewolfyang 发表于 2021-3-25 20:03
Qchem这玩意有问题就写信问开发者吧

我去问一问看看。谢谢

锤子

问题已解决！
The issue comes from the guess generated by promolecule approach.  You can add “CDFTCI_SKIP_PROMOLECULES = 1” to enforce the original constraints.
感谢各位前辈。

动画片

同学，我在算CDFT-CI算electron transfer integral。。不知道如何看输出文件，你的积分是咋读出来的。。谢谢你！