分子内转换速率计算

MOMAP 软件包中分子内转换速率主要由程序 TVCORF_IC [1] （计算关联函数），TVSPEC_IC (计算无辐射谱)完成。此部分计算需要 Dushin计算所得的 *.dus 文件以及 get-nacme 计算所得的 nacme.out 作为输入。编写 momap.inp 文件，并依次执行上述两个程序即可得到结果。实际计算中为了进行并行运算，我们常以 job 文件定义参数生成 momap.inp 文件，并控制程序的运行。内转换速率的计算实例可见 examples/DSB/ic，examples/6T/ic，examples/azulene/ic，在计算实例中运行 命令以及写入 run 文件中，用户仅需 ./run 即可触发程序。

以发光分子 DSB 为例，通过程序 get-nacme 即可读取 Gaussian 计算所得的跃 迁电场以及振动文件而后求得电子耦合项，计算实例请见 examples/DSB/nacme。 为了简化操作步骤，目前 evc 程序也能直接计算 nacme，仅需将 Gaussian 计算所得的跃迁电场文件并在 evc.inp 中开启 nacme 计算即可，计算实例请见 examples/6T/evc。而后仅需修改 job 文件中生成 momap.inp 文件部分即可。job 文件中各参数解释如下：

&control
  DUSHIN      = .t.          # 是否开启 Duschinsky 转动效应的计算
  DUSH12      = 2            # 选择以激发态振动信息进行计算
  D12         = 2            # 选择以激发态振动信息进行计算
  Temp        = 298          # 定义温度
  tmax        = 3000.0 fs    # 定义关联函数积分时间
  dt          = 0.01 fs      # 定义关联函数积分步长
  Ead         = 0.1180424547 # 绝热激发能
  EDMA        = 1 debye      # 吸收(基态)跃迁偶极矩
  EDME        = 1 debye      # 发射(激发态)跃迁偶极矩
  FreqScale   = 1.0          # 频率因子
  logFile     = "log"        # 输出文件文件名
  DSFile      = "evc.cart.dat"       # 定义读取的 evc 文件文件名
  CoulFile    = "nacme.out"  # 定义读取的电子耦合文件文件名
  FtFile      = "ft.dat"     # 定义关联函数文件文件名
/
&gauss
  isgauss      = .f.         # 定义是否开启关联函数展宽
  BroadenType  = lorentzian  # 定义展宽函数类型 (lorentz / gauss)
  Broadenfunc  = time        # 定义展宽物理量 (freq / time)
  lifetime     = 1.0 ps      # 定义展宽大小
  GFile        = "gs.dat"    # 定义展宽输出文件
/
&spectrum
  NScale       = 4
  Emax         = 0.3 au      # 定义光谱频率范围上限
  FoFile       = "fo.dat"    # 定义光谱输出文件
/

计算完成后先确认关联函数是否收敛，即将 ft.dat 画图。确认关联函数收敛后， 可在 log 文件末端读取内转换速率，在 fo.dat 文件中读取内转换谱。

[1]	推荐实际计算中采用并行版TVCORF_IC_para