24小时热门版块排行榜

返回列表

qzhaosdu

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 1401.6
帖子: 1969
在线: 263.5小时
虫号: 493351

[交流] 【讨论】关于G03里guess=only的一个疑惑。

G03里GUESS关键词里说到guess=alter的用二重态NH2举例。
首先，用了guess=only得出初始猜测，说明书里是这样的 Initial guess orbital symmetries.
   Alpha Orbitals:
         Occupied (A1) (A1) (B2) (B1) (A1)
         Virtual  (A1) (B2)
   Beta Orbitals:
         Occupied (A1) (A1) (B2) (A1)
         Virtual  (B1) (A1) (B2)
但是我算了一下却是
Orbital symmetries:
Alpha Orbitals:
   Occupied  (A1) (A1) (B2) (A1) (B1)
   Virtual (A1) (B2)
Beta  Orbitals:
   Occupied  (A1) (A1) (B2) (A1)
   Virtual (B1) (A1) (B2)
也就是说说明书里的占据轨道里4和5是反着的。
但是用了OPT UHF去做的时候，程序在算第一步SCF的时候却给的是和说明书一样的结果。
正常来说，Alpha和Beta的对称性应该是一一对应相同的吧？（我也不知道这样说对不对？）
**************************
另外我刚看了下G09的说明书，GUESS里给的是正常的结果，也就是说G03里写错了？
**************************
另外还有，我算了一个Ni+...C2H2的C2V结构，guess=only给的猜测是2-B2，但是b3lyp/6-31+g*算的时候给的初始猜测却是2-B1.
**************************
所以说，G03里这个guess=only这个关键词是不是有问题啊？（起码第一个SCF是有错误的）

[ Last edited by qzhaosdu on 2010-12-9 at 17:50 ]

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

关于药动学大鼠灌胃的问题已经有5人回复
求几个中国搞纳米材料的牛的重点实验室@一些关于纳米氧化锌的文章（名字即可）！已经有28人回复
关于 MS 5.5并行的问题已经有4人回复
晚辈诚惶诚恐请教：制甲醇，在线色谱产物分析计算方法？已经有7人回复
大侠们请教个问题：离心可否除去粘稠液体里的小分子？已经有10人回复
浙江省自然科学基金委员会关于表彰2010-2011年度先进集体、先进个人的通报已经有6人回复
求助关于酯化反应催化剂的问题已经有5人回复
求大虾指导一个关于atoi函数的问题已经有3人回复
关于二氯甲烷的问题已经有21人回复
关于苯乙炔端位炔基两个碳的带电情况已经有3人回复
关于格氏试剂与羧酸衍生物的反应控制已经有8人回复
求考研调剂指导已经有6人回复
求助一个氨基保护反应已经有4人回复
【讨论】关于反应自发的一个疑惑已经有18人回复
【求助】关于guess=alter的一个不明白的地方。已经有5人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-09 16:56:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ykwang

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 3171.7
帖子: 639
在线: 386.4小时
虫号: 500302

qzhaosdu(金币+10):王老师，第二个我写错了，是alpha和beta轨道，这样的话是一一对应相同的吗？ 2010-12-09 17:51:24

第一个问题：不同版本的程序所用的初始猜测方法可能有所不同(即使方法相同也可能参数或算法有所不同)。不同的方法给出的初始猜测不一定相同，但只要方法合适，优化后的结果应该是一致的。

第二个问题：“正常来说占据轨道和非占据轨道的对称性应该是一一对应相同的吧？”这一理解是错误的！首先因为占据轨道的数目与空轨道的数目一般是不同的，故不可能建立一一对应关系。其次，即使二者的数目相同，也没有任何物理或几何原因能在二者之间建立这种关系。

第三个问题：有些开壳层体系特别是自由基，初始猜测的对称性常常是与基组有关的。原因在于你使用了单Slater行列式近似，而这种近似对于有些开壳层体系并不是好的近似。