电子振动耦合计算

MOMAP软件包含evc程序（electron-vibration coupling，EVC），用于计算Duschinsky转动、振动分析等跃迁信息 [1]。 evc程序 [2] 可读取Gaussian、Turbomole、Chemshell、Dalton、MOLPRO、DFTB、MOPAC程序输出的分子振动频率、力常数矩阵，并同时在内坐标以及直角坐标系下，计算分子跃迁发生初末态间的模式位移、黄昆因子、重整能以及Duschinsky转动矩阵。 最为简单的evc输入文件(evc.inp)可写为 [3]：

&evc
    ffreq(1)    = "Groundstate_vibration.log"
    ffreq(2)    = "Excitedstate_vibration.log"
/

而后运行 bin/evc evc.inp 即可 。 Dushin程序的实例可见examples/DSB/evc，examples/6T/evc 以及 examples/azulene/evc，在计算实例中运行命令以及写入 run 文件中，用户仅需 ./run 即可触发程序。目前evc 程序除其基本功能外，还添加了重整能向内坐标投影，计算同位素效应等功能，详细的高级参数定义请见 examples/6T/evc/evc.inp文件，用户可根据计算需要自行编写evc.inp文件。

evc程序的输出文件主要包括evc.out，evc.cart.dat 以及 evc.dint.dat。其中 evc.cart.dat 为利用直角坐标系计算得到的模式位移、黄昆因子、重整能以及 Duschinsky 转动矩阵的结果；而evc.dint.dat 为内坐标计算模式位移、黄昆因子、重整能，直角坐标计算Duschinsky转动矩阵的结果。请用户比较 evc.cart.dat 以及 evc.dint.dat 文件内的重整能，若重整内相差不大 (< 1000 cm-1) 请用文件 evc.cart.dat 进行后续计算。若重整能相差较大(一般情况 evc.cart.dat 大于 evc.dint.dat)，请用文件 evc.dint.dat 进行后续计算。

[1]	MOMAP 软件包中，还包括 evc 程序，其为读取 Turbomole 的输出并以直角坐标系计算 Duschinsky 转动矩阵以及正则模式位移的程序。

[2]	evc 程序改写自 Jeffrey R. Reimers 教授编写的 dushin 程序，(J. Chem. Phys. 115, 9103-9109, 2001)

[3]	以 Gaussian 软件计算 evc 时需要提供 *.log 以及 *.fchk文件，如/example/DSB/evc/