| 查看: 2543 | 回復(fù): 12 | ||
| 本帖產(chǎn)生 1 個(gè) QC強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | ||
feynman1965金蟲 (正式寫手)
|
[求助]
關(guān)于得到分子波函數(shù)的問題
|
|
|
各位好,請(qǐng)問怎么用高斯得到分子的波函數(shù)的表達(dá)式。其中波函數(shù)用球諧函數(shù)展開(如圖)。大家也可以給出意見或想法,隨便什么,只要有助于問題的解釋!謝謝! 附:高斯中wfn文件怎么看呀!根據(jù)輸出怎么寫出解析的表達(dá)式! wavefunction.png |
木蟲 (著名寫手)
|
高斯中的基函數(shù)可以有直角坐標(biāo)和球諧函數(shù)兩種形式,選擇基組時(shí)可以用5D和7F選項(xiàng)來強(qiáng)制選擇球諧函數(shù),而且大部分基組缺省都是球諧函數(shù)。 你在輸出文件中會(huì)看到P0,D+2,F-1等等符號(hào)來表示基函數(shù),這些符號(hào)就直接表示球諧函數(shù)的l和m下標(biāo)的值。s,p,d,f,...表示l下標(biāo)的值分別是0,1,2,3,...而后面的數(shù)字直接是m的值。這和類輕原子的標(biāo)記方式是基本一致的。 原則上來說,都可能的話都應(yīng)使用球諧函數(shù)基組。直角坐標(biāo)的基組角動(dòng)量不純凈,因此高斯里"純"或pure的基組都是指球諧函數(shù)基組。例如6D就相當(dāng)于5D外加一個(gè)s函數(shù)(dxx+dyy+dzz是一個(gè)s函數(shù))。一般來說使用直角坐標(biāo)基組的唯一意義就是和其他計(jì)算進(jìn)行比對(duì)。 例如d軌道可以用6個(gè)直角坐標(biāo)函數(shù)(d_xx,yy,zz,xy,xz,yz)或者5個(gè)球諧函數(shù)(d+2,+1,0,-1,-2)來表示,所以分別6D和5D表示。而f則有10個(gè)直角坐標(biāo)或者7個(gè)球諧函數(shù)。 也就是說大部分情況下得到分子軌道直接就是球諧函數(shù)線性組合了,不需要任何轉(zhuǎn)化。 但是每個(gè)原子中心不同,分子軌道雖然是以球諧函數(shù)的線性組合表示的,但是這些球諧函數(shù)的中心不同。如果要寫成統(tǒng)一中心的球諧函數(shù),就要進(jìn)行數(shù)值積分得到每一個(gè)Y_lm的徑向函數(shù)。 |
木蟲 (小有名氣)
木蟲 (正式寫手)
| gaussian得到的波函數(shù)是分子的波函數(shù),其中只有S軌道成分的波函數(shù)才能用你所寫的函數(shù)表達(dá)式進(jìn)行書寫(因?yàn)榭梢苑蛛x變量),S軌道成分還必須參考軌道的能級(jí)分布圖才能選擇。gaussian得到的都是解析解,很難寫出分子波函數(shù)的表達(dá)式,甚至完全不可能。 |
金蟲 (正式寫手)
金蟲 (正式寫手)
木蟲 (小有名氣)
金蟲 (正式寫手)
木蟲 (著名寫手)
|
如果想看你的基組具體參數(shù)可以用GFPrint關(guān)鍵字。一個(gè)波函數(shù)涉及幾千幾萬甚至幾百萬個(gè)參數(shù),用眼睛直接看輸出文件不太容易看出端倪,最好還是用波函數(shù)分析軟件,高斯本身會(huì)輸出一些分析,proaim,multiwfn之類的都可以。 可以看一下這一篇,高斯的參考文獻(xiàn) http://www.chem.wayne.edu/schlegel/Pub_folder/174.pdf 重要的內(nèi)容基本上他都引用了 大部分基組都是為球諧函數(shù)設(shè)計(jì)的,一個(gè)特例是pople基組是為直角坐標(biāo)優(yōu)化的 |
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
這兩個(gè)帖子有詳細(xì)討論: 《球諧型與笛卡爾型Gauss函數(shù)的轉(zhuǎn)換關(guān)系》 http://hi.baidu.com/sobereva/item/af178d498db96cea1f19bcdf 《談?wù)?d,6d型d殼層波函數(shù)與它們?cè)诟咚怪械臉?biāo)識(shí)》http://hi.baidu.com/sobereva/item/ae43b6bcd392c8a5ebba93dc |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 362求調(diào)劑 +14 | 西南交材料專碩3 2026-03-31 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +17 | 爭(zhēng)取九點(diǎn)睡 2026-03-30 | 17/850 |
|
|
[基金申請(qǐng)]
請(qǐng)問共同通訊和共同一作的認(rèn)可度問題
10+4
|
psa1234 2026-04-01 | 8/400 |
|
|
[考研] 302求調(diào)劑一志愿華中師范大學(xué) +5 | 小江小江江江 2026-04-02 | 5/250 |
|
|
[考研] 309求調(diào)劑 +8 | 呆菇不是戴夫 2026-04-02 | 8/400 |
|
|
[考研] 342求調(diào)劑 +13 | Mary Keen 2026-03-28 | 14/700 |
|
|
[考研] 262求調(diào)劑 +5 | 勵(lì)志一定發(fā)文章 2026-04-02 | 6/300 |
|
|
[考研] 求調(diào)劑推薦 +3 | 南山南@ 2026-04-01 | 3/150 |
|
|
[考研] 372求調(diào)劑 +3 | jj涌77 2026-04-02 | 3/150 |
|
|
[考研] 一志愿北京科技,085601總分305求調(diào)劑 +9 | 半生瓜! 2026-04-01 | 11/550 |
|
|
[考研] 生物與醫(yī)藥考研調(diào)劑 +5 | 鐵憨憨123425 2026-03-31 | 5/250 |
|
|
[考研] 08工科275求調(diào)劑,可跨考。 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +12 | 26要上岸 2026-03-28 | 12/600 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +10 | @taotao 2026-03-31 | 11/550 |
|
|
[考研] 求調(diào)劑,一志愿北林食品與營(yíng)養(yǎng)095500,301分,已過六級(jí),有科研經(jīng)歷 +4 | 快樂儲(chǔ)蓄罐 2026-03-31 | 4/200 |
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|