计算化学公社

标题: 怎么处理用momap计算固相下共晶磷光无辐射速率报错问题 [打印本页]

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

构建的oniom如下图
(, 下载次数 Times of downloads: 1)

输入文件如下
(, 下载次数 Times of downloads: 0)

其中.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!

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)