计算化学公社

标题: 求助ESD模块计算ISC速率是否需要优化S1和T1结构的问题 [打印本页]

作者
Author: liuxiao.shaw 时间: 2024-6-8 00:40
标题: 求助ESD模块计算ISC速率是否需要优化S1和T1结构的问题
本帖最后由 liuxiao.shaw 于 2024-6-8 10:35 编辑

ORCA使用ESD模块计算ISC速率时，首先需要计算获得分子的S1态和T1态的.hess文件。
有两个问题想请教：
1. 分子的S1态和T1态需要先做td结构优化吗？
2. S1态和T1态的频率分析是在基态做还是在激发态做？

求大神指教！

作者
Author: wzkchem5 时间: 2024-6-8 16:33
1. 必须做结构优化，但不一定要用TDDFT。例如T1态可以用UKS优化
2. 和结构优化所用计算方法必须一致，例如结构优化用UKS，那么频率必须用UKS；结构优化用TDDFT，那么频率必须用TDDFT；结构优化用TDA，那么频率必须用TDA

作者
Author: liuxiao.shaw 时间: 2024-6-8 20:47

wzkchem5 发表于 2024-6-8 16:33
1. 必须做结构优化，但不一定要用TDDFT。例如T1态可以用UKS优化
2. 和结构优化所用计算方法必须一致，例如 ...

好的，谢谢老师！
还有一个问题想要请教，使用ORCA计算TDDFT算了几个小时之后出现error，显示SCF不收敛。
输入文件：
! M062X 6-311+g(d) RIJCOSX tightSCF defgrid3 noautostart miniprint nopop
%maxcore 2000
%pal nprocs 16 end
%tddft
nroots 10
dosoc true
TDA false
printlevel 3
end

分子结构和频率分析是用g16提前做好的，使用的是CAM-B3LYP/6-31+G(d)。

请问这个问题该怎么解决呢？

(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1)

作者
Author: liuxiao.shaw 时间: 2024-6-8 20:56

liuxiao.shaw 发表于 2024-6-8 20:47
好的，谢谢老师！
还有一个问题想要请教，使用ORCA计算TDDFT算了几个小时之后出现error，显示SCF不收敛 ...

另外尝试设置%scf sthresh 1e-6，运行ORCA报错如下。是因为我sthresh设置的格式不对吗？
输入：
! M062X 6-311+g(d) RIJCOSX tightSCF defgrid3 noautostart miniprint nopop
%maxcore  5000
%pal nprocs  32 end
%tddft
nroots 10
dosoc true
TDA false
printlevel 3
%scf
sthresh 1e-6
end

输出错误：
[file orca_tools/Tool-Scanner/qcscan1.cpp, line 88]:
   Unrecognized symbol in CIS block line 9  :
   Last token : %
[file orca_tools/Tool-Scanner/qcscan1.cpp, line 88]:
   Unrecognized symbol in CIS block line 9  :
   Last token : %

作者
Author: wzkchem5 时间: 2024-6-9 16:08

liuxiao.shaw 发表于 2024-6-8 13:56
另外尝试设置%scf sthresh 1e-6，运行ORCA报错如下。是因为我sthresh设置的格式不对吗？
输入：
! M062 ...

%tddft模块漏写了end

作者
Author: wzkchem5 时间: 2024-6-9 16:09

liuxiao.shaw 发表于 2024-6-8 13:47
好的，谢谢老师！
还有一个问题想要请教，使用ORCA计算TDDFT算了几个小时之后出现error，显示SCF不收敛 ...

典型的基组线性相关问题，去掉不必要的弥散函数以后重试

作者
Author: liuxiao.shaw 时间: 2024-6-9 22:22

wzkchem5 发表于 2024-6-9 16:09
典型的基组线性相关问题，去掉不必要的弥散函数以后重试

好的，谢谢老师。这个问题通过调整sthresh解决了~

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