这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:freqofexcitedstates [2016/11/30 14:20] – liu.jun | adf:freqofexcitedstates [2021/03/19 10:56] (当前版本) – [结果查看] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======如何进行激发态的振动频率计算====== | + | ======如何进行单重激发态的振动频率计算====== |
- | + | ||
- | 假定已经完成激发态的几何结构优化,激发态的几何结构优化,参考:单重激发态的几何结构优化。本例以锌酞菁为例说明。 | + | |
=====参数设置===== | =====参数设置===== | ||
- | 这一步需要非常小心,否则会变成基态的频率计算,参数设置如下: | + | 激发态频率计算的参数设置,与激发态结构优化非常相似。 |
- | {{ adf:znpc36.jpg? | + | 激发态结构优化参考: |
+ | * [[adf:geooptofsinglet]] | ||
+ | * [[adf: | ||
- | {{ adf: | + | 激发态的几何结构优化完成后,将最后收敛的结构直接更新到激发态结构优化的Input窗口,然后将Geometry Optimization改为Frequencies即可。如果没有做激发态结构优化,参数设置方面,可以借鉴一下上述激发态结构优化的参数设置,位移区别就是Main中,将Geometry Optimization改为Frequencies: |
- | 本例只计算S1态,因此实际上只有一个激发态,当然这里也可以设置为30~50之间,这样计算的激发态数目增大,精度会提高。但对内存对消耗则随该数值迅速增大。例如设为200,可能单线程消耗内存就可以达到T的量级。 | + | |
- | {{ adf:znpc38.jpg?400 |}} | + | {{ adf:2020npc400.png?650 |}} |
- | 这里我们希望计算A不可约表示的第一个激发态的频率,因此这里输入的是1 A。如果希望计算的是其它不可约表示的激发态,可以类似设置。关于不可约表示,参考:[[adf: | + | ====注意==== |
+ | 激发态频率的计算,实际上是使用数值拟合的方式去计算。也即是说,需要将每个原子逐个从平衡位置做微小偏移动,然后计算这种微小形变之后的激发能,因此对于像这样的大分子,计算激发态频率将是非常耗时的。如果不考虑对称性的话,计算的结构数大约为原子个数的6倍——也即是说需要计算这么多次激发态。 | ||
- | 保存任务,例如名为05Freq_of_Excitation。运行计算。 | + | ADF激发态的并行效率很高,因此高核数并行对于这样的计算很重要,否则几乎不可能完成这样的计算任务。 |
- | + | ||
- | ====注意==== | + | |
- | 激发态频率的计算,实际上是使用数值拟合的方式去计算。也即是说,需要将每个原子逐个从平衡位置做微小偏移动,然后计算这种微小形变之后的激发能,因此对于像这样的大分子,计算激发态频率将是非常耗时的。如果不考虑对称性的话,计算的结构数大约为原子个数的6倍——也即是说需要计算这么多次激发态。ADF激发态的并行效率很高。因此高核数并行对于这样的计算很重要,否则几乎不可能完成这样的计算任务。总体而言ADF的效率不错,本计算在4核台式机上完成,大约用了4天。如果在16核服务器,基本上1天就足够了。 | + | |
=====结果查看===== | =====结果查看===== | ||
- | 同样地,在SCM LOGO > Spectra中可以查看到振动谱: | + | 同样地,在SCM LOGO > Spectra中可以查看到振动谱。默认显示的是激发态的峰,应该选择振动谱(Spectra - Vibration)。振动谱中没有虚频(负数的频率),表示激发态结构优化没有问题。 |
- | + | ||
- | 首先显示的是激发态(我们上面只计算了1 A这一个激发态,因此只有一个峰): | + | |
- | + | ||
- | {{ adf: | + | |
- | + | ||
- | 选择振动谱: | + | |
- | + | ||
- | {{ adf: | + | |
- | 即得到右方的振动谱。振动谱中没有虚频(负数的频率),表示激发态结构优化没有问题。 | + | Spectra窗口显示谱图时,窗口标题栏,会清楚显示打开的是哪个文件,振动谱信息即存储于该文件。 |