24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者ZeroBruce将赠送您 5 个金币

ZeroBruce

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 71.5
帖子: 18
在线: 9.1小时
虫号: 29710479
注册: 2022-06-16

[求助] 【求助】如何做二维材料差分电荷密度图？

各位大佬，我想咨询一下。

我想做单层锗烯和银离子的吸附构型，我下载导入了2D锗烯CIF文件后，首先通过CASTEP进行收敛性测试（启用色散矫正Grimme_DFT-D），得到了合适的Ecut（450eV左右）和Kpoint（8*8*1）。

随后几何优化（启用色散矫正Grimme_DFT-D，勾选了Cell优化晶格）。

得到优化后的构型后分别在H、T、B位添加银离子，又做了一次收敛性测试（启用色散矫正Grimme_DFT-D，此时设置Charge=1），Ecut和Kpoints参数和之前差不多。

随后对每个构型进行了几何优化（启用色散矫正Grimme_DFT-D，此时设置Charge=1）并计算单点能，发现H位的能量最低，结构最稳定。

之后在优化后稳定的构型，我选中了Ag+和几个与银离子成键的锗原子共同设置为1个Set，勾选EDD性质，计算单点能（启用色散矫正Grimme_DFT-D，此时设置Charge=1）。

得到计算结果后，我做出一个电荷差分密度图，发现Ag+周围是负值。感觉与预期不太相似，一般应该阳离子周围是正值才对，比如我看有一篇SA就是Al3+在锗烯附近是吸电子的。

所以我想问问我的建模是否正确？计算步骤有问题吗？谢谢！

回复此楼