辜家芳許可陳文凱

(1福州大學(xué)至誠(chéng)學(xué)院化學(xué)工程系，福州350002)

(2閩江學(xué)院，福州350108)

(3福州大學(xué)化學(xué)系，福州350108)

Cp4An和COT2An(Cp-=C5H5-，COT2-=C8H82-，An=U，Pu)的結(jié)構(gòu)和光譜性質(zhì)的密度泛函理論研究

辜家芳1許可2陳文凱＊,3

(1福州大學(xué)至誠(chéng)學(xué)院化學(xué)工程系，福州350002)

(2閩江學(xué)院，福州350108)

(3福州大學(xué)化學(xué)系，福州350108)

采用密度泛函理論研究氣相和四氫呋喃(THF)溶劑中Cp4An和COT2An(Cp-=C5H5-，COT2-=C8H82-，An=U，Pu)配合物的性質(zhì)。THF溶劑對(duì)配合物的溶劑化效應(yīng)采用類導(dǎo)體極化連續(xù)模型(CPCM)近似計(jì)算。計(jì)算結(jié)果顯示在THF溶液中各配合物結(jié)合能的大小順序?yàn)镃OT2Pu＞COT2U＞Cp4Pu＞Cp4U。溶劑化效應(yīng)降低了該金屬有機(jī)配合物的結(jié)合能。計(jì)算得到的化合物的結(jié)構(gòu)參數(shù)和紅外光譜數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)保持一致。通過(guò)對(duì)Cp4An和COT2An(An=U，Pu)的分子軌道能級(jí)圖分析發(fā)現(xiàn)，采用最高的RSC ECP贗勢(shì)計(jì)算COT2U和Cp4U的基態(tài)分別為三重fφ2和fσ2組態(tài);而COT2Pu和Cp4Pu的基態(tài)分別為五重fσ1fπ1fφ2和fσ3fδ1組態(tài)。

金屬有機(jī)物；鈾；钚；密度泛函理論；溶劑化效應(yīng)

0 引言

錒系元素U和Pu在元素周期表中的地位特殊，U和Pu的外層價(jià)電子組態(tài)分別是5f36d17s2和5f67s2，錒系元素特殊的電子組態(tài)，表現(xiàn)在它們和配體配位時(shí)情況變得比一般的過(guò)渡金屬更加復(fù)雜；其次，由于它們是重原子，核外電子存在相對(duì)論效應(yīng)，量子化學(xué)方法處理起來(lái)也較為困難。在地下水中，由于錒系離子可以與水中的多種無(wú)機(jī)配體形成配合物的形式進(jìn)行遷移而倍受實(shí)驗(yàn)和理論工作者的關(guān)注[1-6]。本課題組先前應(yīng)用密度泛函理論研究了鈾與多種無(wú)機(jī)配體配位的性質(zhì)，并取得一定的成果[7-9]。錒系元素也能與有機(jī)配體形成特殊的結(jié)構(gòu)，例如錒系離子鈾環(huán)辛四烯(COT2-)和環(huán)戊二烯(Cp-)的金屬鈾有機(jī)配合物具有類似二茂鐵的夾心結(jié)構(gòu)。1986年合成出的二環(huán)辛四烯鈾(COT2U)引起理論和實(shí)驗(yàn)科學(xué)家的重視。理論上主要是針對(duì)氣相下該類有機(jī)配合物的結(jié)構(gòu)進(jìn)行研究，而關(guān)于液相下的結(jié)構(gòu)研究很少[10-24]。COT2Pu和COT2U是在四氫呋喃溶液中，通過(guò)將固體[(C2H5)4N2]PuCl5和Py2UCl6加入到K2COT的四氫呋喃溶液中合成出來(lái)[25-27]。Cp4U則是通過(guò)UC14與KC5H5在苯溶液中反應(yīng)制備得到的[28-31]。Fischer等[28]發(fā)現(xiàn)合成得到的紅色晶體狀的Cp4U在空氣中能夠穩(wěn)定存在，但在苯溶液中的結(jié)構(gòu)受空氣影響而不穩(wěn)定。實(shí)驗(yàn)上也得到了Cp4U的X射線晶體結(jié)構(gòu)數(shù)據(jù)[29-30]。因此考慮溶劑化效應(yīng)對(duì)結(jié)構(gòu)穩(wěn)定性的影響具有重要的意義。

