计算化学公社

标题: Qchem SF-TDDFT优化MECP输入文件中的疑惑求助 [打印本页]

作者
Author: 函数与激情 时间: 2020-7-18 14:45
标题: Qchem SF-TDDFT优化MECP输入文件中的疑惑求助
本帖最后由函数与激情于 2020-7-18 14:51 编辑

根据Q-chem 4.3 手册Example 9.6 （P462），作者给出了 SF-TDDFT优化亚乙基 ethylidene 优化S1S0 圆锥交叉的输入文件，我想请教下为何自旋多重度给的是3呢？是因为SF-TDDFT参考态是三重态的原因吗？另外根据手册，在输入的最后几行中MECP_STATE1[i,j] ，i与j分别表示自旋与优化的第j个激发态，j=0 等于SCF基态（此处直译手册），S1S0交叉的两个j不应该是0与1吗？不太了解SF-TDDFT，请教大家，欢迎讨论
The MECP between the S0 and S1 states of ethylidene is optimized using the branching-plane update method:
Example 9.6 MECP optimization for ethylidene using SF-TDDFT
$molecule
0 3
C 0.0446266041 -0.2419241370 0.3571573801
C 0.0089051507 0.6727548956 1.4605006396
H 0.9284257388 -0.1459163900 -0.2720952334
H -0.8310326564 -0.1926895078 -0.2885298629
H -0.0092388670 0.9611331703 2.4799363398
H 0.0683140308 -1.2533580302 0.7788470826
$end
$rem
jobtype opt
mecp_opt true
mecp_methods branching_plane
MECP_PROJ_HESS true !project out y vector from the hessian
GEOM_OPT_COORDS 0 !currently only works for Cartesian coordinate

method bhhlyp
spin_flip true
unrestricted true
basis 6-31G(d,p)
cis_n_roots 4
mecp_state1 [0,1]
mecp_state2 [0,2]
CIS_S2_THRESH 120
$end

作者
Author: myzpku 时间: 2020-7-19 05:32
SF-TDDFT的参考态是high-spin triplet，所以输入文件写的是0 3，TDDFT计算得到的state应该是单重态（alpha->beta激发，会有spin contamination)。

SF-TDDFT从三重态flip down得到的第一激发态对应于能量最低的单重态(S0)，而第二激发态则对应于S1，因此S0-S1 MECP用的state number是1和2

作者
Author: 函数与激情 时间: 2020-7-19 08:11

myzpku 发表于 2020-7-19 05:32
SF-TDDFT的参考态是high-spin triplet，所以输入文件写的是0 3，TDDFT计算得到的state应该是单重态（alpha- ...

谢谢你的答疑！

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