溶劑對鎳連二硫烯與乙烯反應的影響

第57卷第1期2008年1月物理學(xué)報Vol .57 ,No.1 ,January ,20081000-3290/2008/57(01 )/0096-07ACTA PHYSICA SINICA回2008 Chin. Phys .Soe,溶劑對鎳連二硫烯與乙烯反應的影響*韓清珍))2)t耿舂字1)2)趙月紅)戚傳松3)溫浩"1) (中國科學(xué)院過(guò)程工程研究所多相實(shí)驗室，北京10080)2) (中國科學(xué)院研究生院，北京100)3) (北京石油化工學(xué)院化學(xué)工程系,北京102617)(2007年1月9日收到:2007年4月27日收到修改稿)用密度泛函理論在B3LYP/6-31G(d)水平上計算得到了鎳連二硫烯與乙烯反應的勢能面上各駐點(diǎn)(反應物、中間體、產(chǎn)物和兩個(gè)過(guò)渡態(tài))的分子幾何構型、電荷分布和一-些熱力學(xué)參數等,研究了溶劑對鎳連二硫烯與乙烯反應的影響.結果發(fā)現.隨著(zhù)熔劑極性的增強,乙烯和鎳連二碗烯之間的成鍵作用增強,兩個(gè)過(guò)渡態(tài)的前線(xiàn)軌道能量差增大,而產(chǎn)物和中間體的前線(xiàn)軌道能量差卻減小，同時(shí)各駐點(diǎn)的溶劑穩定化能也減小.進(jìn)-步,這表明溶劑極性增強能提高產(chǎn)物的穩定性,有利于反應的進(jìn)行.此外,當溶劑相對介電常數1.00≤e≤7.58時(shí),溶劑效應十分顯著(zhù),而在e>7.58時(shí),落劑效應的影響將達到極限.關(guān)鍵詞:鎳連二硫烯,乙烯,溶劑效應，密度泛函理論PACC: 3130L, 7115M于反應位點(diǎn)是硫原子而不是金屬原子,避免了過(guò)渡1.引 言金屬連二硫烯配合物被裂解氣中H, H2S, CO, CO2和C2H2等組分毒化的問(wèn)題，顯著(zhù)提高了烯烴分離過(guò)烯烴是生產(chǎn)聚烯烴、聚酸、聚酯和聚醚等材料的程的選擇性和經(jīng)濟性.重要原料,我國國民經(jīng)濟的快速發(fā)展使烯烴的需求在理論研究方面,Fan等[61采用密度泛函理論量日益增長(cháng)1,21.以乙烯為代表的烯烴的生產(chǎn)和分離的B3LYP方法07,.8]研究了氣相中乙烯與已經(jīng)成為工業(yè)生產(chǎn)中的關(guān)鍵環(huán)節之一.目前,用于烯Ni(S,CR2)2(R = H, CN, CF,)的反應,發(fā)現該反烴分離的主要方法有深冷分離法和化學(xué)吸附應是一個(gè)兩步的反應過(guò)程,同時(shí)指出,強吸電子基團法[4-8].然而,這些方法存在著(zhù)選擇性低和能耗高的能降低反應的位阻,提高反應中間體和產(chǎn)物的穩定缺點(diǎn).因此,尋找經(jīng)濟高效的烯烴分離方法必將對石性,有利于反應的進(jìn)行.實(shí)際上,過(guò)渡金屬連二硫烯油化工的節能降耗和技術(shù)進(jìn)步產(chǎn)生積極的影響.配合物與烯烴的反應幾乎都是在溶劑中進(jìn)行的,溶20世紀70年代，含雙五元環(huán)結構的過(guò)渡金屬劑將對參與反應的分子的幾何和電子結構產(chǎn)生影連二硫烯配合物因其獨特的幾何結構以及光學(xué)、電響|),進(jìn)而對整個(gè)反應的進(jìn)行產(chǎn)生影響.因此,研究學(xué)和化學(xué)氧化還原特性而受到人們的廣泛關(guān)溶劑效應將對該反應用于烯烴分離的溶劑選擇具有注{9,01 .研究表明11-1)1:含雙五元環(huán)的過(guò)渡金屬絡(luò )合重要的指導意義.據我們所知,目前尚未有過(guò)渡金屬物能在溫和的條件下與降莰二烯和1,3-二甲基-1,連二硫烯配合物與烯烴反應的溶劑效應的報道.3-丁二烯等二烯烴發(fā)生反應,生成1/1的加合物,且反應位點(diǎn)是硫原子.最近, Wang等[.14151利用UV-vis2.計算方法分光光度法對這類(lèi)反應在甲苯中的反應情況進(jìn)行了研究,結果發(fā)現:在甲苯溶液中,常溫常壓下,Ni, Pa,本文首先用密度泛函理論的B3LYP方法對氣Pt等過(guò)渡金屬連二硫烯配合物能選擇性地、可逆地相中中國煤化工了研究.其中,鎳與簡(jiǎn)單烯烴發(fā)生反應,生成1/1的加合物.同時(shí),由原子MYHCNMHG]修正的LANL2DZ.國家自然科學(xué)基金(批準號:20221603.20703047, 10604005)資助的課題.↑E-mail: qphan@ home.ipe. se.cn1期韓清珍等:溶劑對鎳連二硫烯與乙烯反應的影響97(341/341/41)基組,碳、硫、氫等其他原子則采用6-頻率分析確認各駐點(diǎn)的結構.通過(guò)分析不同溶劑中31C(d)基組'0- 21.利用Gaussian03 程序中的線(xiàn)性?xún)确肿拥膸缀螛嬓?、電荷分?前線(xiàn)軌道能級差和溶劑坐標(QST3)法搜索到兩個(gè)過(guò)渡態(tài)[切] ,并對過(guò)渡態(tài)進(jìn)穩定化能,研究了溶劑對鎳連二硫烯與乙烯反應的行了振動(dòng)頻率分析和內稟反應坐標|]8-301計算的驗影響.證.同時(shí)通過(guò)內稟反應坐標計算得到了反應物、中間體和產(chǎn)物的穩定構型,通過(guò)頻率計算的振動(dòng)分析確3.結果與討論認均為勢能面上的極小點(diǎn).進(jìn)一步,用自治反應場(chǎng)"理論中的Onseger溶3.1.反應途徑與幾何構型劑模型(3-在反應溫度T = 298.15K和壓力P =101.3 kPa時(shí),分別在苯(e =2.247)、甲苯(2.279)、四圖1給出了計算得到的鎳連二硫烯與乙烯反應氫呋喃(7.58)、二氯甲烷(8.93).1,2二氯乙烷的反應途徑以及5個(gè)反應駐點(diǎn)的幾何構型.鎳連二(10.36)、丙酮(20.7).乙醇(24.55)、甲醇(32.63)、硫烯與乙烯的反應是分兩步進(jìn)行的.反應的第一步1,2,3-丙三醇(42.50)、二甲亞砜(46.7)和水(78.39)是一個(gè)加成反應,乙烯中的碳原子配位到鎳連二硫共11種非極性和極性溶劑中(351,重新優(yōu)化鎳連二烯的硫原子上.在此過(guò)程中,作為電子受體的硫原子硫烯與乙烯反應勢能面上的所有駐點(diǎn).溶劑效應的計算中,首先根據各個(gè)反應駐點(diǎn)在氣相(ε = 1)中的2p空軌道進(jìn)攻乙烯的兩個(gè)碳原子上的π電子形的優(yōu)化構型來(lái)求得溶質(zhì)分子的球形空穴半徑刈),然成反式結構的中間體配合物.第二步是因該反式結后保持與氣相計算時(shí)相同的計算水平,分別在不同構的中間體不穩定，而迅速轉變?yōu)榉€定的順式結構溶劑中進(jìn)行全優(yōu)化計算,并通過(guò)對所有駐點(diǎn)的振動(dòng)的產(chǎn)物的過(guò)程'06).。C6田54CI日ssC9“由C面C15 g心反應物(R)過(guò)渡有1(TS)中間體(1)C14C產(chǎn)物(P)過(guò)渡毒2(TS2)11反應途徑表1是5個(gè)反應駐點(diǎn)在氣相和幾種典型溶劑中nm)增大了約0.02 nm,這主要是由于在反應物乙烯的主要幾何參數.考慮到計算中采用的是Onsager模中的Cu--C,s為雙鍵,而在產(chǎn)物中電子的重新排布型,且反應物具有D2對稱(chēng)性,因此,反應物的幾何使C鍵杜亦由了的錯守i Tan等的研究結果參數在所有溶劑中均不發(fā)生改變(7,31) ,故表1中僅一致‘中國煤化工?連二硫烯分子中給出了反應物在氣相中的幾何參數.首先,從表中可的Ni.CNMHG..分別增大了約以發(fā)現,氣相中產(chǎn)物的C-C鍵長(cháng)(Ru-1s=0.15330.009 nm和0.003 nm,這是因為隨著(zhù)反應的進(jìn)行,鎳nm)比反應物中相應的C- -C 鍵長(cháng)(Ru-1s =0.1311連二硫烯分子中Ni原子與S原子,以及C原子與S98物理學(xué)報57卷原子之間的相互作用由原來(lái)的多中心大π鍵分別轉順式結構的產(chǎn)物[6),導致處于反應活性中心的C- S變?yōu)榕湮绘I和共價(jià)單鍵,這與Szilagyi等研究的實(shí)驗鍵長(cháng)R_u、鎳連二硫烯的夾角∠S,NiS,以及二面角事實(shí)一-致(川).另外,反應前為平面型的反應物逐漸β(S,NiSS)隨著(zhù)反應的進(jìn)行先增大后減小，而夾角扭曲,形成反式結構的中間體,而后又繼續彎曲形成∠Nis,C,單調減小.表1幾何結構參數駐點(diǎn)溶劑εR_s/m .Rs_g/mmRu- s/omR-ud/m∠SINis/() ∠iSC5/(°) B(S, Ni s,)()反應物0.22060.17050.131188.08124.600.1373過(guò)渡態(tài)11.0000.21830.17270.14110.237796.36107.5522.492.2470.21820.17280.14100.238896.27107.8022.527.5800.17290.14080.240096.15108.020.700.14070.240496.11108.1522.5178.390.21810.240796.08108.2022.50中間體0.22390.17480.15320.186091.35100.7019.050.22370.17510.186191.16101.00.22350.17540. 1531101.5219.040.17550.1531).186190.70101.6819.030.22330.17560.15300.186290.65101.7519.02過(guò)渡態(tài)21000.21490. 17440. 15430.183384.9298.295.4280.21500.17470.15420.183584.7398.695.4880. 15410.18384.7298.715.5440.17520.15410.183884.43).305.5600.21520.17530.15390.183984.3999.38 .5.566產(chǎn)物1.000.22950.17390. 15330.187385.4286.932.7370.22930.17410.187287.402.5887.800.22920.17430. 15320.18785.3687.912.5610.22910.17440.187188.112.5290.22900. 174585.3388.202.031另-方面,溶劑的存在并沒(méi)有導致中間體、過(guò)渡起幾何參 數的顯著(zhù)變化,溶劑效應達到-一個(gè)極限,這態(tài)和產(chǎn)物分子結構的本質(zhì)性變化.因此,溶劑中的反與文獻[40]的結論類(lèi)似.應機理和反應途徑應與氣相中的一致.值得注意的3.2.原子電荷分布是,溶劑效應卻仍然引起了各駐點(diǎn)的鍵長(cháng)和鍵角的些微變化.其中,受溶劑影響較大的是位于反應活性為更好地理解溶劑對鎳連二硫烯與乙烯反應的中心的鎳連二硫烯的C- -S鍵長(cháng)Rg_,,乙烯的C-C影響,表2給出了不同溶劑環(huán)境下分子中的電荷分鍵長(cháng)Ru_-is和連接二者的S-C鍵長(cháng)R_1u,即隨著(zhù)布.在鎳連二硫烯與乙烯的反應中，鎳連二硫烯中S溶劑相對介電常數e的增大,Rs_。和R_.u變長(cháng),而原子的空2p軌道進(jìn)攻乙烯,與乙烯中的π電子配.位,形成新的配合物.以氣相中的反應為例,從反應Ru_1s變短,相應的夾角和二面角變小.這表明隨著(zhù)物到中間體S(3)和C(14)的原子凈電荷數分別從溶劑極性增大,乙烯與鎳連二硫烯間的成鍵作用增0.1981中國煤化工6e和- 0.4416e.強,產(chǎn)物的穩定性提高,從而有利于反應的進(jìn)行.此另外,JM.HCNMH(活性中心位置的外,當溶劑的相對介電常數e約小于7.58時(shí),溶劑時(shí)時(shí)電何雙元增 加后減少,同效應能引起分子幾何參數的明顯變化,但當溶劑的時(shí)與之相連的Ni(1)和C(7)則表現為先減少后增相對介電常數大于7.58時(shí),溶劑的作用基本不再引加,這是因為在反應進(jìn)行中,會(huì )首先形成一個(gè)亞穩態(tài)1期韓清珍等:溶劑對鎳連二硫烯與乙烯反應的影響99的反式中間體,然后再生成穩定的順式產(chǎn)物[6],與反應中心的S(2)原子和C(6)原子在反應過(guò)程中凈前面根據幾何參數的變化得到的結果- .致;而遠離電荷數的變化卻表現出一定的無(wú)序性.表2原子電荷分布(e)駐點(diǎn)溶劑εNi(1)S(2)S(3)C(7)C(6)C(14)反應物- 0.54860. 18970.1981-0.2533-0.2514- 0.2783過(guò)踱態(tài)11.000- 0.50200.1201 .0.2375-0.2367- 0.34562.247-0.50070. 10460.2383- 0.2360- 0.2807-0.34217.580- 0.49920.08920.2391-0.2353- 0.2796-0.338420.70- 0.49870.08340.2394- 0.2351-0.2792-0.336946.70- 0.49850.08130.2395- 0.2350- 0.2791- 0.336478.39 .- 0.49840.0806- 0.2349- 0.2790- 0.3362中間體- 0.40860.05730.3496- 0.2343- 0.3116-0.4416-0.40570.03300.3556-0.2334-0.3088-0.4439- 0.40230.00800.3619-0.2328-0.3059- 0.4463- 0.4011- 0.00130.3641- 0.2325- 0.3048- 0.4472- 0.4006- 0.00460.3649- 0.2324- 0.3045- 0.447578.39- 0.4004- 0.00570.3651- 0.3043- 0.4476過(guò)渡態(tài)21.00- 0.38220.0760).2925- 0.2920.0.2591-0.4240.247-0.38140.04960.295-0.2905- 0.2559-0.4245- 0.38010.02200.2992-0.2889- 0.2527-0.4248- 0.37940.01150.3005- 0.2883- 0.2516- 0.4249- 0.37910.0077.309- 0.2881- 0.2511- 0.37900.00650.3010- 0.28808- 0.2510產(chǎn)物1000- 0.38680.07370.2840- 0.2213. 0.3221- 0.4107- 0.38910.05590.2904- 0.2208.- 0.3199-0.4118 .- 0.39060.03820.2973- 0.2205- 0.3178-0.4128-0.39100.03150.2998-0.2204- 0.3170-0.4132-0.39110.0292 ;0.300-0.2203-0.316-0.4133- 0.3912 ;0.02840. 3010- 0.22030.3166- 0.4134下面進(jìn)一步考慮溶劑極性對電荷分布的影響.(3)原子與C(7)原子之間的相互作用增強.總之,原如表2所示,在不同溶劑中,溶液極性的增強將使中子凈電荷數的變化均說(shuō)明,溶劑極性的增強能加劇間體和過(guò)渡態(tài)中Ni(1)原子的凈電荷數減少,使所乙烯與鎳連二硫烯間的電子傳遞和相互作用,從而有駐點(diǎn)的S(3)原子的凈電荷數增加,這表明Ni(1)促進(jìn)反應的進(jìn)行.原子與S(3)原子之間的相互作用增強.顯然這與前3.3.前線(xiàn)軌道能量差面討論鍵長(cháng)R,_,得到的結論-致.雖然S(3)原子和C(14)原子的凈電荷數隨溶液極性增強均有增加，表3給出了不同溶劑中5個(gè)反應駐點(diǎn)的前線(xiàn)軌但S(3)原子與C(14)原子之間的距離R, -14卻變化道能級差OE,即前線(xiàn)軌道Eowo減去Ewno.在氣相 .很小(見(jiàn)表1),說(shuō)明溶劑極性增強將有利于乙烯與(e=1中國煤化工氣道能級差△E先.鎳連二硫烯之間的電子轉移.與此同時(shí),隨著(zhù)溶劑極從反MHC NM H C后開(kāi)始升高到過(guò)態(tài)1的253.72 kJ,性的增強,鎳連二硫烯中S(2)原子.C(6)原子和C .繼續陶(7)原子的凈電荷數卻均減小,這表明S(2)原子與C渡態(tài)2的248.05 kJ,再繼續升高至反應產(chǎn)物的(6)原子之間的相互作用減弱,而位于反應中心的S272.93kJ.這表明兩個(gè)過(guò)渡態(tài)和中間體的活性較高,100勿理學(xué)指57卷容易發(fā)生電子躍遷;而反應物和產(chǎn)物的活性較低,能表4中數據表明,溶劑相對介電常數對各個(gè)駐點(diǎn)的級差較大,不容易發(fā)生電子躍遷,結構也較穩定.溶劑穩定化能有顯著(zhù)的影響.隨著(zhù)溶劑相對介電常表3前線(xiàn)軌逍能級差(kJ)數從2.247增大到78.39,反應物、過(guò)渡態(tài)1、中間體、溶劑e反應物 過(guò)渡態(tài)1 中間體 過(guò)渡態(tài)2 產(chǎn)物過(guò)渡態(tài)2和產(chǎn)物的溶劑穩定化能則分別從-0.0803,- 3.816, -7.710, -9.482和-5.457 kJ.1.000288.49253.72231.34 248.05 272.93mol~'降低到-0.2274, -9.843,- 19.92, -24.82和2.247245.72230.01248.17 272.78-13.85kJ.mol-'.這表明溶劑化作用將導致鎳連二2.379255.51229.93248.17 272.77硫烯~乙烯分子間相互作用增強,并對反應體系起到7.580257.16 228.56 248.28 272.55了顯著(zhù)的穩定作用.此外,在溶劑相對介電常數約小8.930257.31228.40 248.28 272.54于7.58時(shí),各駐點(diǎn)的溶劑穩定化能顯著(zhù)降低,而在10.36257.42228.32248.27 272.54溶劑介電常數超過(guò)7.58時(shí),溶劑的改變將不再引起20.70257.84227.96248.27 272.46穩定化能的較大變化.24.55257.89227.93 248.34 272.4632.63257.9927.85 248.33 272.45表4靜劑穩定化能(kJ.mol-')42.50258.05227.80248.33 272.44溶劑ε 反應物過(guò)渡態(tài)1中間體過(guò)渡態(tài)2產(chǎn)物46.70258.07227.772.247 - 0.0803-3.816-7.710 -9.482 - 5.45778.39258. 15227.70 248.33 272 .442.379 - 0.0856-4.056-8.198 - 10.09 - 5.7977.580 -0.1749 -7.697 -15.57 -19.31 - 10.89縱向比較表3數據可以發(fā)現,隨著(zhù)溶劑相對介8.930 .-0.1801- 8.025-16.23 - 20.14 -11.34電常數的增大,反應物因具有D2s 的對稱(chēng)性而未使-0.1880-8.285-16.76 -20.81 -11.70OE發(fā)生變化;而過(guò)渡態(tài)1和過(guò)渡態(tài)2的OE則分別20.70 - 0.2090-9.150-18.51 -23.03 - 12.89從253.72 kJ和248.05kJ增大到258.15 kJ和248.3324.55 - 0.2116-9.293. -18.80 -23.40 - 13.09kJ,而產(chǎn)物和中間體的OE則分別從272.93 kJ和32.63 -0.2195 -9.489 -19.20 -23.90 - 13.36.231 34kJ減小到272.44 KJ和227.70 kJ.這表明溶劑-0.2221- 9.628- 19.48-24.26 - 13.55極性的增強對產(chǎn)物和亞穩態(tài)的中間體具有明顯的穩-0.2247- 9.670-19.56 - 24.37 - 13.6定作用，而對過(guò)渡態(tài)則有明顯的降低穩定性和反應- 0.2274- 9.843-19.92 -24.82 - 13.85能壘的作用,從而促進(jìn)反應的發(fā)生和進(jìn)行.此外,各.駐點(diǎn)OE的變化幅度和趨勢與相對介電常數的大小有關(guān).當e≤7.58時(shí),各駐點(diǎn)OE的變化較大;而當e4.結> 7.58時(shí),OEu的變化則趨于平緩,此時(shí)溶劑對OE的影響逐漸減弱,并使OE趨于-一個(gè)穩定值.本文從分子的幾何結構、電荷分布、前線(xiàn)軌道以及溶劑穩定化能等方面研究了溶劑對鎳連二硫烯與3.4.溶劑穩定化能乙烯反應的影響.結果發(fā)現,溶劑將會(huì )引起幾何結表4給出了溶劑相對介電常數對各個(gè)反應駐點(diǎn)構、電荷分布、前線(xiàn)軌道以及溶劑穩定化能等變化，的溶劑穩定化能的影響.從表4中數據可以發(fā)現,在但不會(huì )導致各個(gè)反應駐點(diǎn)的幾何構型和反應機理的i同一種溶劑中(如e=2.247),隨著(zhù)反應的進(jìn)行各駐本質(zhì)變化.溶劑極性增強可以促進(jìn)乙烯與鎳連二硫點(diǎn)的溶劑穩定化能先增大后減小,這可能與5個(gè)駐烯之間的反應,使反應體系更加穩定,表現為溶劑極點(diǎn)的分子偶極矩的大小不同有關(guān).溶劑化作用會(huì )使性的增大能使乙烯與鎳連二硫烯之間的成鍵作用增分子的正、負電荷重心分離得更遠且中心電荷值更強,產(chǎn)物和中間體的前線(xiàn)軌道能量差減小,而過(guò)渡態(tài)大，從而造成偶極矩的增大,因此溶劑化能與反應駐.的前千且各駐點(diǎn)的溶劑點(diǎn)的分子偶極矩常常呈現正相關(guān)的關(guān)系[4).然而，穩定中國煤化工;此外,當溶劑的.溶劑化能增減的幅度并不完全與穩定態(tài)和過(guò)渡態(tài)偶相對介CN MH. G劑效應十分顯.極矩的相對大小成比例.著(zhù),而當ε>7.58時(shí),溶劑效應將達到極限.下面研究溶劑對駐點(diǎn)的溶劑穩定化能的影響.1期韓清珍等:溶劑對鎳連二硫烯與乙烯反應的影響10本工作的部分計算是在中國科學(xué)院超級計算中心深騰6800工作站上完成的,在此表示誠摯的謝意.[1] Sundaram K M, Shrechan M M, Olszewski E F 1995 Encydopedia薛衛東.朱正和、蔣剛、王紅艷 2003物理學(xué)報52 1617]of Chenial Technology (New York: Wily)[22] HayPJ, Wadt w R 1985 J. Chem. Phys .82 270[2] WangK, Siel E12001 Ssience 291 106[23] Coury M, Hall M B 1996 J. Compu. Chem .17 13593] Crantom R L, Royer DJ 1987 UImnn' : Engydopedia of Indwtral[24] Schrdteger P, Dolg M, EugenS w H 1989 J. Chem. Phya.91Chemistry (New York: VCH)1762[4] SukiT, Nobel RD, KovalCA 1997 Inorg. Chem.36 136[25] Andnae D, Hauberman U, Dolg m，Stoll H, Preuss H 19905] TemryP A, Noble R D, Swnnon D, Koval C A 1997 ACHETheor . Chim . Aca.77 123Jounal .43 1709[26] BergnerA, Dolg M. Chle K w, Soll H, Preuss H 1993 Mol.6] Takashi M, Yasushi M 1994 J. Am . Chem . Soc .16 3389Phys .80 1431[7] Shibehan T, Sankane C, Maeyame M, Kobushi H, Yamaimoto T，[27] Frisch MJ, Trucks C w, Schlegel H B 2003 Gausian 03, ReririonWatase T 1996 lnorg . Chim . Aca .251 207B.05; Causian, Inec.: Pitbupb, PA[8] MeKenn M, Wwright L L, Miller DJ, Tarmer L, Haltiwangr RC,[28] Fukui K 1981 Ace. Chem. Rar .14 363DuBois M R 1983 I. Am . Chem . Soe .105 5329[29] Gonzaler C, Schlege! H B 1999 J. Phyr. Chem .94 5523[9] Baker-Hawkes MJ, BlligIE, Cray H B 1966 J. Am. Chem. Soe.[30] Harding L B, Schlegel H B, Krislhnan R, PopleJ A 1980 J. Phys.88 4870Chem .84 3394[10] Krespun C C, Mckueick B C, Caimns TL 1960 J Am. Chem . Soe.[31] Tapia 0, Coscinski 0 1975 Mol . Phys .29 165382 151532] OnsogerL 1938 J. Am. Chem . Soc .58 1486[11] SchruzerCN, HoRKY, Murillo R P 1970 J. Am. Chem. Soe.[33] WibergKB, KeithTA, Frisch MJ, Murcko M A 1995 J. Phys.92 3508Chem .99 9072[12]Schrauner G N, Rabinowitz H N 1968 J. Am. Chem . Soe .90 4297[34]WongMw.WibergKB,FrischMJ1992J.Am.Chem.Soe.11413] Wing RM, Gerald C, Tustin G C, Okarmurn W H 1970 J. Am.23Chem. Sor.92 1935[35] Frisch MJ, TrucksC w, Schlegel H B 20030 Cassian 03 User' ■[14] Wang K, Stiedel E I 2000 United Sates Patent 6120692Referece Gausian, Ine. : Ptsbugh, PA1s] Wang K, Sieel E 12001 Unied States Paent 6296755[36Wong M w, WibengK B, Fisch MJ 1992 J.Am. Chem. Soce .14[16] FunYB, Hall M B 2002 J. Am. Chem. Soe .124 12076164517] Becke A D 1993 J. Chem. Plhys .98 5648[37] Chen zx, XinoJM, XinoH M, ChenY N 999. Phys. Chem.[18] LeeC, Yang w, Par R G 1988 Phys . Rev .37 785A 103 8062I9] Field MJ, Bash P A, Karplus M 1990 J. Comput . Chem.11 700[38wiberg K B, Murcko M A 1987 1. Phys . Chem .91 3616[20] LHH,LiY D, WengCK 2002 Aca Phyx. Sin.51 1239 (in [39] Sailagyi R K, Lim B s, Claser T, Holm R H, Hedmen B, HodgonChinese)[李紅海、李英德、王傳童2002 物理學(xué)報51 1239]KO, Solomon E12003 J. Am. Chem. Soc. 125 9158LiZL, WangCK, Luo Y 2004 Acra Phys. Sin.53 1490(in [40] ShiFQ, XuZC, AnJY, YuJY 2005 Aca Chim. Sin. 4283Chine)[李宗良.王傳董、羅毅 204物理學(xué)報53 1490](in Chine) [史楣強、徐志成、安靜儀俞稼鏞2005 化學(xué)學(xué)報MaY, Zou B, LiZL, WangCK, Luo Y 2006 Acta Phys. Sin.554 283]1974 (in Chinee)[馬勇、鄒斌、李宗良、王傳奎、羅 毅[41] YiPG, HuRD, YuQS, ShangZC 2005 Acta Chim. Sin. 1442006物理學(xué)報55 1974](in Chince) [易平貴、胡瑞定、俞慶森、商志才2005化學(xué)學(xué)報[21] Xioog BT, MengDQ, Xue w D, 2huZH, JingG, Wang HY144]2003 Aca Phy. Sin.52 1617 (in Chinese) [熊必濤、蒙大橋、中國煤化工MYHCNMHG102物理學(xué)報57卷The solvent effects on the complexing reaction of nickel dithioleneand ethylene *Han Qing-hen)2+ Geng Chur-Yul)2) Zhao Yue-Hong”Qi Chuan-Song'”Wen Hao'1)( MuliPhase Reaction Laboratory, Insiue of Proces Enginering. Chinese Academy of Sciences, Beijing 10080, China)2)( Graduste School of Chinese Academy of Seienes, Bejing 1000 China)3)( Chemical Enginering Department, Bejing lninde of Paro cemical Technobey, Bejing 102617, Chin)(Received 9 Janruary 2007; rerised manueip reived 27 April 2007)AbstractTheoretical studies of solvent efects on the reaction of complexing ethylene and nickel dithiolene have been caried out bycalculaing the molecular geometry, eletron disribution, and frequency of all the stagnation points existing in the reactionpolential pofiles by means of density funcional theory methods al the B3LYP/6-31C(d) level. On the basis of frequency analysisand stistial thermodynamic theory, some thermodynamic parameters for the titled compound have been obained for dfferentsolvents. It is shown that the bond energy of ethylene and nickel dithiolene enhances with the increase of the solvent polarity. In .addition, the frotier orbital energy gaps of the two tansition states will increase , while those of the product and the intermediateas well as the coreponding solvation sabiliztion energy will decrease. Furthermore, these resuls demonstrate that in polarsolvents the reaction of complexing ethylene and nickel dithiolene may become easier and faster to occur and yield a more stableproduct. Notably, an obrious solvent efet can be observed when the value of the solvent deletrc constant lie between 1.0and 7.58, while the solvent efet will tend to be maximal when the value of solvent deletriec constant exceeds 7.58.Keywords: mickel dithiolene , ethylene , solvent eft, density-functional theoryPACC: 3130L, 7115M中國煤化工MHCNMHG↑E-mail: qzhan@ home.ipe. se.cn