| 查看: 1443 | 回复: 24 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
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倍。 请问 哪里出了问题？ 多谢您的关注和指点。 |
回复此楼» 猜你喜欢
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有169人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有68人回复
-
申博
已经有2人回复
5楼2012-11-25 10:12:04
★ ★ ★ ★ ★ ★ ★ ★
nodding(金币+1): 谢谢参与
nodding: 金币+5 2012-11-25 12:17:03
gmy1990: 金币+2 2012-11-26 09:41:59
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
11楼2012-11-25 11:06:22
|
首先：多谢 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
简单回复
假大空6楼
2012-11-25 10:18
回复
nodding(金币+1): 谢谢参与
userhung13楼
2012-11-25 11:35
回复
nodding(金币+1): 谢谢参与