请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 121|回复: 1

[Multiwfn使用咨询] bug report

[复制链接]

133

帖子

8

威望

1075

eV
积分
1368

Level 4 (黑子)

发表于 2017-3-16 08:16:47 | 显示全部楼层 |阅读模式
版本3.3.9。如下wfx文件,做电子密度拓扑分析。某临界点的G,K,V,H值:


Lagrangian kinetic energy G(r):  0.3771912731
Hamiltonian kinetic energy K(r):  0.2080016661
Potential energy density V(r): -0.5770091772
Energy density E(r) or H(r): -0.2080016661


发现H 不等于G + V,而是K直接取负值。换其他分子的wfx文件做类似的分析,未发现此错误。

0.wfx

365.27 KB, 下载次数: 1

test

8443

帖子

16

威望

1万

eV
积分
19271

管理员

公社社长

发表于 2017-3-17 00:09:44 | 显示全部楼层
谢谢报告。
导致问题的原因是计算V的时候是根据V=(1/4)der2rho-2G的关系给出的,但是EDF字段只对der2rho有贡献而对G没有贡献。当不读取EDF字段的时候就没这个问题。
等过几天我会在官网上上传新的3.4(dev)版,把V的计算方式改成了V=-K-G,就确保H=V+G=-K了
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://Multiwfn.codeplex.com
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395(2000人已满),思想家公社QQ群2号:466017436(讨论计算化学为主,加入时必须注明研究方向)

Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

Archiver|手机版|小黑屋|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2017-3-28 12:15 , Processed in 0.083838 second(s), 24 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表