24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1281  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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倍。

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

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)


★ ★ ★ ★
nodding(金币+1): 谢谢参与
nodding: 金币+3 2012-11-25 12:16:55
首先不知道楼主如何做的积分,其次27是不是相对质量囧
5楼2012-11-25 10:12:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答
★ ★ ★ ★ ★ ★ ★ ★
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的回帖

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): 谢谢参与
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 材料科学与工程调剂 +8 深V宿舍吧 2026-03-30 8/400 2026-03-30 16:19 by 晶体之美
[考研] 求调剂 +7 研研,接电话 2026-03-24 8/400 2026-03-30 14:14 by 白云朵朵飞
[考研] 337求调剂 +6 《树》 2026-03-29 6/300 2026-03-30 10:15 by herarysara
[考研] 327求调剂 +4 小卡不卡. 2026-03-29 4/200 2026-03-30 06:15 by wxiongid
[考研] 一志愿:西北大学,英一数一408-284分求调剂 +4 12.27 2026-03-27 4/200 2026-03-29 14:40 by zhshch
[考研] 086000生物与医药调剂 +5 Feisty。 2026-03-28 9/450 2026-03-29 12:02 by longlotian
[考研] 0856,材料与化工321分求调剂 +12 大馋小子 2026-03-27 13/650 2026-03-28 10:56 by self2008
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +5 dongfang59 2026-03-27 5/250 2026-03-28 04:53 by wxiongid
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 0703化学求调剂,各位老师看看我!!! +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 308求调剂 +7 墨墨漠 2026-03-25 7/350 2026-03-27 14:47 by 狂炫麦当当
[考研] 调剂推荐 +5 清酒714 2026-03-26 6/300 2026-03-27 11:12 by 不吃魚的貓
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 求调剂 +6 林之夕 2026-03-24 6/300 2026-03-27 08:38 by hypershenger
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
信息提示
请填处理意见