使用Fcclasses3计算IC速率报错

qian.

使用高斯优化au4的s0,s1并进行频率计算，均无虚频；然后对两个fchk文件转化为fcc文件进行IC速率计算。计算fcc.inp文件如下

$$$
PROPERTY     =  IC  ; OPA/EMI/ECD/CPL/RR/TPA/TPCD/MCD/IC/NRSC
MODEL        =   AH   ; AS/ASF/AH/VG/VGF/VH
DIPOLE       =   FC   ; FC/HTi/HTf
TEMP         =   300.00 ; (temperature in K)
BROADFUN     =   GAU  ; GAU/LOR/VOI
HWHM         =   0.06 ; (broadening width in eV)
METHOD       =   TD   ; TI/TD
;VIBRATIONAL ANALYSIS
NORMALMODES  =   COMPUTE   ; COMPUTE/READ/IMPLICIT
COORDS       =   CARTESIAN  ; CARTESIAN/INTERNAL
;INPUT DATA FILES
STATE1_FILE  =   au4-s1-f-1.fcc
STATE2_FILE  =   au4.fcc
ELDIP_FILE   =  nac_au4-s1-f-1_fchk

运行后会报错，报错内容为
NOTE: The calculation for this State will assume non-stationary point since the gradient                  
      is not negligible (maximum atomic gradient:   0.000885 a.u.)                                       
============================================
FATAL ERROR: Non-zero gradient on State2 with Adiabatic model
请问这是由于什么原因引起的呢（之前我算系间穿越速率也出现过这个报错，但是那时候t1有虚频，我以为是虚频的原因，但是现在这个没有虚频呀）

wal

FCclasses3的NAC文件能这么输入吗 我没试过

qian.

wal 发表于 2024-11-25 15:36
FCclasses3的NAC文件能这么输入吗 我没试过

我是看帖子上有人是这么写输入文件的，所以也这么写了

wal

qian. 发表于 2024-11-25 15:57
我是看帖子上有人是这么写输入文件的，所以也这么写了

按手册上的标准写法输入试试

cokie

2L说的对，ELDIP是读跃迁偶极矩的，你要写NAC_FILE = xxxxxxxxxxxxxxxxxxxx

qian.

wal 发表于 2024-11-25 15:59
按手册上的标准写法输入试试

好的好的，这个错误是我个人原因，我没看到这个地方，生成inp模版后只顾着改后面文件名了，不好意思，麻烦了。但是，我改了之后还是会有一样的报错

qian.

cokie 发表于 2024-11-25 16:47
2L说的对，ELDIP是读跃迁偶极矩的，你要写NAC_FILE = xxxxxxxxxxxxxxxxxxxx

好的，我已经改了这个地方，但是还有上面一样的报错

wal

qian. 发表于 2024-11-26 09:27
好的好的，这个错误是我个人原因，我没看到这个地方，生成inp模版后只顾着改后面文件名了，不好意思，麻 ...

那还挺奇怪的，我用的时候没遇到过这个报错，感觉像是没优化到稳定结构

cokie

qian. 发表于 2024-11-26 09:29
好的，我已经改了这个地方，但是还有上面一样的报错

你从什么软件的输出文件获得的au4.fcc这个态文件？

qian.

cokie 发表于 2024-11-26 13:23
你从什么软件的输出文件获得的au4.fcc这个态文件？

我是通过高斯优化得到fchk文件，然后通过gen_fcc_state -i au4.fchk这个命令得到au4.fcc文件

qian.

wal 发表于 2024-11-26 13:04
那还挺奇怪的，我用的时候没遇到过这个报错，感觉像是没优化到稳定结构

可是我高斯上优化是正常结束的，并且我看了没有虚频，所以就搞不懂为啥

cokie

qian. 发表于 2024-11-26 14:41
可是我高斯上优化是正常结束的，并且我看了没有虚频，所以就搞不懂为啥

不知道你这个问题解决了没有。之前解决的报错让我想起了这个帖子。

前些日子也有个朋友有类似的问题，没有虚频但是有类似报错，他上传了出现报错对应的state的计算频率的.log文件。里面虽然没有虚频，但是计算频率的时候只有两个yes，还有两个no，说明他选择的格点精度不够导致虽然符合gaussian的收敛限，但是对于fcclasses收敛限不够，所以报错，你这个情况需要再提升格点精度重新算hessian，再用于fcclasses。