应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 关于VDD中原子的WFN文件
81/1返回列表
查看: 1747  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

172304

铁虫 (小有名气)

[求助] 关于VDD中原子的WFN文件

1.我在分析两个分子间弱相互作用是，用的是wb97xd/6-31g(d),而MULTIWFN中的VDD分析中给出的各个原子的WFN文件好像是HF/6-31G(D)方法的结果，是否会造成较大误差，（我试了一下发现总的VDD值不为零）。
2 是否应该在wb97xd/6-31g(d)方法下自己算一下各原子的WFN。
3 如需重新计算，各个原子的WFN计算时，输入参数（如电荷数，多重度）是指原子状态下的值吗（如对于H,多重度是2），等待各位解惑。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-10-13 07:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
172304: 金币+2, ★★★★★最佳答案 2012-10-14 21:30:50
自带的examples\atomwfn里的原子波函数文件对于主族是用ROHF计算的。VDD计算是基于对变形密度进行积分，这样的积分对于原子的密度在什么级别下产生的不敏感，这并不是导致总和不为0的原因。当然你也可以直接让Multiwfn调用Gaussian在你的计算级别下生成原子波函数文件（只支持对于前四周期元素），但我认为没这个必要。加和不为0是因为数值积分精度问题，VDD这种电荷想要计算精确很难，将settings.ini里的radpot和sphpot设为比默认更大的数值可以增加积分精度，例如分别设为130和974的话精度就很高了。
我个人不建议用VDD电荷。Hirshfeld电荷与之原理相近，物理意义更好，而且容易积分精确（因为Hirshfeld方法定义的原子空间之间是平滑的，而非VDD那样是梯形的，因此便于数值积分）。我也鼓励使用ADCH电荷，是对Hirshfeld进行后校正的结果，结果更合理。
赞一下 回复此楼
2楼2012-10-13 08:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

172304

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by sobereva at 2012-10-13 08:45:10
自带的examples\atomwfn里的原子波函数文件对于主族是用ROHF计算的。VDD计算是基于对变形密度进行积分，这样的积分对于原子的密度在什么级别下产生的不敏感，这并不是导致总和不为0的原因。当然你也可以直接让Multi ...

非常感谢。
回复此楼
3楼2012-10-14 21:30:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

172304

铁虫 (小有名气)
按照二楼的建议，我使用了HIRSHFEL和ADCH 布局分析，可是结果很不对劲，其中Zn上的电荷尽然为-50左右，这太离谱了。什么地方出了问题？
赞一下 回复此楼
4楼2012-10-15 22:39:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
引用回帖:
4楼: Originally posted by 172304 at 2012-10-15 22:39:06
按照二楼的建议，我使用了HIRSHFEL和ADCH 布局分析，可是结果很不对劲，其中Zn上的电荷尽然为-50左右，这太离谱了。什么地方出了问题？

将输入文件和具体操作步骤发到我邮箱sobereva@sina.com我看看
赞一下 回复此楼
5楼2012-10-16 03:38:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

172304

铁虫 (小有名气)
已发到你的sina邮箱，谢谢
赞一下 回复此楼
6楼2012-10-16 09:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

172304

铁虫 (小有名气)
我用2.6版重新计算了一遍，具体结果发到你的sina 邮箱了，Zn的电荷是0.几，但C、N、H的还是分别为4， 3， 1左右，还请你帮忙看看。
赞一下 回复此楼
7楼2012-10-18 20:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ruiying

铜虫 (初入文坛)
引用回帖:
2楼: Originally posted by sobereva at 2012-10-13 08:45:10
自带的examples\atomwfn里的原子波函数文件对于主族是用ROHF计算的。VDD计算是基于对变形密度进行积分，这样的积分对于原子的密度在什么级别下产生的不敏感，这并不是导致总和不为0的原因。当然你也可以直接让Multi ...

请教一下，某个金属算出VDD电荷是正值说明是得到电子啊？还是市区电子？还能说明什么啊？
赞一下 回复此楼
8楼2015-03-11 20:27:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 172304 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +17 花田半亩～白 2026-04-21 17/850 2026-04-23 07:30 by john198197
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[论文投稿] 急需审稿人！！！ +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
关闭
请填处理意见
关闭