怎么处理用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模型计算振动的时候选择范围不合理，所以导致报错。