24小時熱門版塊排行榜

返回列表

feynman1965

金蟲 (正式寫手)

應(yīng)助: 12 (小學生)
金幣: 655.6
散金: 109
紅花: 7
帖子: 374
在線: 391.8小時
蟲號: 1424918
注冊: 2011-10-01
專業(yè): 光學

[求助] 關(guān)于計算改變核間距時兩原子分子HOMO 已有1人參與

各位蟲友好，
我是做實驗的，現(xiàn)在想用高斯做一個這樣的計算，兩原子分子比如H2的某一個排斥態(tài)的HOMO。
這還沒完，分子從基態(tài)躍遷到這個排斥態(tài)肯定核間距不斷增加，直到解離。
怎么計算在這個解離過程中，不同核間距的HOMO?
請大家給個方向性和個別關(guān)鍵細節(jié)指導，詳細的我可以自己學習。
謝謝大家。

回復此樓

» 猜你喜歡

[電子教材]天津大學物理化學（第七版）上、下冊已經(jīng)有3人回復
陰離子交換膜電解二氧化碳還原已經(jīng)有0人回復
物理化學論文潤色/翻譯怎么收費? 已經(jīng)有204人回復
Chemical Bonding at Surfaces and Interfaces,最經(jīng)典的一本表面上化學相互作用教材已經(jīng)有0人回復
九江學院2026年最新高層次人才招聘公告已經(jīng)有0人回復
深圳信息職業(yè)技術(shù)學院2026年教師招聘公告（最新）已經(jīng)有0人回復
誠邀有志之士加入江蘇大學環(huán)境學院合成生物學課題組！(長期接收學碩專碩、碩博連讀) 已經(jīng)有16人回復
研究生招生已經(jīng)有0人回復
一志愿南大理論與計算化學方向（070300）342分，求調(diào)劑已經(jīng)有1人回復
歡迎化學、化工、物理、材料、或計算機背景的考研咨詢已經(jīng)有0人回復
山西大學資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生已經(jīng)有2人回復

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2014-06-22 11:32:11

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

雪狼乖乖

木蟲 (正式寫手)

應(yīng)助: 18 (小學生)
貴賓: 0.005
金幣: 2805.8
散金: 567
紅花: 34
帖子: 901
在線: 402.2小時
蟲號: 532524
注冊: 2008-03-25
性別: GG
專業(yè): 理論和計算化學

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
feynman1965: 金幣+20, ★★★很有幫助 2014-06-23 11:37:35

單重態(tài)H2拉伸過程的電子密度差等值線動畫  轉(zhuǎn)載自Sobereva
http://hi.baidu.com/sobereva/item/c69893149c743a6e70d5e8c3
做電子密度差圖必須先獲得涉及到的各個元素原子的波函數(shù)文件。每次執(zhí)行Multiwfn時程序會從當前目錄下atomwfn目錄里尋找是否已有原子的波函數(shù)文件，沒有則需要調(diào)用Gaussian重算。為了方便、省時，這里先事先算好氫原子的波函數(shù)文件，放在atomwfn目錄下，并命名為H .wfn，注意空格必須有，否則Multiwfn不認�；M建議與計算分子時一致（若基組級別已較高，不一致時問題也不大）。其實計算密度差圖還有個既重要又麻煩的球?qū)ΨQ化原子密度問題，由于氫的密度本來就是球?qū)ΨQ的，這里就不提了。
執(zhí)行壓縮包內(nèi)的H2sin.gjf，這是單重態(tài)H2的scan任務(wù)，從0.2埃開始24步，每步0.1埃。注意盡管是單重態(tài)，但是拉遠后對稱破缺態(tài)是基態(tài)，所以要用非限制性方法，且需要加上guess(always,mix)打破對稱性。計算完畢后將儲存波函數(shù)文件的目錄H2sin拷到Multiwfn所在目錄下。
編寫一個文本文件batchcontourfix.txt放在Multiwfn目錄下，內(nèi)容是：
4
-2
1
2
200,200
6  //自行輸入格點設(shè)定
0.00000  -4.50000  -6.95664  //原點
0.00000 0.04500 0.00000  //平移向量1
0.00000 0.00000 0.06957  //平移向量2
0
[空行]
[空行]
其中-2代表要做變形屬性，由于函數(shù)選了1，計算的就是變形密度（密度差）。第一節(jié)的例子中坐標軸尺寸和刻度在不斷變動，這是為了使當前幀的結(jié)構(gòu)能充滿圖像，然而有些人不喜歡這樣，因此此例中自行設(shè)定原點和平移向量，坐標軸就固定住了。顯然很難直觀判斷數(shù)值應(yīng)該設(shè)多少，我建議用體系涉及空間尺度最廣的一幀作參考獲取這些信息。比如此例H2_0024.wfn對應(yīng)拉得最遠的一幀，就先手動對這個波函數(shù)的任意一個函數(shù)做圖（X=0的YZ平面，延展距離用默認的4.5 bohr），Multiwfn會輸出作圖時所用的格點信息：
Origin is 0.00000  -4.50000  -6.95664 bohr
The x/y/z of transition vector 1 is 0.00000 0.04500 0.00000 bohr
The x/y/z of transition vector 2 is 0.00000 0.00000 0.06957 bohr
令所有幀作圖時都使用這個格點設(shè)定，就能保證每一幀感興趣的區(qū)域都不出界。
然后，將前面batchrun.bat里的H2tri替換為H2sin，batchcolor.txt替換為batchcontourfix.txt，執(zhí)行之，就得到了每幀的.PNG圖像。將圖像復制到/sob/H2sinpic目錄下，并進入/sob目錄，執(zhí)行:
convert -delay 20 -resize 500 -shave '100x0' -colors 128 -monitor -reverse -fill red -annotate +5+395 "Created by Sobereva" H2sinpic/DISLIN.PNG H2sinpic/DISLIN_?.PNG H2sinpic/DISLI_??.PNG H2sinpic/H2sin_all.gif
和第一節(jié)的命令沒什么區(qū)別，只是目錄換了，并且加入了一個新的設(shè)定-shave '100x0'，這是因為這次的圖像的坐標軸橫軸固定了，左右的白邊長度也固定了，將每幀在resize為寬度為500后，左右都截掉100像素可以讓動畫的畫布更緊湊。最終所得動畫H2sin_all.gif如圖所示：

隨距離拉近后，電子密度迅速在原子間聚集，表明共價鍵正在加強。聚集的區(qū)域一開始逐漸變大并逐漸變得扁長，但后來卻開始收縮并變得接近于球型，最后一幀電子密度減少區(qū)域和增加區(qū)域都幾乎變成了球?qū)ΨQ狀，前者包著后者。這是因為原子序數(shù)越大，核吸引勢越強，電子密度分布范圍收縮得越厲害，基態(tài)He的1s電子密度延展范圍小于基態(tài)H的，所以當兩個氫核距離接近到氦核特征明顯出現(xiàn)時內(nèi)部電子密度就增加了，相應(yīng)地外側(cè)減少。

e824b899a9014c08b18022600a7b02087af4f4b8.jpg.gif

贊一下(1人)