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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 量化新手直通车 » 【求助】请教如何理解这句话中的概念
81/1返回列表
查看: 549  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 vigaryang 的 5 个金币

vigaryang

木虫 (正式写手)

[交流] 【求助】请教如何理解这句话中的概念

我看文献时发现有如下一句话：
The value of the electrostatic potential at the donor H nuclear position
(excluding the H‘s own infinite nuclear contribution), EPNUC , gives markedly superior fit statistics to any of the above charge and multipole models.
请问这里所说的EPNUC的含义是什么，它的值是怎样得到的呢？在gview中打开chk文件只能绘制得到一个笼统的静电势分布图，并没有给出空间某点处的静电势数值啊。谢谢！

[ Last edited by vigaryang on 2010-3-7 at 04:05 ]
回复此楼

» 猜你喜欢
1楼 2010-03-07 04:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★
yjcmwgk(金币+3): 2010-03-08 15:56
你需要“给出空间某点处的静电势数值”，我想还是应该通过cubegen命令或cube关键词先得到cube数据文件，这个三维数据文件就是点处的静电势数值，但没有直接给出坐标值，它的前面就是坐标展开的数据，当然读起来不便。

这时可以将这个cube数据文件读入GsGrid程序，它将输出一个为每一个静电势数值配齐坐标的文件：前三列为空间某点处的x、y、z坐标，后一列为该坐标点处的静电势数值。
赞一下(1人) 回复此楼
2楼2010-03-07 14:28:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★ ★
yjcmwgk(金币+2): 2010-03-08 15:56
EPNUC的意思就是指在原子核位置的静电势（不包括原子核对静电势的贡献，否则静电势是无限的），这并不是通用的缩写。

你可以生成静电势格点文件后用gsgrid的功能1提取每个格点位置的静电势然后找到相近的格点（或者直接手动指定网格设置计算某处静电势），也可以用multiwfn计算指定位置的静电势。
赞一下(1人) 回复此楼
3楼2010-03-07 14:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)
引用回帖:
Originally posted by sobereva at 2010-03-07 14:33:05:
EPNUC的意思就是指在原子核位置的静电势（不包括原子核对静电势的贡献，否则静电势是无限的），这并不是通用的缩写。

你可以生成静电势格点文件后用gsgrid的功能1提取每个格点位置的静电势然后找到相近的格点（ ...

非常感谢两位专家的帮助。我想请问一下，对于这个EPNUC，是否就是在gsgrid输出的output文件中寻找与Gaussian log 文件中该原子的笛卡尔坐标相近的一组X、Y、Z、Potential值？我发现并没有和该原子坐标完全吻合的点。如果如您所说进行指定计算某点处的静电势，该如何进行呢？是在gaussian的输入文件中就要指定么？
赞一下 回复此楼
4楼2010-03-08 23:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★ ★ ★
yjcmwgk(金币+3): 2010-03-09 12:48
我在当前gsgrid版本的readme里面最后写了如何自定义网格。文中的EPNUC是指扣除了原子核贡献的静电势（仅由电子云提供），但高斯输出的是完整的静电势，原子核处为无限大。

我可能不久后会在multiwfn 1.3里面加上计算电子云对某点静电势贡献的功能，届时这个问题就很容易了。
赞一下(1人) 回复此楼
5楼2010-03-09 01:20:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)
那再烦问前辈一下，既然gaussian给出的是总静电势，那么在当前条件下，有办法得到这个仅由电子云贡献的静电势值么？在我所看的那篇文献中，他只是说用gaussian 98计算了静电势，并没再描述如何扣掉原子核贡献的。
赞一下 回复此楼
6楼2010-03-09 01:53:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★
yjcmwgk(金币+1): 2010-03-12 13:22
引用回帖:
Originally posted by vigaryang at 2010-03-09 01:53:23:
那再烦问前辈一下，既然gaussian给出的是总静电势，那么在当前条件下，有办法得到这个仅由电子云贡献的静电势值么？在我所看的那篇文献中，他只是说用gaussian 98计算了静电势，并没再描述如何扣掉原子核贡献的。

在源代码里找到计算静电势的部分，把计算核贡献的段落注释掉并重新编译即可。
赞一下(1人) 回复此楼
7楼2010-03-09 19:41:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)
引用回帖:
Originally posted by sobereva at 2010-03-09 19:41:20:


在源代码里找到计算静电势的部分，把计算核贡献的段落注释掉并重新编译即可。

万分感谢！但能否再麻烦前辈告诉一下具体该怎么弄呢？我接触量化才两三个月时间，对很多东西几乎一无所知。这个内容对我很重要，也很急，因此只能恳求前辈都指点指点。
赞一下 回复此楼
8楼2010-03-10 20:13:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vigaryang 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 本人考085602 化学工程 专硕 +7 不知道叫什么！ 2026-03-15 8/400 2026-03-15 20:11 by 棒棒球手
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +3 Liwangman 2026-03-15 3/150 2026-03-15 18:16 by JourneyLucky
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
关闭
请填处理意见
关闭