| 查看: 1219 | 回复: 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倍。 请问 哪里出了问题？ 多谢您的关注和指点。 |
回复此楼» 猜你喜欢
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
PbS量子点如何确定浓度
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有209人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
求一个访问学者邀请函，非常非常感谢
+1/675
-
坐标济南，来碰碰运气
+1/448
-
博士去军队文职怎么样
+5/175
-
上海师范大学生物医用材料方向招收2026级博士研究生
+1/78
-
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+1/77
-
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+1/75
-
深圳大学材料学院黄妍斐教授课题组诚招2026年秋季入学博士生
+1/75
-
南京理工大学曾海波/李伟金 招聘博士后（电磁响应：介电调控等方向）
+1/73
-
北京-89175-事业单位-诚征女友
+1/58
-
华中科技大学郭睿副教授建筑与能源科学课题组诚招博士研究生（2026.1.19截止）
+1/51
-
昆士兰科技大学（QUT）博士招生信息 导师：李志勇教授
+1/43
-
深圳信息职业技术大学-博后招聘（优秀可留校）
+1/29
-
南京邮电大学材料科学与工程学院柔性电子研究所2026年招聘公告
+1/15
-
太原理工大学电工部招聘老师-偏电类专业的博士们快来看啊
+1/11
-
福建师范大学海峡柔性电子学院-2026年9月入学博士-考核制
+1/10
-
湘潭大学沈志刚教授课题组欢迎博士生报考，博士后加盟
+1/9
-
武汉双一流高校干细胞与肿瘤生物学团队招聘2026级申请考核制博士生
+1/7
-
河南师范大学科研助理岗位招聘
+1/6
-
重庆医科大学-药学院-新靶标教育部医药基础研究创新中心-药物化学2026年博士招生
+1/4
-
中科院动物所招收2026年博士生（优先少干专项计划、化学或生命科学背景）
+1/1
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
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+3 2012-11-26 09:43:34
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+3 2012-11-26 09:43:34
|
看来你的命令格式还有问题：cube=(full,cards)。 要算轨道用关键词cube=(Orbitals,cards)；要算电子密度用cube=(density,cards)。 你现在并没有算电子密度，怎么能计算电子数。 还有，算电子密度用cube=(density,cards)，在cube或cubegen计算中，程序默认density=scf，即采用冻结核近似，不计算内层电子，如1s。如计算Al原子，总电荷数给出的是11。只有在指明fdensity=scf，才放弃冻结核近似，作包含全部电子的计算。这时计算Al原子，总电荷数给出的才是13。 另外，计算的电子数会不能完全与实际数值吻合，这是由于格点不够细密引起的误差。格点取到200*200*200就很接近了。 计算电子，已经不必用cube了，用关键词得到wfn文件，在multiwfn程序中更方便。 我已经好久没有用cube了，你的问题又让我把资料找出来复习了一遍。 [ Last edited by zhou2009 on 2012-11-25 at 15:44 ] |
» 本帖已获得的红花（最新10朵）
16楼2012-11-25 14:21:45
送鲜花一朵 |
多谢 zhou版主 再次解答。 我是新手，很多不清楚。多有麻烦。多谢了 目前还有一个疑问，烦请再次指教。 我在 g03 中文手册的 243页 废弃 关键词 cube 看到 如下解释： 立方空间内容选项 Density 只计算电子密度。不能与 Volume 关键字或Cube=Orbitals 选项组合。 FrozenCore 删除 SCF 芯电子密度。这是默认的电子密度，但不适用于势能计算。FC 是FrozenCore 的同义字。 Full 计算包括所有电子的电子密度。 Total 使用总电子密度。这是默认选项。 请问这里 full 选项 不能代替 density 吗，不是指 fdensity 吗？ 另外 这里 说的 total 与 density 选项 有何区别？ 另外：请问 版主 怎么设置 cube fdensity=scf, 是用cubegen 工具么？我设置 cube=(density,Cards) 能运行，而设置 fdenstiy 就不行了，老报错。 还有最后一个问题，怎么得到Al原子核心电子电荷密度分布文件？麻烦您用Al 用multiwfn 或者其他工具 给个指南。 多有麻烦！ 多谢多谢。 |
18楼2012-11-25 19:19:45
19楼2012-11-26 09:48:41
20楼2012-11-26 10:10:52
21楼2012-11-26 10:41:38
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
资料已发。 其实，计算分子的电子密度不是目的，它的图像是圆滚滚地，什么信息也看不见，最终目的是要作电子密度差Δρ。 作电子密度差Δρ： 见Multiwfn作者Sobereva博士写的的帖子“使用Multiwfn作电子密度差图”： http://muchong.com/bbs/viewthread.php?tid=3980314 |
22楼2012-11-26 15:23:54
23楼2012-11-27 09:33:36
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我算了一下你的输入数据。得到的wfn文件进入multiwfn中计算电子密度如下： 取低密度格点（125000）得电子数为25.57； 中密度格点（512000）得电子数为15.77； 高密度格点（1728000）得电子数为13.46。 如果进一步将格点取高，如200^3，结果就是13了。 可见你取的格点密度太稀疏。 为什么会这样，我曾经发帖讨论过类似的问题：“分子片ρ计算的一个值得注意的问题” http://muchong.com/bbs/viewthread.php?tid=2934836 |
24楼2012-11-27 10:50:31
25楼2012-11-27 11:44:46
简单回复
2012-11-24 20:48
回复
nodding(金币+1): 谢谢参与
2012-11-24 20:49
回复
nodding(金币+1): 谢谢参与
2012-11-24 20:49
回复
nodding(金币+1): 谢谢参与
假大空6楼
2012-11-25 10:18
回复
nodding(金币+1): 谢谢参与
neu2347楼
2012-11-25 10:37
回复
nodding(金币+1): 谢谢参与
neu2348楼
2012-11-25 10:38
回复
nodding(金币+1): 谢谢参与
neu2349楼
2012-11-25 10:38
回复
nodding(金币+1): 谢谢参与
userhung13楼
2012-11-25 11:35
回复
nodding(金币+1): 谢谢参与
userhung14楼
2012-11-25 11:35
回复
nodding(金币+1): 谢谢参与
userhung15楼
2012-11-25 11:36
回复
nodding(金币+1): 谢谢参与
轩辕重出江湖17楼
2012-11-25 18:52
回复