这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:nto [2017/03/28 16:36] – liu.jun | adf:nto [2020/12/01 22:54] (当前版本) – [图谱:UV与NTO] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======如何计算紫外-可见吸收谱与自然跃迁轨道NTO====== | + | ======紫外-可见吸收谱与自然跃迁轨道NTO、跃迁偶极矩(非相对论)====== |
+ | =====前言===== | ||
虽然电子的激发实际上比较复杂,电子和电子之间是相关的,而不是彼此孤立的。因此真正严格的说,激发并不是某个电子从某个占据轨道跃迁到另外一个空轨道,而是互相牵连,从一个“多电子”波函数变为另一个“多电子”波函数。 | 虽然电子的激发实际上比较复杂,电子和电子之间是相关的,而不是彼此孤立的。因此真正严格的说,激发并不是某个电子从某个占据轨道跃迁到另外一个空轨道,而是互相牵连,从一个“多电子”波函数变为另一个“多电子”波函数。 | ||
行 6: | 行 6: | ||
NTO就是这样一个直观的图像,列出这些“占据轨道-空轨道”对。 | NTO就是这样一个直观的图像,列出这些“占据轨道-空轨道”对。 | ||
- | + | =====参数设置===== | |
- | 在ADF里面计算NTO非常简单。在ADF2016里面,只要是激发态计算,默认就要计算NTO,所以不需要额外的设置: | + | |
- | + | ||
- | {{: | + | |
- | + | ||
- | 在旧版本里面,只需要在生成的*.run文件中,EXCITATIONS字段中增加一行NTO即可。 | + | |
- | EXCITATIONS | + | |
- | ... | + | |
- | NTO | + | |
- | End | + | |
- | + | ||
- | 关于分子的建模,参考: | + | |
- | + | ||
- | * [[https:// | + | |
- | * [[adf: | + | |
最简单的,计算NTO的过程,图示如下: | 最简单的,计算NTO的过程,图示如下: | ||
- | {{: | + | {{ : |
+ | |||
+ | {{ : | ||
- | {{: | + | 注意,激发态数目实际上是吸收峰的数目,但是计算的时候,是从最长波区域往短波区域计算。因此,激发态数目影响到短波区域的峰数量、峰形。数目本身不太影响计算效率,但是对内存需求剧烈增加,一般较大的分子无法计算超过100个激发态,否则对内存的需求将达到T的量级。 |
- | {{: | + | 另外,这里自动勾选了NTO。 |
- | {{:adf:nto05.png|}} | + | 即完成提交任务。其它提交任务的方式,参考:[[adf:maintance|]] |
- | {{: | + | =====查看结果===== |
+ | ====激发能==== | ||
+ | SCM > Output > Response properties > All Singlet-Singlet Excitation Energies:列出的是,考虑相对论对动能修正的S0→Sn激发能。 | ||
+ | < | ||
+ | All SINGLET-SINGLET excitation energies | ||
- | 提交任务。其它提交任务的方式,参考费米科技WIKI:[[adf: | + | |
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | 10: | ||
+ | 11: | ||
+ | 12: | ||
+ | 13: | ||
+ | 14: | ||
+ | 15: | ||
+ | 16: | ||
+ | 17: | ||
+ | 18: | ||
+ | 19: | ||
+ | 20: | ||
+ | 21: | ||
+ | 22: | ||
+ | 23: | ||
+ | 24: | ||
+ | 25: | ||
+ | 26: | ||
+ | 27: | ||
+ | 28: | ||
+ | 29: | ||
+ | 30: | ||
+ | 31: | ||
+ | 32: | ||
+ | 33: | ||
+ | 34: | ||
+ | 35: | ||
+ | 36: | ||
+ | 37: | ||
+ | 38: | ||
+ | 39: | ||
+ | 40: | ||
+ | </ | ||
+ | 其中tau为激发态的辐射跃迁寿命,其倒数为辐射跃迁速率 | ||
+ | ====跃迁偶极矩==== | ||
+ | 搜索“Transition dipole moments mu”即可看到跃迁偶极矩: | ||
+ | < | ||
+ | | ||
+ | (weak excitations are not printed) | ||
- | 查看结果: | + | |
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | 10 12.588 | ||
+ | 11 12.880 | ||
+ | 12 13.564 | ||
+ | 13 14.471 | ||
+ | 14 14.496 | ||
+ | 15 15.044 | ||
+ | 16 16.030 | ||
+ | 17 16.352 | ||
+ | 18 16.531 | ||
+ | 19 18.016 | ||
+ | 20 18.242 | ||
+ | 21 18.422 | ||
+ | 22 18.711 | ||
+ | 23 19.257 | ||
+ | 24 20.159 | ||
+ | 25 20.302 | ||
+ | 26 20.403 | ||
+ | 27 20.583 | ||
+ | 28 20.803 | ||
+ | 29 21.481 | ||
+ | 30 22.005 | ||
+ | 31 22.020 | ||
+ | 32 22.247 | ||
+ | 33 22.766 | ||
+ | 34 23.029 | ||
+ | 35 23.568 | ||
+ | 36 23.855 | ||
+ | 37 24.311 | ||
+ | 38 24.530 | ||
+ | 39 24.714 | ||
+ | 40 25.156 | ||
+ | </ | ||
+ | ====图谱:UV与NTO==== | ||
- | {{: | + | SCM - Spectra |
- | {{:adf:nto08.png|}} | + | {{ :adf:nto03.png?650 }} |
- | 上图中,下方的窗口按照激发能的从低到高,依次列出每个激发态的所有有贡献的“占据轨道-空轨道”对。 | + | 列表中每一行,对应吸收峰的一个峰,点击将显示该吸收峰的来源,例如上图,表示该吸收峰是8a轨道跃迁到9a轨道。可以通过SCM - level查看能级图对照得到,实际上就是HOMO跃迁到LUMO。 |
- | 例如第一个激发态,激发能是0.310Hartree,是由两对“占据轨道-空轨道”贡献出来的。上图中点击第一对(贡献了99.54%),在上方的窗口就显示出这一对“占据轨道-空轨道”。其中occ表示占据轨道,virt表示空轨道。最左边√则显示该轨道,不√则不显示。上图中,两个轨道都被√了。最后边的数字0.03是等值面的数值,如果知道等值面的意义,则可以自行调整,如果不知道,最好不要动。分别显示occ和virt,可以看到它们处在那个区域上,比如对金属配合物而言,是处于金属上还是处于配体上。 | + | **右下角列表中,一般也会列出NTO信息,蓝色的文字可以点开,即显示NTO轨道。** |
- | 第二个激发态,主要由四对“占据轨道-空轨道”贡献出来的,激发能为0.641Hatree,贡献最大的一对占97.01%。 | + | * 计算有机物的紫外可见吸收谱,往往使用B3LYP能得到很好的结果,但该泛函不适用于多金属中心体系 |
+ | * 选择菜单栏Axes - Molar Adsorption Coefficient,将显示摩尔吸收系数 | ||
+ | * 横坐标单位为Hartree,点击菜单栏Axes - Horizontal Unit - nm可以修改为nm,但是注意横坐标不要出现负值,否则转换的时候会报错 | ||
+ | * 吸收峰的强度只要不为0,往往在实验中就能观察到 |