| 查看: 1457 | 回复: 8 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
yyx19840628木虫 (著名写手)
鹰的眼睛，狼的耳朵，豹的速 ...
|
[交流]
【讨论】电子密度差计算
|
|||
|
这个是我下的教程上的讨论，我实在看不懂。它说Gaussian Test 083有算例，但那个算例中怎么没有cube关键字？能不能给个有注释的说明文件？抱歉，我实在很菜。 下面是教程内容： 以（H2O)2这个典型的分子为例，Gaussian Test 083有算例，并且是进行优化的：opt=tight。 采用关键词cube对（H2O)2算ρ。然后采用关键词cube分别算两个处于（H2O)2中位置独立的H2O的ρ。将它们相差，即得到形成氢键后电子的净变化。请问这步具体的该怎么做呢？ 此时，让氢键O……H－O在XY平面上，先优化出(H2O)2的构型，并同时得到一个电子密度的cube-1，再分别按优化构型的定位算两个独立的水分子，并同时得到两个电子密度的cube-2，cube-3。分别用gsgrid得到三个XY平面的电子密度的截面格点文件：g-1、g-2、g-3，它们都可以在Sigmaplot中分别绘出等值线图来。如果将g-1减去g-2和g-3，就得到二聚水的电子密度与形成氢键之前两个孤立分子的电子密度的差，即电子密度差Δρ。这些表述具体怎么转成输入文件，最后一步步画出图来？有趣的是，在这种相减之中，一些在形成氢键过程中基本没什么变化的内层电子，如1s、2s之类被抵消了，而提供孤对电子形成氢键的水分子的两个O－H键，不在XY平面上，另一个水分子的孤对电子因垂直于XY平面，在XY平面处为节点，都在XY平面上没有值。因此，这个电子密度差Δρ主要只是形成氢键前后的O……H－O电子净变化图像，为我们研究氢键提供信息。特别是相近系列氢键的Δρ变化将一目了然，人们对图形感知比较敏感，些微的差异就会发现，而从量化数据根本无从看起。 按行文习惯将Δρ作成虚实线图，也是很容易的事。 对于量化作图来说，是有章可循的操作，并不困难。倒是对于Δρ如何选择孤立的原子、分子、分子片，以及它们相应MO或ρ需要悟性和技巧，由此而选择设计分子坐标，得到合用的cube、合理的截面，至关重要。 这个是Gaussian Test 083算例#p rhf/6-31g* opt=tight test Gaussian Test Job 83: Water dimer optimization 0,1 O O,1,ROO X,1,1.,2,X3O H,1,RO1H,3,HOX3,2,90.,0 H,1,RO1H,3,HOX3,2,-90.,0 X,2,1.,1,52.5,3,180.,0 H,2,RO2H1,6,H7OX,1,180.,0 H,2,RO2H2,6,H8OX,1,0.,0 ROO=2.98308 RO1H=0.94839 X3O=120.2827 HOX3=52.90868 RO2H1=0.94686 RO2H2=0.95173 H7OX=52.98178 H8OX=51.9632 [ Last edited by yjcmwgk on 2009-11-14 at 13:34 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 理论和计算化学－淘贴专辑 |
» 猜你喜欢
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有136人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
-
求问：LATP煅烧应该使用什么容器？
已经有0人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有0人回复
-
西安交通大学 赵旭课题组 招收研究助理--电催化方向
已经有1人回复
-
西安交通大学 赵旭课题组 2026年招收博士后-电催化方向
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
loovfnd
至尊木虫 (著名写手)
- 应助: 35 (小学生)
- 贵宾: 1.111
- 金币: 14361.4
- 散金: 236
- 红花: 4
- 帖子: 2790
- 在线: 187小时
- 虫号: 219674
- 注册: 2006-03-14
- 专业: 物理有机化学
8楼2009-11-10 15:26:16
yyx19840628
木虫 (著名写手)
鹰的眼睛，狼的耳朵，豹的速 ...
- 应助: 4 (幼儿园)
- 贵宾: 0.41
- 金币: 3521.9
- 散金: 50
- 红花: 8
- 帖子: 2321
- 在线: 129小时
- 虫号: 508615
- 注册: 2008-02-21
- 性别: MM
- 专业: 环境污染化学
2楼2009-11-10 09:43:18
quantum999
金虫 (著名写手)
- 应助: 15 (小学生)
- 金币: 2146.3
- 散金: 39
- 帖子: 1426
- 在线: 198.2小时
- 虫号: 41283
- 注册: 2004-03-05
- 专业: 凝聚态物性 II ：电子结构
3楼2009-11-10 10:15:39
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 11-10 18:51
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 11-10 18:51
|
1、说“Gaussian Test 083有算例”，是说它有一个基本的二聚水坐标输入作单点能计算，不是说它在计算cube和Δρ。改变基组进一步优化之后，才作为计算cube和Δρ的基础。 2、采用关键词cube对（H2O)2和两个水算ρ，是采用关键词cube=(density,cards) 3、这些表述具体怎么转成输入文件，最后一步步画出图来？ GsGrid运算的多种数据，可以归总到Excel中然后去作图。由于我们对作密度差的分子和相应的孤立水分子在分别计算时，采用了同一套格点数、起始点、格点间隔，它们的坐标是完全一致的。因此当我们得到分子整体的GsGrid所取的截面后，在Excel中打开，前三列是x、y、z，第4列是ρ的值。这时可以将另外在Excel中分别打开的、孤立水的ρ值表格，依次将它们的ρ值复制分子整体的Excel中的第5列、第6列，在第7列列式，计算第4列减去第5列、第6列，则第7列就是Δρ的值，就可以作图。 安装sigmaplot后，在Excel的加载项中，就有sigmaplot，单击它就可进入sigmaplot作图。 [ Last edited by zhou2009 on 2009-11-10 at 15:17 ] |
4楼2009-11-10 10:48:02
20