求助cp2k优化激发态结构需要设置那些参数

wmy

请问一下如果优化激发态下的结构，和普通结构优化相比应该设置哪些参数呢？

Eudaimonia

首先确定你的cp2k版本为9.1或者2022.1，低版本用不了

在&DFT部分加入

    &EXCITED_STATES
      STATE N
    &END EXCITED_STATES

N指你需要优化第几个激发态
在&FORCE_EVAL部分加入&TDDFPT部分，注意里面计算的激发态(由NSTATES控制)必须大于你优化的激发态的数，如果是周期性结构还需要在&DFT部分加入

    &AUXILIARY_DENSITY_MATRIX_METHOD
      METHOD BASIS_PROJECTION
      ADMM_PURIFICATION_METHOD NONE
      EXCH_SCALING_MODEL NONE
      EXCH_CORRECTION_FUNC PBEX
    &END AUXILIARY_DENSITY_MATRIX_METHOD

注意这里使用了AUXILIARY_DENSITY_MATRIX_METHOD，所以你的基组需要增加 BASIS_SET AUX_FIT 部分，如果使用的泛函是PBE0，B3LYP和HS06可以由Multiwfn方便的设置，不过默认的AUX_FIT基组为cFIT，对重元素支持不全，官方教程里面使用的是BASIS_ADMM_UZH系列

具体教程见
cp2k手册优化激发态关键词
Google Group里面推荐的教程

wmy

Eudaimonia 发表于 2022-9-19 21:08
首先确定你的cp2k版本为9.1或者2022.1，低版本用不了

在&DFT部分加入

好的，谢谢

wmy

wmy 发表于 2022-9-19 22:19
好的，谢谢

1 优化第几个激发态是什么意思呢？是不是就是优化第几个激发态的结构的意思呢？
2 还有就是任务类型需要改成优化结构吗？还是继续用单点？
3 有没有具体的inp文件呢？看一下应该就懂了

Eudaimonia

1.y
2.GeoOpt
3.教程里面有，自己去翻