请教激发态频率计算

石头

一个纯有机分子共90个原子，S1态优化基于BMK/6-31G(d)基组算的，激发态频率一开始用的是以下方法

%chk=1atd5soptfreq1
%mem=1GB
%nprocshared=12
#p freq td=(singlets,nstates=10,direct) bmk/6-31g(d) scrf=(solvent=toluene) nosymm geom=connectivity

可是算了好久都没算完，每次都是设定的时间用完自己下来了，后来尝试减少NSTATES=3，可是算了好久还是算不完，自己又中断下来。
%mem=1GB
%nprocshared=12
%NProcLinda=1
#p freq td=(singlets,nstates=3,direct) bmk/6-31g(d) scrf=(solvent=toluene) nosymm geom=connectivity

请教各位大神有没有什么好的办法得到这个激发态频率？难道就这样一直算下去？不知道何时能出来结果

小范范1989

为啥要计算激发态的频率呢？
我用高斯计算过，非常慢，当年上研二的时候，半年计算了四个分子。记忆尤深。
不知道高斯下一个新版本能不能改进这一块。
我记得他们说起过别的软件能计算激发态频率。

石头

小范范1989 发表于 2017-1-3 10:50
为啥要计算激发态的频率呢？
我用高斯计算过，非常慢，当年上研二的时候，半年计算了四个分子。记忆尤深。 ...

我们优化完S1态几何构型，不是要计算一下这个构型的频率来确认一下是否存在虚频吗？这个问题好头痛，算了好久了。

yezhonghua

我用g9算过90个原子的激发态频率，24核心的2680V3，32G的内存，算了90个小时才出结果。
但是如果用Turbomore算的化20个小时就可以了。

遗忘爱神

挑几个关键的点看下有无虚频即可，实在不行就不算了

sobereva

g09没有TDDFT解析Hessian，只有解析一阶导数，所以N个原子计算频率的耗时等于6N步几何优化耗时，显然一算就知道耗时得有多高（对比激发态几何优化每一步的耗时）
要么等Gaussian17出来，有解析TDDFT Hessian，要么用Q-Chem，是目前唯一支持TDDFT解析Hessian的程序。
换成其它的程序都没法显著降低耗时，这从原理上一想就明白。

你凭直觉能确信优化到了极小点就没必要再算频率。或者用CIS做opt freq，对于得到的没虚频的结构，如果结构和TDDFT优化的相仿佛，一般也能说明TDDFT优化的没虚频。

石头

yezhonghua 发表于 2017-1-3 12:23
我用g9算过90个原子的激发态频率，24核心的2680V3，32G的内存，算了90个小时才出结果。
但是如果用Turbomo ...

感谢您的回复，我这个分子加上中间间断了几次折腾了有差不多一个月了，还是不行，不过你用的CPU是我的2倍。我再想办法试一试，谢谢！

石头

遗忘爱神 发表于 2017-1-3 14:06
挑几个关键的点看下有无虚频即可，实在不行就不算了

感谢您的回复，请问您说的挑几个关键的点，具体如何做呢？谢谢指点！

石头

sobereva 发表于 2017-1-3 14:20
g09没有TDDFT解析Hessian，只有解析一阶导数，所以N个原子计算频率的耗时等于6N步几何优化耗时，显然一算就 ...

谢谢SOB老师指点，我按照您说的方法再试一下，非常感谢！

yezhonghua

石头 发表于 2017-1-3 14:51
感谢您的回复，我这个分子加上中间间断了几次折腾了有差不多一个月了，还是不行，不过你用的CPU是我的2倍 ...

你实在不行，把文件给我看一下，我寒假设备有空闲的化给你挂上。

石头

yezhonghua 发表于 2017-1-3 15:19
你实在不行，把文件给我看一下，我寒假设备有空闲的化给你挂上。

真的非常感谢您！我先试试SOB老师说的CIS方法，我的分子有好多个，实在是不好意思占用您那么多资源。如果还是算不出来再请您帮助我。再次感谢！！！

遗忘爱神

石头 发表于 2017-1-3 14:51
感谢您的回复，请问您说的挑几个关键的点，具体如何做呢？谢谢指点！

就是你这个反应里比较重要的点，确保他们没有虚频。

石头

遗忘爱神 发表于 2017-1-3 18:43
就是你这个反应里比较重要的点，确保他们没有虚频。

好的，非常感谢您，我再好好想一想。