24小时热门版块排行榜    

查看: 1353  |  回复: 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的回帖

kent1022

金虫 (著名写手)

2楼2022-02-13 16:25:50
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 欢迎211本科同学,过A区国家线,A区非偏远一本,交叉学科课题组 +9 lisimayy 2026-03-04 13/650 2026-03-05 06:53 by chenvision
[考研] 求调剂院校 +6 云朵452 2026-03-02 12/600 2026-03-04 23:17 by 云朵452
[考研] 293求调剂 +3 是乐渝哇 2026-03-04 3/150 2026-03-04 23:12 by wutongshun
[考研] 化学 0703求调剂 总分293 一志愿211 +8 土土小虫 2026-03-03 10/500 2026-03-04 23:04 by 土土小虫
[考研] 320材料与化工,求调剂 +6 鹤遨予卿 2026-03-04 8/400 2026-03-04 20:47 by wutongshun
[考研] 一志愿314求调剂 +7 202111120625 2026-03-03 7/350 2026-03-04 15:56 by zhukairuo
[考研] 347求调剂 +5 啊欧欧欧 2026-03-03 5/250 2026-03-04 14:10 by kakakapanpan
[考研] 325求调剂 +5 学家科 2026-03-04 5/250 2026-03-04 14:04 by kakakapanpan
[考研] 材料325求调剂 30+5 mariusuki 2026-03-02 10/500 2026-03-04 09:24 by njzyff
[论文投稿] EST拒稿重投 5+3 15102603076 2026-03-02 3/150 2026-03-04 00:51 by bobvan
[考研] 293求调剂 +4 是乐渝哇 2026-03-03 4/200 2026-03-03 23:04 by zhukairuo
[基金申请] 请问大家,研究风险与应对措施那里, 大家都怎么写呢 ? +3 cauasen 2026-03-02 3/150 2026-03-03 19:43 by Tide man
[考研] 283求调剂 +9 鹿沫笙 2026-03-02 9/450 2026-03-03 15:58 by tgxtgxtgx9
[考研] 一志愿华中科技大学,化学专业344分,求调剂 +6 邢xing1 2026-03-02 6/300 2026-03-03 15:46 by stewie_jz
[考研] 调剂 +5 13853210211 2026-03-02 7/350 2026-03-03 11:42 by 13853210211
[考研] 288求调剂 +3 少71.8 2026-03-02 5/250 2026-03-03 06:01 by tgxtgxtgx9
[考研] 0856求调剂285 +11 吕仔龙 2026-02-28 11/550 2026-03-02 20:15 by hypershenger
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
信息提示
请填处理意见