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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 关于VDD中原子的WFN文件
81/1返回列表
查看: 1679  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +7 小李0854 2026-03-16 7/350 2026-03-18 14:22 by laoshidan
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +3 1孙悟空 2026-03-17 3/150 2026-03-18 10:28 by macy2011
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 268求调剂 +7 好运连绵不绝 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[考研] 085600材料与化工 +4 安全上岸！ 2026-03-16 4/200 2026-03-17 14:02 by 勇敢太监王公公
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
关闭
请填处理意见
关闭