本文使用AMS2020.101完成计算。本例介绍从三重激发态(T1)回到基态(S0)的发射谱,该发射谱考虑到了从T1回到S0的各个振动态的分布情况。
参考文献:
A.F. Rausch et al. / Chemical Physics Letters 468 (2009) 46–51
该文献中给出了实验测量的20K与300K时的磷光发射谱:
本文演示如何通过计算,得到该谱图。分子结构下载
例子文件名:S0_OptFreq
对基态进行结构优化,并计算频率,参考费米科技维基:如何优化分子的几何结构、如何计算分子的频率、红外、零点能、转动能量、转动惯量、熵、焓、热熔与Gibbs自由能。本文中,使用的参数:
例子中,Task选择Geometry Optimization,同时勾选Frequencies一次性完成结构优化和频率计算。计算完毕,得到S0_OptFreq.results/adf.rkf文件备用。
例子文件名:T1_OptFreq
T1态结构优化后得到的Bond Energy(-219.63454835 eV)与S0的Bond Energy(-222.12126667 eV)相比,大约高2.5eV,也就是大约20000cm-1。这与其他地方得到的实验得到的“T1振动基态”→“S0振动基态”的跃迁符合的比较好。
T1态频率计算,参数设置与S0态的计算一样,唯一区别是Spin polarization为2,但是注意采用T1态优化后的结构。
计算完毕,得到T1_OptFreq.results/adf.rkf文件备用。
用写字板或记事本,创建一个txt文件,内容如下:
$AMSBIN/fcf << eor STATES D:\ADF_DATA\FCF_Pho\fcf\S0_OptFreq.results\adf.rkf D:\ADF_DATA\FCF_Pho\fcf\T1_OptFreq.results\adf.rkf QUANTA 5 0 TRANSLATE ROTATE SPECTRUM -10000 0 1000 eor
创建完毕后,改名为FCF.run,上传到Linux计算。或者直接在Linux创建FCF.run,内容如上。命令行进入FCF所在目录,并赋予可执行权限:chmod 700 FCF.run
注意:
计算完毕,AMSjobs窗口选中FCF任务,点击SCM > Spectra 即可显示振动对发射的一个展宽谱(下图中的峰展宽修改为200了):
可以看到峰的范围是-10000~0之间。
我们前面说了,0 cm-1对应的是“T1振动基态”→“S0振动基态”的跃迁。因此将Spectrum 的谱的0 cm-1平移进行平移,也就是横坐标统一加一个常数:第一步得到的T1态、第二步得到的S0态之间的能量差。即得到与实验图谱相似的结果。我们上面计算得到的(-219.450286+221.95839569) eV = 2.5081 * 8065.5 cm-1 = 20229 cm-1(实验值为21461 cm-1)。
因此,在上图底部的offset填入20229,并将坐标系的横坐标颠倒一下顺序(菜单栏Axes - Flip Herizontal),让它变成从右向左:
这样就可以直接与文章中的实验测量的图(下图中右上角的图)直接比较了:
发射谱的整体误差,是ET1-ES0的误差造成的,而FCF计算看起来是相当准确的,因为谱的形状与实验符合的相当好。理论上讲,温度越低,谱的峰宽越小,在窗口底部的Width处可以调节。
计算文件下载:fcf.rar