分子磷光及系间窜越速率计算

MOMAP 软件包中分子磷光及系间窜越速率目前仅发布了包含一阶项的计算程序,主要由程序 TVCORF_SPEC [1] (计算关联函数),TVSPEC_SPEC (计算光谱以及速率) 完成。此部分计算需要 evc 计算所得的 evc.*.dat 文件作为输入。 编写 input 文件,并依次执行上述两个程序即可得到结果。实际计算中为了进行 并行运算,我们常以 job 文件定义参数生成 input 文件,并控制程序的运行。磷 光及系间窜越速率的计算实例可见 examples/Irppy3,在计算实例中运行命令以及 写入 run 文件中,用户仅需 ./run 即可触发程序。

以磷光光分子 Ir(ppy)3 为例,一般情况下仅需修改 job 文件中生成 input 文件 部分即可。job 文件中各参数解释如下:

&control
  DUSHIN      = .t.          # 是否开启 Duschinsky 转动效应的计算
  HERZ        = .f.          # 是否开启 Herzberg-Teller 效应的计算
  Temp        = 298 K        # 定义温度
  tmax        = 655.36 fs    # 定义关联函数积分时间
  dt          = 0.01 fs      # 定义关联函数积分步长
  Ead         = 0.1008948 au # 绝热激发能
  EDMA        = 1.0 debye    # 吸收(基态)跃迁偶极矩
  EDME        = 0.29215663 debye # 发射(激发态)跃迁偶极矩
  isISC       = .t.              # 是否开启系间窜越速率的计算
  Hso         = 217.825327 cm-1  # 旋轨耦合常数
  FreqScale   = 1.0          # 频率因子
  logFile     = "log"        # 输出文件文件名
  DSFile      = "evc.cart.dat"       # 定义读取的 evc 文件文件名
  FtFile      = "ft.dat"     # 定义关联函数文件文件名
/
&gauss
  isgauss      = .f.         # 定义是否开启关联函数展宽
  BroadenType  = lorentzian  # 定义展宽函数类型 (lorentz / gauss)
  Broadenfunc  = frequency   # 定义展宽物理量 (freq / time)
  FWHM         = 100 cm-1    # 定义展宽大小
  GFile        = "gsfile"    # 定义展宽输出文件
/
&spectrum
  NScale       = 10
  Emin         = -0.3 au     # 定义光谱频率范围下限
  Emax         = 0.3 au      # 定义光谱频率范围上限
  dE           = 0.00001 au  # 定义光谱频率步长
  FoFile       = "fo.dat"    # 定义光谱输出文件
  FoSFile      = "spec.dat"  # 定义光谱输出文件
/
[1]推荐实际计算中采用并行版TVCORF_SPEC_para