α烯烴交替聚合反應(yīng)熱的計算機模擬計算及計算結(jié)果的分析

1、-烯烴交替聚合反應(yīng)熱的計算機模擬計算及計算結(jié)果的分析1 guassian03軟件簡介熱化學(xué)性質(zhì)是根據(jù)正則系綜理想氣體展開得到的。在大話西游論壇上看到過這樣一種說法：氣相下計算的能量（scf能量）是在真空狀態(tài)下計算的體系的能量，而熱化學(xué)性質(zhì)是在298.15k一個大氣壓下條件下計算得到。熱力學(xué)計算中也包括了零點能的輸出，零點能是對分子的電子能量的矯正，計算了在0k溫度下的分子振動能量。當(dāng)比較在0k的能量時，需要在總能量中加上零點能。和頻率一樣，理論模型本身也給零點能計算帶來系統(tǒng)誤差，可以通過矯正因子修正來和實驗值相符。如果沒有設(shè)置freq=readisotopes關(guān)鍵詞，并且設(shè)置矯正因子，那么就需

2、要手工對所計算的能量進行修正。2 乙烯聚合反應(yīng)熱的計算聚合反應(yīng)與一般的小分子反應(yīng)是有差別的?，F(xiàn)已乙烯聚合反應(yīng)為例，說明聚合熱的計算方法。聚合熱計算公式的推導(dǎo)：乙烯聚合反應(yīng)的反應(yīng)式：則聚合反應(yīng)熱（焓）為：用gaussian03計算出當(dāng)n不同時的值（令不同n值時的值為an），結(jié)果如下：當(dāng)n=1時，a1= -79.732696 當(dāng)n=2時，a2= -158.280787 當(dāng)n=3時，a3= -236.828725 當(dāng)n=4時，a4= -315.376650 當(dāng)n=5時，a5= -393.924564 當(dāng)n=6時，a6= -472.472482單體乙烯（）的為a0= -78.516447，氫（）的為由

3、上述數(shù)據(jù)可以得出如下規(guī)律：a2=a1+a1-+0.022569a3=a2+a1-+0.022722a4=a3+a1-+0.022735a5=a4+a1-+0.022746a6=a5+a1-+0.022742以上常數(shù)項比較接近，故令b=0.0227，再把各項代入可得到一個通式為：所以反應(yīng)熱（焓）為：當(dāng)時，則有： 代入數(shù)據(jù)可得:將以上hartree能量單位換算為kj/mol，則有：乙烯聚合熱（氣液聚合熱gc）的實驗值為：-101.5kj/mol，說明這種計算方法的系統(tǒng)誤差為-18.76 kj/mol。3 己烯聚合反應(yīng)熱的計算己烯共聚反應(yīng)的反應(yīng)式為： 己烯聚合物的實際結(jié)構(gòu)如下：n-bu表示正丁基。則

4、己烯聚合反應(yīng)熱（焓）為：用gaussian03計算出當(dāng)n不同時的值（令不同n值時的值為an），結(jié)果如下：當(dāng)n=1時，a1= -236.828725 當(dāng)n=2時，a2= -472.465624 當(dāng)n=3時，a3= -708.103331 當(dāng)n=4時，a4= -943.748826 當(dāng)n=5時，a5= -1179.386868 當(dāng)n=6時，a6= -1415.032157當(dāng)n=7時，a7= -1650.670105當(dāng)n=8時，a8= -1886.315390單體己烯（）的為a0= -235.618500，氫（）的為。從聚己烯的實際結(jié)構(gòu)可以看，主鏈上的正丁基有2種相對位置，即1、4位相隔和1、2位相

5、隔。當(dāng)有偶數(shù)個己烯相聚合時，最末端的2個正丁基是1、4位相隔（間）；當(dāng)有奇數(shù)個己烯相聚合時，最末端的2個正丁基是1、2位相隔（間），這兩種情況是不同的。觀察以上的數(shù)據(jù)（a2 a8），發(fā)現(xiàn)有如下的規(guī)律：偶數(shù)己烯聚合時： 其通式為：所以， 其中，b=0.0203己烯的聚合反應(yīng)熱（焓）為：當(dāng)時，則有：奇數(shù)己烯聚合時： 其通式為： 所以，= 其中己烯的聚合反應(yīng)熱（焓）為：當(dāng)時，則有：因此，己烯聚合熱（液液聚合熱1c）的實驗值為：-83.0kj/mol，說明這種計算方法對于計算聚己烯類型（正丁基1、4位相隔和1、2位相隔交替出現(xiàn)）的聚合熱的系統(tǒng)誤差為-22.1698 kj/mol（約為-22.2 kj/

6、mol）。4 苯乙烯聚合反應(yīng)熱的計算苯乙烯共聚反應(yīng)的反應(yīng)式為：苯乙烯聚合物的實際結(jié)構(gòu)如下：則苯乙烯聚合反應(yīng)熱（焓）為：用gaussian03計算出當(dāng)n不同時的值（令不同n值時的值為an），結(jié)果如下：當(dāng)n=1時，a1= -310.643911 當(dāng)n=2時，a2= -620.103669 當(dāng)n=3時，a3= -929.555628 當(dāng)n=4時，a4= -1239.012667 當(dāng)n=5時，a5= -1548.468890 當(dāng)n=6時，a6= -1857.926043當(dāng)n=7時，a7= -2167.381969當(dāng)n=8時，a8= - 2476.839182單體苯乙烯（）的為a0= -309.4386

7、62，氫（）的為。從聚苯乙烯的實際結(jié)構(gòu)可以看，苯環(huán)有2種相對位置，即1、4位相隔和1、2位相隔。當(dāng)有偶數(shù)個苯乙烯相聚合時，最末端的2個苯環(huán)是1、4位相隔（間）；當(dāng)有奇數(shù)個苯乙烯相聚合時，最末端的2個苯環(huán)是1、2位相隔（間），這兩種情況是不同的。觀察以上的數(shù)據(jù)（a2 a8），發(fā)現(xiàn)有如下的規(guī)律：偶數(shù)苯乙烯聚合時： 其通式為：所以， 其中，b=0.0283苯乙烯的聚合反應(yīng)熱（焓）為：當(dāng)時，則有：奇數(shù)苯乙烯聚合時： 其通式為： 所以，= 其中苯乙烯的聚合反應(yīng)熱（焓）為：當(dāng)時，則有：因此，苯乙烯聚合熱（液液聚合熱1c）的實驗值為：-70.0kj/mol，說明這種計算方法對于計算聚苯乙烯類型（苯基1、4位相隔和1、2位相隔交替出現(xiàn)）的聚合