講課-溶劑效應(yīng)課件

溶劑效應(yīng)的量子化學(xué)研究

化學(xué)是研究分子性質(zhì)和反應(yīng)的科學(xué),有許多化學(xué)過程是在溶液中進(jìn)行的，而氣相反應(yīng)在化學(xué)反應(yīng)中只占很小的部分。只討論分子或離子在氣相環(huán)境中的化學(xué)行為,與溶液中體系還有很大的差別,所以溶劑效應(yīng)是我們必須面對的課題。且大量的研究也已經(jīng)表明,分子在溶液中的性質(zhì),與真空相比,無論是幾何構(gòu)型、構(gòu)象、電子結(jié)構(gòu)還是其穩(wěn)定性和光譜性質(zhì)都會發(fā)生不同程度的變化。另外,溶劑(水)分子還可能通過氫鍵作用直接參與溶質(zhì)與受體之間的結(jié)合，所以對溶劑效應(yīng)的理論研究有著重要的意義,這更能反映出溶液中分子的真實(shí)行為。然而,以往的研究中都是將氣相下的計(jì)算結(jié)果直接應(yīng)用于溶液體系或與溶液中的實(shí)驗(yàn)結(jié)果進(jìn)行比較，這不但無法反映出溶液中分子的真實(shí)行為,而且對一些溶劑敏感的體系來說甚至可能得出錯(cuò)誤的結(jié)論。研究溶劑效應(yīng)的意義溶劑效應(yīng)的發(fā)展

自20世紀(jì)初期schr?dinger,等創(chuàng)立的量子力學(xué)體系以來,科學(xué)家們就開始用量子力學(xué)方法處理化學(xué)問題。到了50年代由于計(jì)算機(jī)的出現(xiàn),為量子化學(xué)提供了有力的工具,分子軌道(MO)理論因易于程序化而蓬勃發(fā)展起來。從70年代開始,MO的從頭計(jì)算研究逐步展開,80年代,量子化學(xué)的研究對象也從中小分子向較大分子、重原子體系發(fā)展,隨之而來了各種計(jì)算方法和程序化,是MO理論和從頭算技術(shù)大發(fā)展時(shí)期。埃爾溫·薛定諤對于溶劑效應(yīng),按溶劑和溶質(zhì)相互作用的性質(zhì)一般可劃分成兩類,一類是把溶劑看成連續(xù)介質(zhì),即非特異性溶劑效應(yīng)，該模型主要從溶劑對溶質(zhì)的極化作用角度來考察溶劑的影響,因而也稱物理溶劑效應(yīng)；另一類是把溶劑看成具體的分子,即特異性溶劑效應(yīng)，該模型主要從溶劑和溶質(zhì)之間的氫鍵和電荷轉(zhuǎn)移相互作用角度來考察溶劑的影響,因而也稱為化學(xué)溶劑效應(yīng)，基于上述兩種模型發(fā)展了各種不同的理論方法，歸結(jié)起來大致有兩類：一類是采用嚴(yán)格的量子化學(xué)計(jì)算方法,如:自洽反應(yīng)場方法(self-ConsistentReactionField,SCRF)、Possion一Boltzmann方法、AMSOL方法、超分子方法(Supermolecule)以及量子力學(xué)與分子力學(xué)相結(jié)合的QM/MM方法等，另一類是采用經(jīng)典力場的動(dòng)力學(xué)模擬方法,如分子動(dòng)力學(xué)(MolecularDynamic,MD)模擬和MonteCarlo模擬等。以上諸方法各有特色和優(yōu)缺點(diǎn),有些方法還可以相互結(jié)合。其中自洽反應(yīng)場方法是目前使用較多的考慮溶劑效應(yīng)的方法,即將靜態(tài)理論方法和連續(xù)介質(zhì)溶劑模型結(jié)合起來研究溶液體系。Poisson-Boltzmann方法就是利用求解PB方程來計(jì)算溶質(zhì)與溶劑分子間的靜電極化作用。AMSOL方法采用Bom模型,進(jìn)行分子軌道方法計(jì)算,來研究溶液體系。超分子方法認(rèn)為溶劑分子可能通過氫鍵或電荷轉(zhuǎn)移作用與溶質(zhì)分子緊密結(jié)合在一起,然后將這些額外原子包含在一個(gè)更大的但仍是孤立的體系中進(jìn)行處理。QM/MM方法是把研究體系設(shè)為幾個(gè)區(qū)域,在中心區(qū)域進(jìn)行高精度的量子化學(xué)計(jì)算(QM),在周邊區(qū)域進(jìn)行半經(jīng)驗(yàn)或分子力學(xué)計(jì)算(MM),該方法也是處理溶液體系的一個(gè)較好的方法，分子動(dòng)力學(xué)或MonteCarlo模擬溶液體系,最適用于研究溶質(zhì)分子周圍溶劑分子的結(jié)構(gòu)性質(zhì),而對溶劑化自由能的計(jì)算則沒有自洽反應(yīng)場方法計(jì)算得到的結(jié)果精確。溶劑效應(yīng)

