計算化學-第一-二章

1、第一章 計算模型兩者的共同點是：1. 計算分子的能量計算分子的能量，分子的性質可以根據能量按照一定的方法得到。2. 進行幾何優(yōu)化進行幾何優(yōu)化，在起始結構的附近尋找具有最低的能量的結構。幾何優(yōu)化是根據能量的一階導數進行的。3. 計算分子內運動的頻率計算分子內運動的頻率。計算依據是能量的二階導數。計算化學的方法主要有計算化學的方法主要有分子力學理論分子力學理論(Molecular Mechanics)和和電子結構理論電子結構理論(Electronic Structure Theory)。1.1 計算化學概述計算化學概述分子力學理論分子力學理論分子理論分子理論采用經典物理對分子進行處理采用經典物理對

2、分子進行處理，可以在，可以在MM3，HyperChem，Quanta，Sybyl，Alchemy 等軟件中看等軟件中看到。根據所采用的力場的不同，分子理論又分為很多種。到。根據所采用的力場的不同，分子理論又分為很多種。分子力學理論方法很便宜分子力學理論方法很便宜(做量化的經常用貴和便宜來描做量化的經常用貴和便宜來描述計算，實際上就是計算時間的長短，因為對于要花錢述計算，實際上就是計算時間的長短，因為對于要花錢上機的而言，時間就是金錢；對于自己有機器的，要想上機的而言，時間就是金錢；對于自己有機器的，要想算的快，也要多在機器上花錢算的快，也要多在機器上花錢)，可以計算多達幾千個原，可以計算多達幾

3、千個原子的體系。其缺點是子的體系。其缺點是:1. 每一系列參數都是針對特定原子得出的。沒有對于原每一系列參數都是針對特定原子得出的。沒有對于原子各個狀態(tài)的統一參數。子各個狀態(tài)的統一參數。2. 計算中忽略了電子，只考慮鍵和原子，自然就不能處計算中忽略了電子，只考慮鍵和原子，自然就不能處理有很強電子效應的體系，比如不能描述鍵的斷裂。理有很強電子效應的體系，比如不能描述鍵的斷裂。電子結構理論電子結構理論這一理論基于薛定鄂方程薛定鄂方程，采用量子化學方法對分子進行處理。主要有兩類：1. 半經驗方法半經驗方法 包括AM1，MINDO/3，PM3，常見的軟件包有MOPAC，AMPAC，HyperChem，

4、以及Gaussian。半經驗方法采用了一半經驗方法采用了一些實驗得來的參數些實驗得來的參數，來幫助對薛定鄂方程的求解。2. 從頭算從頭算 從頭算，在解薛定鄂方程的過程中，只采用了幾個物理常數，包括光速，電子和核的質量，普朗克常數，在求解薛定鄂方程的過程中采用一系列的數學近似數學近似，不同的近似也就導致了不同的方法。最經典的是Hartree-Fock 方法，縮寫為HF。從頭算能夠在很廣泛的領域提供比較精確的信息，當然計算量要比前面講的方法大的多，就是貴得多了。密度泛函密度泛函(Density Functional Methods)密度泛函是最近幾年興起的第三類電子結構理論第三類電子結構理論方法方

5、法。它采用泛函采用泛函(以函數為變量的函數以函數為變量的函數)對薛定對薛定鄂方程進行求解鄂方程進行求解，由于密度泛函包涵了電子相關，它的計算結果要比HF 方法好，計算速度也快。1.2 化學模型化學模型(Model Chemistries)Gaussian 認為所謂理論是，一個理論模型，必須適用一個理論模型，必須適用于任何種類和大小體系，它的應用限制只應該來自于計于任何種類和大小體系，它的應用限制只應該來自于計算條件的限制。算條件的限制。這里包括兩點：1. 一個理論模型應該對于任何給定的核和電子有唯一的定義，就是說，對于解薛定鄂方程來講，分子結構本身就可以提供充分的信息。2. 一個理論模型是沒有