理論上應(yīng)用相對(duì)論lcao Hartree-Fock-Slater (HFS)方法計(jì)算COT2M(M=Th，Pa，U，Pu)結(jié)構(gòu)中金屬與環(huán)的鍵距離分別為0.208、0.202、0.198、0.197和0.196 nm[32]。采用多組態(tài)自洽場(chǎng)MCSCF方法計(jì)算COT2U得到鈾與環(huán)的鍵距離為0.207 nm[33]。Moritiz和Dlog[18]采用將5f電子凍芯的贗勢(shì)energy-consistent PPs和核極化贗勢(shì)CPPs研究了COT2An(An=Th-Pu)的結(jié)構(gòu)和結(jié)合能等性質(zhì)，并比較MP2，CCDC以及密度泛函理論B3LYP和PW91不同方法優(yōu)化得到的結(jié)構(gòu)參數(shù)。該計(jì)算采用的贗勢(shì)將5f電子當(dāng)做內(nèi)層電子進(jìn)行凍芯處理，并且B3LYP計(jì)算得到鈾與環(huán)的距離為0.202 5 nm[18]。Kerridge等[19]也通過(guò)CASPT2，在固定環(huán)結(jié)構(gòu)參數(shù)情況下局部?jī)?yōu)化COT2U(和COT2Pu(結(jié)構(gòu)得到金屬到COT環(huán)中心的距離分別為0.194 4和0.189 8 nm。以上方法得到COT2U的鈾與環(huán)的鍵距離都比實(shí)驗(yàn)值[34-35]0.192 4 nm大，并且都只針對(duì)氣相下的結(jié)構(gòu)進(jìn)行計(jì)算，并未考慮溶劑化效應(yīng)。理論上對(duì)Cp4U的研究主要是在20世紀(jì)70年代和80年代時(shí)期，采用的方法多為半經(jīng)驗(yàn)方法和早期的密度泛函理論[36]。隨后也很少通過(guò)更現(xiàn)代的方法進(jìn)行探討。

本文應(yīng)用密度泛函理論的B3LYP方法，通過(guò)采用不同的有效核勢(shì)(ECP)計(jì)算研究環(huán)辛四烯(COT2-)和環(huán)戊二烯(Cp-)的金屬鈾有機(jī)配合物(Cp4U和COT2U)和具有類似結(jié)構(gòu)的钚金屬有機(jī)化合物(Cp4Pu和COT2Pu)在氣相和四氫呋喃溶劑中的結(jié)構(gòu)和化學(xué)性質(zhì)。

1 計(jì)算方法

標(biāo)量相對(duì)論效應(yīng)方法主要是在非相對(duì)論量子化學(xué)的波函數(shù)的哈密頓算符中增加考慮相對(duì)論mass-velosity和Darwin兩項(xiàng)[36]。包含相對(duì)論效應(yīng)計(jì)算的高斯軟件主要是通過(guò)有效核勢(shì)(ECP)考慮相對(duì)論。用標(biāo)量相對(duì)論計(jì)算含鈾體系的基態(tài)性質(zhì)是行之有效的手段。本文采用Gaussian09軟件[37]的B3LYP泛函[38-39]計(jì)算鈾體系的基態(tài)性質(zhì)。EMSL數(shù)據(jù)庫(kù)[40]中有CRENBL ECP，“大核芯”有效贗勢(shì)(RLC ECP)和“小核芯”有效贗勢(shì)(RSC ECP)的贗勢(shì)基組[41-43]。對(duì)U和Pu采用3種不同的贗勢(shì)基組，其它原子采用6-31+g(d)全電子基組進(jìn)行結(jié)構(gòu)優(yōu)化。CRENBL ECP和RLC ECP對(duì)U和Pu外層14或16個(gè)電子當(dāng)成價(jià)電子處理，把內(nèi)層78個(gè)電子凍芯處理。RSC-ECP則對(duì)對(duì)U和Pu外層32或34個(gè)電子當(dāng)成價(jià)電子處理，把內(nèi)層60個(gè)電子凍芯處理。計(jì)算中溶劑化效應(yīng)采用CPCM模型。采用聯(lián)合原子拓?fù)淠Ｐ?UATM)構(gòu)造空穴，原子半徑采用UA0(simple united atom topological model)。四氫呋喃(THF)的介電常數(shù)采用ε=7.42。

