 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2857 | 回復(fù): 9 | |||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)5次，作者beefly增加金幣 4 個(gè) | |||||
| 本帖產(chǎn)生 1 個(gè) QC強(qiáng)帖 ，點(diǎn)擊這里進(jìn)行查看 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
beefly專家顧問 (職業(yè)作家)
|
[資源]
在Gaussian09計(jì)算中考慮標(biāo)量相對(duì)論和旋軌耦合
|
||||
|
在Gaussian09中，加入了一種考慮旋軌耦合相對(duì)論效應(yīng)的方法DKHSO。原則上它支持所有的化學(xué)元素，只要有合適的基組。DKHSO目前僅能做Hatree-Fock級(jí)別的旋軌耦合計(jì)算。旋軌耦合DFT計(jì)算尚未充分測試，需要修改G09源代碼才能執(zhí)行。 下面簡單說說DKHSO的執(zhí)行方法。 1，標(biāo)量相對(duì)論計(jì)算 默認(rèn)情況下，DKHSO計(jì)算四階Douglas-Kroll-Hess單電子標(biāo)量相對(duì)論的哈密頓，文獻(xiàn)中一般叫做DKH4，選項(xiàng)為Integral=DKHSO。標(biāo)量DKH4可以與任何全電子的從頭算和DFT方法結(jié)合使用，但不能與半經(jīng)驗(yàn)方法和ECP混用。注意，使用DKH4的方法都沒有解析梯度和解析頻率，做結(jié)構(gòu)優(yōu)化和振動(dòng)頻率計(jì)算可能會(huì)非常耗時(shí)。 DKH4需要使用非收縮基組（計(jì)算非常耗時(shí)，甚至超過Gaussian程序的處理能力），或者用全電子相對(duì)論方法收縮的基組（通常用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級(jí)別做過1-103號(hào)元素的最小基，但是沒有應(yīng)用價(jià)值。 如果不追求高精度，18號(hào)以下的元素（定性計(jì)算甚至可以放寬到Kr元素）也可以用一般的非相對(duì)論基組。 注意：Gaussian的DKHSO選項(xiàng)強(qiáng)制使用有限尺寸原子核模型，而以上DKH基組都是用點(diǎn)原子核模型收縮的，會(huì)帶來微小的誤差。 例：Ar原子的標(biāo)量DKH4相對(duì)論Hartree-Fock計(jì)算。作為演示，這里使用非相對(duì)論基組3-21G。輸入： #p hf/3-21g int=dkhso test 0 1 Ar 計(jì)算得到的能量為-525.857790223 a.u. 2，旋軌耦合（也叫二分量）相對(duì)論計(jì)算 在標(biāo)量DKH4基礎(chǔ)上考慮旋軌耦合，需要做復(fù)數(shù)軌道Hartree-Fock計(jì)算，選項(xiàng)為GHF和Int=DKHSO。 如果有Gaussian09源代碼，還可以做復(fù)數(shù)軌道的 DFT計(jì)算（GB3LYP，等），只需把L502的一處出錯(cuò)語句注釋掉。 很多年前人們就發(fā)現(xiàn)，單電子旋軌耦合積分和雙電子旋軌耦合積分的貢獻(xiàn)在同一數(shù)量級(jí)，缺一不可（J. Chem. Phys. 82, 5577, 1985）。但是雙電子旋軌耦合積分會(huì)急劇增加計(jì)算量。Gaussian的DKHSO是一個(gè)純單電子近似，結(jié)果會(huì)不會(huì)毫無意義呢？不用擔(dān)心，這是因?yàn)椋珼KHSO使用了Boettger提出的單電子旋軌耦合積分屏蔽公式，可以得到與嚴(yán)格旋軌耦合積分計(jì)算非常接近的結(jié)果。細(xì)節(jié)參見J. Chem. Phys. 120, 5875, 2004。 例：Ar原子的二分量DKH4相對(duì)論Hartree-Fock計(jì)算。 基組同上。輸入： #p ghf/3-21g int=dkhso test 0 1 Ar 考慮旋軌耦合以后，不再使用軌道的概念，而是使用旋量。旋量的滿占據(jù)數(shù)是1。由于GHF不使用對(duì)稱性，我們只能看到旋量的能量，看不到對(duì)稱符號(hào)。不過，對(duì)于原子，旋量對(duì)稱性（J）與旋量簡并性（D）的關(guān)系為D=2J+1。根據(jù)簡并性，并結(jié)合前面標(biāo)量DKH4計(jì)算輸出的軌道，不難猜出這些旋量的對(duì)稱性依次為：1s_1/2（兩個(gè)），2s_1/2（兩個(gè)），2p_3/2（四個(gè)），2p_1/2（兩個(gè)），等等。 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 開殼層原子計(jì)算會(huì)破壞簡并性，無法得到旋量對(duì)稱性，只能去分析旋量的成分和空間分布。分子計(jì)算的情況類似。 這個(gè)計(jì)算得到的能量為-525.857816992 a.u.。和前一個(gè)計(jì)算結(jié)果比較，旋軌耦合使Ar原子的總能量下降約6個(gè)波數(shù)。 3，注意事項(xiàng) 以上旋軌耦合計(jì)算方法適合處理閉殼層的體系，以及非簡并的開殼層體系。如果是高對(duì)稱體系的簡并態(tài)，以上方法得到的結(jié)果可能很差。例如碘原子。這是HF單參考波函自身的缺陷造成的。 不要把“旋軌耦合效應(yīng)”與“旋軌耦合分裂”這兩個(gè)概念混淆。后者是在前者的基礎(chǔ)上，計(jì)算多個(gè)旋量態(tài)之間的能量差。DKHSO目前只能得到一個(gè)旋量態(tài)（通常是最低的），因此不能計(jì)算旋軌耦合分裂。 [ Last edited by beefly on 2012-10-6 at 09:42 ] 注意：Gaussian09所有的版本在計(jì)算DKH4旋軌耦合積分的步驟中都使用了錯(cuò)誤的換算因子，導(dǎo)致旋軌耦合效應(yīng)被低估，無法重復(fù)文獻(xiàn)中的結(jié)果。需要修改源代碼才能正確執(zhí)行DKH4旋軌耦合計(jì)算。 [ Last edited by beefly on 2012-11-15 at 10:47 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化知識(shí) |
» 本帖已獲得的紅花（最新10朵）
» 猜你喜歡
研究發(fā)現(xiàn)一個(gè)新的臨界系數(shù)，獨(dú)立于臨界壓縮因子
已經(jīng)有7人回復(fù)
-
200年來人們一直無法破解的難題： 氣液臨界點(diǎn)狀態(tài)方程
已經(jīng)有7人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有210人回復(fù)
-
細(xì)胞培養(yǎng)，這22個(gè)細(xì)節(jié)一定要注意�。ㄒ唬�
已經(jīng)有0人回復(fù)
-
關(guān)于引力延遲下雙星系統(tǒng)角動(dòng)量不守恒的疑惑---求高手解答！
已經(jīng)有10人回復(fù)
-
CSC訪學(xué)博后項(xiàng)目獲批，外方學(xué)校暫停合作該怎么辦？
已經(jīng)有51人回復(fù)
-
半透性容器進(jìn)行高濕或者低濕研究
已經(jīng)有0人回復(fù)
-
分子動(dòng)力學(xué)模擬合作需求
已經(jīng)有2人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 用Gaussian09執(zhí)行計(jì)算出現(xiàn)如圖的問題，怎么解決？
已經(jīng)有5人回復(fù)
- 請(qǐng)問大俠，Gaussian09中計(jì)算出來的能量單位 a.u 如何轉(zhuǎn)換成焦耳？
已經(jīng)有3人回復(fù)
- molpro里MRCI的能量到底該用哪個(gè)
已經(jīng)有6人回復(fù)
- 用MOLPRO構(gòu)造原子自然軌道（ANO）基函數(shù)
已經(jīng)有3人回復(fù)
- 怎么補(bǔ)算電聲子耦合？
已經(jīng)有3人回復(fù)
- 3萬以內(nèi)配運(yùn)算Gaussian09的服務(wù)器，大伙給個(gè)配置唄！
已經(jīng)有4人回復(fù)
- 標(biāo)量相對(duì)論效應(yīng)是什么意思捏
已經(jīng)有3人回復(fù)
- 關(guān)于CASSCF和CASPT2計(jì)算時(shí)間
已經(jīng)有6人回復(fù)
- 關(guān)于gaussian09中BSSE校正
已經(jīng)有10人回復(fù)
- MolCAS計(jì)算的一些問題
已經(jīng)有8人回復(fù)
- Gaussian如何計(jì)算旋軌耦合常數(shù)
已經(jīng)有6人回復(fù)
- 學(xué)校買了gaussian09, 我若用g03的發(fā)文章，算侵權(quán)么？
已經(jīng)有7人回復(fù)
- 如何在計(jì)算中考慮范德華力？
已經(jīng)有6人回復(fù)
- 關(guān)于H4分子軌道問題
已經(jīng)有12人回復(fù)
- casscf計(jì)算空間占用
已經(jīng)有4人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 高斯計(jì)算中zORA方法是否可以使用
已經(jīng)有9人回復(fù)
- 請(qǐng)問MS里怎么考慮自旋軌道耦合作用？
已經(jīng)有9人回復(fù)
- 【求助】基組疊加誤差 （BSSE）
已經(jīng)有13人回復(fù)
- 【求助】td計(jì)算三重態(tài)問題
已經(jīng)有8人回復(fù)
- 【討論】大家計(jì)算重元素的性質(zhì)時(shí)有沒有考慮相對(duì)論效應(yīng)呢？
已經(jīng)有7人回復(fù)
4樓2012-10-06 11:17:34
luolun2008
至尊木蟲 (著名寫手)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 20299.7
- 帖子: 2601
- 在線: 1106.9小時(shí)
- 蟲號(hào): 579917
3樓2012-10-06 09:23:18
beefly
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +458 - QC強(qiáng)帖: 12
- 應(yīng)助: 88 (初中生)
- 金幣: 13024.2
- 帖子: 4673
- 在線: 927.1小時(shí)
- 蟲號(hào): 408372
|
旋軌耦合常數(shù)全稱是“旋軌耦合分裂常數(shù)”，上面說過，G09的GHF+DKHSO只能算單根，計(jì)算單個(gè)態(tài)（通常是基態(tài)）在旋軌耦合的影響下的總能量，因此不能得到旋軌耦合常數(shù)。如果復(fù)軌道方法將來能推廣到CIS，TDDFT，CASSCF，等等，就有可能算旋軌耦合常數(shù)。 旋軌耦合效應(yīng)受動(dòng)態(tài)相關(guān)能、基組級(jí)別的影響不大，只要是適合HF描述的單參考體系就能計(jì)算。為了解決HF的總能量精度不高的問題，可以在GHF+DKHSO和HF+DKHSO級(jí)別分別計(jì)算總能量，得到旋軌耦合對(duì)總能量的校正項(xiàng)，疊加到其他方法的總能量上，例如CCSD/ECP。 這種方法類似于NWChem和Turbomole中的SODFT方法，不過后者是基于DFT+旋軌耦合ECP，計(jì)算量小很多。所有這些方法都不能算旋軌耦合常數(shù)。 |
5樓2012-10-06 21:53:32
送鮮花一朵 | 為什么Gaussian09計(jì)算旋軌耦合效應(yīng)一定要先進(jìn)行標(biāo)量相對(duì)論計(jì)算，我用Molcas的計(jì)算完全不采用RCC基組，就是ANO-L基組，不過原子都很小，第二周期的，這樣誤差會(huì)有多大呢。本身HF的精度就不高。molcas的casscf級(jí)別計(jì)算還是不錯(cuò)的，怪不得我在caspt2基礎(chǔ)上計(jì)算旋軌耦合出來的數(shù)值和不包含caspt2模塊的一樣，不知道是本身誤差就不大，還是沒有開發(fā)到caspt2的層次。說明書上說好像目前只能計(jì)算casscf 水平的。另外有限尺寸原子核模型和點(diǎn)原子核模型是怎么回事呢。哪個(gè)更接近實(shí)際情況，可能原子太小沒有影響，只會(huì)造成一點(diǎn)點(diǎn)數(shù)據(jù)上的誤差。 |
7樓2012-10-07 11:50:35
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[論文投稿]
chinese chemical letters英文版投稿求助
 130+4
|
Yishengeryi 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +7 | 雙馬尾痞老板2 2026-03-31 | 7/350 |
|
|
|
[考研] 材料調(diào)劑 +9 | Eujd1 2026-03-31 | 10/500 |
|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +12 | 小黃鴨寶 2026-03-30 | 12/600 |
|
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑，總分281 +4 | 橙子（勝意） 2026-03-30 | 4/200 |
|
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來財(cái)aa 2026-03-25 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|