6、偏見的，指不依靠于任何的化學結構和化學過程。這樣的理論可以被認為是化學理論模型(theoretical-model chemistry)，簡稱化學模型(modelchemistry)Gaussian 包含多種化學模型(計算方法)，比如Gaussian 關鍵詞關鍵詞 方法HF Hartree-Fock 自恰場模型B3LYP Becke型3 參數密度泛函模型，采用Lee-Yang-Parr 泛函MP2 二級Moller-Plesset 微擾理論MP4 四級Moller-Plesset 微擾理論QCISD(T) 二次CI基組基組基組是分子軌道的數學表達，具體見第七章開殼層，閉殼層開殼層，閉殼層指電子

7、的自旋狀態(tài)，對于閉殼層閉殼層，采用限制性計算方法，在方法關鍵詞前面加R；對于開殼層開殼層，采用非限制性計算方法，在方法關鍵詞前面加U。比如開殼層的HF 就是UHF。對于不加的，程序默認為是閉殼層。一般采用開殼層的可能性是1. 存在奇數個電子，如自由基，一些離子2. 激發(fā)態(tài)3. 有多個單電子的體系4. 描述鍵的分裂過程模型的組合模型的組合高精度的計算往往要幾種模型進行組合，比如用中等算法進行結構優(yōu)化，然后用高精度算法計算能量?；瘜W模型的應用化學模型的應用1. 模型的建立模型的建立2. 分子的穩(wěn)定構形分子的穩(wěn)定構形3. 分子的各種譜圖分子的各種譜圖4. 化學反應機理化學反應機理5. 分子中的化學鍵

8、分子中的化學鍵A.化學模型的可計算性化學模型的可計算性最大原子數可計算量分子力學1000-100萬粗略的幾何結構半經驗200-1000幾何結構(有機分子)HF(DFT)50-200能量(含過渡金屬)MP225-50能量(弱, 氫鍵)CCSD(T)8-12精確能量(弱作用)CASPT210磁性(多個多重度)模型建立的步驟模型建立的步驟1.找出活性中心: 配位中心, 反應中心2.劃定最感興趣的區(qū)域: 主要配體及官能團3.根據化學知識進行剪切斷鍵,用類似原子或官能團補全: 用氫原子代替甲基等4.對削減太多的模型, 至少建立兩個層級的模型, 互相比較結果以確定能反映本質的最小模型5.利用最小模型進行詳

9、細研究, 解釋實驗結果和現象: 結構穩(wěn)定性, 最穩(wěn)構型等 幾何結構來自同一體系能量差值， 容易比較， 因此幾何結構的優(yōu)化一般都比較準確B. 分子的穩(wěn)定構形分子的穩(wěn)定構形 分子的穩(wěn)定構形分子的穩(wěn)定構形 理論計算的第一步一般是先優(yōu)化結構， 得到理論上本模型的最穩(wěn)定結構！ 確定分子的穩(wěn)定結構就是找到分子的極小能量點, 比較不同構形的能量差，確定最穩(wěn)定構形。結構決定性能！極小點最小點計算穩(wěn)定構形的熱化學數據計算穩(wěn)定構形的熱化學數據 體系總能量(電子能量+核排斥能) 零點振動能(從振動頻率計算出) 從0K到298K的焓, 熵(按理想氣體計算)C. 分子的各種譜圖分子的各種譜圖 在得到分子的穩(wěn)定結構后,

10、可以從其波函數得到相關的譜圖, 并與實驗進行對比 分子簡正振動頻率(紅外和拉曼) (Freq) 核磁共振屏蔽張量和磁化系數 (NMR) 振動圓二色性(Freq=VCD) 旋光性(Polar=OptRot)用用NICS度量環(huán)的芳香性度量環(huán)的芳香性 NICS（Nuclear-Independent Chemical Shift）與核無關的化學位移 芳香性是化學分子穩(wěn)定的一個重要方面 芳香性來源于環(huán)上的大p鍵 芳香性可以用位于環(huán)中心處虛 原子的化學位移上（NMR) 可以解釋很多環(huán)化合物的穩(wěn)定性差別D. 化學反應機理化學反應機理 化學反應機理的研究主要在于確定其單步反應的過渡態(tài), 即勢能面上的鞍點。反

