这是本文档旧的修订版!
假定已经完成激发态的几何结构优化,激发态的几何结构优化,参考:单重激发态的几何结构优化。本例以锌酞菁为例说明。
这一步需要非常小心,否则会变成基态的频率计算,参数设置如下:
本例只计算S1态,因此实际上只有一个激发态,当然这里也可以设置为30~50之间,这样计算的激发态数目增大,精度会提高。但对内存对消耗则随该数值迅速增大。例如设为200,可能单线程消耗内存就可以达到T的量级。
这里我们希望计算A不可约表示的第一个激发态的频率,因此这里输入的是1 A。如果希望计算的是其它不可约表示的激发态,可以类似设置。关于不可约表示,参考:分子对称性、点群与不可约表示、轨道对称性的通俗理解
保存任务,例如名为05Freq_of_Excitation。运行计算。
激发态频率的计算,实际上是使用数值拟合的方式去计算。也即是说,需要将每个原子逐个从平衡位置做微小偏移动,然后计算这种微小形变之后的激发能,因此对于像这样的大分子,计算激发态频率将是非常耗时的。如果不考虑对称性的话,计算的结构数大约为原子个数的6倍——也即是说需要计算这么多次激发态。ADF激发态的并行效率很高。因此高核数并行对于这样的计算很重要,否则几乎不可能完成这样的计算任务。总体而言ADF的效率不错,本计算在4核台式机上完成,大约用了4天。如果在16核服务器,基本上1天就足够了。