2 結(jié)果與討論

2.1 COT2An和Cp4An(An=U，Pu)基態(tài)的確定

如圖1，環(huán)辛四烯(COT2-)和環(huán)戊二烯(Cp-)的金屬鈾和钚有機(jī)配合物的結(jié)構(gòu)中金屬與配體成鍵形成夾心結(jié)構(gòu)。Hayes和Edelstein[44]通過(guò)基礎(chǔ)分子軌道計(jì)算預(yù)測(cè)COT2U和COT2Pu的基態(tài)結(jié)構(gòu)得到|Jz|= 4，0，并認(rèn)為COT2U是磁性化合物而COT2Pu為反磁性化合物。為確定化合物的基態(tài)，計(jì)算考慮了多種自旋態(tài)的優(yōu)化，結(jié)果如表1所示。COT2U和COT2Pu配合物中金屬的化合價(jià)為+4價(jià)，因?yàn)殁欕x子電子結(jié)構(gòu)有2個(gè)5f電子，而钚離子電子結(jié)構(gòu)中有4個(gè)5f電子。因此計(jì)算分別考慮了單重態(tài)、三重態(tài)的COT2U和單重態(tài)、三重態(tài)、五重態(tài)的COT2Pu。根據(jù)計(jì)算得到的體系能量大小可知，三重態(tài)COT2U和五重態(tài)COT2Pu屬于基態(tài)結(jié)構(gòu)。能量上，其它態(tài)的體系能量比基態(tài)體系能量高出140～500 kJ·mol-1。Kerridge等[19]也通過(guò)CASPT2優(yōu)化金屬-COT環(huán)中心的距離計(jì)算表明旋-軌耦合對(duì)COT2U和COT2Pu結(jié)構(gòu)影響很小。由于COT2Pu的理論計(jì)算表明基態(tài)為五重態(tài)，而實(shí)際上卻是反磁性物質(zhì)的特征，在理論研究中也倍受爭(zhēng)議[26,44]。因此在計(jì)算中有必要考慮單重態(tài)的COT2Pu的結(jié)構(gòu)。

表1 氣相下B3LYP/CRENBL ECP方法計(jì)算不同自旋態(tài)的COT2An和Cp4An(An=U,Pu)結(jié)構(gòu)的優(yōu)化結(jié)果(nm)和相對(duì)能量(kJ·mol-1)Table1 Optimized geometry parameters(nm)and relative energies(kJ·mol-1)of COT2An and Cp4An with An=U,Pufor different spin states by B3LYP/CRENBL ECP in gas phase

表1中，與相關(guān)的實(shí)驗(yàn)數(shù)據(jù)比較，三重態(tài)的COT2U和Cp4U中U-ligand的距離比實(shí)驗(yàn)值分別長(zhǎng)0.013 6和0.016 6 nm。低自旋態(tài)的結(jié)構(gòu)參數(shù)卻與實(shí)驗(yàn)值更為接近，這可能是由于CRENBL ECP贗勢(shì)高估了鍵長(zhǎng)。因此有必要進(jìn)行其它贗勢(shì)基組計(jì)算。

2.2 不同的相對(duì)論贗勢(shì)基組對(duì)幾何結(jié)構(gòu)的影響

圖1 環(huán)辛四烯(COT2-)和環(huán)戊二烯(Cp-)的金屬鈾和钚有機(jī)化合物的結(jié)構(gòu)(Cp4An和COT2An)Fig.1 Structures of cyclopentadienyl(Cp-)and annulene(COT2-)sandwiched uranium and plutonium compounds(Cp4An and COT2An)