11、應物產物與確定分子的最穩(wěn)定構型不同, 要研究反應機理, 必須找到勢能面上的鞍點, 在此鞍點處, 振動頻率有且只有一個頻率是虛數.計算過渡態(tài)注意事項計算過渡態(tài)注意事項 尋找過渡態(tài)需要的時間要比尋找分子穩(wěn)定結構要多很多 過渡態(tài)必須進行頻率驗證，確保振動頻率中只有一個簡正振動頻率為虛數。E. 分子中的化學鍵分子中的化學鍵 從化學家的角度解釋各種分子體系的現象, 比如鍵型, 鍵的強弱, 電荷分布, 反應活性位點等化學鍵性質的主要方面化學鍵性質的主要方面分子是由原子通過化學鍵結合而成的量子化學語言: LCAOMO LCAO 各原子上的電荷：Mulliken布居分析， Mayer鍵級, 自然布居分析(NB

12、O) HOMO和LUMO成份前線軌道（前線軌道（HOMO, LUMO) 前線軌道理論確定反應活性點前線軌道理論確定反應活性點 HOMO成份多的原子是親核活性點成份多的原子是親核活性點 LUMO成份多的原子是親電活性點成份多的原子是親電活性點第二章 GaussView 使用簡介一一 簡介簡介GaussView 是一個專門設計于高斯配套使用的軟件，其主要用途有兩個1 構建高斯的輸入文件構建高斯的輸入文件2 以圖的形式顯示高斯計算的結果以圖的形式顯示高斯計算的結果除了可以自己構建輸入文件外，GVIEW 還可讀入CHEM3D，HYPERCHEM 和晶體數據等諸多格式的文件。從而使其可以于諸多圖形軟件連

13、用，大大拓寬了使用范圍（詳見下圖）GVIEW 的界面及主要功能鍵的介紹界面簡介第一行為菜單欄:File:主要功能是建立，打開，保存和打印當前的文件Save Inage 將當前文件保存為圖片格式Preferences?？梢栽诶锩娓淖僄VIEW 默認的各種顯示的設置。Edit: 在這里可以完成對分子的剪且，拷貝，刪除 和抓圖等。Atom List,顯示當前分子的內坐標，笛卡兒坐標，分數坐標等。Point Group 可以顯示當前分子的點群及可能有的點群。PBC 顯示晶體文件（可以將CIF 文件轉換為圖形，在點PBC 按鈕后所給并的對話框中根據選項調節(jié)具體顯示的格式。Mos 用于顯示分子軌道（只有檢

14、查點文件，此選項才能給出分子軌道圖。Symmetrize，對當前體系進行對稱性控制。View 這里面的選項都是于分子的顯示有關的，如顯示氫原子，顯示鍵，顯示元素符號，顯示啞原子顯示坐標軸等Calculalte：可從GVIEW 中直接向高斯提交計算。這是GVIEW 作為高斯軟件配套功能的重要體現。從所給的對話框中可以選擇工作類型Job Type（如優(yōu)化，能量或頻率等）;計算方法Method（如半經驗方法，HF 方法，DFT 方法，MP 方法等，還可以選定基組）;Title(對所要做的計算給一個說明，以備以后的查看) Link 0（給檢查點文件命名，還可以在此用RWF 命令設置臨時數據交換文件的大

