计算化学公社

标题: 密度差局部积分曲线绘制疑问 [打印本页]

作者
Author: 古与春分 时间: 2025-9-1 10:18
标题: 密度差局部积分曲线绘制疑问
本帖最后由古与春分于 2025-9-1 14:48 编辑

老师，我用cp2k在revPBE/DZVP-MOLOPT-SR-GTH做能量计算得到四个cube文件，按照博文http://sobereva.com/638绘制密度差局部积分曲线，计算公式是ρ=ρ_total(r)-ρ_(sur+water)(r)-ρ_(slab+water)+ρ_(water)。然后发现绘制出的密度差局部积分曲线在slab区域有波动，请问老师这个合理吗？相关inp和pdb文件已上传，cube文件太大没有上传。

作者
Author: sobereva 时间: 2025-9-1 11:52
重新编辑帖子，把内容整理好

作者
Author: 古与春分 时间: 2025-9-1 14:48

sobereva 发表于 2025-9-1 11:52
重新编辑帖子，把内容整理好

老师，已经重新编辑

作者
Author: sobereva 时间: 2025-9-2 09:37
我完全不清楚你算的到底是什么密度差，算这种密度差要做什么
密度差求完后先在Multiwfn里观看密度差的等值面图，确保合理、有意义，之后再说z方向积分的事

作者
Author: 古与春分 时间: 2025-9-2 10:51

sobereva 发表于 2025-9-2 09:37
我完全不清楚你算的到底是什么密度差，算这种密度差要做什么
密度差求完后先在Multiwfn里观看密度差的等值 ...

老师，是密度差局部积分曲线，考察固液界面吸附导致的电子转移情况。等值面图是这样的，cubiso是0.004。

作者
Author: sobereva 时间: 2025-9-2 10:55

古与春分发表于 2025-9-2 10:51
老师，是密度差局部积分曲线，考察固液界面吸附导致的电子转移情况。等值面图是这样的，cubiso是0.004。

既然等值面没问题，根据局部积分曲线的定义，结合等值面图该怎么解释曲线怎么解释。取消显示分子结构，减小等值面数值，更有助于通过观看等值面理解微小的波动的来由。

作者
Author: 古与春分 时间: 2025-9-2 11:08

sobereva 发表于 2025-9-2 10:55
既然等值面没问题，根据局部积分曲线的定义，结合等值面图该怎么解释曲线怎么解释。取消显示分子结构，减 ...

好的老师，我再研究研究

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