24小时热门版块排行榜    

查看: 2171  |  回复: 5

tephoon78

木虫 (正式写手)

[交流] 【求助】关于guess=alter的一个不明白的地方。 已有5人参与

在gaussian说明书guess关键字中有一个例子,关于amino 自由基的计算。
第一次计算使用 guess=only,产生的轨道如下

注意右边的那一列,也就是β自旋轨道4和5,4是a1对称性的,5是b1对称性的。这时候5(b1)的轨道能量比4(a1)的高。



第二次进行计算,使用guess=alter,将4和5的轨道交换,按照书上说的是2A1的激发态。轨道如下:


我的问题是为什么轨道5的能量比4要高?(这里的轨道5相当于第一次计算的轨道4)


个人理解:guess=alter调节轨道,换成了计算激发态,那么电子应该占据能量高一点的轨道,而空出一个能量低一点的轨道。我们一般好像接受的激发态就是这样的一个概念,为什么计算结果会这样呢?是我计算错了吗?
回复此楼

» 猜你喜欢

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

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

lingznpmg

铁虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhou2009(金币+2): 2010-11-08 16:50:28
楼主的问题很好!
我记得有朋友说过,高斯在计算中只有占据轨道的电子参加SCF自洽,空轨道不会参与。空轨道的轨道能好像是通过计算一个电子的激发能得到的,并不是很准。

我自己不是很明白,随后也好好研究一下!
2楼2010-11-07 23:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)

楼上说的很有道理
莫强求!
3楼2010-11-08 12:19:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhou2009(金币+2): 2010-11-08 16:46:27
个人的理解:这里说的激发态并不是激发电子,而是同一个电子结构下的两种态。
产生这两种态的原因并不是由于电子激发,而是由于轨道对称性匹配的变化而产生的。
氨基自由基的基态是2B1,这也是默认guess得到的组态。但是如果想要得到2A1组态,那么必须要改变初始猜测,也就是guess=alter,就是说,楼主你的第二张图已经和前一张没关系了,是在算氨基自由基的另一个态,所以轨道的能量顺序还是遵循既定原则。
Anewday,anewhour,anewminute,anewpeople.
4楼2010-11-08 16:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhou2009(金币+5): 2010-11-09 14:10:10
LS说得很对。这种利用对称性不同通过SCF过程计算激发态的方法也叫ΔSCF。只有用guess=(alter,only)的时候才会出现你预期的结果,也就是能量低的轨道没有被占据,形成“空穴”态。
当涉及到迭代过程,在第一次迭代后占据的βb1会立刻降低,βa1空轨道能量会立刻上升并高于βb1,然后程序自动根据aufbau定律按能量由低到高重新排布电子,βb1就是占据的,而βa1就是空的。
这是由于每一步的单电子有效势算符是由上一步的占据轨道构建的,这算符的本征函数中与上一步占据轨道相似者能量会比较低。在HF中,收敛的占据轨道电子感受到的是n-1个电子的互斥势,而空轨道感受到的是全部n个电子互斥势,明显后者更大,故空轨道能量更高。
5楼2010-11-08 22:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tephoon78

木虫 (正式写手)

感谢楼上的回复!!很详细!
6楼2010-11-09 13:52:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tephoon78 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +4 小贾同学123 2026-03-15 6/300 2026-03-17 07:06 by laoshidan
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 4/200 2026-03-16 23:52 by zbzihdhd
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
信息提示
请填处理意见