15、小）； General, Guess,（這兩個選項主要是給出體系中各原子的連接關系及如何給出初始猜測）；NBO（可在此設定NBO 計算）,PBC（可在此設定晶體的有關計算）,Solvation（可在此設定溶液中的計算，除了選擇溶劑外，還要選擇模擬溶劑的理論模型）Result ：顯示計算的結果，包括電荷，靜電勢表示的表面，振動頻率，核磁，勢能面掃描，優(yōu)化等。注意有些結果只能用檢查點文件才能顯示。構建一個間氟苯乙烷分子并從GVIEW 里遞交計算(1)雙擊GVIEW 圖標，打開GVIEW ，下圖就是GVIVE 打開后的窗口(2)雙擊窗口中的 圖標，得到如下窗口里面有常用的環(huán)狀官能團。選中苯環(huán)（單擊即

16、可選中）(3).在當前工作窗口（打開GIVEW 時程序自動打開一個工作窗口，如下圖）也可通過File-new 路徑新建一個工作窗口在這個窗口中點鼠標左鍵窗口中就會出現苯分子，見下圖將鼠標放在分子上，按左鍵左右或前后移動，可以調節(jié)分子的角度將鼠標放在分子上，前后移動，可以將分子放大或縮小Shift+Alt+鼠標左鍵組合可以在窗口內平移分子當工作窗口內有多個分子時【在構建大的分子時，這種情況很容易出現】?？梢杂肅trl+Alt+鼠標左鍵組合調節(jié)其中一個分子的角度，以調節(jié)各個分子間的角度。(4)雙擊GVIEW界面上的 圖標。出現窗口:點擊氟的元素符號“F”，就選中氟原子。(5).回到工作窗口在苯環(huán)上

17、的任一個H上單擊左鍵，將H置換成F。但要注意此時僅是元素符號發(fā)生了改變，C-F間的距離仍是C-H鍵的距離。需要用鍵長工具進行調整。單擊GVIEW界面 圖標。 然后再點擊工作窗口中的C原子和F原子，看到被選中的兩個原子于周圍的原子再亮度上有差異 此時會出現右邊的窗口根據C-F鍵的長度在 0.675和2.700之間的方框內進行C-F鍵的調整。完畢后點擊OK即可。(6)雙擊GVIEW界面上的 圖標，出現以下窗口這是GVIEW里內置的鏈烴庫，選中乙烷(7).在工作窗口內空處點左鍵用我們前面講過的命令調節(jié)苯和乙烷間的距離和角度將鼠標放在分子上，按左鍵左右或前后移動，可以調節(jié)分子的角度；將鼠標放在分子上，

18、移動滾輪，可以將分子放大或縮小。Shift+Alt+鼠標左鍵組合可以在窗口內平移分子當工作窗口內有多個分子時【在構建大的分子時，這種情況很容易出現】這時可用以下命令可以用Shift+Alt+鼠標左鍵組合移動想要移動的分子，以調節(jié)各個分子間的距離。可以用Ctrl+Alt+鼠標左鍵組合調節(jié)其中一個分子的角度，以調節(jié)各個分子間的角度Ctrl+Alt+鼠標左鍵這個組合常和將鼠標左鍵放在分子上，左右或前后移動，可以調節(jié)分子的角度這個功能連用。具體的好處需要自己去試。（8）單擊GVIEW界面上 圖標。然后點擊與乙烷正對的苯上的C原子和H原子選擇“None”，乙烷上的C-H鍵作類似處理然后單擊GVEW上的圖

19、標，點選工作窗口中的兩個氫原子，即將它們刪去再單擊GVEW界面上 圖標，根據實際鍵長將兩個C原子聯起來。（9）查看分子的對稱性分子的對稱性從Edit-Point groupEdit-Point group路徑可以查看所構建分子的點群。點擊Point groupPoint group后，出現如下窗口：為C1點群，其下拉菜單中的為可能的點群（改變Tolerance,也可幫助我們判斷所構建體系可能有的點群）(10)查看分子坐標。單擊GVIEW界面上的 圖標。出現下面窗口 Z表示時內坐標，C表示直角坐標?？梢栽诶锩鎸ψ鴺俗鲞m當的調整（11）向GAUSS遞交計算。點GVIEW界面上Calculation

