如何用Gaussian掃描探討鍵角、鍵長與鍵能的關係

作者：
我們都學過CO2是直鏈型分子，C以sp軌域和O形成鍵結，剩下兩對LP。
那我們要怎麼確認CO2在直鏈排列下或是水在真空下以109.5度展開時，能量是最低的呢？

其實Gaussian具有掃描的功能。以下是路谷兮自己研究出的教學。

  • 該如何掃描鍵角呢？
1.    創建分子，此處以CO2為範例。
2.    此處我使用STO-3G估算，以加速計算。先將物質結構OPT，確保鍵長、鍵角是在最低能的位置。
3.    開啟log檔，確認已優化完成

4.    調整鍵角，選取鍵角調整工具
最後點OK

5.    回到GaussView視窗，選擇Redundant coordinate editor，上面有R的圖案。

6.    點擊Add，添加新的參數。然後選擇Angle，輸入指定的鍵角，Gaussian讀取鍵角方式aussian讀取鍵角方式，同數學讀取方式，即中間那位數字是中間的原子，必須要看GaussView上面的Label來輸入。

輸入要掃描幾次，每次增加幾度。然後點選OK
(圖中起始角度應為160度)

6.    ctrl+G，開啟計算編輯，此處選擇Scan

5.    選擇Redundant coord.

7.    Submit完，靜待計算，之後選擇顯示log檔，就可以看到囉！
點選左上角的綠點，可以播放動畫。

#鍵角為180度時，G09座標會有問題，所以若要算180度的能量，請先直接用OPT完的能量。

  • 那該如何掃描鍵長與能量的關係呢？
          很簡單只要改變調整鍵長，輸入指定的原子就好了！
接下來就剩等待電腦計算囉！



若是您正在研究量子化學或是有興趣的話，歡迎繼續點閱後續的文章。
敬請不吝指教！！


參考：
1.Gaussian官網：
2.成大化學系鄭沐政老師-分子與材料模擬技術課程

留言

張貼留言

這個網誌中的熱門文章

如何使用ImageJ分析金相圖相分率

作者：
圖片
材料學中，判斷金相圖是所有材料人必經的過程。但是，一樣的晶相，每個人說的結論都不太相同，所以同樣一張照片，用電腦判斷金相圖有其必要性。 本次，路谷兮將以ImageJ作為操作軟體，並以錫鉛合金作為範例，教你如何判斷相分率。 第一步：開啟ImageJ，並載入金相圖 第二步：調整圖片成8-bit 第三步：點擊Threshold 此時，畫面會長的像這樣，若是前面未調整金相圖Type，將會分析不了照片強度。 此處可以調整要選取的值上下限、顯示顏色。 此處以預設為例。 第四步：點擊分析粒子數 此時會跳出此視窗，此處已預設為例。 第五步：分析完成，Summary中的%Area就是相分率 因為此處以錫鉛合金為例，觀察原金相圖，可以得知原本白色，現在是紅色被選取的部分是Beta相。雖然說本次分析的是截面下的相分率，但是我們假設整個材料是均勻的，所以可以將此值當作體積分率。若是要精準量測相分率，請用其他軟體測量。 若是您正在對ImageJ有興趣的話，歡迎繼續點閱後續的文章。 敬請不吝指教！！ 感謝成大材料系材料實驗課-顏紹宇助教的指導。
閱讀完整內容

如何用Gaussian看波函數

作者：
圖片
Gaussian可以拿來觀看分子的波函數，可以選擇要看的是LUMO，HOMO，或是總電子雲密度。 1.     創建分子團 首先，用GaussView創建新的分子團，角度跟方向先不用刻意處理，只要分子相對位置及鍵結差不多即可！ 此處以乙酸為例 2.     ctrl+s存檔，並以記事本打開此檔案，即可檢視 。 圖中會顯示chk檔位置，使用方法，電荷，多重度，座標，鍵結級數 3.     回到GaussView視窗，點擊Gaussian calculation setup(快捷ctrl+G) 4.      選取能量，然後點選submit，後面參數先不用調整 5.     等待G09計算，直到執行完成，即可先關閉G09， 此視窗請點選.chk，不然看不到波函數 計算完成的分子，乍看下是一樣的，但是上面附檔名會是.chk 6.     回到GaussView視窗，點選Result>Surfaces/Contours 在此畫面中，選取New Cube 這裡可以選取顯示的有，MO、電子雲密度 在MO中，可以選取LUMO、HOMO，或是一次選取 生成後，選取Cube，再點選New Surface，就可以看到選定的MO了 乙酸HOMO 乙酸LUMO 要看總電子密度也是可以的 選取Mapped Surface可以看到密度分布 選取只顯示表層，然後選ESP 表層電子雲密度分布就出現了 選取Hide，可以隱藏掉不要顯示的部分喔！ 回到GaussView，選取Save image file或是Windows+S，就能保存當前2D畫面囉！ 那要怎麼保存3D圖呢？那請把chk檔給想要看的人吧！ 注意！此分子的波函數不會是現實中分子的波函數，因為他可能結構有異，或是因為在現實中他是在水溶液環境中，(Gaussian預設為在真空中氣體分子的互動)。 那該如何調整或優化結構，或是調整環境呢？請待後續文章揭曉。 若是您正在研究量子化學或是有興趣的話，歡迎繼續點閱後續的文章。 敬請不吝指教！！ 參考： 1.Gaussian官網： http://gaussian.com/ 2.成大化學系 鄭沐政老師 - 分子與材料模擬技術課程
閱讀完整內容

ICSD基礎教學

作者：
圖片
Gaussian能計算晶體結構，但是你必須要先有晶體結構參數才能計算。 此處，路谷兮利用ICSD的晶體結構數據做分析，ICSD全名是無機結晶結構資料庫(Inorganic Crystal Structure Database)，成功大學有購買權限，只要你是用學校IP位置讀取此網站，那就能無限自由搜索，且不會過期。但若是你們學校沒有購買權限，那你仍然可以申請一個30天的DEMO版本，以下是創立教學： 1.   ICSD創立教學：  到 ICSD官網 創立DEMO帳號 此處選擇Get Demo Account 輸入個人資料，輸入完，回到信箱 用下方的Username跟Password登入就可以囉！ 2.     ICSD操作教學： ICSD可以列很多搜尋條件，我們可以利用結晶學的知識，如Space group、元素、cell parameters等等，代入參數搜尋，在 此處我們想查詢C的同素異形體(Allotropes)，我們可以在composition中輸入C，然後在後方不同元素個數輸入1，即可搜尋到所有碳的同素異形體。 下方的Count basic search可計算符合條件的資料數量 右方會顯示數據 點擊Run Query進入數據庫 因為我們使用Demo版本，所以只會顯示graphite跟Diamond的數據，不過這樣也夠了。 可以點擊進入，看更詳細的數據 選擇Export cif，將參數儲存以利後續分析。 ICSD輸出的cif檔，可以利用Gaussian或是VESTA開啟，那你想要知道分析教學嗎？敬啟鎖定後續文章。 參考： 1.成大化學系 鄭沐政老師 - 分子與材料模擬技術課程
閱讀完整內容