有机分子吸附在金表面，CP2K用杂化泛函PBE0结合ADMM计算PDOS的问题

lycheeho

大家好，目前体系为有机分子吸附在金表面，有250个原子左右，想对比PBE与PBE0计算的PDOS差异。考虑大体系PBE0计算需要ADMM加速，OT要打开，请问这种情况该如何绘制PBE0产生molden文件的PDOS呢？谢谢大家！

UW_0728.

可以用目前GitHub上最新的开发版本，已经能够在OT情况下正确往molden文件写入空轨道的波函数（轨道系数）信息了。在&MO_MOLDEN模块下加一个子关键词OT_NLUMO n（你要的空轨道数目）就可以了。

northern_mar

UW_0728. 发表于 2026-3-17 19:56
可以用目前GitHub上最新的开发版本，已经能够在OT情况下正确往molden文件写入空轨道的波函数（轨道系数）信 ...

您好！请问您说的开发板是指2026.1版嘛？还是用于测试的版本呢？

UW_0728.

northern_mar 发表于 2026-3-17 21:57
您好！请问您说的开发板是指2026.1版嘛？还是用于测试的版本呢？

https://github.com/cp2k/cp2k/tree/master 这个开发分支。2026.1是release（已发布版本）

或者直接 git clone https://github.com/cp2k/cp2k.git cp2k 也可以

Chao666

UW_0728. 发表于 2026-3-17 19:56
可以用目前GitHub上最新的开发版本，已经能够在OT情况下正确往molden文件写入空轨道的波函数（轨道系数）信 ...

OT的短板也一个个补上来了，以后会越来越方便了

UW_0728.

Chao666 发表于 2026-3-17 22:21
OT的短板也一个个补上来了，以后会越来越方便了

其实OT情况下的空轨道信息一直就能产生，在SCF收敛后执行一次eigensolver额外计算一下就可以了；CP2K培训班里介绍的MO_INDEX_RANGE和MO_CUBES之所以能给出空轨道能量就是这样算出来的。上面所说的其实就是稍微修改了一下Molden输出时的行为，使得空轨道信息能够被输出到molden里面。

其实OT最大的一个硬伤还是无法考虑k点；但这个真的没办法，太难了，Hutter教授在相关issue的一个回复中说他们试了几次都失败了。

northern_mar

UW_0728. 发表于 2026-3-17 21:59
https://github.com/cp2k/cp2k/tree/master 这个开发分支。2026.1是release（已发布版本）

或者直接 g ...

谢谢您

lycheeho

UW_0728. 发表于 2026-3-17 19:56
可以用目前GitHub上最新的开发版本，已经能够在OT情况下正确往molden文件写入空轨道的波函数（轨道系数）信 ...

非常感谢！这样就能杂化泛函开OT算完，一并输出含有空轨道信息的molden文件了

lycheeho

UW_0728. 发表于 2026-3-17 22:59
其实OT情况下的空轨道信息一直就能产生，在SCF收敛后执行一次eigensolver额外计算一下就可以了；CP2K培训 ...

你好，还有一个小问题，这样开OT的情况下也是无法开ADDED_MOS求解空轨道，包括smear，只不过是把算出来的空轨道信息输出了出来？

Stardust0831

UW_0728. 发表于 2026-3-17 19:56
可以用目前GitHub上最新的开发版本，已经能够在OT情况下正确往molden文件写入空轨道的波函数（轨道系数）信 ...

我有点好奇，金属Au作为导体，是否适用于OT，计算结果是否真的是可靠的？

UW_0728.

lycheeho 发表于 2026-3-18 15:31
你好，还有一个小问题，这样开OT的情况下也是无法开ADDED_MOS求解空轨道，包括smear，只不过是把算出来的 ...

都不能，因为OT本身就设定所求解的所有轨道都是整数占据的。eigensolver实质上是对OT已收敛的结果做了一次对角化

UW_0728.

Stardust0831 发表于 2026-3-18 16:06
我有点好奇，金属Au作为导体，是否适用于OT，计算结果是否真的是可靠的？

不适合，金属用OT通常也很难SCF收敛

lycheeho

UW_0728. 发表于 2026-3-18 17:41
都不能，因为OT本身就设定所求解的所有轨道都是整数占据的。eigensolver实质上是对OT已收敛的结果做了一 ...

明白了，感谢！