怎么处理用momap计算固相下共晶磷光无辐射速率报错问题

Yangd

共晶分子式为C18H14Br2N2O6·CHCl3，其中溶剂分子是三氯甲烷，是一个无序性溶剂分子，在构建oniom模型时已经将无序性处理好。溶剂分子和分子的数量比为1：1。在构建oniom模型时，我把一个溶剂分子和另一个分子同时选为high层。输入文件我使用的是GAUSSIAN计算的S0和T1的.log和.fch文件，使用m062x泛函和def2svp基组，无虚频。并且我还计算了把共晶中溶剂分子去除的固相的磷光无辐射速率，无溶剂分子晶体的固相下的磷光无辐射速率在momap的计算中无报错。

构建的oniom如下图


输入文件如下


其中.inp文件如下
do_evc            = 1
do_spec_tvcf_ft   = 1
do_spec_tvcf_spec = 1
do_ic_tvcf_ft     = 0
do_ic_tvcf_spec   = 0
do_isc_tvcf_ft    = 1
do_isc_tvcf_spec  = 1
do_spec_sums      = 0

&evc
  ffreq(1)      = "s0.log"
  ffreq(2)      = "t1.log"
  sort_mode     = 1
  proj_reorg    = .t.
/

&spec_tvcf
  DUSHIN        = .f.
  HERZ          = .f.
  mode          = 0
  Temp          = 300 K
  tmax          = 1000 fs
  dt            = 0.01     fs  
  Ead           = 0.11348254 au
  EDMA          = 1 debye
  EDME          = 0.024734636 debye
  FreqScale     = 1.0
  logFile       = "spec.tvcf.log"
  DSFile        = "evc.dint.dat"
  FtFile        = "spec.tvcf.ft.dat"
  DDplFile      = "evc.dint.dip"
  debug         = .t.
  isgauss       = .f.
  BroadenType   = "lorentzian"
  Broadenfunc   = "frequency"
  #Broadenfunc  = "time"
  FWHM          = 100 cm-1
  GFile         = "spec.tvcf.gauss.dat"
  NScale        = 10
  Emin          = -0.3 au
  Emax          =  0.3 au
  dE            = 0.00001 au
  FoFile        = "spec.tvcf.fo.dat"
  FoSFile       = "spec.tvcf.spec.dat"
  spectra0      = .f.
  IntEmin       = 0.0 au
  IntEmax       = 0.09 au
/



&isc_tvcf
  DUSHIN        = .f.
  mode          = 0
  Temp          = 300 K
  tmax          = 1500 fs
  dt            = 0.01    fs  
  Ead           = 0.11348254 au
  isISC         = .true.
  Hso           = 1.523253518 cm-1
  FreqScale     = 1.0
  logFile       = "isc.tvcf.log"
  DSFile        = "evc.dint.dat"
  FtFile        = "isc.tvcf.ft.dat"
  debug         = .f.
  isgauss       = .f.
  BroadenType   = "lorentzian"
  Broadenfunc   = "frequency"
  #Broadenfunc  = "time"
  FWHM          = 100 cm-1
  GFile         = "isc.tvcf.gauss.dat"
  NScale        = 10
  Emin          = -0.3 au
  Emax          =  0.3 au
  dE            = 0.00001 au
  FoFile        = "isc.tvcf.fo.dat"
  FoSFile       = "isc.tvcf.spec.dat"
  spectra0      = .f.
/

evc.out文件报错如下
Warning! nv is incorrect
  Try to set set_cart = t
  or set
  check_gmatrix_stop = f
Stop!

cokie

报错已经把解决方案给你了，写set_cart = t或者check_gmatrix_stop = f。你这个是振动自由度不匹配导致的，自由的情况自然是匹配的，但也许你在构造ONIOM模型计算振动的时候选择范围不合理，所以导致报错。

Yangd

cokie 发表于 2026-5-19 20:29
报错已经把解决方案给你了，写set_cart = t或者check_gmatrix_stop = f。你这个是振动自由度不匹配导致的， ...

我把这两个命令行加入了.inp文件里面
&evc
  ffreq(1)      = "s0.log"
  ffreq(2)      = "t1.log"
  sort_mode     = 1
  proj_reorg    = .t.
  check_gmatrix_stop = f
  set_cart = t
/

momap.out文件出现了新的报错
Error terminated for spec_ft calculation
Stop
cokie老师我应该重新计算或者重新oniom吗，如果重新计算，振动的范围应该怎么选择

cokie

你可以再试试 check_gmatrix_stop = f 这个命令，是跳过G矩阵检查的，强行计算的。

如果还不行，你可以检查一下你的T1,S0两个态的文件的原子顺序是否一致，以及振动模式数量是否匹配，以及是否有虚频。

此外我不知道你ONIOM做振动分析是怎么做的，是不是把Low layer 也加进去了。我一般做的时候，会做完自由的 ONIOM 之后，再把 lowlayer 冻结再做OPT+FREQ，lowlayer 通常坐标变化很大，如果你的 lowlayer 是自由的并且参与了频率计算，那得到的结果算磷光一方面各个坐标对不上，此外就算对上了结果也大概率用不了。