• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    CH2BrCl在265 nm光解產(chǎn)生的氯甲基自由基的振動(dòng)內(nèi)能

    2012-11-30 10:48:48唐碧峰
    物理化學(xué)學(xué)報(bào) 2012年5期
    關(guān)鍵詞:內(nèi)能角動(dòng)量母體

    張 秀 吳 東 唐碧峰,*

    (1湖北工程學(xué)院物理系,湖北孝感432000;2華中師范大學(xué)物理科學(xué)與技術(shù)學(xué)院,武漢430079)

    CH2BrCl在265 nm光解產(chǎn)生的氯甲基自由基的振動(dòng)內(nèi)能

    張 秀1吳 東2唐碧峰1,*

    (1湖北工程學(xué)院物理系,湖北孝感432000;2華中師范大學(xué)物理科學(xué)與技術(shù)學(xué)院,武漢430079)

    利用離子速度影像技術(shù)研究了CH2BrCl在265 nm附近的激光光解.利用2+1共振增強(qiáng)多光子電離分別獲得光解產(chǎn)物Br(2P1/2)和Br(2P3/2)的離子速度圖像,從而得出Br(2P1/2)和Br(2P3/2)的速度分布,以及光解碎片的總平動(dòng)能分布.據(jù)此,運(yùn)用角動(dòng)量守恒碰撞模型獲得了解離氯甲基自由基(·CH2Cl)的振動(dòng)內(nèi)能分布.研究結(jié)果表明:CH2BrCl+hv→Br(2P1/2)+CH2Cl通道產(chǎn)生的氯甲基自由基中被激發(fā)的振動(dòng)模主要是v4、v3+v4、v2+v4和v2+v6;CH2BrCl+hv→Br(2P3/2)+CH2Cl通道產(chǎn)生的氯甲基自由基中被激發(fā)的振動(dòng)模主要是v2+v6、v1+v3、v2+v5、v2+v3+v5和v1+v5;母體分子CH2BrCl在吸收光解光子后除有v5(CBr stretch)振動(dòng)模被激發(fā)外,還有v7(CH2a-stretch)等其它振動(dòng)模也被激發(fā).

    CH2BrCl;光解;離子速度影像;角動(dòng)量守恒碰撞模型;振動(dòng)內(nèi)能分布

    1 引言

    溴烷烴是一類具有特殊氣味并對(duì)人體有害的物質(zhì),許多水源受到化工、醫(yī)藥等廢水污染會(huì)產(chǎn)生溴烷烴.溴烷烴及其它含溴化合物的光解受到廣泛關(guān)注,這一方面是因?yàn)槲覀円呀?jīng)了解到溴原子對(duì)地球大氣臭氧層有很大的破壞力.1-3另一方面是因?yàn)殇逋闊N的光解對(duì)光解動(dòng)力學(xué)研究具有重要的理論意義.4-13

    近年來隨著各種實(shí)驗(yàn)技術(shù)的發(fā)展,光解動(dòng)力學(xué)的研究越來越關(guān)注態(tài)-態(tài)層面上的細(xì)節(jié).1997年P(guān)arker等14發(fā)展了離子速度影像(ion velocity map imaging)技術(shù),不僅使速度分辨率得到極大的提高,而且提高了離子通過率.Liu15和Suits16等發(fā)展了時(shí)間切片離子速度成像技術(shù)(time-sliced velocity-mapped ion imaging),使得在振動(dòng)態(tài)層面上研究分子光解動(dòng)力學(xué)成為可能.應(yīng)用這些新技術(shù)所進(jìn)行的研究17-19表明,解離自由基的振動(dòng)激發(fā)及振動(dòng)內(nèi)能的分布反映著解離機(jī)理的細(xì)節(jié).所以,解離自由基的振動(dòng)激發(fā)和振動(dòng)內(nèi)能分布研究對(duì)深入研究光解動(dòng)力學(xué)有重要意義.

    Butler研究小組20利用離子速度影像方法研究了2-溴乙醇(BrCH2CH2OH)在193 nm的光解.2-溴乙醇在193 nm的光解通道為BrCH2CH2OH→C2H4OH+Br(2P3/2)或C2H4OH+Br(2P1/2),而初生的2-羥乙基自由基(·C2H4OH)可能進(jìn)一步發(fā)生分解反應(yīng)C2H4OH→C2H4+OH,產(chǎn)生乙烯(C2H4)和羥自由基(·OH).因此,2-溴乙醇光解產(chǎn)生的2-羥乙基自由基是不穩(wěn)定的,只有振動(dòng)內(nèi)能較小、不能克服分解反應(yīng)勢(shì)壘的2-羥乙基自由基才能最后形成穩(wěn)定的產(chǎn)物.Butler研究小組20-23提出了一個(gè)精確的碰撞模型,該碰撞模型應(yīng)用了角動(dòng)量守恒、動(dòng)量守恒和能量守恒定律,我們稱之為角動(dòng)量守恒碰撞模型(impulsive model invoking angular momentum conservation). Butler研究小組20利用角動(dòng)量守恒碰撞模型得出了穩(wěn)定2-羥乙基自由基的平動(dòng)能分布,與實(shí)驗(yàn)測量的結(jié)果一致.能得到正確結(jié)果的關(guān)鍵因素是角動(dòng)量守恒碰撞模型準(zhǔn)確計(jì)算出了2-羥乙基自由基的振動(dòng)內(nèi)能分布.由此可見,如果光解產(chǎn)生的自由基是中間體,則其振動(dòng)內(nèi)能分布決定著其下一步發(fā)生化學(xué)反應(yīng)的方向和限度.

    CH2BrCl擁有兩個(gè)不同的鹵素原子.相比只有單個(gè)鹵素原子的鹵代甲烷,它有更加復(fù)雜的紫外(UV)吸收譜和更多的光解通道.其光解通道CH2BrCl+ hv→CH2Cl+Br(2P3/2,1/2)(下面以Br代表Br(2P3/2),以Br*代表Br(2P1/2))產(chǎn)生的氯甲基自由基(·CH2Cl)進(jìn)一步解離的勢(shì)壘較高,相對(duì)穩(wěn)定,是研究自由基碎片振動(dòng)內(nèi)能的較好范例,從而引起人們的興趣.Tzeng等24利用平動(dòng)能光譜學(xué)(PTS)研究了氯溴甲烷(CH2BrCl)分子在248和193 nm的光解.發(fā)現(xiàn)在這兩波長下有四個(gè)可能的光解通道起作用:

    他們發(fā)現(xiàn)在193 nm波長下C-Br鍵解離和C-Cl鍵解離的分支比是4.5:1;在248 nm波長下,只有C-Br鍵的解離.通過與CH3Br分子的吸收譜比較, CH2BrCl分子的A吸收帶中心位于203 nm處,范圍大約為3.0×104-5.2×104cm-1(333-192 nm),被歸屬到n(Br)→σ*(C-Br)躍遷.CH2BrCl分子的n(Cl)→σ*(C-Cl)躍遷對(duì)應(yīng)的B吸收帶其中心預(yù)期在164 nm,25因?yàn)镃F2BrCl分子的n(Cl)→σ*(C-Cl)躍遷對(duì)應(yīng)的吸收帶中心就在這一波長.26B吸收帶帶頭能量大約為5.1×104cm-1(196 nm),25所以,在論文描述的實(shí)驗(yàn)中不太可能發(fā)生反應(yīng)式(3)和(4)所示的反應(yīng)產(chǎn)生氯原子.

    North等27,28利用共振增強(qiáng)多光子電離(REMPI)和飛行時(shí)間質(zhì)譜儀研究了CH2BrCl從193到268 nm的光解,這一波段覆蓋了n(Br)→σ*(C-Br)躍遷對(duì)應(yīng)的整個(gè)UV吸收譜.他們發(fā)現(xiàn)與CH3Br的光解相比,非絕熱效應(yīng)在CH2BrCl的光解過程中起很小的作用.他們認(rèn)為引起這個(gè)差別的原因可能是因?yàn)镃H2BrCl屬于Cs對(duì)稱性,且CH2BrCl解離后產(chǎn)生的氯甲基自由基的速度較低.他們得到的Br*和Br通道的總平動(dòng)能分布都可用單高斯函數(shù)擬合.Jung等29利用分辨率較高的離子成像技術(shù)研究了CH2BrCl在234 nm的光解.他們發(fā)現(xiàn)Br通道的總平動(dòng)能分布出現(xiàn)雙峰,反映該通道的氯甲基自由基碎片的內(nèi)能有兩個(gè)值,并認(rèn)為是由σ*(C-Br)勢(shì)能面和σ*(C-Cl)勢(shì)能面相互作用產(chǎn)生的.他們認(rèn)為快速成分對(duì)應(yīng)無振動(dòng)激發(fā)的硬碰撞模型,慢速成分對(duì)應(yīng)有振動(dòng)激發(fā)的軟碰撞模型,兩部分的能量差對(duì)應(yīng)CH2Cl的傘形振動(dòng)模.Ng等30利用時(shí)間切片離子速度影像技術(shù)研究了CH2BrCl在193.3-266.6 nm的光解,他們同樣發(fā)現(xiàn)Br*通道的總平動(dòng)能分布都可用單高斯函數(shù)擬合,而Br通道的總平動(dòng)能分布都需用兩個(gè)高斯函數(shù)擬合.他們認(rèn)為不同的高斯成分來源于不同的解離態(tài)勢(shì)能面,而沒有討論對(duì)應(yīng)的氯甲基自由基的振動(dòng)激發(fā).

    我們利用離子速度影像技術(shù),研究了CH2BrCl在265 nm的光解.通過分析Br*和Br的離子速度影像,我們得到Br*和Br通道對(duì)應(yīng)的總平動(dòng)能分布以及各向異性化參數(shù)(anisotropy parameter).然后,運(yùn)用角動(dòng)量守恒碰撞模型分別推測出Br*和Br通道對(duì)應(yīng)的氯甲基自由基的振動(dòng)內(nèi)能分布.并探究兩通道對(duì)應(yīng)的氯甲基自由基振動(dòng)模的激發(fā),分析引起不同振動(dòng)模激發(fā)的機(jī)理.

    2 實(shí)驗(yàn)

    實(shí)驗(yàn)是在一個(gè)經(jīng)過改裝的自制的飛行時(shí)間質(zhì)譜儀(TOF)上完成的.具體改裝情況我們以前報(bào)道過,4,7,8主要是加裝一個(gè)靜電離子透鏡.

    使用的氯溴甲烷(CH2BrCl)樣品是從上海試劑公司購買,純度為99.9%,未經(jīng)進(jìn)一步提純.以1.0× 105Pa的氦(He)作為種子氣體與氯溴甲烷(CH2BrCl)液體樣品蒸汽混合,通過噴嘴直徑為0.6 mm的脈沖閥(park,general valve)噴入真空室形成超聲分子束,脈沖閥溫度保持在30°C,由脈沖閥驅(qū)動(dòng)器(general valve,IOTA-1)驅(qū)動(dòng).超聲分子束運(yùn)動(dòng)大約6 cm后經(jīng)孔徑為1.0 mm的漏勺(skimmer)準(zhǔn)直后進(jìn)入反應(yīng)作用區(qū).在實(shí)驗(yàn)過程中,進(jìn)樣前束源室真空保持為4.0× 10-6Pa,電離室真空保持在5.0×10-6Pa以上;開脈沖閥進(jìn)樣后束源室真空介于1.0×10-4到3.1×10-4Pa之間,電離室的真空度高于1.0×10-4Pa.為了減少團(tuán)簇生成,我們對(duì)脈沖閥的脈寬、載氣壓強(qiáng)、樣品濃度以及激光與超聲分子束作用的位置等實(shí)驗(yàn)條件進(jìn)行了優(yōu)化.最后在質(zhì)譜中沒發(fā)現(xiàn)團(tuán)簇,實(shí)驗(yàn)就在此優(yōu)化條件下進(jìn)行.

    實(shí)驗(yàn)中使用的離子透鏡是在用SIMION 3D軟件31模擬優(yōu)化的基礎(chǔ)上設(shè)計(jì)的,其構(gòu)造有利于離子速度影像的生成.

    激光是由Nd:YAG(Quantel,YG980)激光器的三倍頻激光(355 nm)泵浦染料激光器(Lambda Physik,ScanMate 2E OG)產(chǎn)生.染料激光器輸出的基頻光再經(jīng)BBO晶體倍頻后輸出UV光.得到的UV激光強(qiáng)度大約為每脈沖400 μJ,脈寬為10 ns,倍頻前線寬為0.03 cm-1.染料激光輸出波長經(jīng)由Ne光電池的電流(Optogalvanic,OG)信號(hào)校準(zhǔn),也就是經(jīng)由已知波長校準(zhǔn),因?yàn)镺G信號(hào)反映的是氖的吸收譜,而這些譜線的波長是準(zhǔn)確已知的.在離子速度影像實(shí)驗(yàn)中要求有非常好的線偏振激光,且要求線偏振方向平行于探測器平面,實(shí)驗(yàn)中通過使用可變波片(Soleil-Bebubat)來實(shí)現(xiàn).UV光通過焦距f=20 cm的聚焦透鏡聚焦后與分子束垂直交叉.我們的實(shí)驗(yàn)是單色光實(shí)驗(yàn),就是用同一束激光解離分子,并通過(2+1)共振增強(qiáng)多光子電離(REMPI)來電離碎片, CH2BrCl→CH2Cl+Br*通道對(duì)應(yīng)的波長為264.86 nm,CH2BrCl→CH2Cl+Br通道對(duì)應(yīng)的波長為264.77 nm.我們通過降低光強(qiáng)盡量壓制高階多光子過程、飽和效應(yīng)和空間電荷效應(yīng).

    電離后的光解碎片被排斥場(排斥極與拉出極之間的電場)和拉出場(拉出極和地之間的電場)加速,并通過極板中心的圓孔拉出光解作用區(qū),然后,經(jīng)50 cm的自由飛行后到達(dá)一個(gè)直徑為40 mm的二維位置敏感探測器.位置敏感探測器由一個(gè)雙層的微通導(dǎo)板(MCP)和一個(gè)反應(yīng)速度很快的熒光屏(PS)組成.帶電碎片的質(zhì)量大小可以由其飛行到探測器的時(shí)間來確定.為了只對(duì)選定的光解碎片成像,我們?cè)谠摴饨馑槠竭_(dá)MCP的同時(shí)在MCP上加一脈沖電壓——質(zhì)量門.脈沖電壓加上前,MCP上的直流電壓低于其工作電壓;脈沖電壓加上后,MCP上的瞬間電壓達(dá)到工作電壓,對(duì)信號(hào)進(jìn)行倍增放大,產(chǎn)生電子流轟擊熒光屏,形成所選離子的速度影像.我們利用光電倍增管(PMT)接收來自熒光屏上的光信號(hào),將其送到數(shù)字化示波器(Tektronix,TDS 2012)顯示離子的飛行時(shí)間質(zhì)譜.熒光屏上離子速度影像的形狀和大小可以通過調(diào)節(jié)排斥極和拉出極上的電壓(VR,VE)來優(yōu)化.實(shí)驗(yàn)中的優(yōu)化條件是VE/ VR=0.71(VR~2200 V).我們使用一個(gè)市售的CCD (MS-2821E,Mintron)攝像機(jī)采集熒光屏上圖像,將圖像信號(hào)通過同軸電纜線輸入到一個(gè)配置了圖像采集卡的計(jì)算機(jī)上,并用圖像采集軟件(CODA32)進(jìn)行累加平均,累加平均次數(shù)對(duì)應(yīng)30000個(gè)激光脈沖.考慮到多普勒(Doppler)效應(yīng),在累加的過程中激光波長在±0.8 cm-1范圍內(nèi)掃描,以覆蓋離子碎片的所有速度分量.激光脈沖、分子束和探測系統(tǒng)的工作頻率為10 Hz,整個(gè)時(shí)序控制由一個(gè)多通道的時(shí)序控制器(DG535)來實(shí)現(xiàn).

    由于離子透鏡形成的電場不是平行電場,離子在探測器上形成弧形的半徑r不再是簡單的由碎片速度v和飛行時(shí)間t決定,即r=vt,而是存在一個(gè)放大系數(shù)N,即r=Nvt.在實(shí)驗(yàn)過程中,需要對(duì)N值進(jìn)行校準(zhǔn).我們是用Br2在265 nm32光解的碎片Br的速度影像對(duì)N進(jìn)行校準(zhǔn)的,得到N=1.306.為了保證放大系數(shù)N一定,在實(shí)驗(yàn)中保持裝置條件不變.

    3 實(shí)驗(yàn)結(jié)果

    圖1是氯溴甲烷(CH2BrCl)在265 nm附近光解后生成的Br*和Br的原始二維離子速度影像及其反Abel變換圖,每個(gè)原始二維離子速度影像都是對(duì)應(yīng)光解碎片的三維速度分布在二維平面的投影.光解碎片的速度分布是圓柱對(duì)稱的,激光的偏振方向就是對(duì)稱軸,如圖1中箭頭所示.圖像的形狀由碎片的速度分布和角度分布決定.

    對(duì)二維離子速度影像進(jìn)行反Abel變換可重構(gòu)光解碎片的三維速度分布.對(duì)重構(gòu)的三維分布在一定的角度進(jìn)行速率積分,可得到光解碎片的角度分布P(θ),用如下公式(5)去擬合光解碎片的角度分布P(θ),可得到各向異性參數(shù)β的值.

    其中θ表示光解碎片反沖速度和光解激光偏振方向的夾角.P(θ)和P2(cosθ)分別表示離子信號(hào)強(qiáng)度和二階勒讓德多項(xiàng)式.β在2(平行躍遷)到-1(垂直躍遷)之間取值.在265 nm附近光解時(shí)β的取值是β(Br*)= 1.67±0.16和β(Br)=0.94±0.09.先前實(shí)驗(yàn)27,28,30得到的結(jié)果為:β(Br*)在(1.4±0.2)與(1.51±0.05)之間,β(Br)在(0.7±0.1)與(0.8±0.1)之間,與我們的結(jié)果在實(shí)驗(yàn)誤差范圍內(nèi)是一致的.從氯溴甲烷的吸引光譜28可以看出,在265 nm附近3Q1[1A?,2A?]、3Q0[3A?]和1Q1[2A?,4A?]態(tài)的吸收截面在同一數(shù)量級(jí).可見,Br通道源于這三個(gè)態(tài),其中,3Q1[1A?,2A?]和1Q1[2A?,4A?]態(tài)(都是垂直躍遷)是絕熱解離,3Q0[3A?]態(tài)(平行躍遷)經(jīng)非絕熱解離產(chǎn)生大部分的Br.Br*通道主要源于3Q0[3A?]態(tài),1Q1[2A?,4A?]態(tài)中的4A?態(tài)也通過非絕熱交叉產(chǎn)生Br*.這與文獻(xiàn)27,28,30報(bào)道的實(shí)驗(yàn)結(jié)果一致.

    圖1 CH2BrCl在265 nm附近光解產(chǎn)生的Br*和Br碎片的原始離子影像(a,b)及其反Abel變換圖(c,d)Fig.1 Raw ion images(a,b)and inverseAbel images (c,d)of Br*and Br fragments from the photolysis of CH2BrCl at near 265 nmThe arrow indicates the laser polarization direction. Br represents Br(2P3/2),Br*represents Br(2P1/2).

    圖2 在265 nm附近光解碎片Br*(a)和Br(b)的速度分布Fig.2 Velocity distributions for the Br*(a)and Br(b) fragments at near 265 nm

    對(duì)重構(gòu)的三維分布在一定的速度下進(jìn)行角度積分,可以得到光解碎片的速度分布P(v).圖2給出CH2BrCl在265 nm附近光解時(shí)產(chǎn)生的碎片Br*和Br的速度分布.用雙光子共振增強(qiáng)電離掃描了對(duì)應(yīng)HBr的X1Σ+→g3Σ-(v=0)躍遷的吸收帶(雙光子能量:74000-85000 cm-1),33發(fā)現(xiàn)有少量的HBr存在.我們推測:樣品中可能存在雜質(zhì)HBr,經(jīng)光解HBr+ hv→H+Br后產(chǎn)生速度較低的Br*和Br.圖2中的高斯背景可能主要源于此.我們?cè)谒俣确植糚(v)中減去高斯背景,然后利用公式(6)和(7)可將Br*和Br的速度分布轉(zhuǎn)換成總的質(zhì)心平動(dòng)能分布P(ET),

    Br*和Br通道的總平動(dòng)能分布P(ET)如圖3所示.

    4 分析與討論

    氯溴甲烷(CH2BrCl)在265 nm附近光解時(shí)滿足的能量關(guān)系為:

    圖3 在265 nm附近Br*(a)和Br(b)通道的總平動(dòng)能分布Fig.3 Total translational energy distributions for Br*(a) and Br(b)channels near 265 nmThe arrows indicate the available energies,Eavl.

    上式中光解光子能Ehv為451.8519(Br*)或452.0189 kJ·mol-1(Br),C-Br鍵的解離能D0(C-Br)為(266.40±0.80)kJ·mol-1.34Eint(Br[2PJ])是Br原子的電子能級(jí),基態(tài)Br的能量值定為0,激發(fā)態(tài)Br*的能量為44.12 kJ·mol-1.35Eint(CH2BrCl)是母體分子CH2BrCl的轉(zhuǎn)動(dòng)內(nèi)能和振動(dòng)內(nèi)能之和.在超聲分子束中,可以認(rèn)為母體分子的轉(zhuǎn)動(dòng)溫度接近絕對(duì)零度,而振動(dòng)溫度與脈沖閥的溫度相同,為30°C.傳統(tǒng)上認(rèn)為超聲分子束中母體分子的轉(zhuǎn)動(dòng)溫度和振動(dòng)溫度都接近絕對(duì)零度,但未給出直接的實(shí)驗(yàn)證據(jù).Butler等20,21研究表明母體分子振動(dòng)溫度與脈沖閥的溫度相同.后面我們將進(jìn)一步討論這個(gè)問題.用Gaussian 03軟件36在B3LYP/6-311++G(3df,3pd)水平上計(jì)算母體分子CH2BrCl的振動(dòng)頻率,并用修正因子0.9854修正后計(jì)算得到母體分子CH2BrCl的平均振動(dòng)內(nèi)能,即Eint(CH2BrCl)=2.507 kJ·mol-1.Eint(CH2Cl)是解離氯甲基自由基的內(nèi)能,等于轉(zhuǎn)動(dòng)內(nèi)能Erot(CH2Cl)與振動(dòng)內(nèi)能Evib(CH2Cl)之和.

    對(duì)于每個(gè)測得的光解碎片的總平動(dòng)動(dòng)能ET,由角動(dòng)量守恒碰撞模型可得到對(duì)應(yīng)的氯甲基自由基的轉(zhuǎn)動(dòng)動(dòng)能.光子的角動(dòng)量可忽略不計(jì),因?yàn)檫@部分角動(dòng)量一般耦合到光電子的角動(dòng)量中去了.由角動(dòng)量守恒定律可知:反射碎片Br的經(jīng)典軌道角動(dòng)量與氯甲基自由基的轉(zhuǎn)動(dòng)角動(dòng)量大小相等、方向相反.用公式表示為:

    其中μ是CH2Cl+Br體系的約合質(zhì)量,是Br和 CH2Cl的相對(duì)速度vrel的大小,b是碰撞參數(shù),ω是氯甲基自由基的轉(zhuǎn)動(dòng)角速度,I是氯甲基自由基的轉(zhuǎn)動(dòng)慣量.假定解離過程中Br和CH2Cl之間的相互作用力沿C-Br鍵,則氯甲基自由基的轉(zhuǎn)軸通過CH2Cl的質(zhì)心,且垂直于CH2Cl的質(zhì)心和C-Br鍵構(gòu)成的平面.由公式(9)可得到解離氯甲基自由基的轉(zhuǎn)動(dòng)內(nèi)能Erot(CH2Cl)與光解碎片的總平動(dòng)動(dòng)能ET的關(guān)系:

    將(10)式代入(8)式,整理后可得:

    在應(yīng)用公式(11)計(jì)算解離氯甲基自由基的振動(dòng)內(nèi)能Evib(CH2Cl)時(shí),需要考慮解離氯甲基自由基的幾何構(gòu)型與CH2BrCl分子內(nèi)CH2Cl部分的幾何構(gòu)型可能有差異.由于CH2BrCl分子的A態(tài)解離是在Franck-Condon區(qū)域完成的,這種差異很小,可以忽略.另外,解離氯甲基自由基的振動(dòng)會(huì)改變碰撞參數(shù)b和轉(zhuǎn)動(dòng)慣量I,這樣在角動(dòng)量不變的情況下會(huì)改變解離氯甲基自由基的轉(zhuǎn)動(dòng)內(nèi)能.我們先利用Gaussian 03軟件在B3LYP/6-311++G(3df,3pd)水平上計(jì)算CH2BrCl分子的振動(dòng)模,然后沿各個(gè)振動(dòng)模的振動(dòng)方向(由Gaussian 03軟件的輸出文件給出)位移至振幅處,計(jì)算各坐標(biāo)點(diǎn)的μb2/I因數(shù).發(fā)現(xiàn)Cl-C-Br剪刀振動(dòng)模(CBrCl scissors,v6=229 cm-1)對(duì)μb2/I因數(shù)變化的影響最顯著,而只和H原子運(yùn)動(dòng)相關(guān)的振動(dòng)模對(duì)μb2/I因數(shù)的影響很小.在后面的計(jì)算中,我們計(jì)算沿v6振動(dòng)模位移矢量方向各坐標(biāo)點(diǎn)的μb2/I因數(shù),并以振動(dòng)模的簡正振動(dòng)波函數(shù)的平方作為對(duì)應(yīng)的概率密度.得到每個(gè)測得的光解碎片的總平動(dòng)動(dòng)能ET對(duì)應(yīng)的解離氯甲基自由基振動(dòng)內(nèi)能的子分布.最后用這個(gè)振動(dòng)內(nèi)能子分布卷積光解碎片的總平動(dòng)動(dòng)能分布P(ET),得到解離氯甲基自由基的振動(dòng)內(nèi)能分布P(Evib).20,21

    由下列公式(12)可得到兩個(gè)通道的可資用能Eavl,

    由兩個(gè)通道的總平動(dòng)能分布和下列公式(13)可得到兩個(gè)通道產(chǎn)生的氯甲基自由基的內(nèi)能分布,

    圖4和圖5分別是得到的源于Br*和Br通道的氯甲基自由基的內(nèi)能分布和振動(dòng)內(nèi)能分布.

    由圖4和圖5可見,氯甲基自由基的內(nèi)能較低時(shí),振動(dòng)內(nèi)能接近零,內(nèi)能主要是轉(zhuǎn)動(dòng)能,與硬碰撞模型假設(shè)基本相符.當(dāng)內(nèi)能增加時(shí),振動(dòng)內(nèi)能所占比重增加,直到內(nèi)能較高時(shí),轉(zhuǎn)動(dòng)內(nèi)能接近零,內(nèi)能主要是振動(dòng)能.所以,硬碰撞模型和軟碰撞模型都不能很好地模擬分子的光解過程.

    由圖4(b)和圖5(b)還可以看到,氯甲基自由基的振動(dòng)內(nèi)能出現(xiàn)負(fù)值.這是由于母體分子CH2BrCl的振動(dòng)內(nèi)能轉(zhuǎn)化成碎片的平動(dòng)能和轉(zhuǎn)動(dòng)能造成的.在公式(11)中Eint(CH2BrCl)用的是平均振動(dòng)內(nèi)能,當(dāng)高于平均振動(dòng)內(nèi)能的振動(dòng)能轉(zhuǎn)化成碎片的平動(dòng)能和轉(zhuǎn)動(dòng)能時(shí),碎片的平動(dòng)能和轉(zhuǎn)動(dòng)能之和就大于可資用能,計(jì)算得到的氯甲基自由基的振動(dòng)內(nèi)能就為負(fù).這部分的振動(dòng)內(nèi)能分布反映了吸收光解光子前母體分子CH2BrCl的振動(dòng)內(nèi)能分布.我們用它們來分析超聲分子束中母體分子CH2BrCl的振動(dòng)溫度,圖6是Br*通道的數(shù)據(jù)分析示意圖.

    圖4 源于CH2Cl+Br*通道的氯甲基自由基的內(nèi)能分布(a)和振動(dòng)內(nèi)能分布(b)Fig.4 Internal energy distributions(a)and vibrational internal energy distributions(b)of the chloromethyl radical from CH2Cl+Br*channel

    圖5 源于CH2Cl+Br通道的氯甲基自由基的內(nèi)能分布(a)和振動(dòng)內(nèi)能分布(b)Fig.5 Internal energy distributions(a)and vibrational internal energy distributions(b)of the chloromethyl radical from CH2Cl+Br channel

    如果采用Boltamann分布擬合,由Br*通道的數(shù)據(jù)得到母體分子CH2BrCl的振動(dòng)溫度為1425 K,由Br通道的數(shù)據(jù)得到振動(dòng)溫度為1277 K.如果采用Bose分布擬合,振動(dòng)溫度也都在1000 K以上.這個(gè)溫度顯然過高,可能是母體分子CH2BrCl高振動(dòng)能級(jí)的弛豫幾率較高所致.這也能說明超聲分子束中母體分子CH2BrCl的振動(dòng)溫度不是絕對(duì)零度.Butler等20-23設(shè)定母體分子振動(dòng)溫度與脈沖閥的溫度相同,得出的穩(wěn)定自由基的平動(dòng)能分布計(jì)算結(jié)果與實(shí)驗(yàn)測量結(jié)果一致,可以說這個(gè)認(rèn)定是合理的.

    由圖4(b)可見,Br*通道的氯甲基自由基的振動(dòng)內(nèi)能分布呈現(xiàn)4個(gè)峰值.由圖5(b)可見,Br通道的氯甲基自由基的振動(dòng)內(nèi)能分布呈現(xiàn)5個(gè)峰值.我們?cè)诒?中給出了各個(gè)峰值對(duì)應(yīng)的振動(dòng)能,表1中的誤差是各峰的半寬度.對(duì)比氯甲基自由基振動(dòng)模的振動(dòng)頻率,37,38對(duì)各個(gè)振動(dòng)內(nèi)能分布峰進(jìn)行了標(biāo)識(shí),如表1所示.

    圖6 超聲分子束中母體分子(CH2BrCl)的振動(dòng)溫度分析示意圖Fig.6 Schematic of vibrational temperature analysis of parent molecules(CH2BrCl)in the supersonic molecular beam

    表1 氯甲基自由基的振動(dòng)模與振動(dòng)內(nèi)能分布各峰值能量的比較Table 1 Comparison between the vibrational modes and the peak energies of the vibrational internal energy distributions of the chloromethyl radical

    從表1可以看到,兩個(gè)通道的氯甲基自由基的振動(dòng)峰多數(shù)為和頻的激發(fā).Jung等29研究CH2BrCl在234 nm的光解時(shí),發(fā)現(xiàn)Br通道的總平動(dòng)能分布出現(xiàn)雙峰,能量間隔為25 kJ·mol-1.他們把它歸結(jié)為氯甲基自由基的v2(CH2scissors)振動(dòng)模的激發(fā).但是,CH2Cl的v2(CH2scissors)振動(dòng)模的能量約為17 kJ·mol-1,這樣標(biāo)識(shí)不十分合理.由于內(nèi)能分布中包含有振動(dòng)能和轉(zhuǎn)動(dòng)能,且峰寬過寬(約10 kJ·mol-1),確實(shí)給標(biāo)識(shí)帶來困難.角動(dòng)量守恒碰撞模型則可以解決這些困難.

    圖7 CH2Cl自由基的v4振動(dòng)模(a)和CH2BrCl的v5振動(dòng)模(b)的比較Fig.7 Comparison between the vibrational mode v4of CH2Cl radical(a)and the vibrational mode v5of CH2BrCl(b)

    由圖4(b)和圖5(b)可見,兩個(gè)通道的氯甲基自由基的振動(dòng)內(nèi)能分布都不是熱平衡分布.Br通道的氯甲基自由基的振動(dòng)內(nèi)能分布幾率最大的峰出現(xiàn)在振動(dòng)內(nèi)能較大的位置,與Br*通道的氯甲基自由基的振動(dòng)內(nèi)能分布有明顯的不同.其原因可能不僅是Br通道的可資用能多出44.12 kJ·mol-1.比較圖4 (b)、圖5(b)和表1中的數(shù)據(jù)后,我們發(fā)現(xiàn)Br*通道中氯甲基自由基的ν4(CH2wag)振動(dòng)模的激發(fā)幾率最大,而Br通道中氯甲基自由基的v5(CH2a-stretch)振動(dòng)模的激發(fā)幾率最大.對(duì)比CH2BrCl的ν5(CBr stretch)振動(dòng)模39和氯甲基自由基的ν4(CH2wag)振動(dòng)模,如圖7所示,我們發(fā)現(xiàn)這兩個(gè)振動(dòng)模的CH2搖擺(CH2wag)振動(dòng)位移部分是基本一致的.可見CH2Cl-Br的解離振動(dòng)(CBr stretch)可以激發(fā)氯甲基自由基的ν4(CH2wag)振動(dòng).CH2BrCl的ν5(CBr stretch)振動(dòng)模中沒有與氯甲基自由基的ν5(CH2a-stretch)振動(dòng)模相近的位移,但是CH2BrCl的ν7(CH2a-stretch)振動(dòng)模39和氯甲基自由基的ν5(CH2a-stretch)振動(dòng)模的位移基本相同,見圖8.可見氯甲基自由基的ν5(CH2a-stretch)振動(dòng)模的激發(fā)可能是由CH2BrCl的ν7(CH2a-stretch)振動(dòng)模弛豫產(chǎn)生的.這意味著在吸收光解光子后母體分子CH2BrCl的ν7(CH2a-stretch)振動(dòng)模被激發(fā)了.氯甲基自由基的其它振動(dòng)模的激發(fā)也可作相似的分析.用Gaussian 03軟件,36在B3LYP/6-311++G(3df,3pd)水平上優(yōu)化了CH2BrCl基態(tài)的結(jié)構(gòu),在UCIS/6-311++G(3df,3pd)水平上優(yōu)化了3Q1、3Q0和1Q1態(tài)的結(jié)構(gòu).發(fā)現(xiàn)相對(duì)于基態(tài),三個(gè)激發(fā)態(tài)的C-Br鍵長顯著增長,C-H鍵、C-Cl鍵的鍵長變短,∠H-C-H、∠H-C-Cl等鍵角變大.這可能就是CH2BrCl母體分子n(Br)→σ* (C—Br)躍遷后剩余能量向除ν5(CBr stretch)振動(dòng)模外的其它振動(dòng)模弛豫的原因.Br*通道中氯甲基自由基的ν4(CH2wag)振動(dòng)模激發(fā)幾率最大,可能是因?yàn)槟阁w分子CH2BrCl中ν5(CBr stretch)振動(dòng)模的激發(fā)幾率最大.Br通道中氯甲基自由基的ν5(CH2a-stretch)振動(dòng)模激發(fā)幾率最大.我們分析可能的原因是:由于Br主要源于3Q0態(tài)與1Q1態(tài)的非絕熱交叉,而結(jié)構(gòu)優(yōu)化的結(jié)果表明這兩個(gè)態(tài)的C—H鍵長、H-C-Cl鍵角都有不同,所以這兩個(gè)態(tài)的相互作用增加了母體分子CH2BrCl的ν7(CH2a-stretch)振動(dòng)模被激發(fā)的幾率.總之,在吸收光解光子后母體分子CH2BrCl除ν5(CBr stretch)振動(dòng)模被激發(fā)外,還存在其它振動(dòng)模被激發(fā).當(dāng)然我們的結(jié)論還需要更加精細(xì)的實(shí)驗(yàn)和解離氯甲基自由基振動(dòng)激發(fā)的直接測量加以進(jìn)一步驗(yàn)證,這是我們今后準(zhǔn)備做的工作.

    圖8 CH2Cl自由基的v5振動(dòng)模(a)和CH2BrCl的v7振動(dòng)模(b)的比較Fig.8 Comparison between the vibrational mode v5of CH2Cl radical(a)and the vibrational mode v7of CH2BrCl(b)

    5 結(jié)論

    依據(jù)CH2BrCl在265 nm附近光解的實(shí)驗(yàn)結(jié)果,和角動(dòng)量守恒碰撞模型的分析,得出以下結(jié)論:

    (1)CH2BrCl+hv→Br*+CH2Cl通道產(chǎn)生的氯甲基自由基中被激發(fā)的振動(dòng)能級(jí)主要是v4、v3+v4、v2+v4和v2+v6振動(dòng)模.

    (2)CH2BrCl+hv→Br+CH2Cl通道產(chǎn)生的氯甲基自由基中被激發(fā)的振動(dòng)能級(jí)主要是v2+v6、v1+v3、v2+ v5、v2+v3+v5和v1+v5振動(dòng)模.

    (3)母體分子CH2BrCl在吸收光解光子后除有v5(CBr stretch)振動(dòng)模激發(fā)外,還有其它振動(dòng)模的激發(fā),如v7(CH2a-stretch)振動(dòng)模被激發(fā)了.

    致謝: 感謝美國University of Chicago的BUTLER J.Laurie教授和WOMACK C.Caroline博士為我們提供了用于角動(dòng)量守恒碰撞模型計(jì)算的基本C語言程序源代碼Rotprog2.c Version 1.0!

    (1) Barrie,L.A.;Bottenheim,J.W.;Schnell,R.C.;Crutzen,P.J.; Rasamussen,R.A.Nature 1988,334,138.

    (2)Anderson,J.G.;Toohey,D.W.;Brune,W.H.Science 1991, 251,39.

    (3) Foster,K.L.;Plastridge,R.A.;Bottenheim,J.W.;Shepso,P.B.; Finlayson-Pitts,B.J.;Spicer,C.W.Science 2001,291,471.

    (4) Zhu,R.;Tang,B.;Zhang,X.;Zhang,B.J.Phys.Chem.A 2010, 114(21),6188.

    (5) Niu,M.L.;Song,L.;Tang,X.F.;Zhou,X.G.;Liu,S.L.;Liu, F.Y.;Shan,X.B.;Sheng,L.S.Acta Phys.-Chim.Sin.2011,27, 1797.[牛銘理,宋 磊,唐小鋒,周曉國,劉世林,劉付軼,單曉斌,盛六四.物理化學(xué)學(xué)報(bào),2011,27,1797.]

    (6) Gougousi,T.;Samartzis,P.C.;Kitsopoulos,T.N.J.Chem. Phys.1998,108,5742.

    (7) Liu,Y.Z.;Qin,C.C.;Zhang,S.;Wang,Y.M.;Zhang,B.Acta Phys.-Chim.Sin.2011,27,965. [劉玉柱,秦朝朝,張 嵩,王艷梅,張 冰.物理化學(xué)學(xué)報(bào),2011,27,965.]

    (8)Zhang,F.;Wang,Y.M.;Zhang,B.;Feng,W.L.ActaPhys.-Chim. Sin.2010,26,1903.[張 鋒,王艷梅,張 冰,馮文林.物理化學(xué)學(xué)報(bào),2010,26,1903.]

    (9) Zhang,C.H.;Zhang,Y.;Zhang,S.;Zhang,B.Acta Phys.-Chim. Sin.2009,25,1708.[張昌華,張 延,張 嵩,張 冰.物理化學(xué)學(xué)報(bào),2009,25,1708.]

    (10) Cao,Z.Z.;Zhang,C.H.;Wang,Y.M.;Zhang,F.;Hua,L.Q.; Zhang,B.Acta Phys.-Chim.Sin.2009,25,423.[曹振洲,張昌華,王艷梅,張 鋒,華林強(qiáng),張 冰.物理化學(xué)學(xué)報(bào),2009,25, 423.]

    (11) Rubio-Lago,L.;Rodríguez,J.D.;García-Vela,A.;González, M.G.;Amaral,G.A.;Ba?ares,L.Phys.Chem.Chem.Phys. 2011,13,8186.

    (12) Sage,A,G.;Oliver,T.A.A.;Murdock,D.;Crow,M.B.; Ritchie,G.A.D.;Harvey,J.N.;Ashfold,M.N.R.Phys.Chem. Chem.Phys.2011,13,8075.

    (13) Song,L.;Yu,F.;Wu,L.X.;Zhou,X.G.;Liu,S.L.Acta Phys.-Chim.Sin.2011,27,801. [宋 磊,于 鋒,吳琍霞,周曉國,劉世林.物理化學(xué)學(xué)報(bào),2011,27,801.]

    (14) Eppink,A.T.J.B.;Parker,D.H.Rev.Sci.Instrum.1997,68(9), 3477.

    (15) Lin,J.J.;Zhou,J.;Shiu,W,;Liu,K.Rev.Sci.Instrum.2003,74 (4),2495.

    (16) Townsend,D.;Minitti,M.P.;Suits,A.G.Rev.Sci.Instrum. 2003,74(4),2530.

    (17) Chen,Z.;Liu,F.;Jiang,B.;Yang,X.;Parker,D.H.J.Phys. Chem.Lett.2010,1,1861.

    (18) Li,J.;Zhang,C.;Zhang,Q.;Chen,Y.;Huang,C.;Yang,X. J.Chem.Phys.2011,135,116102.

    (19) Tang,X.;Zhou,X.;Wu,M.;Liu,S.;Liu,F.;Shan,X.;Sheng,L. J.Chem.Phys.2012,136,034304.

    (20) Ratliff,B.J.;Womack,C.C.;Tang,X.N.;Landau,W.L.; Butler,L.J.;Szpunar,D.E.J.Phys.Chem.A 2010,114,4934.

    (21) Womack,C.C.;Fang,W.H.;Straus,D.B.;Butler,L.J. J.Phys.Chem.A 2010,114,13005.

    (22) Ratliff,B.J.;Alligood,B.W.;Butler,L.J.;Lee,S.H.;Lin,J.J. M.J.Phys.Chem.A 2011,115,9097.

    (23)Womack,C.C.;Booth,R.S.;Brynteson,M.D.;Butler,L.J.; Szpunar,D.J.Phys.Chem.A 2011,115,14559.

    (24) Tzeng,W.B.;Lee,Y.R.;Lin,S.M.Chem.Phys.Lett.1994, 227,467.

    (25) Cadman,P.;Simsons,J.P.Trans.Faraday Soc.1966,62,631.

    (26) Doucet,J.;Gilbert,R.;Sauvageau,P.;Sandorfy,C.J.Chem. Phys.1975,62,366.

    (27) McGivern,W.S.;Li,T.;Zou,P.;North,S.W.J.Chem.Phys. 1999,111,5771.

    (28)Zou,P.;McGivern,W.S.;North,S.W.Phys.Chem.Chem. Phys.2000,2,3785.

    (29) Lee,S.H.;Jung,Y.J.;Jung,K.H.Chem.Phys.2000,260,143.

    (30) Zhou,J.G.;Lau,K.C.;Hassanein,E.;Xu,H.F.;Tian,S.X.; Jones,B.;Ng,C.Y.J.Chem.Phys.2006,124,34309.

    (31) Dahl,D.A.;Appelhans,A.D.SIMION 3D,Version 7.0; Scientific Instrument Services,Inc.:Ringoes,USA,2004.

    (32) Jee,Y.J.;Jung,Y.J.;Jung,K.H.J.Chem.Phys.2006,115, 9739.

    (33) Callaghan,R.;Gordon,R.J.J.Chem.Phys.1990,93,4624.

    (34) Li,J.;Yang,J.;Mo,Y.;Lau,K.C.;Qian,X.M.;Song,Y.;Liu, J.;Ng,C.Y.J.Chem.Phys.2007,126,184304.

    (35) Arepalli,S.;Presser,N.;Robie,D.;Gordon,R.J.Chem.Phys. Lett.1985,117,64.

    (36) Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 03, RevisionA.01;Gaussian Inc.:Pittsburgh,PA,2003.

    (37) Whitney,E.S.;Haeber,T.;Schuder,M.D.;Blair,A.C.;Nesbitt, D.J.J.Chem.Phys.2006,125,054303.

    (38) Fridgen,T.D.;Zhang,X.K.;Parnis,J.M.;March,R.E. J.Phys.Chem.A 2000,104,3487.

    (39) Shimanouchi,T.Tables of Molecular Vibrational Frequencies Consolidated;National Bureau of Standards:Gaithersburg,MD, 1972;Vol.I,pp 1-160.

    December 21,2011;Revised:March 5,2012;Published on Web:March 6,2012.

    Vibrational Internal Energy of Chloromethyl Radical Formed by the Photodissociation of CH2BrCl at 265 nm

    ZHANG Xiu1WU Dong2TANG Bi-Feng1,*
    (1Department of Physics,Hubei Engineering University,Xiaogan 432000,Hubei Province,P.R.China;2School of Physical Science and Technology,Huazhong Normal University,Wuhan 430079,P.R.China)

    Photodissociation of CH2BrCl was investigated around 265 nm using resonance-enhanced multiphoton ionization technique combined with velocity map ion-imaging detection.The ion images of Br (2P1/2)and Br(2P3/2)were analyzed to obtain the corresponding velocity distributions and total translational energy distributions.Using an impulsive model invoking angular momentum conservation,the vibrational internal energy distributions of chloromethyl radical(·CH2Cl)formed by the photodissociation of CH2BrCl, were derived from the total translational energy distributions.In the CH2BrCl+hv→Br(2P1/2)+CH2Cl channel, v4,v3+v4,v2+v4and v2+v6vibrational modes were found to be excited in the radical;while in the CH2BrCl+hv→Br(2P3/2)+CH2Cl channel,the excited vibrational modes were v2+v6,v1+v3,v2+v5,v2+v3+v5,and v1+v5.The results further implied that,following absorption of one photon by the parent molecule CH2BrCl,other vibrational modes besides v5(CBr stretch)mode,such as v7(CH2a-stretch)mode,are excited in the parent molecule.

    CH2BrCl;Photodissociation;Velocity map ion-imaging;Impulsive model invoking angular momentum conservation;Vibrational internal energy distribution

    10.3866/PKU.WHXB201203061

    O643

    ?Corresponding author.Email:tangbf@wipm.ac.cn.Tel:+86-13986465587.

    The project was supported by the Science Foundation of the Education Department of Hubei Province,China(Z20082601).

    湖北省教育廳科學(xué)技術(shù)研究項(xiàng)目(Z20082601)

    猜你喜歡
    內(nèi)能角動(dòng)量母體
    “內(nèi)能”“內(nèi)能的利用”綜合測試題
    對(duì)經(jīng)典力學(xué)中的軌道角動(dòng)量和自轉(zhuǎn)角動(dòng)量的探討
    蒲公英
    遼河(2021年10期)2021-11-12 04:53:58
    基于角動(dòng)量模型的流場渦旋提取方法
    “內(nèi)能”“內(nèi)能的利用”綜合測試題
    用角動(dòng)量的方法解決并推廣一個(gè)功能關(guān)系問題
    夏季角動(dòng)量輸送變化與中國東部降水的關(guān)系
    “內(nèi)能和內(nèi)能的利用”易錯(cuò)點(diǎn)剖析
    “內(nèi)能”“內(nèi)能的利用”綜合測試題
    多胎妊娠發(fā)生的原因及母體并發(fā)癥處理分析
    午夜福利视频1000在线观看| 三级国产精品欧美在线观看| 少妇高潮的动态图| 精品少妇黑人巨大在线播放 | 欧美激情在线99| 中国国产av一级| 真实男女啪啪啪动态图| 色视频www国产| a级毛色黄片| 少妇猛男粗大的猛烈进出视频 | 97超视频在线观看视频| 99久久久亚洲精品蜜臀av| 尾随美女入室| 一卡2卡三卡四卡精品乱码亚洲| 亚洲第一电影网av| 狂野欧美白嫩少妇大欣赏| 人妻夜夜爽99麻豆av| 免费无遮挡裸体视频| 国产精品不卡视频一区二区| 国产精品一区二区三区四区免费观看| 插逼视频在线观看| 日韩av不卡免费在线播放| 黄片无遮挡物在线观看| 欧美不卡视频在线免费观看| 日韩视频在线欧美| 国产成人午夜福利电影在线观看| 亚洲欧美日韩东京热| 九九在线视频观看精品| 三级国产精品欧美在线观看| 99热网站在线观看| 国产亚洲精品久久久com| 亚洲国产精品国产精品| 亚洲第一区二区三区不卡| 亚洲av男天堂| 久久久久九九精品影院| 免费人成在线观看视频色| 九九在线视频观看精品| 亚洲av免费在线观看| 99riav亚洲国产免费| 国产久久久一区二区三区| 婷婷亚洲欧美| 美女国产视频在线观看| 久久午夜福利片| 人妻少妇偷人精品九色| 国产真实乱freesex| 午夜免费激情av| 51国产日韩欧美| 中出人妻视频一区二区| 99热这里只有是精品50| 日韩在线高清观看一区二区三区| 边亲边吃奶的免费视频| 久久久久久久久大av| av国产免费在线观看| 99热这里只有是精品50| 狂野欧美白嫩少妇大欣赏| 国产69精品久久久久777片| 亚洲精品粉嫩美女一区| 夜夜爽天天搞| 韩国av在线不卡| 国产 一区 欧美 日韩| 国内精品一区二区在线观看| 日韩一本色道免费dvd| 日韩欧美三级三区| 成人漫画全彩无遮挡| 美女cb高潮喷水在线观看| 极品教师在线视频| 十八禁国产超污无遮挡网站| 中文资源天堂在线| 我要看日韩黄色一级片| 国产av不卡久久| 两个人视频免费观看高清| 免费在线观看成人毛片| 国产又黄又爽又无遮挡在线| 欧美人与善性xxx| 一区二区三区高清视频在线| 久久婷婷人人爽人人干人人爱| 校园人妻丝袜中文字幕| 日韩 亚洲 欧美在线| 99久久人妻综合| 搞女人的毛片| 2022亚洲国产成人精品| 99久久精品热视频| 少妇裸体淫交视频免费看高清| 亚洲五月天丁香| 精品人妻熟女av久视频| 国产乱人偷精品视频| 亚洲精华国产精华液的使用体验 | 听说在线观看完整版免费高清| 超碰av人人做人人爽久久| av卡一久久| 亚洲国产精品成人久久小说 | 草草在线视频免费看| 我要看日韩黄色一级片| 色播亚洲综合网| 欧美色视频一区免费| 成熟少妇高潮喷水视频| 欧美潮喷喷水| 网址你懂的国产日韩在线| 熟女电影av网| 日本五十路高清| 美女黄网站色视频| 欧美zozozo另类| 日本黄色片子视频| 免费一级毛片在线播放高清视频| 在线a可以看的网站| 亚洲欧美精品专区久久| 久久人人爽人人爽人人片va| 亚洲精品粉嫩美女一区| videossex国产| 久久精品国产99精品国产亚洲性色| 中文字幕制服av| 天堂av国产一区二区熟女人妻| 久久久久久久久久成人| 国产一级毛片七仙女欲春2| av国产免费在线观看| 色吧在线观看| 久久精品影院6| 在线观看午夜福利视频| 色视频www国产| 91精品国产九色| 欧美潮喷喷水| 寂寞人妻少妇视频99o| 日日撸夜夜添| 女人十人毛片免费观看3o分钟| 在现免费观看毛片| 欧美色视频一区免费| 91精品国产九色| 国产伦精品一区二区三区四那| 中文亚洲av片在线观看爽| 国产伦精品一区二区三区视频9| 免费看美女性在线毛片视频| 久99久视频精品免费| 国产色爽女视频免费观看| av天堂中文字幕网| 精品国产三级普通话版| 亚洲av熟女| 久久久久性生活片| 中文在线观看免费www的网站| 国产精品一区二区三区四区免费观看| 成人三级黄色视频| 美女xxoo啪啪120秒动态图| 人妻少妇偷人精品九色| 天堂网av新在线| 国产高清不卡午夜福利| 少妇丰满av| 我要搜黄色片| 可以在线观看的亚洲视频| 精品久久久久久久久久免费视频| 高清午夜精品一区二区三区 | 国产v大片淫在线免费观看| 亚洲人成网站在线观看播放| 国产成人精品久久久久久| 少妇人妻一区二区三区视频| 国产精品爽爽va在线观看网站| 亚洲成人久久爱视频| 99久久精品一区二区三区| 久久久午夜欧美精品| 国产精品日韩av在线免费观看| 插阴视频在线观看视频| 在线播放国产精品三级| av天堂中文字幕网| 美女国产视频在线观看| videossex国产| 午夜久久久久精精品| 国产激情偷乱视频一区二区| 天天躁夜夜躁狠狠久久av| a级毛片免费高清观看在线播放| 简卡轻食公司| 99久久精品一区二区三区| 一进一出抽搐gif免费好疼| 丰满人妻一区二区三区视频av| 精品免费久久久久久久清纯| 免费在线观看成人毛片| 日韩欧美国产在线观看| 看十八女毛片水多多多| av免费观看日本| 综合色av麻豆| 国产真实乱freesex| 内地一区二区视频在线| 亚洲成av人片在线播放无| 日本免费一区二区三区高清不卡| 国产一区亚洲一区在线观看| 一本久久精品| 中文字幕人妻熟人妻熟丝袜美| 日韩国内少妇激情av| 免费av毛片视频| 乱系列少妇在线播放| 91狼人影院| 欧美bdsm另类| 非洲黑人性xxxx精品又粗又长| 看十八女毛片水多多多| 一级二级三级毛片免费看| 久久精品国产亚洲av天美| 成人毛片60女人毛片免费| 亚洲av第一区精品v没综合| 丝袜喷水一区| 一区福利在线观看| 中文字幕精品亚洲无线码一区| 黄色配什么色好看| 99在线人妻在线中文字幕| 99久国产av精品| 不卡一级毛片| а√天堂www在线а√下载| av在线老鸭窝| 亚洲精品成人久久久久久| 午夜精品在线福利| 一个人看的www免费观看视频| 国产一区二区激情短视频| 99热这里只有精品一区| 亚洲乱码一区二区免费版| 亚洲高清免费不卡视频| 亚洲不卡免费看| 成人综合一区亚洲| 国产v大片淫在线免费观看| 少妇的逼水好多| 亚洲欧美精品自产自拍| 色噜噜av男人的天堂激情| 中国国产av一级| 啦啦啦啦在线视频资源| 在线播放无遮挡| 波多野结衣高清无吗| or卡值多少钱| 最近的中文字幕免费完整| 亚洲性久久影院| 深夜精品福利| 亚洲第一电影网av| www日本黄色视频网| 久久久久网色| 日本成人三级电影网站| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 亚洲aⅴ乱码一区二区在线播放| 两个人视频免费观看高清| 女人十人毛片免费观看3o分钟| 久久婷婷人人爽人人干人人爱| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品电影小说 | 晚上一个人看的免费电影| 搞女人的毛片| 亚洲成人精品中文字幕电影| 在线观看美女被高潮喷水网站| 国产精品麻豆人妻色哟哟久久 | 1000部很黄的大片| 99riav亚洲国产免费| 麻豆成人av视频| АⅤ资源中文在线天堂| 亚洲美女搞黄在线观看| 深夜a级毛片| 国产片特级美女逼逼视频| www.av在线官网国产| 亚洲国产欧美人成| 日韩欧美国产在线观看| 久久这里只有精品中国| 特大巨黑吊av在线直播| 日本免费一区二区三区高清不卡| 亚洲四区av| 日本熟妇午夜| 欧美性猛交╳xxx乱大交人| 国产色婷婷99| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产视频首页在线观看| 欧美bdsm另类| 国产一区二区三区在线臀色熟女| 青春草国产在线视频 | 欧美激情国产日韩精品一区| 蜜臀久久99精品久久宅男| 九九热线精品视视频播放| 欧美xxxx性猛交bbbb| 少妇的逼水好多| 免费看a级黄色片| 日日干狠狠操夜夜爽| 亚洲不卡免费看| 欧美三级亚洲精品| 亚洲熟妇中文字幕五十中出| 亚洲成人久久性| 国产亚洲91精品色在线| 国产精品嫩草影院av在线观看| 成人欧美大片| 村上凉子中文字幕在线| 欧美日韩乱码在线| 99热这里只有精品一区| 午夜视频国产福利| 免费av观看视频| 黄片wwwwww| av卡一久久| 观看免费一级毛片| 国产免费一级a男人的天堂| 91麻豆精品激情在线观看国产| 人妻系列 视频| h日本视频在线播放| 国产一区二区激情短视频| 麻豆乱淫一区二区| 亚洲色图av天堂| 少妇人妻一区二区三区视频| 午夜老司机福利剧场| 亚洲国产精品国产精品| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站| 亚洲久久久久久中文字幕| av.在线天堂| 一级毛片电影观看 | 激情 狠狠 欧美| 亚洲久久久久久中文字幕| 国产成人精品久久久久久| 亚洲一级一片aⅴ在线观看| 村上凉子中文字幕在线| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 久久久午夜欧美精品| 国产精品久久久久久精品电影| 波多野结衣高清无吗| 亚洲乱码一区二区免费版| 免费看美女性在线毛片视频| 久久精品久久久久久久性| av在线老鸭窝| 国产精品久久久久久亚洲av鲁大| 精品免费久久久久久久清纯| 最近手机中文字幕大全| 午夜老司机福利剧场| 成人鲁丝片一二三区免费| 成人高潮视频无遮挡免费网站| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说 | av免费在线看不卡| 精品无人区乱码1区二区| 久久久久久久久久久免费av| www.色视频.com| 亚洲精品久久国产高清桃花| 国产探花极品一区二区| 99热只有精品国产| 国产单亲对白刺激| 国产熟女欧美一区二区| 亚洲国产欧美人成| 非洲黑人性xxxx精品又粗又长| 亚洲在久久综合| 26uuu在线亚洲综合色| 国产午夜精品久久久久久一区二区三区| 最近中文字幕高清免费大全6| 久久久久久久亚洲中文字幕| 日产精品乱码卡一卡2卡三| 一边摸一边抽搐一进一小说| 亚洲在线自拍视频| 中文欧美无线码| 欧美成人一区二区免费高清观看| 亚洲第一电影网av| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 99久久人妻综合| 亚洲欧美精品自产自拍| 激情 狠狠 欧美| 日本三级黄在线观看| 国内精品美女久久久久久| 偷拍熟女少妇极品色| 国产综合懂色| 国产av在哪里看| 亚洲五月天丁香| 国产毛片a区久久久久| 欧美一区二区亚洲| 国产成人freesex在线| 深爱激情五月婷婷| 久久精品国产亚洲av涩爱 | 青春草视频在线免费观看| 黑人高潮一二区| 老女人水多毛片| 国产麻豆成人av免费视频| 99久久九九国产精品国产免费| 久久精品影院6| 久久精品久久久久久久性| 亚洲中文字幕一区二区三区有码在线看| 中文字幕熟女人妻在线| 亚洲美女搞黄在线观看| 99国产极品粉嫩在线观看| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲熟妇中文字幕五十中出| 亚洲一区高清亚洲精品| 成人特级黄色片久久久久久久| 18禁在线播放成人免费| 精品人妻一区二区三区麻豆| 日产精品乱码卡一卡2卡三| 亚洲精品国产av成人精品| 麻豆成人av视频| 国产精品人妻久久久影院| 天堂影院成人在线观看| 日韩欧美三级三区| 亚洲av成人精品一区久久| 亚洲欧美成人综合另类久久久 | 欧美激情久久久久久爽电影| 日本爱情动作片www.在线观看| 日韩国内少妇激情av| 国产黄片视频在线免费观看| 亚洲av二区三区四区| 一边摸一边抽搐一进一小说| 天堂中文最新版在线下载 | .国产精品久久| 18禁在线无遮挡免费观看视频| 免费av不卡在线播放| videossex国产| 99视频精品全部免费 在线| 午夜老司机福利剧场| 全区人妻精品视频| 国产老妇伦熟女老妇高清| 两个人的视频大全免费| 秋霞在线观看毛片| 久久国产乱子免费精品| 久久精品国产亚洲av天美| 婷婷六月久久综合丁香| 久久久a久久爽久久v久久| 日韩欧美三级三区| 国产欧美日韩精品一区二区| a级毛片a级免费在线| kizo精华| 日本爱情动作片www.在线观看| 黄色一级大片看看| 97热精品久久久久久| 熟女电影av网| 国产人妻一区二区三区在| 91午夜精品亚洲一区二区三区| 韩国av在线不卡| 日韩欧美精品v在线| 亚洲av电影不卡..在线观看| 色吧在线观看| 免费一级毛片在线播放高清视频| 18禁在线播放成人免费| 亚洲色图av天堂| 18禁裸乳无遮挡免费网站照片| 精品熟女少妇av免费看| 国产精品美女特级片免费视频播放器| 国产av在哪里看| 久久草成人影院| 亚洲欧美成人综合另类久久久 | 亚洲欧美日韩东京热| 99国产精品一区二区蜜桃av| 熟女电影av网| 国产麻豆成人av免费视频| 免费看a级黄色片| 女的被弄到高潮叫床怎么办| 国产蜜桃级精品一区二区三区| 国产大屁股一区二区在线视频| 亚洲18禁久久av| 亚洲av中文字字幕乱码综合| 搡老妇女老女人老熟妇| 夫妻性生交免费视频一级片| 亚洲国产精品成人久久小说 | 久久久久久久久久久丰满| 久久久久久久久久黄片| 午夜免费男女啪啪视频观看| 日本av手机在线免费观看| 老熟妇乱子伦视频在线观看| 午夜福利视频1000在线观看| 91午夜精品亚洲一区二区三区| 欧美3d第一页| 免费av毛片视频| 日韩视频在线欧美| 一级毛片久久久久久久久女| 亚洲精品久久久久久婷婷小说 | 久99久视频精品免费| 最好的美女福利视频网| 亚洲美女搞黄在线观看| 国产亚洲91精品色在线| 国产伦一二天堂av在线观看| 国产成人a∨麻豆精品| 日本-黄色视频高清免费观看| 国产精品一区www在线观看| 自拍偷自拍亚洲精品老妇| 老司机福利观看| 丰满的人妻完整版| 99热精品在线国产| 超碰av人人做人人爽久久| 亚洲av二区三区四区| 国内精品宾馆在线| 2021天堂中文幕一二区在线观| 内地一区二区视频在线| АⅤ资源中文在线天堂| 搞女人的毛片| 一级毛片我不卡| 国产精品综合久久久久久久免费| 女人被狂操c到高潮| 久久久成人免费电影| 精品久久久噜噜| 我要看日韩黄色一级片| 国产精品一区二区在线观看99 | 男女边吃奶边做爰视频| 精品少妇黑人巨大在线播放 | av在线老鸭窝| 久久久精品欧美日韩精品| 伦理电影大哥的女人| 久久九九热精品免费| 午夜福利在线在线| 日本在线视频免费播放| 婷婷六月久久综合丁香| 国内久久婷婷六月综合欲色啪| 日韩欧美一区二区三区在线观看| 最近手机中文字幕大全| 亚洲在线观看片| 久久久久久伊人网av| 国产v大片淫在线免费观看| 五月伊人婷婷丁香| 内地一区二区视频在线| www.色视频.com| 精品人妻一区二区三区麻豆| 爱豆传媒免费全集在线观看| 日本一本二区三区精品| 黑人高潮一二区| 成人一区二区视频在线观看| 国产国拍精品亚洲av在线观看| 99热精品在线国产| 一边亲一边摸免费视频| 联通29元200g的流量卡| 久久人人爽人人爽人人片va| 久久久a久久爽久久v久久| 高清日韩中文字幕在线| 久久精品国产鲁丝片午夜精品| 久久精品夜夜夜夜夜久久蜜豆| 蜜桃久久精品国产亚洲av| 久久精品国产自在天天线| 欧美成人精品欧美一级黄| 日本色播在线视频| 亚洲精品456在线播放app| 亚洲精品成人久久久久久| 午夜福利成人在线免费观看| 美女黄网站色视频| 国产精品不卡视频一区二区| 欧美成人a在线观看| 国产女主播在线喷水免费视频网站 | 国产精品人妻久久久久久| 免费人成视频x8x8入口观看| 成年女人看的毛片在线观看| 日本撒尿小便嘘嘘汇集6| 热99在线观看视频| 国产av麻豆久久久久久久| 狠狠狠狠99中文字幕| av天堂中文字幕网| 三级毛片av免费| 久久精品久久久久久噜噜老黄 | 少妇熟女aⅴ在线视频| 亚洲五月天丁香| 亚洲一区二区三区色噜噜| 桃色一区二区三区在线观看| 搡女人真爽免费视频火全软件| 亚洲欧美日韩高清在线视频| 久久综合国产亚洲精品| 日本成人三级电影网站| 51国产日韩欧美| 亚洲精品成人久久久久久| 男女那种视频在线观看| 免费大片18禁| 高清午夜精品一区二区三区 | 深夜a级毛片| 国产久久久一区二区三区| 久久精品国产清高在天天线| 国产成人午夜福利电影在线观看| videossex国产| 男人舔奶头视频| 别揉我奶头 嗯啊视频| 九色成人免费人妻av| 国产精品乱码一区二三区的特点| 天天躁日日操中文字幕| 男的添女的下面高潮视频| 特大巨黑吊av在线直播| 国产av不卡久久| 久久久久久久久久成人| 美女cb高潮喷水在线观看| 中文欧美无线码| 国产亚洲av片在线观看秒播厂 | 夜夜夜夜夜久久久久| 国产精品一及| 91精品一卡2卡3卡4卡| 好男人视频免费观看在线| 男女啪啪激烈高潮av片| 亚洲精品乱码久久久v下载方式| 亚洲内射少妇av| 国产毛片a区久久久久| 麻豆国产av国片精品| 波多野结衣高清无吗| 精品久久久久久久久久久久久| or卡值多少钱| 插逼视频在线观看| 国产精品一二三区在线看| 国产精品,欧美在线| 尾随美女入室| 国产亚洲精品久久久com| 亚洲自拍偷在线| 精品久久久噜噜| 日本爱情动作片www.在线观看| 麻豆久久精品国产亚洲av| 亚洲av不卡在线观看| 国产一区二区激情短视频| 亚洲av二区三区四区| 精品无人区乱码1区二区| 久久久久久久久大av| 色5月婷婷丁香| 久久午夜亚洲精品久久| 久久精品夜夜夜夜夜久久蜜豆| 免费av观看视频| 国内精品宾馆在线| 日韩精品有码人妻一区| 亚洲国产日韩欧美精品在线观看| 久久九九热精品免费| 国内揄拍国产精品人妻在线| 日韩在线高清观看一区二区三区| 亚洲四区av| 男女下面进入的视频免费午夜| 午夜免费男女啪啪视频观看| 国产亚洲av嫩草精品影院| 少妇丰满av| 日韩成人伦理影院| 久久久国产成人免费| 国产成人91sexporn| 久久精品国产亚洲av涩爱 | 国产不卡一卡二| 亚洲成av人片在线播放无| 一本一本综合久久| 久久99热这里只有精品18| 国产精品99久久久久久久久|