計算化學第一二章

計算化學第一二章第1頁，共54頁，2023年，2月20日，星期二第一章計算模型第2頁，共54頁，2023年，2月20日，星期二兩者的共同點是：1.

計算分子的能量，分子的性質可以根據能量按照一定的方法得到。2.

進行幾何優(yōu)化，在起始結構的附近尋找具有最低的能量的結構。幾何優(yōu)化是根據能量的一階導數(shù)進行的。3.

計算分子內運動的頻率。計算依據是能量的二階導數(shù)。計算化學的方法主要有分子力學理論(MolecularMechanics)和電子結構理論(ElectronicStructureTheory)。第3頁，共54頁，2023年，2月20日，星期二1.1計算化學概述分子力學理論分子理論采用經典物理對分子進行處理，可以在MM3，HyperChem，Quanta，Sybyl，Alchemy等軟件中看到。根據所采用的力場的不同，分子理論又分為很多種。分子力學理論方法很便宜(做量化的經常用貴和便宜來描述計算，實際上就是計算時間的長短，因為對于要花錢上機的而言，時間就是金錢；對于自己有機器的，要想算的快，也要多在機器上花錢)，可以計算多達幾千個原子的體系。其缺點是:1.每一系列參數(shù)都是針對特定原子得出的。沒有對于原子各個狀態(tài)的統(tǒng)一參數(shù)。2.計算中忽略了電子，只考慮鍵和原子，自然就不能處理有很強電子效應的體系，比如不能描述鍵的斷裂。第4頁，共54頁，2023年，2月20日，星期二電子結構理論這一理論基于薛定鄂方程，采用量子化學方法對分子進行處理。主要有兩類：1.半經驗方法包括AM1，MINDO/3，PM3，常見的軟件包有MOPAC，AMPAC，HyperChem，以及Gaussian。半經驗方法采用了一些實驗得來的參數(shù)，來幫助對薛定鄂方程的求解。2.從頭算從頭算，在解薛定鄂方程的過程中，只采用了幾個物理常數(shù)，包括光速，電子和核的質量，普朗克常數(shù)，在求解薛定鄂方程的過程中采用一系列的數(shù)學近似，不同的近似也就導致了不同的方法。最經典的是Hartree-Fock方法，縮寫為HF。從頭算能夠在很廣泛的領域提供比較精確的信息，當然計算量要比前面講的方法大的多，就是貴得多了。第5頁，共54頁，2023年，2月20日，星期二密度泛函(DensityFunctionalMethods)密度泛函是最近幾年興起的第三類電子結構理論方法。它采用泛函(以函數(shù)為變量的函數(shù))對薛定鄂方程進行求解，由于密度泛函包涵了電子相關，它的計算結果要比HF方法好，計算速度也快。第6頁，共54頁，2023年，2月20日，星期二1.2化學模型(ModelChemistries)Gaussian認為所謂理論是，一個理論模型，必須適用于任何種類和大小體系，它的應用限制只應該來自于計算條件的限制。這里包括兩點：1.一個理論模型應該對于任何給定的核和電子有唯一的定義，就是說，對于解薛定鄂方程來講，分子結構本身就可以提供充分的信息。2.一個理論模型是沒有偏見的，指不依靠于任何的化學結構和化學過程。這樣的理論可以被認為是化學理論模型(theoretical-modelchemistry)，簡稱化學模型(modelchemistry)第7頁，共54頁，2023年，2月20日，星期二Gaussian包含多種化學模型(計算方法)，比如Gaussian關鍵詞方法HFHartree-Fock自恰場模型B3LYPBecke型3參數(shù)密度泛函模型，采用Lee-Yang-Parr泛函MP2

二級Moller-Plesset微擾理論MP4