表2 氣相和THF液相下不同贗勢(shì)基組計(jì)算COT2An(和Cp4An(An=U,Pu)基態(tài)結(jié)構(gòu)的優(yōu)化結(jié)果(nm)Table2 Optimized geometry parameters(nm)for ground states of COT2Anand Cp4An with An=U,Puby different ECPs in gas and THF solution phase

表2中CRENBL ECP，RLC ECP和RSC ECP贗勢(shì)基組的優(yōu)化結(jié)果存在一定程度的差異。CRENBL和RLC ECPs計(jì)算得到的鍵長(zhǎng)比RSC ECP計(jì)算的鍵長(zhǎng)結(jié)果長(zhǎng)。對(duì)含鈾和钚的體系采用小核芯的RSC ECP贗勢(shì)基組可以得到與實(shí)驗(yàn)值[34-35]更為接近的結(jié)果。含钚體系中CRENBL和RLC ECPs計(jì)算得到的結(jié)構(gòu)參數(shù)中鍵長(zhǎng)比RSC ECP的計(jì)算結(jié)果約長(zhǎng)0.01～0.02 nm。采用RSC ECP贗勢(shì)基組計(jì)算得到COT2U中金屬到環(huán)中心距離U-ligand為0.194 3 nm與實(shí)驗(yàn)數(shù)據(jù)0.192 4 nm一致[34-35]。同時(shí)采用RSC ECP贗勢(shì)基組計(jì)算得到Cp4U中金屬到環(huán)中心距離U-ligand為0.260 4 nm只比實(shí)驗(yàn)值[30]長(zhǎng)0.007 6 nm。關(guān)于COT2Pu確切的結(jié)構(gòu)參數(shù)未見(jiàn)實(shí)驗(yàn)報(bào)道，但有報(bào)道通過(guò)已知的錒系金屬環(huán)狀有機(jī)配合物的鍵長(zhǎng)結(jié)果推算出Pu-ligand的數(shù)值為0.183nm[19,45]。不同自旋態(tài)COT2Pu的結(jié)構(gòu)參數(shù)如表2所示，采用RSC ECP贗勢(shì)基組得到相對(duì)穩(wěn)定的五重態(tài)的Pu-ligand的數(shù)值為0.197 4 nm，而不穩(wěn)定的單重態(tài)的Pu-ligand的數(shù)值為0.193 4 nm。

在THF溶液中，An-ligand和An-C鍵受到溶劑化效應(yīng)的影響都比氣相中的結(jié)果長(zhǎng)。其中COT2U的U-ligand增長(zhǎng)了0.009 2 nm，五重態(tài)COT2Pu的Pu-ligand增長(zhǎng)了0.022 9 nm，但是單重態(tài)COT2Pu在兩相中相應(yīng)鍵長(zhǎng)相似。環(huán)戊二烯化合物Cp4U(和Cp4Pu的結(jié)構(gòu)參數(shù)受THF溶劑的影響也很小，但溶劑化效應(yīng)對(duì)COT2U和COT2Pu的影響很大。Berthet等[46]合成了彎曲的環(huán)辛四烯鈾的衍生物說(shuō)明，被2個(gè)COT2-環(huán)夾心的鈾在溶液中可能受到環(huán)境中其它配體的作用而使兩環(huán)成一定夾角而彎曲。

2.3 結(jié)合能計(jì)算

COT2An和Cp4An(An=U，Pu)中金屬與環(huán)之間的結(jié)合能如表3所示。結(jié)合能采用下面的公式計(jì)算：

其中Eb、Ecomplex、Eactinideion和Ecyclo-ligand分別表示結(jié)合能、金屬有機(jī)化合物體系的能量、錒系金屬離子的能量和環(huán)狀有機(jī)配體的能量。所有能量值都通過(guò)零點(diǎn)能(ZPE)校正。n表示環(huán)狀有機(jī)配體的數(shù)目。

表3中所采用的不同ECPs贗勢(shì)基組都表明COT2An和Cp4An(An=U，Pu)中金屬與環(huán)之間的結(jié)合能很大，進(jìn)而說(shuō)明在氣相和THF液相中環(huán)狀有機(jī)配體與錒系金屬離子之間有很強(qiáng)的成鍵作用。含COT2-配合物的結(jié)合能比相應(yīng)的含Cp-配合物的能量高。Moritiz和Dlog[18]研究發(fā)現(xiàn)氣相中COT2An錒系金屬有機(jī)配合物的結(jié)合能隨中心金屬核電荷數(shù)的增加而增大。表3中含钚的金屬有機(jī)配合物也比相應(yīng)的含鈾有機(jī)配合物的結(jié)合能大。CRENBL和RLC ECPs計(jì)算的結(jié)果表明，氣相中金屬钚有機(jī)配合物比類似結(jié)構(gòu)的金屬鈾有機(jī)配合物穩(wěn)定，其結(jié)合能大1 eV左右。金屬钚有機(jī)配合物的穩(wěn)定性在THF溶液中尤為突出，其結(jié)合能比類似結(jié)構(gòu)的金屬鈾配合物大4 eV左右。但是采用RSC-ECP的計(jì)算卻表明不管在氣相還是THF溶劑中金屬钚有機(jī)配合物的結(jié)合能比類似結(jié)構(gòu)的金屬鈾有機(jī)配合物的結(jié)合能大4 eV。溶劑化效應(yīng)使得配合物的結(jié)合能相對(duì)氣相中的結(jié)合能下降了接近70%。實(shí)驗(yàn)上合成得到的紅色晶體狀的Cp4U[28]在空氣中能夠穩(wěn)定存在，但在苯溶液中的結(jié)構(gòu)卻不穩(wěn)定，也說(shuō)明了溶劑化對(duì)結(jié)構(gòu)穩(wěn)定性的影響很大。從計(jì)算給出的結(jié)果可以說(shuō)明，該類化合物在氣相當(dāng)中比較穩(wěn)定，保存的環(huán)境考慮為氣相。使用不同的相對(duì)論贗勢(shì)對(duì)結(jié)合能的數(shù)值影響比較大。氣相中RSC ECP的結(jié)果顯示COT2Pu五重態(tài)的結(jié)合能為95.9 eV，而其單重態(tài)的結(jié)合能為96.3 eV。而考慮THF溶劑化效應(yīng)后，COT2Pu五重態(tài)的結(jié)合能反而比其單重態(tài)的結(jié)合能高。Cp4Pu配合物的結(jié)合能也有同樣的趨勢(shì)。這說(shuō)明溶劑可以使Cp4Pu和COT2Pu的基態(tài)發(fā)生變化。

表3 氣相和THF液相下不同贗勢(shì)基組計(jì)算得到的COT2An和Cp4An(An=U,Pu)結(jié)合能(eV)Table3 Actinide-ring binding energies(eV)for Cp4An and COT2Anwith An=U,Puby different ECPs in gas and THF solution phase

2.4 紅外光譜性質(zhì)

采用最佳的RSC ECP贗勢(shì)基組得到的紅外光譜模擬圖如圖2所示。通過(guò)計(jì)算得到氣相中的COT2U的紅外吸收鋒有724、923 cm-1(medium)和3 183 cm-1(strong)；THF液相中其紅外吸收峰有712 cm-1(medium)、918 cm-1(medium)和3 177 cm-1(medium)。而在氣相和THF液相中COT2U的理論紅外吸收分別在1 556 cm-1(strong)和1 518 cm-1(strong)處有峰，該峰是C8H8環(huán)中C-H的上下振動(dòng)頻率。氣相中COT2U、Cp4U、COT2Pu和Cp4Pu分別在724 cm-1(strong)、788 cm-1(medium)、702 cm-1(strong)和772 cm-1處的紅外吸收光譜屬于有機(jī)環(huán)的伸縮振動(dòng)特征光譜。而結(jié)構(gòu)中環(huán)中C-H鍵紅外特征吸收光譜出現(xiàn)在1 000～1 500 cm-1波數(shù)范圍，同時(shí)在THF溶劑中配合物的紅外光譜吸收峰具有更高的強(qiáng)度。由此可見(jiàn)，理論計(jì)算預(yù)算得到的環(huán)辛四烯(COT2-)和環(huán)戊二烯(Cp-)配體在在800 cm-1波數(shù)以下的指紋區(qū)域有強(qiáng)的有機(jī)環(huán)的伸縮振動(dòng)特征吸收峰。實(shí)驗(yàn)測(cè)得COT2U的紅外光譜吸收峰有：594 cm-1(strong)、741 cm-1(medium)、772 cm-1、787 cm-1、900 cm-1(medium)和3 000 cm-1[47]。結(jié)果表明理論紅外光譜計(jì)算結(jié)果與實(shí)驗(yàn)側(cè)得的數(shù)據(jù)吻合，也進(jìn)一步說(shuō)明該化合物的結(jié)構(gòu)優(yōu)化結(jié)果合理。

2.5 分子軌道

圖2 氣相和THF液相下COT2U,COT2Pu,Cp4U和Cp4Pu的紅外光譜模擬Fig.2 IR absorption spectra of COT2U,COT2Pu,Cp4U,and Cp4Pu in gas phase and THF solution

RSC ECP計(jì)算得到An(Cp)4和An(COT)2(An= U，Pu)的分子級(jí)軌道能級(jí)圖(圖3)中能量低于-8 eV能級(jí)的電子受自旋影響小，可明顯觀測(cè)到能量相等α電子與β電子配對(duì)存在于分子軌道能級(jí)中。由于錒系元素電子排布的復(fù)雜性，很難準(zhǔn)確確定其基態(tài)電子結(jié)構(gòu)。通過(guò)采用RSC贗勢(shì)的計(jì)算進(jìn)行預(yù)測(cè)，配合物An(COT)2的分子軌道能級(jí)圖顯示含5fφ電子的軌道不與COT環(huán)的軌道作用，5f電子主要分布在能量范圍在-4～-8 eV區(qū)域的最高占據(jù)軌道上。配合物U(COT)2的H-0和H-1軌道屬于5fφ電子軌道，而配合物Pu(COT)2中H-3，H-4，H-5，H-6軌道屬于5fφ，5fφ，5fπ和5fσ電子軌道。由此可見(jiàn)U4+和Pu4+離子中5f電子在形成配合物時(shí)部分5f電子不與環(huán)配體電子作用，而單獨(dú)占據(jù)電子軌道。含Pu有機(jī)配合物的結(jié)合能大小可能與軌道中的5f單電子有關(guān)。同時(shí)從圖3中也可以清晰地看到配合物An(COT)2的成鍵軌道H-2，H-3(U(COT)2)和H-0(Pu (COT)2)存在環(huán)配體與金屬5f電子相互作用形成的成鍵軌道。該成鍵軌道是通過(guò)配體提供電子對(duì)，金屬供空的fδ軌道而形成的。與2個(gè)COT環(huán)形成成鍵軌道還包含了部分空的6dδ軌道，比如H-4(U(COT)2)和H-1，H-2(Pu(COT)2)軌道。金屬U或Pu與COT環(huán)的成鍵作用，比其與Cp環(huán)的成鍵作用強(qiáng)。根據(jù)圖3的分子軌道能級(jí)圖可判斷U(COT)2是fφ2組成的三重態(tài)結(jié)構(gòu)，而Pu(COT)2則是fφ2fπ1fσ1組成的五重態(tài)結(jié)構(gòu)。King[48]應(yīng)用拓?fù)滠壍览碚擃A(yù)測(cè)An(Cp)4和An(COT)2(An=Th，Pa，U，Np)金屬有機(jī)化合物的5f價(jià)電子首先占據(jù)5fφ軌道。本章節(jié)對(duì)U(COT)2的軌道能級(jí)計(jì)算與King的結(jié)論一致[48]。但對(duì)于Pu(COT)2化合物f電子只有2個(gè)電子占據(jù)5fφ軌道。由于Pu電子結(jié)構(gòu)的復(fù)雜性，其配合物的基態(tài)結(jié)構(gòu)還有待實(shí)驗(yàn)上進(jìn)一步驗(yàn)證。配合物An(Cp)4中，配位鍵主要是由5fφ空軌道與Cp環(huán)的2p軌道作用形成。U(Cp)4是三重態(tài)的fσ2結(jié)構(gòu)，Pu(Cp)4則為五重態(tài)的fσ3fδ1結(jié)構(gòu)。

圖3 RSC ECP計(jì)算得到的An(Cp)4和An(COT)2(An=U,Pu)分子軌道能級(jí)圖(部分前線軌道圖顯示在能級(jí)圖上,其中H代表HOMO軌道,L代表LUMO軌道)Fig.3 MO energy level diagrams for An(Cp)4and An(COT)2with An=U,Pucalculated with RSC ECP(Partial frontier orbital are shown in energy level diagrams with H and L shorted for the highest occupied and lowest unoccupied orbital respectively)

3 結(jié)論

本文采用Gaussian09中的密度泛函中的B3LYP泛函研究配合物COT2U、COT2Pu、Cp4U和Cp4Pu的幾何結(jié)構(gòu)和電子結(jié)構(gòu)。優(yōu)化結(jié)果表明對(duì)含鈾和含钚的體系進(jìn)行計(jì)算需采用小核凍芯的RSC相對(duì)論有效贗勢(shì)。計(jì)算得到該類有機(jī)物的IR光譜與實(shí)驗(yàn)值一致，其中紅外特征吸收峰在500～1 000 cm-1區(qū)域。結(jié)合能計(jì)算表明配合物的穩(wěn)定性呈COT2Pu＞COT2U＞Cp4Pu＞Cp4U的趨勢(shì)，因此含U和Pu的金屬有機(jī)配合物主要是通過(guò)配位鍵而穩(wěn)定存在的。配位鍵主要是由5f空軌道與COT或Cp環(huán)的2p占據(jù)軌道作用形成。分子軌道能級(jí)圖可判斷U (COT)2是組成的三重態(tài)結(jié)構(gòu)，而Pu(COT)2則是組成的五重態(tài)結(jié)構(gòu)，U(Cp)4是三重態(tài)的結(jié)構(gòu)，Pu(Cp)4則為五重態(tài)的結(jié)構(gòu)。

[1]Clark D L,Hobart D E,Neu M P.Chem.Rev.,1995,95(1): 25-48

[2]Craw J S,Vincent M A,Hillier I H.J.Phys.Chem.A,1995, 99:10181-10185

[3]Hemmingsen L,Amara P,Ansoborlo E,et al.J.Phys.Chem. A,2000,104(17):4095-4101

[4]Privalov T,Schimmelpfennig B,Wahlgren U,et al.J.Phys. Chem.A,2002,106(46):11277-11282

[5]de Jong W A,Aprà E,Windus T L,et al.J.Phys.Chem.A, 2005,109(50):11568-11577

[6]Denecke M A.Coord.Chem.Rev.,2006,250(7/8):730-754

[7]GU Jia-Fang(辜家芳),LU Chun-Hai(陸春海),CHEN Wen-Kai(陳文凱),et al.Acta Phys.-Chim.Sin.(物理化學(xué)學(xué)報(bào)), 2009,25(4):655-660

[8]GU Jia-Fang(辜家芳),LU Chun-Hai(陸春海),CHEN Wen-Kai(陳文凱),et al.Acta Phys.-Chim.Sin.(物理化學(xué)學(xué)報(bào)), 2012,28(4):792-798

[9]GU Jia-Fang(辜家芳),MAN Mei-Ling(滿梅玲),LU Chun-Hai(陸春海),et al.Chinese J.Inorg.Chem.(無(wú)機(jī)化學(xué)學(xué)報(bào)), 2012,28(7):1324-1332

[10]Tatsumi K,Nakamura A.J.Organomet.Chem.,1984,272(2): 141-154

[11]Tatsumi K,Nakamura A.J.Am.Chem.Soc.,1987,109(11): 3195-3206

[12]Chang A H H,Pitzer R M.J.Am.Chem.Soc.,1989,111(7): 2500-2507

[13]Gourier D,Caurant D,Berthet J C,et al.Inorg.Chem.,1997, 36(25):5931-5936

[14]Li J,Bruce E B.J.Am.Chem.Soc.,1997,119:9021-9032

[15]Li J,Bursten B E.J.Am.Chem.Soc.,1998,120(44):11456-11466

[16]Hager J S,Zahardis J,Pagni R M,et al.J.Chem.Phys., 2004,120(6):2708-2718

[17]BenYahia M,Belkhiri L,Boucekkine A.J.Mol.Struct., 2006,777(1/2/3):61-73

[18]Moritz A,Dolg M.Chem.Phys.,2007,337(1/2/3):48-54

[19]Kerridge A,Kaltsoyannis N.J.Phys.Chem.A,2009,113 (30):8737-8745

[20]Streitwieser A.J.Org.Chem.,2009,74(12):4433-4446

[21]Elkechai A,Meskaldji S,Boucekkine A,et al.J.Mol.Struct., 2010,954(1/2/3):115-123

[22]Yahia A,Castro L,Maron L.Chem.Eur.J.,2010,16(19): 5564-5567

[23]Zhou J,Sonnenberg J L,Schlegel H B.Inorg.Chem.,2010, 49(14):6545-6551

[24]Castro L,Yahia A,Maron L.C.R.Chim.,2010,13(6/7):870-875

[25]Streitwieser A,Mueller-Westerhoff U.J.Am.Chem.Soc., 1968,90(26):7364-7364

[26]Karraker D G,Stone J A,Jones E R,et al.J.Am.Chem. Soc.,1970,92(16):4841-4845

[27]Streitwieser J A,Muller-Westerhoff U,Sonnichsen G,et al. J.Am.Chem.Soc.,1972,95(26):8644-8649

[28]FischerVEO,HristiduY.Z.Naturforsch.,1962,17b:275-276 [29]Burns J H.J.Organomet.Chem.,1974,69(2):225-233

[30]Burns J H.J.Am.Chem.Soc.,1973,95:3815-3817

[31]Kanellakopulos B,Maier R,Heuser J.J.Alloys Compd., 1991,176(1):89-96

[32]Boerrigter P M,Baerends E J,Snijders J G.Chem.Phys., 1988,122(3):357-374

[33]Liu W,Dolg M,Fulde P.J.Chem.Phys.,1997,107(9):3584-3591

[34]Zalkin A,Raymond K N.J.Am.Chem.Soc.,1969,91:5667-5668

[35]Avdeef A,Raymond K N,Hodgson K O,et al.Inorg.Chem., 1972,11:1083-1088

[36]Kaltsoyannis N,Hay P J,Li J,et al.Chemistry of the Actinide and Transactinide Elements.Morss L R,Edelstein N M, Fuger J Eds.,Netherlands:Springer,2006:1893-2012

[37]Frisch M J,Trucks G W,Schlegel H B,et al.Gaussian 09, Wallingford CT:Gaussian,Inc.,2009.

[38]Becke A D.Phys.Rev.A,1988,38:3098-3100

[39]Bardin N,Rubini P,Madic C.Radiochim.Acta,1998,83: 189-194

[40]Gaussian Basis Set Order Form:http://www.emsl.pnl.gov/ forms/basisform.html

[41]Hay P J.J.Chem.Phys.,1983,79:5469-5482

[42]Kahn L R,Hay P J,Cowan R D.J.Chem.Phys.,1978,68: 2386-2397

[43]Feller D.J.Comput.Chem.,1996,17(13):1571-1586

[44]Hays R G,Edelstein N M.J.Am.Chem.Soc.,1972,94(25): 8688-8691

[45]Shannon R.Acta Crystallogr.Sect.A,1976,32(5):751-767

[46]Berthet J C,Thuéry P,Ephritikhine M.Organometallics, 2008,27(8):1664-1666

[47]Weakliem P C,Carter E A.J.Chem.Phys.,1992,96:3240-3250

[48]King R B.Inorg.Chem.,1992,31(10):1978-1980

DFT Studies on Structure and Spectral Properties of Organoactinide Complexes of Cp4An and COT2An(Cp-=C5H5-,COT2-=C8H82-,An=U,Pu)

GU Jia-Fang1XU Ke2CHEN Wen-Kai＊,3
(1Department of Chemical Engineering,Fuzhou University Zhicheng College,Fuzhou 350002,China)
(2Minjiang College,Fuzhou 350108,China)
(3Department of Chemistry,Fuzhou University,Fuzhou 350108,China)

Properties of Cp4An and COT2An(Cp-=C5H5-,COT2-=C8H82-,An=U,Pu)species both in gas phase and tetrahydrofuran(THF)solution are systematically studied by DFT(density functional theory).The solvent effect of THF is simulated by a conductor-like polarized continuum model(CPCM).Theoretical calculations show that the binding energies of the complexes are in the order of COT2Pu＞COT2U＞Cp4Pu＞Cp4U.Solvent effects decreased binding energies of the structures.The optimized geometry structures and the IR spectra for the complexes in gas and aqueous phase are in accordance to the available experimental data very well.Detail studies on MO energy levels for Cp4An and COT2An with An=U,Pucalculated with rsc ecp indicate that COT2U and CP4U favors triplets fφ2and fσ2configuration respectively,meanwhile,COT2Pu and Cp4Pu favor quintets fσ1fπ1fφ2and fσ3fδ1configuration respectively.

organometallic;uranium;plutonium;density functional theory;solvent effect

O641

A

1001-4861(2017)09-1579-08

10.11862/CJIC.2017.179

2017-03-21。收修改稿日期：2017-06-14。

國(guó)家自然科學(xué)基金(No.10676007)資助項(xiàng)目和福建省教育廳科研基金(No.JB14222)資助項(xiàng)目。＊

。E-mail：qc2008@fzu.edu.cn，Tel：+86591-22866162