北京石油化工学院2026年研究生招生接收调剂公告
查看: 1376  |  回复: 2
【悬赏金币】回答本帖问题,作者kent1022将赠送您 66 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kent1022

金虫 (著名写手)

[求助] 计算原子能级怎样设置wf卡的对称性 已有1人参与

计算了一分子体系的能量,含有一个重原子As。想要验证选用的基组aug-cc-pwCV5Z-PP合不合适,特别是相对论效应考虑得够不够。
所以想到用aug-cc-pwCV5Z-PP基组计算As原子的基态及第一、第二激发态的能级能量,谱项分别是4S3/2,2D3/2,2D5/2,2P1/2,2P3/2
用了以下命令
geometry={As}
basis=aug-cc-pwcV5Z-PP
rhf
multi
ci
发现对单个原子的计算molpro用的D2h群,默认的Full valence occupancy是4   2   2   1   2   1   1   0
在wf中怎样设置4S3/2,2D3/2,2D5/2,2P1/2,2P3/2的对称性,后续还要进行SOC计算的。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)

【答案】应助回帖

在multi中做state averaged计算,molpro一般用D2h点群,所以只有1 - 8 个irrep, 其中 1 , 4 , 6 , 7 是g, 2, 3 ,5, 8是u。
计算As的基态,configuration为4s^2 4p^3 , s为g, p为u, 则 4s^2 4p^3  为u, 在2,3,5,8的irrep下做计算。

举例 如果active space只包括 4s4p, 则 ( 根据你的full valence occupancy 推断  是10-e 的PP?)
{multi; occ, 4   2   2   1   2   1   1   0; closed, 3   1   1   1   1   1   1   0;! open up 4s 4p
wf, 23, 8, 3; state, 1; !4S
}

关于第一第二激发态2D跟2P我有点忘了每个对称性下具体的态的数目了,但是你可以这么做来自己推断:
我们知道2D不算spin带来的degeneracy,那么degenerate的态的数目是 2*2+1 = 5个, 2P是 1*2+1 = 3个
加起来一共有8个态,那么就每个irrep下算尽可能多的态,然后自己看输出文件的能量对称性推断state的分布应该是怎么样的,
比如以下我每个irrep下算3个态(可能会多于最大态数量报错,可以自己调整)
{multi; occ, 4   2   2   1   2   1   1   0; closed, 3   1   1   1   1   1   1   0;! open up 4s 4p
wf, 23, 2, 1; state, 3;
wf, 23, 3, 1; state, 3;
wf, 23, 5, 1; state, 3;
wf, 23, 8, 1; state, 3;
}
计算得到每个irrep下面的态能量分布,3-fold degenerate就是2P,5-fold degenerate就是2D。
然后按照相应的irrep下面的态的数量填上正确的state数目,就可以算出来2P跟2D的能量了。
3楼2022-05-25 23:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 0703化学求调剂,各位老师看看我!!! +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 295求调剂 +4 1428151015 2026-03-27 5/250 2026-03-27 19:53 by JourneyLucky
[考研] 求调剂推荐 材料 304 +14 荷包蛋hyj 2026-03-26 14/700 2026-03-27 17:49 by kiokin
[考研] 329求调剂 +4 星野? 2026-03-26 4/200 2026-03-27 12:00 by 不吃魚的貓
[考研] 333求调剂 +3 question挽风 2026-03-23 3/150 2026-03-27 11:29 by 不吃魚的貓
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 085601求调剂总分293英一数二 +4 钢铁大炮 2026-03-24 4/200 2026-03-26 16:28 by dick_runner
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 299求调剂 +4 15188958825 2026-03-25 4/200 2026-03-25 22:56 by 418490947
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 各位老师您好:本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 求调剂 +6 研研,接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
信息提示
请填处理意见