四級Moller-Plesset微擾理論QCISD(T)二次CI第8頁，共54頁，2023年，2月20日，星期二基組基組是分子軌道的數(shù)學表達，具體見第七章開殼層，閉殼層指電子的自旋狀態(tài)，對于閉殼層，采用限制性計算方法，在方法關鍵詞前面加R；對于開殼層，采用非限制性計算方法，在方法關鍵詞前面加U。比如開殼層的HF就是UHF。對于不加的，程序默認為是閉殼層。一般采用開殼層的可能性是1.存在奇數(shù)個電子，如自由基，一些離子2.激發(fā)態(tài)3.有多個單電子的體系4.描述鍵的分裂過程模型的組合高精度的計算往往要幾種模型進行組合，比如用中等算法進行結構優(yōu)化，然后用高精度算法計算能量。第9頁，共54頁，2023年，2月20日，星期二化學模型的應用1.模型的建立2.分子的穩(wěn)定構形3.分子的各種譜圖4.化學反應機理5.分子中的化學鍵第10頁，共54頁，2023年，2月20日，星期二A.化學模型的可計算性最大原子數(shù)可計算量分子力學1000-100萬粗略的幾何結構半經驗200-1000幾何結構(有機分子)HF(DFT)50-200能量(含過渡金屬)MP225-50能量(弱,氫鍵)CCSD(T)8-12精確能量(弱作用)CASPT2<10磁性(多個多重度)第11頁，共54頁，2023年，2月20日，星期二模型建立的步驟找出活性中心:配位中心,反應中心劃定最感興趣的區(qū)域:主要配體及官能團根據化學知識進行剪切斷鍵,用類似原子或官能團補全:用氫原子代替甲基等對削減太多的模型,至少建立兩個層級的模型,互相比較結果以確定能反映本質的最小模型利用最小模型進行詳細研究,解釋實驗結果和現(xiàn)象:結構穩(wěn)定性,最穩(wěn)構型等第12頁，共54頁，2023年，2月20日，星期二幾何結構來自同一體系能量差值，容易比較，因此幾何結構的優(yōu)化一般都比較準確B.分子的穩(wěn)定構形第13頁，共54頁，2023年，2月20日，星期二

分子的穩(wěn)定構形理論計算的第一步一般是先優(yōu)化結構，得到理論上本模型的最穩(wěn)定結構！

確定分子的穩(wěn)定結構就是找到分子的極小能量點,比較不同構形的能量差，確定最穩(wěn)定構形。結構決定性能！極小點最小點第14頁，共54頁，2023年，2月20日，星期二計算穩(wěn)定構形的熱化學數(shù)據體系總能量(電子能量+核排斥能)零點振動能(從振動頻率計算出)從0K到298K的焓,熵(按理想氣體計算)第15頁，共54頁，2023年，2月20日，星期二C.分子的各種譜圖在得到分子的穩(wěn)定結構后,可以從其波函數(shù)得到相關的譜圖,并與實驗進行對比分子簡正振動頻率(紅外和拉曼)(Freq)核磁共振屏蔽張量和磁化系數(shù)(NMR)振動圓二色性(Freq=VCD)旋光性(Polar=OptRot)第16頁，共54頁，2023年，2月20日，星期二用NICS度量環(huán)的芳香性NICS（Nuclear-IndependentChemicalShift）與核無關的化學位移芳香性是化學分子穩(wěn)定的一個重要方面芳香性來源于環(huán)上的大p鍵芳香性可以用位于環(huán)中心處虛原子的化學位移上（NMR)可以解釋很多環(huán)化合物的穩(wěn)定性差別第17頁，共54頁，2023年，2月20日，星期二D.化學反應機理化學反應機理的研究主要在于確定其單步反應的過渡態(tài),即勢能面上的鞍點。反應物產物與確定分子的最穩(wěn)定構型不同,要研究反應機理,必須找到勢能面上的鞍點,在此鞍點處,振動頻率有且只有一個頻率是虛數(shù).第18頁，共54頁，2023年，2月20日，星期二計算過渡態(tài)注意事項尋找過渡態(tài)需要的時間要比尋找分子穩(wěn)定結構要多很多過渡態(tài)必須進行頻率驗證，確保振動頻率中只有一個簡正振動頻率為虛數(shù)。第19頁，共54頁，2023年，2月20日，星期二E.分子中的化學鍵從化學家的角度解釋各種分子體系的現(xiàn)象,比如鍵型,鍵的強弱,電荷分布,反應活性位點等化學鍵性質的主要方面分子是由原子通過化學鍵結合而成的量子化學語言:LCAO→MO→LCAO各原子上的電荷：Mulliken布居分析，Mayer鍵級,自然布居分析(NBO)HOMO和LUMO成份第20頁，共54頁，2023年，2月20日，星期二前線軌道（HOMO,LUMO)前線軌道理論確定反應活性點HOMO成份多的原子是親核活性點LUMO成份多的原子是親電活性點第21頁，共54頁，2023年，2月20日，星期二第二章GaussView

