“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 89|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] Pt配合物ORCA与FCclasses3的ISC速率计算问题

[复制链接 Copy URL]

40

帖子

0

威望

235

eV
积分
275

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师,求助Pt配合物ISC速率计算问题,应该如何解决
体系90原子S1,T1结构无虚频,SOC大概30cm-1,ORCA计算速率只有10的2次方数量级,严重偏小,同体系另两个分子计算得到的速率为8-10次方数量级
使用ORCA计算,尝试过增加NPOINTS,LINEW试过10,20,25,100,默认条件计算是负值,随着改变增加到2次方便不再有变化
inp文件
! ESD(ISC) noiter
%ESD
  ISCISHessian "N1-s1.hess"
  ISCFSHessian "N1-t1.hess"
  UseJ true
  Dele 1685.46
  SOCME 0.00013139
  Temp 77
  LINEW 25
  NPOINTS 1048576
end
* xyzfile 0 1 N1-t1.xyz

out文件
The calculated ISC rate constant is                1.144603e+02 s-1



换用FCclasses3计算,将ORCA得到的hess文件转换成fcc文件,运行fcclasses3 fcc.inp后出现报错

$$$
PROPERTY     =   NR0  ; OPA/EMI/ECD/CPL/RR/TPA/MCD/IC
MODEL        =   AH   ; AS/ASF/AH/VG/VGF/VH
DE           =   0.2090
NR0_COUPL    =   1.3139E-4 ; SOC
TEMP         =   298.15 ; (temperature in K)
BROADFUN     =   GAU  ; GAU/LOR/VOI
HWHM         =   0.06 ; (broadening width in eV)
METHOD       =   TD   ; TI/TD
ROT          =   1    ; Rotate state 1
;VIBRATIONAL ANALYSIS
NORMALMODES  =   COMPUTE   ; COMPUTE/READ/IMPLICIT
COORDS       =   CARTESIAN ; CARTESIAN/INTERNAL
;INPUT DATA FILES
STATE1_FILE  =   N1-s1.fcc
STATE2_FILE  =   N1-t1.fcc


报错信息
ERROR in diagonalization. (LAPACK error code: 269)

Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-4-13 19:01 , Processed in 0.152588 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list