请教CP2K怎么正确计算赝氢

djjj148

准备用赝氢钝化表面模型，同一个练手的小模型，用VASP的赝氢可以把表面的悬挂键全部钝化(算出来的PDOS的带隙和形状几乎和bulk的完全一致，费米能级附近也没有中间态)，但是用CP2K总是捣鼓不出来，由于实际计算的模型较大，需要用CP2K。以下是自己做的一些调研和尝试：


http://bbs.keinsci.com/thread-11093-1-1.html
八楼有提到用CORE_CORRECTION 关键词可以设置，但有时候仍会出现中间态，我在练手的时候就出现了，故根据开发者的指示优化赝势，但是目前没有找到优化赝氢/其他元素同位素的赝势的输入文件，只能根据现有的资料和tests/ATOM/regtest-pseudo中的例子自己修改，尝试了好几天，算出来的DOS仍然会出现中间态。


故向大家求助~