使用簡介第22頁，共54頁，2023年，2月20日，星期二一簡介GaussView是一個專門設計于高斯配套使用的軟件，其主要用途有兩個1．構建高斯的輸入文件2．以圖的形式顯示高斯計算的結果第23頁，共54頁，2023年，2月20日，星期二除了可以自己構建輸入文件外，GVIEW還可讀入CHEM3D，HYPERCHEM和晶體數(shù)據等諸多格式的文件。從而使其可以于諸多圖形軟件連用，大大拓寬了使用范圍（詳見下圖）第24頁，共54頁，2023年，2月20日，星期二GVIEW的界面及主要功能鍵的介紹界面簡介第25頁，共54頁，2023年，2月20日，星期二第一行為菜單欄:File:主要功能是建立，打開，保存和打印當前的文件SaveInage將當前文件保存為圖片格式Preferences?？梢栽诶锩娓淖僄VIEW默認的各種顯示的設置。第26頁，共54頁，2023年，2月20日，星期二Edit:在這里可以完成對分子的剪且，拷貝，刪除和抓圖等。AtomList,顯示當前分子的內坐標，笛卡兒坐標，分數(shù)坐標等。PointGroup可以顯示當前分子的點群及可能有的點群。PBC

顯示晶體文件（可以將CIF文件轉換為圖形，在點PBC按鈕后所給并的對話框中根據選項調節(jié)具體顯示的格式。Mos

用于顯示分子軌道（只有檢查點文件，此選項才能給出分子軌道圖。Symmetrize，對當前體系進行對稱性控制。第27頁，共54頁，2023年，2月20日，星期二View

