| 查看: 1307 | 回复: 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倍。 请问 哪里出了问题？ 多谢您的关注和指点。 |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有65人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
湖南大学刘巧玲课题组2026年博士研究生招生信息
已经有28人回复
-
南方科技大学表面催化反应动力学研究课题组招收2026届博士生
已经有0人回复
» 抢金币啦！回帖就可以得到:
-
湖南大学化工院招收2026年入学博士
+1/86
-
有经费能报销的老师合作发表SCI开源期刊吗？
+1/82
-
太行山徒步之凌水河（一）
+1/79
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士（5月15号截止）
+1/39
-
上海大学 “生物有机电子材料及器件”团队副教授、博士后招聘
+1/37
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收博士生
+1/33
-
山东中医药大学2026年招聘有机化学、药物化学等相关方向师资博士后
+1/32
-
谷歌学术引用量
+1/32
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士（5月15号截止）
+1/29
-
又是一年基金申报，焦虑，某些学术大佬能不能少点干预
+1/16
-
类器官/器官芯片 有编制研究员
+1/11
-
211高校2026年入学博士名额
+1/11
-
2026年西安交通大学博士招生！
+1/10
-
中国科学院深圳先进院成会明院士—王锋项目高级工程师诚聘博士后(医用分离膜方向)
+1/10
-
北航国家优青万思杰团队招收材料、化学博士一名
+1/10
-
哈尔滨工业大学（深圳）——急招能动/材料/机械——26年9月份博士
+1/9
-
东北林业大学（211、双一流）招收高分子、化学、林业工程、生物质材料等方向博士生
+1/6
-
昆明理工大学交通载运车辆方向招收2026年秋季入学博士生2名
+4/4
-
南京大学医学院医学数字孪生方向博士后招聘
+1/2
-
博士生急招！哈尔滨工业大学（深圳）绿色化学团队招聘26秋和27年入学的博士生
+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
回复