24小时热门版块排行榜    

查看: 263  |  回复: 1
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 wandong_1984 的 8 个金币

wandong_1984

木虫 (正式写手)

[交流] 【求助】小分子双键电荷密度的模拟

现在有几个结构类似的都带双键的化合物,想从电子结构这块加以解释。
不知道能否用模拟的方法计算一下电荷密度?

万分感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wandong_1984(金币+2,VIP+0): 1-7 19:07
yjcmwgk(金币+8,VIP+0): 1-8 13:33
当然可以用模拟的方法计算一下电荷密度。如果要我作的话,我会进行如下操作:

1、如果几个结构类似的分子是A=B,可分别对它们进行优化。
然后将优化结构算单点能,并将双键相连的两个原子放在坐标的X轴方向,以便作图比较。
这时可以查看计算输出的out文件,找到第n个MO是要观察的π轨道。由于对称性不同、MO之间正交,π轨道是一个独立的轨道。如果分子中没有孤对电子,这π可能就是HOMO,如果有孤对电子,这π应在孤对电子之下。

2、算单点能时可用关键词cube=( Orbitals,cards),指定算这第n个MO。关于它的使用本论坛多次讨论过,我还编过一个手册,可找来一读。
这样就会得到三维空间π-MO格点文件cube。

将这cube输入GsGrid程序中,对cube进行平方,再次输入GsGrid程序中乘以2(MO是单电子的,乘以2是这个π-MO有两个电子)。这样就得到了新的π的三维空间电子密度格点文件cube。

将这新的cube输入GsGrid程序中,对双键作截面,就可以得到这双键的电子密度等值线图了。设A=B所在的平面为xy,可以作出它的双键相连的两个原子的Pz轨道融合成的MO的电子密度ρ,并选择xz平面(y=0)作为截面,可作出π电子密度ρ的等值线图。

几个结构类似的分子可以比较它们的电子密度ρ的不同。
如果不易看出区别,还要作下面的工作。

3、进一步还可作出π电子密度差Δρ的等值线图。甚至对Δρ进一步定量看有多少电荷量。这样更有利于相互之间比较。

[ Last edited by zhou2009 on 2010-1-7 at 17:50 ]
2楼2010-01-07 17:29:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wandong_1984 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 WqtG0Aem8Ufl 2026-10-24 6/300 2026-10-26 08:26 by 13760108691
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +8 7xHjyB3bAiwg 2026-10-24 8/400 2026-10-26 08:02 by Equinoxhua
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 4/200 2026-10-26 03:04 by cLOojpOxkX0G
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 6/300 2026-10-26 02:42 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3xN9SualEUPt 2026-10-24 6/300 2026-10-26 02:42 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:38 by cLOojpOxkX0G
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:29 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:25 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:12 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:05 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-26 01:50 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:39 by cLOojpOxkX0G
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:33 by cLOojpOxkX0G
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-26 01:24 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:45 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-26 00:28 by cLOojpOxkX0G
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 GVGRDgNfv7P4 2026-10-23 5/250 2026-10-26 00:08 by cLOojpOxkX0G
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
信息提示
请填处理意见