求助：按照multiwfn手册4.13.6计算电荷位移曲线（CDC），图像右边不为零

wdg741

整个体系是苯甲醇和水，想知道问题出在哪

Uus/pMeC6H4-/キ

geom=connectivity，pressure和temperature都对当前单点能计算没有作用，建议去掉并就后两者阅读http://sobereva.com/423一文。
电场严重极化电子结构的情形下，基组不带弥散函数应该会使描述定性不准。
field=z+300对应的电场强度是300*0.0001=0.03 a.u.，换算成V/m达到10^10数量级了，确认现实中的苯甲醇和水体系不会被这么强的电场击穿么？确认体系此时没有电离、波函数稳定么？体系是孤立团簇还是周期性溶液，在算这个单点能之前有无强电场下的几何优化？

sobereva

细节根本没什么都没清楚，算的是谁和谁的密度差也不提、完整的输入文件也没给，具体怎么产生的密度差格点数据都只字未提，不予回复。提问时候别敷衍了事，认真把此文看了：
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

除了手册外，下文该说的都说了
使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线
http://sobereva.com/638（http://bbs.keinsci.com/thread-28225-1-1.html）

wdg741

sobereva 发表于 2025-4-13 01:49
细节根本没什么都没清楚，算的是谁和谁的密度差也不提、完整的输入文件也没给，具体怎么产生的密度差格点数 ...

抱歉，描述不清楚。我这个回复加了附件，用multiwfn主功能5之后“z+300.wfn”减去“bjc-h2o.wfn”计算电子密度差（中等质量），之后返回主菜单，用主功能13→18并选择Z轴得到电荷转移曲线与局部积分曲线。

wdg741

Uus/pMeC6H4-/キ 发表于 2025-4-12 23:56
geom=connectivity，pressure和temperature都对当前单点能计算没有作用，建议去掉并就后两者阅读http://sob ...

谢谢，我再去试一试

sobereva

你这两个wfn文件都有问题，Multiwfn主功能100的子功能4对电子密度全空间积分都明显偏离整数。应当直接用formchk转出来的fch当Multiwfn输入文件，要么用out=wfn/wfx直接由Gaussian产生的wfn/wfx文件。不要用非主流、奇怪的方式产生wfn文件，否则自找麻烦。

认真阅读
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）

wdg741

sobereva 发表于 2025-4-13 14:53
你这两个wfn文件都有问题，Multiwfn主功能100的子功能4对电子密度全空间积分都明显偏离整数。应当直接用for ...

我是用multiwfn主功能100转化的文件，不过听您的直接用fch文件计算已经可以了，谢谢老师。

sobereva

wdg741 发表于 2025-4-13 20:08
我是用multiwfn主功能100转化的文件，不过听您的直接用fch文件计算已经可以了，谢谢老师。

你之前绝对不是这么转换的
你的wfn文件的第一行的文字暴露了真相，Multiwfn转出来的文件根本不会有那一行。Multiwfn产生出的wfn第一行写的是 Generated by Multiwfn

wdg741

sobereva 发表于 2025-4-14 01:35
你之前绝对不是这么转换的
你的wfn文件的第一行的文字暴露了真相，Multiwfn转出来的文件根本不会有那一 ...

抱歉，检查了一下，用来计算的文件，当时嫌麻烦是用的“fch2wfn”。