如何用Gaussian看化學反應生成熱

作者：
        在化學反應中，我們學過生成物的能量扣反應物的能量，就是反應生成熱。(Product - Reactant = delta internal energy)
        在Gaussian中，我們可以輕易計算生成熱。以下是教學：
此處以CO&CO2反應作為範例：
CO + 1/2 O2 = CO2

1.    將CO和CO2分別作OPT，確保二者位於最穩定結構。氧氣不需要，因為氧氣生成熱為0，所以不影響計算。
打開LOG檔，確認結構無誤。

2.    找到檔案位置，然後以記事本等文書軟體打開檔案。找尋優化後，分子結構的能量。秘訣：ctrl+F，搜尋E(，不要點選Wrap around，然後Enter鍵按住，最後一個搜尋到的值就是優化下的能量。
由圖中可得，優化後，
CO能量為-111.225449517 hartree
CO2能量為-185.068390601 hartree
將以上二者相減就是反應熱
一個Hartree等於627.51 kcal/mol，加上習慣上kcal/mol只留小數點下一位估計值。
故，本次化學反應生成熱為-46377.2 kcal/mol

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


參考：
1.Gaussian官網：
2.成大化學系鄭沐政老師-分子與材料模擬技術課程
位置： 402台灣台中市南區興大路145號

留言

張貼留言

這個網誌中的熱門文章

如何使用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.成大化學系 鄭沐政老師 - 分子與材料模擬技術課程
閱讀完整內容