Excitation analyze based on molecular fragments in TDDFT

有时，我们希望将某个TDDFT激发态到底由分子片内或分子片段之间的激发所贡献，因此，可以利用通过定义分子片段，将分子轨道（CMO，LMO或FLMO）指定到分子片段，再分析TDDFT的激发态。需要注意的是，CMO通常比较弥散，所以分子片段激发分析是近似的。

BDF的输入域&DATABASE中可以定义分子片段, 例如下面的例子，对于分子C8H10我们定义了四个分子片段，分子片段可以重叠，也可以不重叠，这里我们使用的有重叠的分子片段。我们希望使用CMO做分子片段激发分析

$COMPASS 
Title
 CH2 Molecule test run, cc-pvqz 
Basis
 6-31G
Geometry
  c   0.5833330000  0.0   0.0000000000
  c   1.9203330000  0.0   0.0000000000
  h   0.0250410000  0.0  -0.9477920000
  h   0.0250620000  0.0   0.9477570000
  h   2.4703130000  0.0  -0.9525920000
  c   2.6718330000  0.0   1.3016360000
  c   4.0088330000  0.0   1.3016360000
  h   2.1218540000  0.0   2.2542280000
  h   4.5588130000  0.0   0.3490440000
  c   4.7603330000  0.0   2.6032720000
  c   6.0973330000  0.0   2.6032720000
  h   4.2103540000  0.0   3.5558650000
  h   6.6473130000  0.0   1.6506800000
  c   6.8488330000  0.0   3.9049090000
  c   8.1858330000  0.0   3.9049090000
  h   6.2988540000  0.0   4.8575010000
  h   8.7441260000  0.0   4.8527010000
  h   8.7441050000  0.0   2.9571520000
End geometry
Check
Skeleton
Group
 C(1)
$END

$XUANYUAN
Direct
Schwarz
$END

$SCF
RHF
iprtmo
 2
threshconverg
 1.d-8 1.d-6
$END

# 如果需要利用LMO做分子片段激发分析，打开这下面两个注释行。
#$localmo  
#$end


$TDDFT
imethod
1
isf
0
itda
1
idiag
1
nexi
10
AOKXC
crit_e
 1.d-5
fraganalyze
$END

&DATABASE
nfragment 4
fragment 1  9
 1 2 3 4 5 6 7 8 9
fragment 2 12
 1 2 4 5 6 7 8 9 10 11 12 13
fragment 3 12
 6 7 8 9 10 11 12 13 14 15 16 18 
fragment 4 9
 10 11 12 13 14 15 16 17 18 
&END

输出结果如下

 No.     1    w=      5.5290 eV     -308.3680814226 a.u.  f= 2.2716   D<Pab>= 0.0000   Ova= 0.8714
      CV(0):    A(  28 )->   A(  31 )  c_i:  0.2111  Per:  4.5%  Fragment (   3) to Fragment (   3)  Oai: 0.9039
      CV(0):    A(  29 )->   A(  30 )  c_i: -0.9577  Per: 91.7%  Fragment (   2) to Fragment (   2)  Oai: 0.8783

 No.     2    w=      7.5222 eV     -308.2948336397 a.u.  f= 0.0000   D<Pab>= 0.0000   Ova= 0.7865
      CV(0):    A(  27 )->   A(  31 )  c_i: -0.1334  Per:  1.8%  Fragment (   3) to Fragment (   3)  Oai: 0.8432
      CV(0):    A(  28 )->   A(  30 )  c_i: -0.7328  Per: 53.7%  Fragment (   3) to Fragment (   2)  Oai: 0.8135
      CV(0):    A(  29 )->   A(  31 )  c_i:  0.6100  Per: 37.2%  Fragment (   2) to Fragment (   3)  Oai: 0.7726