這里面的選項都是于分子的顯示有關的，如顯示氫原子，顯示鍵，顯示元素符號，顯示啞原子顯示坐標軸等第28頁，共54頁，2023年，2月20日，星期二Calculalte：可從GVIEW中直接向高斯提交計算。這是GVIEW作為高斯軟件配套功能的重要體現(xiàn)。從所給的對話框中可以選擇工作類型JobType（如優(yōu)化，能量或頻率等）;計算方法Method（如半經驗方法，HF方法，DFT方法，MP方法等，還可以選定基組）;Title(對所要做的計算給一個說明，以備以后的查看)Link0（給檢查點文件命名，還可以在此用RWF命令設置臨時數(shù)據交換文件的大?。籊eneral,Guess,（這兩個選項主要是給出體系中各原子的連接關系及如何給出初始猜測）；NBO（可在此設定NBO計算）,PBC（可在此設定晶體的有關計算）,Solvation（可在此設定溶液中的計算，除了選擇溶劑外，還要選擇模擬溶劑的理論模型）第29頁，共54頁，2023年，2月20日，星期二Result：顯示計算的結果，包括電荷，靜電勢表示的表面，振動頻率，核磁，勢能面掃描，優(yōu)化等。注意有些結果只能用檢查點文件才能顯示。第30頁，共54頁，2023年，2月20日，星期二構建一個間氟苯乙烷分子并從GVIEW里遞交計算(1)雙擊GVIEW圖標，打開GVIEW，下圖就是GVIVE打開后的窗口第31頁，共54頁，2023年，2月20日，星期二(2)．雙擊窗口中的圖標，得到如下窗口里面有常用的環(huán)狀官能團。選中苯環(huán)（單擊即可選中）第32頁，共54頁，2023年，2月20日，星期二(3).在當前工作窗口（打開GIVEW時程序自動打開一個工作窗口，如下圖）也可通過File-new路徑新建一個工作窗口在這個窗口中點鼠標左鍵窗口中就會出現(xiàn)苯分子，見下圖將鼠標放在分子上，按左鍵左右或前后移動，可以調節(jié)分子的角度將鼠標放在分子上，前后移動，可以將分子放大或縮小Shift+Alt+鼠標左鍵組合可以在窗口內平移分子當工作窗口內有多個分子時【在構建大的分子時，這種情況很容易出現(xiàn)】?？梢杂肅trl+Alt+鼠標左鍵組合調節(jié)其中一個分子的角度，以調節(jié)各個分子間的角度。第33頁，共54頁，2023年，2月20日，星期二(4)．雙擊GVIEW界面上的圖標。出現(xiàn)窗口:點擊氟的元素符號“F”，就選中氟原子。(5).回到工作窗口在苯環(huán)上的任一個H上單擊左鍵，將H置換成F。第34頁，共54頁，2023年，2月20日，星期二但要注意此時僅是元素符號發(fā)生了改變，C-F間的距離仍是C-H鍵的距離。需要用鍵長工具進行調整。單擊GVIEW界面圖標。然后再點擊工作窗口中的C原子和F原子，看到被選中的兩個原子于周圍的原子再亮度上有差異

此時會出現(xiàn)右邊的窗口根據C-F鍵的長度在0.675和2.700之間的方框內進行C-F鍵的調整。完畢后點擊OK即可。第35頁，共54頁，2023年，2月20日，星期二(6)．雙擊GVIEW界面上的圖標，出現(xiàn)以下窗口這是GVIEW里內置的鏈烴庫，選中乙烷第36頁，共54頁，2023年，2月20日，星期二(7).在工作窗口內空處點左鍵用我們前面講過的命令調節(jié)苯和乙烷間的距離和角度將鼠標放在分子上，按左鍵左右或前后移動，可以調節(jié)分子的角度；將鼠標放在分子上，移動滾輪，可以將分子放大或縮小。Shift+Alt+鼠標左鍵組合可以在窗口內平移分子當工作窗口內有多個分子時【在構建大的分子時，這種情況很容易出現(xiàn)】這時可用以下命令可以用Shift+Alt+鼠標左鍵組合移動想要移動的分子，以調節(jié)各個分子間的距離。可以用Ctrl+Alt+鼠標左鍵組合調節(jié)其中一個分子的角度，以調節(jié)各個分子間的角度Ctrl+Alt+鼠標左鍵這個組合常和[將鼠標左鍵放在分子上，左右或前后移動，可以調節(jié)分子的角度]這個功能連用。具體的好處需要自己去試。第37頁，共54頁，2023年，2月20日，星期二（8）．單擊GVIEW界面上圖標。然后點擊與乙烷正對的苯上的C原子和H原子選擇“None”，乙烷上的C-H鍵作類似處理然后單擊GVEW上的 圖標，點選工作窗口中的兩個氫原子，即將它們刪去第38頁，共54頁，2023年，2月20日，星期二再單擊GVEW界面上圖標，根據實際鍵長將兩個C原子聯(lián)起來。第39頁，共54頁，2023年，2月20日，星期二（9）．查看分子的對稱性從Edit-Pointgroup路徑可以查看所構建分子的點群。點擊Pointgroup后，出現(xiàn)如下窗口：為C1點群，其下拉菜單中的為可能的點群（改變Tolerance,也可幫助我們判斷所構建體系可能有的點群）第40頁，共54頁，2023年，2月20日，星期二(10)．查看分子坐標。單擊GVIEW界面上的圖標。出現(xiàn)下面窗口

