计算化学公社

标题: 使用Fcclasses3计算IC速率报错 [打印本页]

作者
Author: qian. 时间: 2024-11-25 15:27
标题: 使用Fcclasses3计算IC速率报错
使用高斯优化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有虚频，我以为是虚频的原因，但是现在这个没有虚频呀）

作者
Author: wal 时间: 2024-11-25 15:36
FCclasses3的NAC文件能这么输入吗我没试过

作者
Author: qian. 时间: 2024-11-25 15:57

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

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

作者
Author: wal 时间: 2024-11-25 15:59

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

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

作者
Author: cokie 时间: 2024-11-25 16:47
2L说的对，ELDIP是读跃迁偶极矩的，你要写NAC_FILE = xxxxxxxxxxxxxxxxxxxx

作者
Author: qian. 时间: 2024-11-26 09:27

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

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

作者
Author: qian. 时间: 2024-11-26 09:29

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

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

作者
Author: wal 时间: 2024-11-26 13:04

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

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

作者
Author: cokie 时间: 2024-11-26 13:23

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

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

作者
Author: qian. 时间: 2024-11-26 14:38

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

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

作者
Author: qian. 时间: 2024-11-26 14:41

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

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

作者
Author: cokie 时间: 2024-12-17 13:03

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

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

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

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