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

计算化学公社

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

[Multiwfn使用咨询] bug report

[复制链接]

202

帖子

12

威望

1558

eV
积分
2000

Level 5 (御坂)

发表于 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

1万

帖子

16

威望

1万

eV
积分
24760

管理员

公社社长

发表于 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了
北京科音自然科学研究中心:http://www.keinsci.com  不定期开办各层次量子化学、分子动力学、Multiwfn程序培训
思想家公社的门口Blog:http://sobereva.com
Multiwfn量子化学波函数分析程序主页:http://sobereva.com/multiwfn
计算化学公社论坛:http://bbs.keinsci.com
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2017-10-24 12:03 , Processed in 0.093182 second(s), 24 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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