| 查看: 1378 | 回复: 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计算的。 |
回复此楼» 猜你喜欢
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有236人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
-
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有17人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
070300总分298求调剂
已经有2人回复
-
26博士申请
已经有0人回复
-
一志愿厦门大学化学学硕307求调剂
已经有8人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有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
