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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】如何做电荷密度图和差分电荷密度图?
252/3返回列表上一页123下一页
查看: 5064  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[文学芳草园] 伙伴们，祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 一志愿中科院，化学方向，295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。，。，。，。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
关闭
请填处理意见
关闭