请问CP2K可以准确计算带电荷的体系的能量吗？

xexlalalan

准备从VASP转到CP2K。VASP计算带电荷体系的结果不可信，不知道CP2K是否存在这个问题。

biogon

算算不就知道了

jiangning198511

你所指的带电荷体系具体指什么体系，一般bulk带电荷体系VASP是可以的，好像表面带电荷就有问题，不知道CP2K是如何解决这个问题的@卡开发发  

wolfli369

表面带电体系也能处理，例如厦门大学田中群院士的文章https://www.nature.com/articles/s41563-019-0356-x。但是这类模拟，对初始结构要求较高

卡开发发

jiangning198511 发表于 2021-4-17 10:15
你所指的带电荷体系具体指什么体系，一般bulk带电荷体系VASP是可以的，好像表面带电荷就有问题，不知道CP2K ...

听说vasp的新版本能处理charged slab。也许更改CP2K_INPUT/FORCE_EVAL/DFT/POISSON的PBC（如XY）是可行的，不过我不确定具体行不行，可以试着做一些测试。

jiangning198511

wolfli369 发表于 2021-4-17 10:40
表面带电体系也能处理，例如厦门大学田中群院士的文章https://www.nature.com/articles/s41563-019-0356-x ...

文章中使用了CP2K，但体系整体为电中性，通过不同阳离子来调节电极表面，这样的体系VASP也可以处理

wolfli369

jiangning198511 发表于 2021-4-17 11:15
文章中使用了CP2K，但体系整体为电中性，通过不同阳离子来调节电极表面，这样的体系VASP也可以处理

这种处理对阳离子种类有要求吗？和实际研究体系关系大吗？

jiangning198511

wolfli369 发表于 2021-4-17 21:33
这种处理对阳离子种类有要求吗？和实际研究体系关系大吗？

这个具体没有研究过，你可以查查那篇文章(nature)相关作者的文章 看看

wolfli369

jiangning198511 发表于 2021-4-19 09:13
这个具体没有研究过，你可以查查那篇文章(nature)相关作者的文章 看看

好的，谢谢你的指导

九月九

卡开发发 发表于 2021-4-17 11:12
听说vasp的新版本能处理charged slab。也许更改CP2K_INPUT/FORCE_EVAL/DFT/POISSON的PBC（如XY）是可行的 ...

找到了这篇文章，用VASP算的：PHYSICAL REVIEW LETTERS 126, 076401 (2021)摘要：Supercell models are often used to calculate the electronic structure of local deviations from the idealperiodicity in the bulk or on the surface of a crystal or in wires. When the defect or adsorbent is charged, ajellium counter charge is applied to maintain overall neutrality, but the interaction of the artificially repeatedcharges has to be corrected, both in the total energy and in the one-electron eigenvalues and eigenstates.This becomes paramount in slab or wire calculations, where the jellium counter charge may inducespurious states in the vacuum. We present here a self-consistent potential correction scheme and providesuccessful tests of it for bulk and slab calculations

Zhangpc

jiangning198511 发表于 2021-4-17 10:15
你所指的带电荷体系具体指什么体系，一般bulk带电荷体系VASP是可以的，好像表面带电荷就有问题，不知道CP2K ...

请问是不是在水盒子bulk里面加一个OH-离子，VASP就能算呢？

jiangning198511

Zhangpc 发表于 2022-9-27 11:42
请问是不是在水盒子bulk里面加一个OH-离子，VASP就能算呢？

bulk应该没问题