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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gauss cube=(full,Cards) 的得到电子总数 不等于 原子总电子数目。
51/1返回列表
查看: 1446  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

nodding

木虫 (小有名气)

[交流] gauss cube=(full,Cards) 的得到电子总数 不等于 原子总电子数目。

我计算单个原子 Al 的电子分布情况，
输入文件为：
# b3lyp/3-21g cube=(full,cards)

Title Card Required

0 2
Al 0. 0. 0.

full.cube
-1 -4.0495 -4.0495 -4.0495
64 0.130629 0.0 0.0
64 0.0 0.130629 0.0
64 0.0 0.0 0.130629

结果把 full.cube 中所用格点的电子密度积分后，发现电子总数 为27 ，比Al 原子电子总数 高出约 2倍。

请问 哪里出了问题？ 多谢您的关注和指点。
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-11-24 20:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★
nodding(金币+1): 谢谢参与
nodding: 金币+5 2012-11-25 12:17:03
gmy1990: 金币+2 2012-11-26 09:41:59
你如何把 full.cube 中所用格点的电子密度积分的？
cube已经是离散的数字了，需要把每一个数字乘以微体积元，再将它们加和起来，就得到电子数。这在excel中很容易操作。
因为cube是用原子单位表示的，微体积元要用波尔为单位的数据。
如果用埃为单位的，数值就会差约一倍。你是否误用了埃为单位的微体积元？
在multiwfn的处理结果中，用波尔为单位的数据和用埃为单位的都列出了。
赞一下 回复此楼
10楼2012-11-25 11:03:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

lishijunzong

铁杆木虫 (正式写手)
★ ★ ★ ★
nodding(金币+1): 谢谢参与
nodding: 金币+3 2012-11-25 12:16:55
首先不知道楼主如何做的积分，其次27是不是相对质量囧
赞一下 回复此楼
5楼2012-11-25 10:12:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nodding

木虫 (小有名气)
引用回帖:
5楼: Originally posted by lishijunzong at 2012-11-25 10:12:04
首先不知道楼主如何做的积分，其次27是不是相对质量囧

我把 格点电荷密度分布* 每个格点体积 然后求和  得到 总的电子数。

27不是相对质量。
赞一下 回复此楼
11楼2012-11-25 11:06:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nodding

木虫 (小有名气)
引用回帖:
10楼: Originally posted by zhou2009 at 2012-11-25 11:03:34
你如何把 full.cube 中所用格点的电子密度积分的？
cube已经是离散的数字了，需要把每一个数字乘以微体积元，再将它们加和起来，就得到电子数。这在excel中很容易操作。
因为cube是用原子单位表示的，微体积元要用 ...

首先：多谢 zhou 版主 光顾 指点。

列出full.cbue 几行数据, 便于 高人 诊断 指点。

Title Card Required
SCF Total Density
    1   -7.652446   -7.652446   -7.652446
   64    0.246853    0.000000    0.000000
   64    0.000000    0.246853    0.000000
   64    0.000000    0.000000    0.246853
   13   13.000000    0.000000    0.000000    0.000000
  9.58947E-12  1.54167E-11  2.44023E-11  3.80283E-11  5.83476E-11  8.81410E-11
  1.31091E-10  1.91958E-10  2.76746E-10  3.92820E-10  5.48967E-10  7.55332E-10
  1.02322E-09  1.36471E-09  1.79205E-09  2.31686E-09  2.94909E-09  3.69587E-09

full.cube 输出按照 原子单位给出的。对应的 电荷密度 e/Bohr^3, 格点间隔为 bohr。 没有 转化为 埃 为单位。

我 按照 版主的 方法 积分的。首先用 full.cube 格点电荷密度 * 每格点的体积 即 0.246853^3， 然后 再求和，结果 为27.207.

这个数值还是不对啊。要不我附上 我的 full.cube 文件。看您方便处理一下吗？
full.cube(3.34MB)
http://kuai.xunlei.com/d/GCVZLZQVYXAY?p=130497

[ Last edited by nodding on 2012-11-25 at 12:19 ]
赞一下 回复此楼
12楼2012-11-25 11:27:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
假大空6楼
2012-11-25 10:18   回复  
nodding(金币+1): 谢谢参与
userhung13楼
2012-11-25 11:35   回复  
nodding(金币+1): 谢谢参与
查看全部 25 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