| 3 | 1/1 | 返回列表 |
| 查看: 1350 | 回復(fù): 2 | ||
| 【懸賞金幣】回答本帖問題,作者kent1022將贈送您 66 個金幣 | ||
[求助]
計算原子能級怎樣設(shè)置wf卡的對稱性 已有1人參與
|
||
|
計算了一分子體系的能量,含有一個重原子As。想要驗證選用的基組aug-cc-pwCV5Z-PP合不合適,特別是相對論效應(yīng)考慮得夠不夠。 所以想到用aug-cc-pwCV5Z-PP基組計算As原子的基態(tài)及第一、第二激發(fā)態(tài)的能級能量,譜項分別是4S3/2,2D3/2,2D5/2,2P1/2,2P3/2 用了以下命令 geometry={As} basis=aug-cc-pwcV5Z-PP rhf multi ci 發(fā)現(xiàn)對單個原子的計算molpro用的D2h群,默認(rèn)的Full valence occupancy是4 2 2 1 2 1 1 0 在wf中怎樣設(shè)置4S3/2,2D3/2,2D5/2,2P1/2,2P3/2的對稱性,后續(xù)還要進(jìn)行SOC計算的。 |
|
|
在multi中做state averaged計算,molpro一般用D2h點(diǎn)群,所以只有1 - 8 個irrep, 其中 1 , 4 , 6 , 7 是g, 2, 3 ,5, 8是u。 計算As的基態(tài),configuration為4s^2 4p^3 , s為g, p為u, 則 4s^2 4p^3 為u, 在2,3,5,8的irrep下做計算。 舉例 如果active space只包括 4s4p, 則 ( 根據(jù)你的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 } 關(guān)于第一第二激發(fā)態(tài)2D跟2P我有點(diǎn)忘了每個對稱性下具體的態(tài)的數(shù)目了,但是你可以這么做來自己推斷: 我們知道2D不算spin帶來的degeneracy,那么degenerate的態(tài)的數(shù)目是 2*2+1 = 5個, 2P是 1*2+1 = 3個 加起來一共有8個態(tài),那么就每個irrep下算盡可能多的態(tài),然后自己看輸出文件的能量對稱性推斷state的分布應(yīng)該是怎么樣的, 比如以下我每個irrep下算3個態(tài)(可能會多于最大態(tài)數(shù)量報錯,可以自己調(diào)整) {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下面的態(tài)能量分布,3-fold degenerate就是2P,5-fold degenerate就是2D。 然后按照相應(yīng)的irrep下面的態(tài)的數(shù)量填上正確的state數(shù)目,就可以算出來2P跟2D的能量了。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | 博斯特525 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +5 | 等春來, 2026-03-04 | 5/250 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 306求調(diào)劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|
|
[考研] 理學(xué),工學(xué),農(nóng)學(xué)調(diào)劑,少走彎路,這里歡迎您! +8 | likeihood 2026-03-02 | 11/550 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 環(huán)境調(diào)劑 +5 | 柒槿levana 2026-03-01 | 5/250 |
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 0856調(diào)劑 +10 | 劉夢微 2026-02-28 | 10/500 |
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 一志愿中石油(華東)本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考博] 26申博 +4 | 想申博! 2026-02-26 | 6/300 |
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|