 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 2854 | 回復(fù): 9 | |||||
| 【獎勵】 本帖被評價5次，作者beefly增加金幣 4 個 | |||||
| 本帖產(chǎn)生 1 個 QC強帖 ，點擊這里進(jìn)行查看 | |||||
beefly專家顧問 (職業(yè)作家)
|
[資源]
在Gaussian09計算中考慮標(biāo)量相對論和旋軌耦合
|
||||
|
在Gaussian09中，加入了一種考慮旋軌耦合相對論效應(yīng)的方法DKHSO。原則上它支持所有的化學(xué)元素，只要有合適的基組。DKHSO目前僅能做Hatree-Fock級別的旋軌耦合計算。旋軌耦合DFT計算尚未充分測試，需要修改G09源代碼才能執(zhí)行。 下面簡單說說DKHSO的執(zhí)行方法。 1，標(biāo)量相對論計算 默認(rèn)情況下，DKHSO計算四階Douglas-Kroll-Hess單電子標(biāo)量相對論的哈密頓，文獻(xiàn)中一般叫做DKH4，選項為Integral=DKHSO。標(biāo)量DKH4可以與任何全電子的從頭算和DFT方法結(jié)合使用，但不能與半經(jīng)驗方法和ECP混用。注意，使用DKH4的方法都沒有解析梯度和解析頻率，做結(jié)構(gòu)優(yōu)化和振動頻率計算可能會非常耗時。 DKH4需要使用非收縮基組（計算非常耗時，甚至超過Gaussian程序的處理能力），或者用全電子相對論方法收縮的基組（通常用DKH2或DKH3收縮；DKH4收縮更好，但是很少見；ZORA收縮的芯軌道不佳，不要使用），用Gen關(guān)鍵詞讀入外部基組。這類基組可以在以下網(wǎng)站找到： https://bse.pnl.gov/bse/portal ANO-RCC基組，SARC-DKH，以及帶DK后綴的cc-pVnZ基組 http://tyr0.chem.wsu.edu/~kipeters/basis.html 帶DK后綴的基組 http://setani.sci.hokudai.ac.jp/sapporo/Welcome.do Sapporo-DK系列和DK3-Gen-Tk系列基組 此外，Hirao等人還在DKH3-HF級別做過1-103號元素的最小基，但是沒有應(yīng)用價值。 如果不追求高精度，18號以下的元素（定性計算甚至可以放寬到Kr元素）也可以用一般的非相對論基組。 注意：Gaussian的DKHSO選項強制使用有限尺寸原子核模型，而以上DKH基組都是用點原子核模型收縮的，會帶來微小的誤差。 例：Ar原子的標(biāo)量DKH4相對論Hartree-Fock計算。作為演示，這里使用非相對論基組3-21G。輸入： #p hf/3-21g int=dkhso test 0 1 Ar 計算得到的能量為-525.857790223 a.u. 2，旋軌耦合（也叫二分量）相對論計算 在標(biāo)量DKH4基礎(chǔ)上考慮旋軌耦合，需要做復(fù)數(shù)軌道Hartree-Fock計算，選項為GHF和Int=DKHSO。 如果有Gaussian09源代碼，還可以做復(fù)數(shù)軌道的 DFT計算（GB3LYP，等），只需把L502的一處出錯語句注釋掉。 很多年前人們就發(fā)現(xiàn)，單電子旋軌耦合積分和雙電子旋軌耦合積分的貢獻(xiàn)在同一數(shù)量級，缺一不可（J. Chem. Phys. 82, 5577, 1985）。但是雙電子旋軌耦合積分會急劇增加計算量。Gaussian的DKHSO是一個純單電子近似，結(jié)果會不會毫無意義呢？不用擔(dān)心，這是因為，DKHSO使用了Boettger提出的單電子旋軌耦合積分屏蔽公式，可以得到與嚴(yán)格旋軌耦合積分計算非常接近的結(jié)果。細(xì)節(jié)參見J. Chem. Phys. 120, 5875, 2004。 例：Ar原子的二分量DKH4相對論Hartree-Fock計算。 基組同上。輸入： #p ghf/3-21g int=dkhso test 0 1 Ar 考慮旋軌耦合以后，不再使用軌道的概念，而是使用旋量。旋量的滿占據(jù)數(shù)是1。由于GHF不使用對稱性，我們只能看到旋量的能量，看不到對稱符號。不過，對于原子，旋量對稱性（J）與旋量簡并性（D）的關(guān)系為D=2J+1。根據(jù)簡并性，并結(jié)合前面標(biāo)量DKH4計算輸出的軌道，不難猜出這些旋量的對稱性依次為：1s_1/2（兩個），2s_1/2（兩個），2p_3/2（四個），2p_1/2（兩個），等等。 occ. eigenvalues -- -118.26821-118.26821 -12.34229 -12.34229 -9.52602 occ. eigenvalues -- -9.52602 -9.52602 -9.52602 -9.48972 -9.48972 occ. eigenvalues -- -1.27958 -1.27958 -0.58938 -0.58938 -0.58938 occ. eigenvalues -- -0.58938 -0.58603 -0.58603 virt. eigenvalues -- 0.72886 0.72886 0.86389 0.86389 0.86389 virt. eigenvalues -- 0.86389 0.86772 0.86772 開殼層原子計算會破壞簡并性，無法得到旋量對稱性，只能去分析旋量的成分和空間分布。分子計算的情況類似。 這個計算得到的能量為-525.857816992 a.u.。和前一個計算結(jié)果比較，旋軌耦合使Ar原子的總能量下降約6個波數(shù)。 3，注意事項 以上旋軌耦合計算方法適合處理閉殼層的體系，以及非簡并的開殼層體系。如果是高對稱體系的簡并態(tài)，以上方法得到的結(jié)果可能很差。例如碘原子。這是HF單參考波函自身的缺陷造成的。 不要把“旋軌耦合效應(yīng)”與“旋軌耦合分裂”這兩個概念混淆。后者是在前者的基礎(chǔ)上，計算多個旋量態(tài)之間的能量差。DKHSO目前只能得到一個旋量態(tài)（通常是最低的），因此不能計算旋軌耦合分裂。 [ Last edited by beefly on 2012-10-6 at 09:42 ] 注意：Gaussian09所有的版本在計算DKH4旋軌耦合積分的步驟中都使用了錯誤的換算因子，導(dǎo)致旋軌耦合效應(yīng)被低估，無法重復(fù)文獻(xiàn)中的結(jié)果。需要修改源代碼才能正確執(zhí)行DKH4旋軌耦合計算。 [ Last edited by beefly on 2012-11-15 at 10:47 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化知識 |
» 本帖已獲得的紅花（最新10朵）
» 猜你喜歡
研究發(fā)現(xiàn)一個新的臨界系數(shù)，獨立于臨界壓縮因子
已經(jīng)有7人回復(fù)
-
200年來人們一直無法破解的難題： 氣液臨界點狀態(tài)方程
已經(jīng)有7人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有76人回復(fù)
-
細(xì)胞培養(yǎng)，這22個細(xì)節(jié)一定要注意�。ㄒ唬�
已經(jīng)有0人回復(fù)
-
關(guān)于引力延遲下雙星系統(tǒng)角動量不守恒的疑惑---求高手解答！
已經(jīng)有10人回復(fù)
-
CSC訪學(xué)博后項目獲批，外方學(xué)校暫停合作該怎么辦？
已經(jīng)有51人回復(fù)
-
分子動力學(xué)模擬合作需求
已經(jīng)有2人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 用Gaussian09執(zhí)行計算出現(xiàn)如圖的問題，怎么解決？
已經(jīng)有5人回復(fù)
- 請問大俠，Gaussian09中計算出來的能量單位 a.u 如何轉(zhuǎn)換成焦耳？
已經(jīng)有3人回復(fù)
- molpro里MRCI的能量到底該用哪個
已經(jīng)有6人回復(fù)
- 用MOLPRO構(gòu)造原子自然軌道（ANO）基函數(shù)
已經(jīng)有3人回復(fù)
- 怎么補算電聲子耦合？
已經(jīng)有3人回復(fù)
- 3萬以內(nèi)配運算Gaussian09的服務(wù)器，大伙給個配置唄！
已經(jīng)有4人回復(fù)
- 標(biāo)量相對論效應(yīng)是什么意思捏
已經(jīng)有3人回復(fù)
- 關(guān)于CASSCF和CASPT2計算時間
已經(jīng)有6人回復(fù)
- 關(guān)于gaussian09中BSSE校正
已經(jīng)有10人回復(fù)
- MolCAS計算的一些問題
已經(jīng)有8人回復(fù)
- Gaussian如何計算旋軌耦合常數(shù)
已經(jīng)有6人回復(fù)
- 學(xué)校買了gaussian09, 我若用g03的發(fā)文章，算侵權(quán)么？
已經(jīng)有7人回復(fù)
- 如何在計算中考慮范德華力？
已經(jīng)有6人回復(fù)
- 關(guān)于H4分子軌道問題
已經(jīng)有12人回復(fù)
- casscf計算空間占用
已經(jīng)有4人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 高斯計算中zORA方法是否可以使用
已經(jīng)有9人回復(fù)
- 請問MS里怎么考慮自旋軌道耦合作用？
已經(jīng)有9人回復(fù)
- 【求助】基組疊加誤差 （BSSE）
已經(jīng)有13人回復(fù)
- 【求助】td計算三重態(tài)問題
已經(jīng)有8人回復(fù)
- 【討論】大家計算重元素的性質(zhì)時有沒有考慮相對論效應(yīng)呢？
已經(jīng)有7人回復(fù)
2樓2012-10-06 08:31:28
3樓2012-10-06 09:23:18
4樓2012-10-06 11:17:34
beefly
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +458 - QC強帖: 12
- 應(yīng)助: 88 (初中生)
- 金幣: 13024.2
- 帖子: 4673
- 在線: 927.1小時
- 蟲號: 408372
|
旋軌耦合常數(shù)全稱是“旋軌耦合分裂常數(shù)”，上面說過，G09的GHF+DKHSO只能算單根，計算單個態(tài)（通常是基態(tài)）在旋軌耦合的影響下的總能量，因此不能得到旋軌耦合常數(shù)。如果復(fù)軌道方法將來能推廣到CIS，TDDFT，CASSCF，等等，就有可能算旋軌耦合常數(shù)。 旋軌耦合效應(yīng)受動態(tài)相關(guān)能、基組級別的影響不大，只要是適合HF描述的單參考體系就能計算。為了解決HF的總能量精度不高的問題，可以在GHF+DKHSO和HF+DKHSO級別分別計算總能量，得到旋軌耦合對總能量的校正項，疊加到其他方法的總能量上，例如CCSD/ECP。 這種方法類似于NWChem和Turbomole中的SODFT方法，不過后者是基于DFT+旋軌耦合ECP，計算量小很多。所有這些方法都不能算旋軌耦合常數(shù)。 |
5樓2012-10-06 21:53:32
6樓2012-10-07 00:10:51
送鮮花一朵 | 為什么Gaussian09計算旋軌耦合效應(yīng)一定要先進(jìn)行標(biāo)量相對論計算，我用Molcas的計算完全不采用RCC基組，就是ANO-L基組，不過原子都很小，第二周期的，這樣誤差會有多大呢。本身HF的精度就不高。molcas的casscf級別計算還是不錯的，怪不得我在caspt2基礎(chǔ)上計算旋軌耦合出來的數(shù)值和不包含caspt2模塊的一樣，不知道是本身誤差就不大，還是沒有開發(fā)到caspt2的層次。說明書上說好像目前只能計算casscf 水平的。另外有限尺寸原子核模型和點原子核模型是怎么回事呢。哪個更接近實際情況，可能原子太小沒有影響，只會造成一點點數(shù)據(jù)上的誤差。 |
7樓2012-10-07 11:50:35
beefly
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +458 - QC強帖: 12
- 應(yīng)助: 88 (初中生)
- 金幣: 13024.2
- 帖子: 4673
- 在線: 927.1小時
- 蟲號: 408372
|
1，有些比較輕的元素，標(biāo)量相對論效應(yīng)和一階旋軌耦合都不大。這種情況一般可以忽略全部相對論效應(yīng)。 例子有H2、HF的離解能。 2，有些比較輕的元素，標(biāo)量相對論效應(yīng)不大，但是一階旋軌耦合比較大。這種情況可以只考慮后者。 例子有HBr的離解能。Br的旋軌耦合使離解能下降3 kcal/mol。 3，有些比較重的元素，標(biāo)量相對論效應(yīng)很大，但是一階旋軌耦合很小或是0。這種情況可以只考慮前者。 例子有Hg2的離解能。 4，有些比較重的元素，標(biāo)量相對論效應(yīng)和一階旋軌耦合都很大。這種情況二者都要考慮。 例子有I2的離解能。 5，Rn之后的重元素，即便一階旋軌耦合為零，二階旋軌耦合通常也很大，因此二者都要考慮。 |
8樓2012-10-07 23:28:52
beefly
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +458 - QC強帖: 12
- 應(yīng)助: 88 (初中生)
- 金幣: 13024.2
- 帖子: 4673
- 在線: 927.1小時
- 蟲號: 408372
|
很少用molcas算旋軌耦合，因為感覺不如molpro方便。molpro早先允許用caspt2波函計算旋軌耦合，但是不太穩(wěn)定，不存在入侵態(tài)問題的情況下，有時候會得到奇怪的結(jié)果，所以我一般都用casscf或mrci波函旋軌耦合。后來的某些版本干脆禁用caspt2波函算旋軌耦合，不知道現(xiàn)在解決了沒有。 真實的原子核是有尺寸的，因此有限尺寸原子核更接近真實情況。用了有限尺寸原子核還有其他一些好處。首先，Dirac方程不再有原子序數(shù)小于137的限制。此外，有限尺寸原子核可以減小基函數(shù)的指數(shù)，即使構(gòu)造超重元素的高級別基組，指數(shù)一般也不會超過10^9。但是如果用點核模型構(gòu)造同樣級別的基組，指數(shù)往往會達(dá)到10^11以上，在做雙精度計算的時候容易出現(xiàn)數(shù)值不穩(wěn)定的問題。 |
9樓2012-10-07 23:51:27
簡單回復(fù)
83644936610樓
2013-06-05 17:15
回復(fù)
五星好評 
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 beefly 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分調(diào)劑 +10 | Faune 2026-03-30 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
|
[考研]
 304材料求調(diào)劑
+5
|
鐘llll 2026-03-26 | 5/250 |
|
|
|
[有機交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研]
一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+10
|
y7czhao 2026-03-26 | 10/500 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|