【求助】怎样计算一个原子的轨道能级 - 量子化学 - 量化新手直通车

11楼2009-12-25 15:54:59

dy_414

木虫 (正式写手)

应助: 20 (小学生)
金币: 2316.5
散金: 26
红花: 4
帖子: 390
在线: 29.8小时
虫号: 830319
注册: 2009-08-18
性别: GG
专业: 基因表达调控与表观遗传学

不错不错，顶一下！

我坚信，rp是攒出来的

12楼2009-12-25 15:57:00

gnli

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 298.1
散金: 20
帖子: 273
在线: 24.1小时
虫号: 294810
注册: 2006-11-11
性别: MM
专业: 理论和计算化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

S ( 1)
   3S : 4.49    4S :12.72    5S :42.66    6S :10.08    0PZ: 1.10
C ( 2)
   2S : 1.82    2PZ: 1.07    3S :13.70    3PZ: 2.51    4S : 3.46
   5S : 1.46

请问小卒兄这些 S P 轨道的数据是咋分析出来的？
还有个问题是：
CS只有两个原子，分析起来简单一些，直接看分子轨道图就能判断是sigma键还是pi键，以及每个原子贡献的S or P 轨道
如果大一些的有机分子，只用gview看分子轨道图感觉很晕，怎么知道某个分子轨道到底是哪些原子的那种轨道的贡献呢？
不知道我表述清楚了没有

争取早日成为合格的菜鸟

13楼2009-12-29 04:31:35

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 553.8
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by gnli at 2009-12-29 04:31:
S ( 1)
   3S : 4.49    4S :12.72    5S :42.66    6S :10.08    0PZ: 1.10
C ( 2)
   2S : 1.82    2PZ: 1.07    3S :13.70    3PZ: 2.51    4S : 3.46
   5S : 1.46

请问小卒兄这些 S P ...

你添加iop(6/7=3)就行啦，呵呵

14楼2010-01-04 15:50:55

dym851206

应助: (幼儿园)
在线:
虫号: 771862

长见识了

15楼2010-01-13 09:29:10

imunbg

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 4469.5
散金: 1
帖子: 159
在线: 35.9小时
虫号: 879742
注册: 2009-10-21
性别: GG
专业: 凝聚态物性I:结构、力学和

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

本人也关注相同的问题，请问你们用的是什么软件？

学习学习再学习

16楼2010-01-13 12:37:15

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 553.8
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by imunbg at 2010-1-13 12:37:
本人也关注相同的问题，请问你们用的是什么软件？

gaussian

17楼2010-01-13 12:46:50

zhangmt

至尊木虫 (著名写手)

我叫MT

QC强帖: 5
应助: 99 (初中生)
金币: 6961.8
散金: 10406
红花: 49
帖子: 1761
在线: 763.2小时
虫号: 880392
注册: 2009-10-22
性别: GG
专业: 理论和计算化学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by yjcmwgk at 2010-1-4 15:50:

你添加iop(6/7=3)就行啦，呵呵

这个iop(6/7)=3代表什么意思？用它就可以直接给出各个MO中对应的每个AO的百分比么？如果是的话，那可太省事了啊。请小卒版主指教。

一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以

18楼2010-01-13 13:22:55

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 553.8
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by zhangmt at 2010-1-13 13:22:
这个iop(6/7)=3代表什么意思？用它就可以直接给出各个MO中对应的每个AO的百分比么？如果是的话，那可太省事了啊。请小卒版主指教。

晕啊，不是iop(6/7)=3，而是iop(6/7=3)。

这个命令不是给出各个MO对应的每个AO的百分比，而是给出各个MO对应的每个波函的组合系数。（请zhangmt先生注意，组合系数并不是百分比，要对组合系数进行平方归一化才能得到百分比；而且，波函并不对应真实的AO，真实的AO对应的是波函的线性组合）
现在解释iop(6/7=3)是什么意思
Overlap6中，对IOP(7)的解释是MOLECULAR ORBITAL COEFFICIENTS. 即iop(6/7)指的是对“组成MO的各个波函”的组合系数的打印控制
当这个值为3的时候，输出全部组合系数。如果这个值为1（默认值），则输出从HOMO-4到LUMO+4这10个轨道的组合系数
我看您还是下载一本Gaussian iop大全吧，至少作为工具书查阅呀，呵呵

什么是iop？iop是高斯任务的内部控制字符。比如下面这一段源码（从L602.F中复制过来的,对应的是IOp(6/14=?)的程序执行过程）：

引用回帖:

   If(IOp(14).eq.-1) then
      MaxDer = -1
   else if(IOp(14).eq.0) then
      If(DoChrg) then
      MaxDer = 0
      else
      MaxDer = 2
      endIf
   else if(IOp(14).eq.1) then
      MaxDer = 2
   else if(IOp(14).eq.2) then
      MaxDer = 1
   else if(IOp(14).eq.3) then
      MaxDer = 0
   else if(IOp(14).eq.4) then
      MaxDer = -1
   else
      MaxDer = IOp(14)
      endIf
   If(MaxDer.lt.-1.or.MaxDer.gt.2) then
      Write(IOut,1030) MaxDer
      Call Lnk1E(0)
      endIf

你可以看到，全部命令均由iop来控制，事实上我们敲的那些命令行都被翻译成了iop才能被高斯执行（把高斯关键字翻译成iop，这就是L1.exe做的事情）

[ Last edited by yjcmwgk on 2010-1-13 at 14:01 ]

19楼2010-01-13 13:52:22

zhangmt

至尊木虫 (著名写手)

我叫MT

QC强帖: 5
应助: 99 (初中生)
金币: 6961.8
散金: 10406
红花: 49
帖子: 1761
在线: 763.2小时
虫号: 880392
注册: 2009-10-22
性别: GG
专业: 理论和计算化学

★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+6,VIP+0):呵呵，有一段时间高斯公司隐藏了iop(6/7)，因为用户抱怨说用了iop(6/7=3)之后，输出文件奇大无比（这也确实是事实），结果高斯公司真的就隐藏了它了。哈哈 1-14 16:16
yjcmwgk(金币+0,VIP+0):有时候高斯公司也够恶搞的。把iop(6/7)弄得跟游戏里的隐藏任务似的。哈哈。zhangmt先生 1-14 16:18

iop控制我知道，手边有gaussian iop详解，不过没有找到6/7的说明。
应该先找源码去看看再问，呵呵，有点懒了。
当然应该用系数平方和归一化之后再算百分比了——显然你写出来的是计算后的百分比——他提问的也是针对百分比的数据怎么出来的——我以为用6/7就可以直接输出百分比而不用自己算了呢，呵呵。想偷懒而已。
我的gaussian iop详解是03年4月版的，里面6开头的iop是从15-83的，没有iop(6/7)，可能是版本太低了……小卒有新版的没？共享一下。
刚查了601的源码，找到了iop(6/7)的说明，呵呵。。
刚去gaussian网站看了，gaussian09的iop里确实有iop(6/7)，俺old了。
再次确认，06年9月版的gaussian03的iop里面也没有提供iop(6/7)的介绍。
不过在早期B05的源码中就有iop(6/7)的解释了，看来这就是传说中的隐藏参数吧。

[ Last edited by zhangmt on 2010-1-13 at 15:51 ]

赞一下(21人)

20楼2010-01-13 14:34:11