计算化学公社

标题: 有机分子吸附在金表面，CP2K用杂化泛函PBE0结合ADMM计算PDOS的问题 [打印本页]

作者
Author: lycheeho 时间: 2026-3-17 19:51
标题: 有机分子吸附在金表面，CP2K用杂化泛函PBE0结合ADMM计算PDOS的问题
大家好，目前体系为有机分子吸附在金表面，有250个原子左右，想对比PBE与PBE0计算的PDOS差异。考虑大体系PBE0计算需要ADMM加速，OT要打开，请问这种情况该如何绘制PBE0产生molden文件的PDOS呢？谢谢大家！

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

作者
Author: northern_mar 时间: 2026-3-17 21:57

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

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

作者
Author: UW_0728. 时间: 2026-3-17 21:59

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 也可以

作者
Author: Chao666 时间: 2026-3-17 22:21

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

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

作者
Author: UW_0728. 时间: 2026-3-17 22:59

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

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

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

作者
Author: northern_mar 时间: 2026-3-18 09:38

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

或者直接 g ...

谢谢您

作者
Author: lycheeho 时间: 2026-3-18 15:12

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

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

作者
Author: lycheeho 时间: 2026-3-18 15:31

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

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

作者
Author: Stardust0831 时间: 2026-3-18 16:06

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

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

作者
Author: UW_0728. 时间: 2026-3-18 17:41
本帖最后由 UW_0728. 于 2026-3-18 17:54 编辑

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

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

作者
Author: UW_0728. 时间: 2026-3-18 18:01

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

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

作者
Author: lycheeho 时间: 2026-3-18 18:51

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

明白了，感谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)