24小时热门版块排行榜    

查看: 2773  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:02 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见