Z表示時內坐標，C表示直角坐標?？梢栽诶锩鎸ψ鴺俗鲞m當?shù)恼{整第41頁，共54頁，2023年，2月20日，星期二（11）．向GAUSS遞交計算。點GVIEW界面上Calculation中的Gaussian第42頁，共54頁，2023年，2月20日，星期二會出來一個遞交計算的對話框。從所給的對話框中可以選擇工作類型JobType（如優(yōu)化，能量或頻率等）;計算方法Method（如半經驗方法，HF方法，DFT方法，MP方法等，還可以選定基組;Title(對所要做的計算給一個說明，以備以后的查看)Link0（給檢查點文件命名，還可以在此用RWF命令設置臨時數(shù)據交換文件的大?。籊eneral,Guess,（這兩個選項主要是給出體系中各原子的連接關系及如何給出初始猜測）；NBO（可在此設定NBO計算）,PBC（可在此設定晶體的有關計算）,Solvation（可在此設定溶液中的計算，除了選擇溶劑外，還要選擇模擬溶劑的理論模型）選擇完畢后，點Submit即可遞交計算。有時由于安裝的原因GVIEW無法與GAUSS建立關聯(lián)，就不能直接從GVIEW里遞交計算。這時可以在GVIEW里保存用于GAUSS計算的輸入文件，然后從GAUSS里調出文件進行計算。第43頁，共54頁，2023年，2月20日，星期二3．構建分子中的注意事項在構建分子前，應對分子的鍵長，鍵角和點群等有一個詳細的了解。這樣在構建過程中才能更好的把握細節(jié)，如構建大分子時，各個集團間的鍵長，鍵角和二面角一定要給的準確。以我的經驗而言，相當一部分的錯誤（如SCF不收斂，1502報錯等）都是分子建模不合理造成的。分子模型不合理導致程序根據分子模型所得的薛定鄂方程不合理，也就導致了對方程的解不能順利完成，所以報錯。對于大分子，不會出現(xiàn)分子直線排步的情況，鏈烴通常為鋸齒型的排列，環(huán)之間總是有一定的夾角，共面情況雖有但不多見，這些都是在構建分子模型的過程中需要注意的。再者，分子建模是量化計算的一個難點，特別是對金屬配合物，塊狀體系和納米體系的構建。只有給出合理的結構，程序才能給出合理的解。我是做有機體系的計算，對以上三種不熟，還請各位同學見諒。如果不能體現(xiàn)出分子的對稱性，則程序不能正確判斷分子軌道的對稱性，這在涉計到反應機理和軌道相互作用的研究中可能會導致錯誤的判斷。第44頁，共54頁，2023年，2月20日，星期二練習練習2.1利用GVIEW構建一個二茂鐵分子練習2.2利用GVIEW構建一個2-苯基吲哚分子練習2.3利用GVIEW構建一個三苯胺分子練習2.4利用GVIEW構建一個醋酸鈀分子第45頁，共54頁，2023年，2月20日，星期二第三章單點能計算第46頁，共54頁，2023年，2月20日，星期二單點能計算是指對給定幾何構性的分子的能量以及性質進行計算，由于分子的幾何構型是固定不變的，只是"一個點"，所以叫單點能計算。單點能計算可以用于：1.計算分子的基本信息2.可以作為分子構型優(yōu)化前對分子的檢查3.在由較低等級計算得到的優(yōu)化結果上進行高精度的計算4.在計算條件下，體系只能進行單點計算5.單點能的計算可以在不同理論等級，采用不同基組進行第47頁，共54頁，2023年，2月20日，星期二例

甲醛的單點能計算第48頁，共54頁，2023年，2月20日，星期二輸出文件中的信息能量找到SCFDone：E(RHF)=-113.863697598A.U.after6cycles這里的數(shù)值就是能量，單位是hartree。在一些高等級計算中，往往有不止一個能量值，比如下一行：E2=-0.3029540001D+00EUMP2=-0.1141