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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】如何做电荷密度图和差分电荷密度图?
252/3返回列表上一页123下一页
查看: 5460  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

zhou2009

版主 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yuhuobuku(金币+2,VIP+0):谢谢交流 7-7 13:56
生成电荷密度时，用SCF。

在控制行中加关键词cube=(Orbitals,cards)

或者调用cubegen命令来作。

[ Last edited by zhou2009 on 2009-7-4 at 14:59 ]
赞一下(11人) 回复此楼
11楼2009-07-04 08:46:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
也就是说直接优化结构，然后生成FCHECK文件，再生成 CUBE文件就可以了？可是有大侠讨论过，说SCF生成的电荷分布不太合理，甚至是有些离谱，建议用NBO布居分析的结果来代替，那么，用这样的数据来讨论是否会有问题呢？
赞一下(1人) 回复此楼
12楼2009-07-04 09:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+5,VIP+0):感谢交流 7-4 20:21
是的，直接优化结构，然后生成FCH文件，再生成 CUBE文件就可以了。
但是如果是作Δρ，直接用关键词CUBE，就直接生成CUBE文件了。

说SCF生成的电荷分布不太合理，甚至是有些离谱，建议用NBO布居分析的结果来代替，那么，用这样的数据来讨论是否会有问题呢？

这是两回事。

SCF生成的电荷分布不太合理，甚至是有些离谱，不是SCF的问题，而是将SCF的数据此后进一步分析为集居数产生的问题，如将成键原子中间的电子均分给原子双方、更严重的是将σ、π、pi、δ等电子的变化头发胡子一把抓地归结到一个原子上形成净电荷造成的，实与SCF无关。

NBO布居分析也是为了解决将电荷划归原子而作出的后处理，也还远不能令人满意，同样会生成的电荷分布不太合理，甚至是有些离谱。并没有从根本上解决问题。

这些，都是习惯将电荷人为地划分给原子惹的祸。

而ρ、Δρ它直接从SCF来，它是空间分布在原子周围的，ρ是真实的形象表达，Δρ是如实的科学抽象的图形表达，原子上的净电荷还原为电荷在空间的净变化了，都没有什么人为造成的错处。
它既不将成键原子中间的电子均分给原子双方，它就是一个空间如实的分布。还会将σ、π、pi、δ等电子的变化分别去作出图来研究、考察。

为此，我正在写一篇这方面的文章，专门讲电荷。

[ Last edited by zhou2009 on 2009-7-4 at 10:19 ]
赞一下(11人) 回复此楼
13楼2009-07-04 09:47:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
希望能早日看到你的大作。想学习理论计算，看来要学的东西还很多。不过目前就你看来，对于原子电荷，哪种理论方法更合理一些？毕竟时常有人借原子电荷的分布来讨论反应的亲电性或亲核性。
赞一下(1人) 回复此楼
14楼2009-07-04 09:55:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+3,VIP+0):感谢交流 7-4 20:21
我也一直在寻找，但我认为，目前把分子中电荷归结为原子的净电荷方法尚没有合理的、合用的。
也许传统的一定要把电荷人为划归原子本来就不是合理的，不反映客观存在。

所以可以用ρ图或Δρ图来看。

[ Last edited by zhou2009 on 2009-7-4 at 15:24 ]
赞一下(11人) 回复此楼
15楼2009-07-04 10:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+3,VIP+0):感谢交流 7-5 17:47
引用回帖:
Originally posted by zhou2009 at 2009-7-4 10:25:
我也一直在寻找，但我认为，目前把分子中电荷归结为原子的净电荷方法尚没有合理的、合用的。
也许传统的一定要把电荷人为划归原子本来就不是合理的，不反映客观存在。

所以可以用ρ图或Δρ图来看。

[ La ...

菜鸟的问题：
NBO分析可同时产生AOs和MOs（同时有PNAOs和PNMOs）能否利用这些轨道来计算电子密度图和差图？这样可以省略单算AOs的过程，同时保证所有的数据均在同等环境下产生。
下图是用NBO的一个π轨道经GSGRID平方乘2生成的图形。请大侠指教。
赞一下(11人) 回复此楼
16楼2009-07-05 11:10:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+3,VIP+0):感谢交流 7-5 17:47
这种方式我没有试过。

我想，关键是看能否进入cubegen中，另外计算出在空间格点的每一个点的电子密度值来，有了这些格点值，才能选取截面作出等值线来。

仅仅是给出一个值、或球棒式的表达，是不能计算出在空间格点的每一个点的电子密度值的。

这种截面或说剖面，比球棒式的表达更加深入内部。而且才能作出密度差来。

[ Last edited by zhou2009 on 2009-7-5 at 14:18 ]
赞一下(11人) 回复此楼
17楼2009-07-05 14:14:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+3,VIP+0):感谢交流 7-5 17:47
我是这样考虑的：
1. 在CHEMCRAFT中，可以选择任意的轨道生成CUB文件（选择PNAOs和PNMOs），这种文件可以被GV识别并作图，但不知这是不是标准的或格式化的CUB文件。
2. 利用GSGIRD将轨道文件平方（成键轨道再乘以2），变成电子密度文件（？）
3. 利用GV可以查看3D或指定平面的等高线图（如上帖之图）
4. 利用GSGRID提取指定平面的电子密度数据（成键轨道^2及相应的原子轨道^2）
5. 利用其它数据处理软件作等高线密度差图。
不知这样是否合理？有没有确定的意义？
赞一下(11人) 回复此楼
18楼2009-07-05 15:38:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+3,VIP+0):感谢交流 7-5 17:47
关键一步是：选择任意的轨道生成CUB文件（选择PNAOs和PNMOs）。

我想如果运行cubegen，它自动去找指定SCF的MO，不会去找PNAOs和PNMOs。PNAOs和PNMOs可能是cubegen不能识别的。

我还认为SCF的MO比PNMOs要好。PNMOs只是经过处理使原子净电荷可能好一点。真正的最原始的信息还是SCF的MO。
赞一下(11人) 回复此楼
19楼2009-07-05 16:48:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★
erylingjet(金币+3,VIP+0):感谢交流 7-5 17:47
引用回帖:
Originally posted by yytsnake at 2009-7-5 15:38:
我是这样考虑的：
1. 在CHEMCRAFT中，可以选择任意的轨道生成CUB文件（选择PNAOs和PNMOs），这种文件可以被GV识别并作图，但不知这是不是标准的或格式化的CUB文件。
2. 利用GSGIRD将轨道文件平方（成键轨道再乘以2），变成电子密度文件（？）
3. 利用GV可以查看3D或指定平面的等高线图（如上帖之图）
4. 利用GSGRID提取指定平面的电子密度数据（成键轨道^2及相应的原子轨道^2）
5. 利用其它数据处理软件作等高线密度差图。
不知这样是否合理？有没有确定的意义？

如果是我，我会这样做，yytsnake不要见笑：
1. 在Gaussian中，可以选择任意的轨道，利用cubegen生成CUBE文件（MO）
2. 在Gaussian中，可以选择density，利用cubegen生成CUBE文件（电子密度图）
3. 利用GV可以查看3D或指定平面的等高线图
4. 利用GV提取指定平面的电子密度数据（成键轨道^2及相应的原子轨道^2）
5. 利用GV作等高线密度差图。
赞一下(10人) 回复此楼
20楼2009-07-05 16:49:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 erylingjet 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖？ +4 CrisMessi 2026-08-18 4/200 2026-08-19 00:53 by gatelove
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了，能看出什么问题？ +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息，中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
关闭
请填处理意见
关闭