| 查看: 1299 | 回复: 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计算的。 |
回复此楼» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有278人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有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