计算化学公社

标题: q-chem中CDFT-CI输入文件$cdft部分写法问题 [打印本页]

作者
Author: 锤子 时间: 2021-3-24 17:34
标题: 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

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

作者
Author: 喵星大佬 时间: 2021-3-24 19:55
我觉得这种东西并不应该用CDFT算

作者
Author: 锤子 时间: 2021-3-25 18:58

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

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

作者
Author: alonewolfyang 时间: 2021-3-25 20:03
Qchem这玩意有问题就写信问开发者吧

作者
Author: 锤子 时间: 2021-3-25 20:23

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

我去问一问看看。谢谢

作者
Author: 锤子 时间: 2021-4-2 10:46
问题已解决！
The issue comes from the guess generated by promolecule approach. You can add “CDFTCI_SKIP_PROMOLECULES = 1” to enforce the original constraints.
感谢各位前辈。

作者
Author: 动画片 时间: 2022-9-19 22:56
同学，我在算CDFT-CI算electron transfer integral。。不知道如何看输出文件，你的积分是咋读出来的。。谢谢你！

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