| 查看: 1394 | 回复: 2 | |||
| 【悬赏金币】回答本帖问题，作者kent1022将赠送您 66 个金币 | |||
[求助]
计算原子能级怎样设置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计算的。 |
回复此楼» 猜你喜欢
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有93人回复
-
研究生招生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
-
中国电机工程学报的投稿系统
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
 |
2楼2022-02-13 16:25:50
【答案】应助回帖
|
在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
20