第一步,优化基态的分子结构(详情请点击);如果计算垂直激发,则不需要优化激发态,否则还需要优化激发态(不支持基态方法优化T1的情况),参考:
任务类型必须选择为Geometry Optimization,而且后面还必须设置优化0步:
如果计算单重态可以选择Allowed Only或Singlet Only,否则选择Triplet Only:
这里设置需要计算哪个激发态的CDD与电子密度,n A表示第n个激发态,通过下方选择Singlet/Triplet来决定是计算Sn还是Tn。这里设置与上图的设置一致。
注意:LDA、GGA低估了交换作用,因此会产生一种不符合事实的后果 —— 对于能量很低的激发态,也出现了电荷的转移。一个比较粗暴的解决方式就是替换为杂化泛函,例如B3LYP。
保存并运行任务。计算完成的时候,会显示:
ERROR: GEOMETRY DID NOT CONVERGE
不必理会。
SCM > View > Add > Isosurface:With Phase > Select Field > Excited State > Difference Density
将右下角的数值调整(该值为等值面的数值,不同的数值,则窗口显示该值的等值面),例如此例改为0.02。这就是该几何结构下,该激发态(激发态几何结构优化那一步指定的那个激发态)与基态的电子密度之差的空间分布。其中红色区域为负值、蓝色为正值。
基态电子密度+CDD=激发态电子密度。因此执行如下操作:Field — Calculate,之后最下方下方选择如图所示
左下角的编号“C-”1,这是计算出来的“CDD+基态密度=激发态密度”的数据编号。
显示该密度(C-1这一项):
Add - ISOSurface > Select Field > Others > C-1:
即显示激发态电子密度的等值面。
设置图的透明度,参考:如何设置View中各种空间分布图的透明度
AMS软件提供免费试用,试用申请方式参考:AMS免费试用**