所謂溶劑就是能夠溶解其他物質(zhì)的液體。水是一種特殊的溶劑,有著極高的相對介電常數(shù),既有很弱的酸性,又有很弱的堿性,大多數(shù)的化學(xué)反應(yīng)和所有的生物學(xué)過程都能在水這種溶劑中進(jìn)行。而水作為溶劑,對這些過程有著深刻的影響:一種是溶質(zhì)與溶劑分子發(fā)生相互作用,如形成水合離子、氫鍵等,溶質(zhì)分子的幾何構(gòu)型、振動(dòng)頻率、總能、電子光譜等均受溶劑分子的影響;另外一種是溶劑分子不直接與溶質(zhì)分子發(fā)生作用,而是作為環(huán)境強(qiáng)烈地影響溶質(zhì)的行為,即溶劑可作為某種場來處理。因此,在理論上理解和精確計(jì)算溶劑效應(yīng)是至關(guān)重要的。

溶劑化過程可被定義為在常溫下,溶質(zhì)分子從真空的一個(gè)固定位置移至溶劑中的一個(gè)固定位置的過程。更精微地說,溶劑效應(yīng)是溶質(zhì)和溶劑分子相互吸引的結(jié)果,即溶劑分子通過它們和溶質(zhì)的相互作用,累積在溶質(zhì)周圍的過程。這個(gè)相互作用的性質(zhì)和溶質(zhì)與溶劑的本質(zhì)有關(guān)。描述溶劑效應(yīng)的一個(gè)關(guān)鍵參數(shù)就是溶劑化自由能ΔGsol(solvationfreeenergy),它是將溶質(zhì)分子從真空移至溶劑中自由能的變化值。一般來說,溶劑效應(yīng)的計(jì)算由三部分組成：靜電勢、范德華力和孔洞能。靜電效應(yīng)(ElectrostaticEffects)靜電力在分子間相互作用中處于首要地位,它的強(qiáng)度大而且是長程力。在分子中,原子核周圍的電子分布產(chǎn)生了一個(gè)靜電場,它和其他核相互作用。溶質(zhì)分子和溶劑分子的電荷分布在溶劑化過程中扮演一個(gè)根本的角色。當(dāng)在溶液中，溶質(zhì)分子由于極化而出現(xiàn)的偶極矩μ作用于溶劑分子,使得溶劑分子產(chǎn)生誘導(dǎo)偶極矩,而溶劑分子的誘導(dǎo)偶極矩反過來作用于溶質(zhì)分子,從而對溶質(zhì)分子產(chǎn)生一個(gè)額外的穩(wěn)定作用,即在溶質(zhì)分子周圍產(chǎn)生電場。這就是所謂的反應(yīng)場，它作用在溶質(zhì)分子上,改變了溶質(zhì)分子的自身的能量。而且由于溶劑分子的出現(xiàn),溶質(zhì)分子的分子內(nèi)庫侖相互作用也被屏蔽了(圖2一1)?？锥茨苁侨苜|(zhì)分子要在溶劑分子中形成孔洞所需要的能量。本質(zhì)上,這一項(xiàng)屬于熵效應(yīng)。它說明了由于在非極性溶質(zhì)分子周圍水分子的重新分布而使得熵減少了。尤其對于水,熵減少是因?yàn)橐粋€(gè)非氫鍵溶質(zhì)的出現(xiàn)而使有利于形成氫鍵的方法減少了，孔洞能還包括溶劑分子與溶劑分子間吸引-排斥相互作用的改變,而這個(gè)變化是由于孔洞內(nèi)溶劑分子的遺漏以及局部的溶劑結(jié)構(gòu)的改變。因?yàn)殪販p少,孔洞能對于溶劑化是不利的。分子的運(yùn)動(dòng)對溶劑化自由能也有貢獻(xiàn),它包括零點(diǎn)能(ZEP),熱振動(dòng)、轉(zhuǎn)動(dòng)、平動(dòng)。這些貢獻(xiàn)一般不被明確地考慮,如果計(jì)算程序(例如,G98)提供振動(dòng)頻率、標(biāo)準(zhǔn)統(tǒng)計(jì)力學(xué)計(jì)算,那么這些項(xiàng)的計(jì)算就被增加到其他自由能的項(xiàng)中一起計(jì)算。正確地對待溶劑化自由能對于理解溶液化學(xué)是至關(guān)重要的。對于化學(xué)反應(yīng),反應(yīng)物和生成物的自由能可以確定化學(xué)平衡,過渡態(tài)和生成物的自由能的差額能夠控制反應(yīng)速率因子，溶液結(jié)構(gòu),例如徑向分布函數(shù)RDF(r),也能夠從溶劑化自由能中得到。但溶劑化自由能是很復(fù)雜的,為了使得計(jì)算易處理,需做許多假設(shè)和近似,一個(gè)關(guān)鍵的假設(shè)就是我們能夠把溶劑化自由能分成不同的部分,短程范德華作用能,熵效應(yīng)和長程靜電作用能。溶劑化模型處理溶液體系的模型有很多種,但大致上可以分為兩類:不連續(xù)模型(discretemodels)和連續(xù)模型(continummodels),這兩種模型各有千秋,可互相補(bǔ)充。下面我們從不同的方面來介紹這兩種模型的特點(diǎn)。不連續(xù)模型對于不連續(xù)模型,溶劑分子被明確的考慮,所選的許多溶劑分子的構(gòu)型形式被包括在系統(tǒng)的量子力學(xué)描述中。這樣在計(jì)算中,溶劑分子的微觀結(jié)構(gòu)和溶質(zhì)-溶劑分子的相互作用能夠被清楚地描述。但是這種模型也是有缺點(diǎn)的,大量的溶劑分子需要模擬稀溶液中的溶劑分子,而且無論是QM方法還是MM計(jì)算方法都會限制不連續(xù)模型的應(yīng)用。不連續(xù)模型通常也不考慮溶質(zhì)的電子極化效應(yīng),而且溶質(zhì)和溶劑之間的靜電作用也依賴于所選力場的電荷參數(shù)。另外這種模型的計(jì)算效率不高,因?yàn)樵谟?jì)算中每個(gè)原子都被清楚地對待,就要增加三個(gè)自由度。另一方面,不連續(xù)計(jì)算提出將整個(gè)波函數(shù)分成兩個(gè)獨(dú)立的部分,分別來描述溶劑和溶質(zhì)分子,溶劑-溶劑和溶質(zhì)-溶劑的相互作用用半經(jīng)典的方法來處理。因此,不連續(xù)模型仍被廣泛地應(yīng)用于一些量子力學(xué)或分子力學(xué)方法中,并能得出好的結(jié)果。連續(xù)模型1.SASA方法緊密表面層溶劑法的模型第一次被Eisenberg在1986年提出。在這個(gè)模型里,溶質(zhì)分子由最靠近它的與它作用最大的一層溶劑分子環(huán)繞著,用測定每個(gè)原子或基團(tuán)與溶劑接觸的表面積計(jì)算溶液的自由能ΔGsol式中:N表示一個(gè)分子中有N個(gè)原子,σi是溶質(zhì)區(qū)域的表面矢量,Ai是該區(qū)域的表面積。該方法不區(qū)分各種能量的貢獻(xiàn),表面矢