这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
adf:cddofexcitedstate [2017/12/15 11:44] – liu.jun | adf:cddofexcitedstate [2019/12/08 17:33] – [如何计算激发态CDD(激发态与基态电子密度之差)、激发态电子密度] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ====== | + | ====== 激发态CDD(激发态与基态电子密度之差)、激发态电子密度====== |
+ | =====第一步:优化结构===== | ||
- | **本例适用于对激发态发生明显电荷转移的情况。** | + | 第一步,[[adf: |
+ | * [[adf: | ||
+ | * [[adf: | ||
+ | * [[adf: | ||
+ | | ||
+ | | ||
- | 第一步,[[adf: | + | =====第二步:计算CDD、激发态电子密度===== |
+ | 任务类型必须选择为Geometry Optimization,而且后面还必须设置优化0步: | ||
- | 第二步,优化激发态几何结构(如果关心的是垂直激发的CDD,那么设置Details > Geometry Convergence > Maximum number of iterations: 0即可),参数设置如下(以CO分子为例): | + | {{ :adf: |
- | {{ : | + | 如果计算单重态可以选择Allowed Only或Singlet Only,否则选择Triplet Only: |
- | {{ : | + | {{ : |
- | {{ : | + | 这里设置需要计算哪个激发态的CDD与电子密度,n A表示第n个激发态,通过下方选择Singlet/ |
- | {{ : | + | {{ : |
- | 注意:LDA、GGA低估了交换作用,因此会产生一种不符合事实的后果 —— 对于能量很低的激发态,也出现了电荷的转移。一个比较粗暴的解决方式就是替换为杂化泛函,例如B3LYP。激发态几何优化取消对称性的原因是:我们并不清楚激发态分子结构对称性是不是会降低(当然对于这个例子里面,CO无论怎么变,都是C< | + | {{ : |
- | 如果是计算垂直激发的CDD,设置同上,唯独增加一项: | + | {{ : |
- | Details | + | 注意:LDA、GGA低估了交换作用,因此会产生一种不符合事实的后果 |
+ | |||
+ | 保存并运行任务。计算完成的时候,会显示: | ||
+ | < | ||
ERROR: GEOMETRY DID NOT CONVERGE | ERROR: GEOMETRY DID NOT CONVERGE | ||
- | | + | </ |
不必理会。 | 不必理会。 | ||
+ | =====第三步,结果分析===== | ||
+ | ====CDD==== | ||
+ | SCM > View > Add > Isosurface: | ||
- | 保存并运行任务。 | + | {{ : |
- | + | ||
- | 第三步,结果分析: | + | |
- | + | ||
- | 在ADFinput窗口,SCM LOGO > View: | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | 将右下角的数值调整(该值为等值面的数值,不同的数值,则窗口显示该值的等值面),例如此例改为0.003 | + | |
- | + | ||
- | {{ : | + | |
- | 这就是该几何结构下,该激发态(激发态几何结构优化那一步指定的那个激发态)与基态的电子密度之差的空间分布。 | + | 将右下角的数值调整(该值为等值面的数值,不同的数值,则窗口显示该值的等值面),例如此例改为0.02。这就是该几何结构下,该激发态(激发态几何结构优化那一步指定的那个激发态)与基态的电子密度之差的空间分布。其中红色区域为负值、蓝色为正值。 |
- | **补充内容:如何要得到总的激发态电子密度?** | + | =====激发态电子密度===== |
- | 1,Field — Calculate | + | 基态电子密度+CDD=激发态电子密度。因此执行如下操作:Field — Calculate,之后最下方下方选择如图所示 |
- | 之后最下方下方选择如图所示 | + | {{ : |
- | {{ : | + | 左下角的编号“C-”1,这是计算出来的“CDD+基态密度=激发态密度”的数据编号。 |
- | 注意右下角红色圆圈里面的C-1,这是计算出来的“基态密度+激发态基态密度差=激发态密度”的数据编号。 | + | 显示该密度(**C-1**这一项): |
- | 2,显示该密度(**C-1**这一项): | + | Add - ISOSurface > Select Field > Others > C-1: |
- | Add- ISOSurface Double(+/ | + | {{ : |
- | {{ : | + | 即显示激发态电子密度的等值面。 |
- | 即显示激发态电子密度的等值面。注意倒数第二行编号为**C-1**,因此,最后一行选择Other — 1 — C-1,即显示激发态的电子密度空间分布。 | + | 设置图的透明度,参考:[[adf: |
- | ADF软件提供**免费试用**(一般为一个月),试用申请方式参见**费米科技维基百科:[[adf: | + | AMS软件提供**免费试用**(一般为一个月),试用申请方式参考:[[adf: |