20、 中的Gaussian會出來一個遞交計算的對話框。從所給的對話框中可以選擇工作類型Job Type（如優(yōu)化，能量或頻率等）;計算方法Method（如半經驗方法，HF方法，DFT方法，MP方法等，還可以選定基組;Title(對所要做的計算給一個說明，以備以后的查看) Link 0（給檢查點文件命名，還可以在此用RWF命令設置臨時數據交換文件的大?。籊eneral,Guess,（這兩個選項主要是給出體系中各原子的連接關系及如何給出初始猜測）；NBO（可在此設定NBO計算）,PBC（可在此設定晶體的有關計算）, Solvation（可在此設定溶液中的計算，除了選擇溶劑外，還要選擇模擬溶劑的理論模型

21、）選擇完畢后，點Submit即可遞交計算。有時由于安裝的原因GVIEW無法與GAUSS建立關聯，就不能直接從GVIEW里遞交計算。這時可以在GVIEW里保存用于GAUSS計算的輸入文件，然后從GAUSS里調出文件進行計算。3 3構建分子中的注意事項構建分子中的注意事項在構建分子前，應對分子的鍵長，鍵角和點群等有一個詳細的在構建分子前，應對分子的鍵長，鍵角和點群等有一個詳細的了解了解。這樣在構建過程中才能更好的把握細節(jié)，如構建大分子時，各個集團間的鍵長，鍵角和二面角一定要給的準確。以我的經驗而言，相當一部分的錯誤（如SCF不收斂，1502報錯等）都是分子建模不合理造成的。分子模型不合理導致程序根

22、據分子模型所得的薛定鄂方程不合理，也就導致了對方程的解不能順利完成，所以報錯。對于大分子，不會出現分子直線排步的情況，鏈烴通常為鋸齒型的排列，環(huán)之間總是有一定的夾角，共面情況雖有但不多見，這些都是在構建分子模型的過程中需要注意的。再者，分子建模是量化計算的一個難點，特別是對金屬配合物，塊狀體系和納米體系的構建。只有給出合理的結構，程序才能給出合理的解。我是做有機體系的計算，對以上三種不熟，還請各位同學見諒。如果不能體現出分子的對稱性，則程序不能正確判斷分子軌道的對稱性，這在涉計到反應機理和軌道相互作用的研究中可能會導致錯誤的判斷。練習練習練習2.1 利用GVIEW構建一個二茂鐵分子練習2.2

23、利用GVIEW構建一個2-苯基吲哚分子練習2.3 利用GVIEW構建一個三苯胺分子練習2.4 利用GVIEW構建一個醋酸鈀分子第三章 單點能計算單點能計算是指對單點能計算是指對給定幾何構性的分子的能量以給定幾何構性的分子的能量以及性質進行計算及性質進行計算，由于分子的幾何構型是固定不由于分子的幾何構型是固定不變的，只是變的，只是一個點一個點，所以叫單點能計算，所以叫單點能計算。單點。單點能計算可以用于：能計算可以用于：1. 計算分子的基本信息計算分子的基本信息2. 可以作為分子構型優(yōu)化前對分子的檢查可以作為分子構型優(yōu)化前對分子的檢查3. 在由較低等級計算得到的優(yōu)化結果上進行高精在由較低等級計算得到的優(yōu)化結果上進行高精度的計算度的計算4. 在計算條件下，體系只能進行單點計算在計算條件下，體系只能進行單點計算5. 單點能的計算可以在不同理論等級，采用不同單點能的計算可以在不同理論等級，采用不同基組進行基組進行例例 甲醛的單點能計算甲醛的單點能計算輸出文件中的信息輸出文件中的信息能量能量找到SCF Done：E(RHF)= -113.863697598 A.U. after 6 cycles這里的數值就是能量，單位是hartree。在一些高等級計算中，往往有不止一個能