用户工具

站点工具

本页面的其他翻译:
  • zh

adf:emissionofpho2023

如何计算磷光发射谱(计算、实验对比)

前言

本文使用AMS2020.101完成计算。本例介绍从三重激发态(T1)回到基态(S0)的发射谱,该发射谱考虑到了从T1回到S0的各个振动态的分布情况。

参考文献:

A.F. Rausch et al. / Chemical Physics Letters 468 (2009) 46–51

该文献中给出了实验测量的20K与300K时的磷光发射谱:

本文演示如何通过计算,得到该谱图。分子结构下载

一,计算S0态的频率

例子文件名:S0_OptFreq

对基态进行结构优化,并计算频率,参考费米科技维基:如何优化分子的几何结构如何计算分子的频率、红外、零点能、转动能量、转动惯量、熵、焓、热熔与Gibbs自由能。本文中,使用的参数:

  • 泛函:BP86
  • 相对论:Scalar
  • 基组:DZP(Pt会自动选为TZP)
  • 冻芯:small
  • 精度:Good

例子中,Task选择Geometry Optimization,同时勾选Frequencies一次性完成结构优化和频率计算。计算完毕,得到S0_OptFreq.results/adf.rkf文件备用。

二,计算T1态的频率

例子文件名: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文件备用。

三,计算Franck-Condon Factors

参数设置: Details → File(Restart)读取末态(即S0态)的*.results文件夹 打开FCF开关,并选择参考态(即初态T1),类别选择Emission:

保存作业并提交运行。

四,得到发射谱

计算完毕,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

adf/emissionofpho2023.txt · 最后更改: 2023/04/27 12:56 由 liu.jun

© 2014-2022 费米科技(京ICP备14023855号