用户工具

站点工具


adf:geooptoftri

差别

这里会显示出您选择的修订版和当前版本之间的差别。

到此差别页面的链接

两侧同时换到之前的修订记录前一修订版
后一修订版
前一修订版
adf:geooptoftri [2018/03/04 09:15] – [优化Tn] liu.junadf:geooptoftri [2020/12/01 19:23] (当前版本) – 移除 liu.jun
行 1: 行 1:
-======如何进行三重高激发态(Tn,N>1)的几何结构优化====== 
  
-有两种方式:1)直接优化Tn态(适用于基态为S0的体系);2)先优化T1态,然后优化Tn态(适用于基态为T1的体系,例如<chem>O2</chem> 
- 
-=====一、直接优化Tn态(适用于基态为S0的体系)===== 
- 
-以<chem>CH2O</chem>为例说明: 
- 
-====参数设置==== 
-注意此时的计算中,基态是单重态,激发态才是三重态,我们是希望通过计算激发态来优化<color green>(如果要考虑相对论,但不考虑旋轨耦合,这里可以设置Relativity (ZORA) : Scalar)</color> 
- 
-{{ :adf:t1-opt04.png?600 |}} 
- 
-选择只计算三重激发态: 
- 
-{{ :adf:t1-opt05.png?600 |}} 
- 
-选择要优化哪个激发态,n A表示优化Tn态,例如下图表示优化T3态: 
- 
-{{ :adf:t1-opt10.png?600 |}} 
- 
-因为我们预先并不知道<chem>CH2O</chem>的T3态有没有对称性,因此和前一种方法一样,取消掉对称性,让程序在没有对称性限制的条件下去优化: 
- 
-{{ :adf:t1-opt09.png?600 |}} 
- 
-保存任务,并提交。 
- 
-====查看优化结果==== 
- 
-点击ADF LOGO > Movie 可以查看优化完成之后的结构。 
- 
- 
-=====二、先优化T1态,然后优化Tn态(适用于基态为T1的体系)===== 
-以<chem>O2</chem>为例进行说明,<chem>O2</chem>的基态为三重态。 
- 
-====优化T1==== 
-将体系设置为三重态,进行普通的结构优化,这实际上就是在做T1的优化,因为能量最低的三重态就是T1,参考:[[adf:t1opt]]中的第一种方法。 
- 
-==== 优化Tn==== 
- 
-基于优化之后的分子结构,计算更高的三重激发态(三重态的高激发态T2、T3、T4……)。 
- 
-===参数设置=== 
- 
-{{ :adf:tnopt02.png |}} 
- 
-说明:<color green>(如果要考虑相对论,但不考虑旋轨耦合,这里可以设置Relativity (ZORA) : Scalar)</color> 
- 
-{{ :adf:tnopt03.png |}} 
- 
-{{ :adf:tnopt04.png |}} 
- 
-=== 结果查看=== 
- 
-计算完成之后,检查该T2激发态是否合理: 
- 
-查看*.out文件中电子激发的信息 
- 
-{{ :adf:tnopt05.png |}} 
- 
-对应的轨道,可以在SCM LOGO > level中分别查看 
- 
-{{ :adf:tnopt06.png |}} 
- 
-优化所需的三重态高激发态(T2,T3,T4……),例如T2(三重态的第二个态)。 
- 
- 
-{{ :adf:tnopt07.png |}} 
- 
-{{ :adf:tnopt08.png |}} 
- 
-{{ :adf:tnopt09.png |}} 
- 
-**注意:如果优化T3,T4……,则将此处设置为2A、3A……即可。此处spin type不需要修改,这是目前版本的一个bug,保留原样即可,不必改为Triplet。** 
- 
-{{ :adf:tnopt04.png |}} 
- 
-计算结束,收敛正常(得到5个T,尤其是第二个T最重要): 
- 
-{{ :adf:tnopt10.png |}} 
- 
-查看激发能的变化情况: 
- 
-{{ :adf:tnopt11.png |}} 
- 
-最后一个能量,就是优化完成的激发能。 
- 
-优化经历了6步。最后用第6步的结构,重新计算了激发能。 
- 
-SCM LOGO > Spectra,显示的就是最后这个几何结构对应的激发能。要注意,T2的总能量应该是该结构的“T1态能量”加上“激发能”,也就是: 
- 
-“logfile最后显示的Bonding Energy”加上“激发能”。 
- 
-本例中logfile最后显示的Bonding Energy: 
- 
-{{ :adf:tnopt12.png |}} 
- 
-因此T2的总能量就是:-7.87621845 eV+3.15671 eV = -4.71950845 eV。在新版中logfile会列出 current energy 实际上就是该激发态的能量,也就是当前结构的基态能量+激发能。 
- 
-至此完成。如果要考虑旋轨耦合,请参考: 
- 
-费米维基:[[adf: socmatrix]] 
- 
-费米维基:[[adf:nonrels-t]] 
- 
-费米维基:[[s-trelative]] 
- 
-费米维基:[[adf:splitoftriplet|三重态的劈裂]] 
adf/geooptoftri.1520126119.txt.gz · 最后更改: 2018/03/04 09:15 由 liu.jun

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