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

    N-(4-磺酸丁基)吡啶硫酸氫鹽離子液體的結(jié)構(gòu)分析

    2017-01-07 03:16:13張文林錢浩智李俊飛李春利
    關(guān)鍵詞:兩性離子吡啶氫鍵

    張文林,錢浩智,陳 瑤,李俊飛,李春利

    (河北工業(yè)大學(xué) 化工學(xué)院,天津 300130)

    N-(4-磺酸丁基)吡啶硫酸氫鹽離子液體的結(jié)構(gòu)分析

    張文林,錢浩智,陳 瑤,李俊飛,李春利

    (河北工業(yè)大學(xué) 化工學(xué)院,天津 300130)

    采用密度泛函理論(DFT),在基組6-31+G*的作用下,對(duì) N-(4-磺酸丁基)吡啶硫酸氫鹽離子液體的兩性化合物 [C4SO3Py]±、陽(yáng)離子 [C4SO3HPy]+、離子液體 [C4SO3HPy]HSO4進(jìn)行了結(jié)構(gòu)優(yōu)化、頻率分析及能量計(jì)算,并進(jìn)行了相互作用能計(jì)算和電荷分析.研究發(fā)現(xiàn):在兩性離子、陽(yáng)離子、離子液體的穩(wěn)定構(gòu)型中,分子間和分子內(nèi)單一或同時(shí)存在著氫鍵.離子液體結(jié)構(gòu)最穩(wěn)定時(shí),陰陽(yáng)離子在功能基團(tuán)附近相互作用形成氫鍵,HSO3中的氫與HSO4中的氧結(jié)合,最后以兩性離子與 H2SO4相互作用的形式存在.電荷分析顯示離子液體中陰陽(yáng)離子功能基團(tuán)上的氫具酸性,吡啶環(huán)中的氫原子也顯一定酸性.

    離子液體;結(jié)構(gòu)分析;密度泛函;氫鍵;算法

    0 引言

    離子液體是由有機(jī)陽(yáng)離子和無(wú)機(jī)或有機(jī)陰離子組成,在室溫或室溫附近的溫度下呈液體狀態(tài)的鹽類[1-2].由于其優(yōu)于常規(guī)有機(jī)溶劑及催化劑的特性,逐漸成為國(guó)際綠色高新科技前沿研究熱點(diǎn),離子液體構(gòu)型關(guān)系與分子設(shè)計(jì)是當(dāng)今離子液體亟待解決的關(guān)鍵技術(shù)之一[3].大量有關(guān)離子液體的研究?jī)H集中在催化反應(yīng)[4-10]方面,理論研究較少.Elizabeth等[11]運(yùn)用abinitio法研究了二烷基咪唑離子液體的構(gòu)型及陰陽(yáng)離子間的相互作用關(guān)系.李浩然等[12]通過密度泛函理論DFT對(duì)離子液體與水間的作用模型進(jìn)行了系統(tǒng)的研究.張鎖江等[13]用DFT對(duì)不同陰離子的咪唑類離子液體結(jié)構(gòu)、陰陽(yáng)離子構(gòu)型、氫鍵作用及氫鍵網(wǎng)絡(luò)進(jìn)行研究分析.

    目前對(duì)離子液體的理論研究主要集中在咪唑類離子液體[14-16],對(duì)其他類型的離子液體及功能化離子液體的研究較少.只有從分子層面[17]上了解離子液體的微觀結(jié)構(gòu)及相互作用關(guān)系,才能深入理解其優(yōu)越特性的深層次原因,具有重要的研究?jī)r(jià)值與科學(xué)意義.

    本文以N-(4-磺酸丁基)吡啶硫酸氫鹽 [C4SO3HPy]HSO4為研究對(duì)象,采用DFT對(duì)離子液體的兩性離子 [C4SO3Py]±、陽(yáng)離子 [C4SO3HPy]+及離子對(duì)[C4SO3HPy]HSO4的結(jié)構(gòu)進(jìn)行全優(yōu)化,并對(duì)其進(jìn)行頻率分析、相互作用能計(jì)算及電荷分析,探討功能基團(tuán)對(duì)離子液體構(gòu)型的影響及功能化離子液體相互作用特點(diǎn).

    1 計(jì)算方法

    在離子液體的量化計(jì)算中,密度泛函理論(DFT)方法是基于電荷密度自洽的方法,此方法考慮了電子之間的相關(guān)作用,因此得到的能量要較哈特里—??朔匠蹋℉artree-Fock,HF)方法來(lái)得精確,它是目前最常用的量子化學(xué)計(jì)算方法.其中,DFT理論中的B3LYP(hyPrid Pecke3-Lee-Yang-Parr)方法,既保持了從頭算MO法的很多優(yōu)點(diǎn),又節(jié)省了機(jī)時(shí).考慮到 [C4SO3HPy]HSO4離子液體的環(huán)狀結(jié)構(gòu)及離子對(duì)的相互作用,選用的基組為6-31+G*[18-21].

    通過對(duì)所有可能構(gòu)型進(jìn)行全優(yōu)化,在優(yōu)化結(jié)構(gòu)的基礎(chǔ)上對(duì)其進(jìn)行振動(dòng)頻率分析,沒有虛頻則為穩(wěn)定結(jié)構(gòu).使用經(jīng)過零點(diǎn)能(Zero pointenergy,ZPE)校正后的能量比較各構(gòu)型間的穩(wěn)定性順序.對(duì)于同一物質(zhì)的不同構(gòu)型,能量越低結(jié)構(gòu)的穩(wěn)定性越高,也就是說總能量最低的構(gòu)型為最優(yōu)構(gòu)型.

    在計(jì)算離子對(duì)間的相互作用能時(shí),考慮了基組重疊誤差(Pasissetsuperpositionerror,PSSE)的影響,基組疊加誤差采用Poys和Pornardi提出的Counterpoise(CP)方法消除.對(duì)于AB作用體系,其相互作用能可用式 (1)表示

    在計(jì)算相互作用能時(shí),為了使計(jì)算結(jié)果更準(zhǔn)確,所用的計(jì)算能量均為經(jīng)過零點(diǎn)能校正后的能量.

    本文以N-(4-磺酸丁基)吡啶硫酸氫鹽 [C4SO3HPy]HSO4為研究對(duì)象,采用DFT對(duì)離子液體的兩性離子 [C4SO3Py]±、陽(yáng)離子 [C4SO3HPy]+及離子液體 [C4SO3HPy]HSO4的結(jié)構(gòu)進(jìn)行了全優(yōu)化及能量計(jì)算,并對(duì)其進(jìn)行頻率分析、相互作用能計(jì)算及電荷分析,確定了功能基團(tuán)對(duì)離子液體構(gòu)型的影響及功能化離子液體相互作用特點(diǎn).

    2 結(jié)果與討論

    2.1 兩性離子液體的結(jié)構(gòu)分析

    首先對(duì)兩性離子液體的結(jié)構(gòu)進(jìn)行考察.為了得到離子液體[C4SO3HPy]HSO4前體化合物的穩(wěn)定結(jié)構(gòu),參考前人的計(jì)算方法[11],在吡啶環(huán)的穩(wěn)定構(gòu)型上依次引入甲基,最后優(yōu)化得到離子液體陽(yáng)離子 [C5Py]+的所有穩(wěn)定結(jié)構(gòu),然后將 [C5Py]+陽(yáng)離子碳鏈末端的甲基換成 SO3,并對(duì)得到的所有構(gòu)型進(jìn)行優(yōu)化計(jì)算,通過頻率分析,無(wú)虛頻確定為穩(wěn)定結(jié)構(gòu).圖1為兩性離子液體的最穩(wěn)定構(gòu)型.對(duì)所得穩(wěn)定構(gòu)型進(jìn)行分析,表1列出其主要的結(jié)構(gòu)參數(shù)及相對(duì)能量.

    從圖1可以看出,由于靜電作用的影響,功能基團(tuán)上帶負(fù)電荷的SO3具有較強(qiáng)的形成分子內(nèi)氫鍵的能力.對(duì)于形成分子內(nèi)氫鍵的兩性離子液體,同一個(gè)氧原子可以同時(shí)與2個(gè)不同的氫原子形成氫鍵(a1中的25O、26O;a2中的26O),而同一個(gè)氫原子也可以同時(shí)是2個(gè)氧原子的氫鍵給予體(a1中的11H).

    通過對(duì)比所形成的兩性離子液體穩(wěn)定構(gòu)型的能量發(fā)現(xiàn)分子內(nèi)形成的氫鍵數(shù)目越多,體系能量越低,結(jié)構(gòu)越穩(wěn)定.不含分子內(nèi)氫鍵的構(gòu)型經(jīng)頻率計(jì)算也無(wú)虛頻,確定為穩(wěn)定結(jié)構(gòu)(如a3).通過對(duì)計(jì)算結(jié)果,如鍵長(zhǎng)(O26…H11)、鍵角(∠C5-H11…O26)以及二面角(DN6-C12-C15-C18)的比較,可知不同穩(wěn)定結(jié)構(gòu)的鍵長(zhǎng)鍵角是不同的,主要是因?yàn)榉肿觾?nèi)氫鍵的影響,而吡啶側(cè)鏈上的氫影響更大,由此可以推知,兩性離子液體不同穩(wěn)定結(jié)構(gòu)的形成主要是受吡啶側(cè)鏈烷基空間位置的影響.

    圖1 [C4SO3Py]±兩性離子的最穩(wěn)定構(gòu)型Fig.1 Themoststable configurationsof the zw itterion[C4SO3Py]±

    表1 [C4SO3Py]±兩性離子的最穩(wěn)定構(gòu)型結(jié)構(gòu)參數(shù)及能量差Tab.1 The parameterand energy differencesof themoststable configurationsof the zw itterion

    2.2 陽(yáng)離子的結(jié)構(gòu)分析

    對(duì)離子液體的陽(yáng)離子 [C4SO3Py]+進(jìn)行結(jié)構(gòu)全優(yōu)化計(jì)算,得到含有分子內(nèi)氫鍵及不含分子內(nèi)氫鍵的穩(wěn)定結(jié)構(gòu),如圖2所示.陽(yáng)離子最穩(wěn)定構(gòu)型的主要結(jié)構(gòu)參數(shù)及相對(duì)能量如表2所示.

    計(jì)算結(jié)果表明,當(dāng)SO3獲得一個(gè)氫原子變?yōu)?SO3H后,整個(gè)體系的能量降低,且體系內(nèi)的氫鍵無(wú)論從鍵數(shù)和鍵能上都削弱了.b1構(gòu)型中C5-H11…O25的氫鍵和C12-H14…O25的氫鍵距離都比相應(yīng)的兩性離子長(zhǎng),說明所形成氫鍵的鍵能減弱了.陽(yáng)離子 [C4SO3Py]+只形成分子內(nèi)三氫鍵而沒有四氫鍵的穩(wěn)定構(gòu)型.通過對(duì)比所有陽(yáng)離子穩(wěn)定結(jié)構(gòu)的能量發(fā)現(xiàn),陽(yáng)離子內(nèi)氫鍵數(shù)目越多,總能量越低,結(jié)構(gòu)越穩(wěn)定.

    圖2 [C4SO3Py]+陽(yáng)離子的最穩(wěn)定構(gòu)型Fig.2 Themoststable configurationsof the cations[C4SO3Py]+

    表2 [C4SO3HPy]+陽(yáng)離子的最穩(wěn)定構(gòu)型結(jié)構(gòu)參數(shù)及能量差Tab.2 The parameterand energy differencesof themoststable configurationsof the cations

    2.3 離子液體的結(jié)構(gòu)分析

    將陰離子HSO4放在陽(yáng)離子附近不同位置,對(duì)其可能的初始結(jié)構(gòu)進(jìn)行優(yōu)化,得到代表性的幾種構(gòu)型如圖3所示.

    由圖3可以看出,離子液體中陰離子能從吡啶環(huán)陽(yáng)離子附近不同位置與其碳原子及硫原子上的氫相互作用形成氫鍵,在所有構(gòu)型中c1的能量最低,所形成的氫鍵最強(qiáng)(25O…34H的鍵長(zhǎng)僅為1.659?),構(gòu)型最穩(wěn)定.從 c1結(jié)構(gòu)中可以看出,除陽(yáng)離子本身形成氫鍵外,離子對(duì)間也形成了三氫鍵的穩(wěn)定構(gòu)型,在模擬計(jì)算過程中,陽(yáng)離子功能基團(tuán)HSO3上的28H會(huì)與27O發(fā)生斷裂,而與陰離子HSO4上的32O結(jié)合,在此期間,鍵長(zhǎng)由27O-28H的0.979?拉長(zhǎng)到27O…28H的1.738?.與此同時(shí)28H與陰離子HSO4上的32O形成鍵長(zhǎng)為1.004?的28H-32O鍵.陰離子HSO4上原有的32O-34H之間的鍵長(zhǎng)由0.972?拉長(zhǎng)到1.009?.此時(shí)離子液體體系中硫酸的H-O鍵長(zhǎng)(1.004?、1.009?)較純硫酸液體中鍵長(zhǎng)(0.977?、0.977?)增加,鍵能降低,但由于其與兩性離子液體 [C4SO3Py]±之間形成氫鍵作用而致使氫原子不易解離下來(lái).從計(jì)算結(jié)果可以看出,當(dāng)陰離子與陽(yáng)離子相互作用時(shí),其主要作用位置在陰離子HSO4與陽(yáng)離子功能基團(tuán)上的HSO3之間,而當(dāng)它們之間形成氫鍵時(shí),HSO3中的氫更易與HSO4中的氧結(jié)合,最后以兩性離子與H2SO4相互作用的形式存在.

    2.4 相互作用能計(jì)算

    對(duì)離子液體 [C4SO3HPy]HSO4進(jìn)行相互作用能的計(jì)算時(shí),采用counterpoise(CP)理論對(duì)其進(jìn)行了基組重疊誤差的校正,通過公式 (1)計(jì)算求得.EAB、EA、EB分別為離子液體 [C4SO3HPy]HSO4、陽(yáng)離子[C4SO3HPy]+、陰離子穩(wěn)定構(gòu)型的能量.對(duì)所得上述8種離子液體的穩(wěn)定構(gòu)型進(jìn)行相互作用能計(jì)算,結(jié)果列于表3.

    圖3 離子液體 [C4SO3HPy]HSO4的穩(wěn)定構(gòu)型Fig.3 Themoststable configurationsof the ionic liquid[C4SO3HPy]HSO4

    表3 離子液體 [C4SO3HPy]HSO4不同構(gòu)型的總能量,相對(duì)能量及相互作用能Tab.3 The totalenergies EZPE,the relative energies,the interaction energies

    表3 離子液體 [C4SO3HPy]HSO4不同構(gòu)型的總能量,相對(duì)能量及相互作用能Tab.3 The totalenergies EZPE,the relative energies,the interaction energies

    構(gòu)型 總能量EZPE( a u) 相對(duì)能量EZPE* / ( k J / m o l ) ECP/ ( k J / m o l ) c 1 1 7 2 9 . 3 5 0 5 7 8 0 3 8 1 . 8 6 0 6 0 1 c 2 1 7 2 9 . 3 5 0 0 4 2 1 . 4 0 7 2 6 8 3 8 0 . 4 5 3 3 3 1 c 3 1 7 2 9 . 3 3 2 7 1 8 4 6 . 8 9 1 4 3 1 3 6 0 . 0 1 3 8 0 9 c 4 1 7 2 9 . 3 3 2 6 4 7 4 7 . 0 7 7 8 4 1 3 3 4 . 7 8 2 7 6 1 c 5 1 7 2 9 . 3 2 8 4 3 4 5 8 . 1 3 9 0 7 2 3 2 3 . 7 2 1 5 2 1 c 6 1 7 2 9 . 3 4 8 3 1 4 5 . 9 4 4 1 3 2 3 7 5 . 9 1 6 4 6 1 c 7 1 7 2 9 . 3 3 5 7 5 9 3 8 . 9 0 7 2 8 4 3 4 2 . 9 5 3 3 0 9 c 8 1 7 2 9 . 3 1 5 7 6 3 9 1 . 4 0 6 7 8 2 2 9 0 . 4 5 3 8 0 9 c 9 1 7 2 9 . 3 3 1 0 1 7 5 1 . 3 5 7 4 0 6 3 3 0 . 5 0 3 1 9 1 c 1 0 1 7 2 9 . 3 3 4 3 4 5 4 2 . 6 1 9 7 4 2 3 3 9 . 2 4 0 8 5 1 c 1 1 1 7 2 9 . 3 2 8 6 4 7 5 7 . 5 7 9 8 4 1 3 4 9 . 3 2 5 4 1 2 c 1 2 1 7 2 9 . 3 3 3 8 2 8 4 3 . 9 7 7 1 2 5 3 3 7 . 8 8 3 4 7 1

    從表3可以看出,在所有穩(wěn)定構(gòu)型中,c1的總能量最低,相互作用能量最大,證明c1的穩(wěn)定性最好.

    2.5 電荷分析

    在離子液體穩(wěn)定構(gòu)型基礎(chǔ)上,對(duì)其進(jìn)行頻率分析,得到[C4SO3HPy]HSO4中所有原子電荷,如表4所示.

    從表4可以看出,陽(yáng)離子功能基團(tuán)磺酸基中氧原子(25O、26O、27O)上帶有部分負(fù)電荷,而與之相連的硫原子上帶有部分正電荷,離子液體陰陽(yáng)離子間發(fā)生相互作用,從計(jì)算得到的各原子的電荷布局分析可以預(yù)測(cè)陽(yáng)離子磺酸根上的氧原子在與其它物質(zhì)作用時(shí)更容易參與配位;對(duì)于氫原子除了陰離子HSO4中34H及陽(yáng)離子功能基團(tuán)中的28H外,吡啶環(huán)中的11H也有較高的電荷,有較強(qiáng)形成分子內(nèi)氫鍵的趨勢(shì),同樣表現(xiàn)出一定的酸性.由酸量測(cè)定結(jié)果可知,由于11H呈現(xiàn)酸性,故多消耗了NaOH.

    從模擬結(jié)果結(jié)合電荷分析可以看出,離子液體陰陽(yáng)離子間發(fā)生相互作用,由功能基團(tuán)及陰離子中的O原子提供孤對(duì)電子與H原子形成氫鍵.氫鍵越強(qiáng),氫鍵構(gòu)型越穩(wěn)定,此時(shí)形成氫鍵的H原子的電荷密度較高(c1中28H、34H的電荷密度較其它氫原子的大).從而進(jìn)一步說明了對(duì)于不同類型離子液體氫原子酸性的強(qiáng)弱與其穩(wěn)定性和電子密度有關(guān),其穩(wěn)定性和電子密度越低時(shí),酸性越強(qiáng).對(duì)于一類離子液體,當(dāng)離子液體中所形成的氫鍵被削弱時(shí),H原子的酸性則會(huì)增強(qiáng).因此可以通過削弱離子對(duì)中所形成的氫鍵,使氫更容易解離下來(lái),從而增強(qiáng)離子液體的酸性.

    3 結(jié)論

    1)對(duì)N-(4-磺酸丁基)吡啶硫酸氫鹽 [C4SO3HPy]HSO4的理論分析表明,以兩性離子a1,陽(yáng)離子b1,離子液體c1結(jié)構(gòu)最穩(wěn)定,能量最低,其能量分別為 1029.123306au、 1029.506576au、 1029.350578au.

    在審美由城市主導(dǎo)、偽文藝開始崩塌的時(shí)代,滋生于民間的土味文化,其野性、真誠(chéng)和賣力正在刺激大批城市人口以旁觀者的獵奇心理加入其中,逐步形成新一輪的審丑文化潮流。

    2)在兩性離子中,形成了含有分子內(nèi)四氫鍵、三氫鍵及不含分子內(nèi)氫鍵的穩(wěn)定結(jié)構(gòu),當(dāng)兩性離子中的SO3接受一個(gè)氫變?yōu)殛?yáng)離子之后,形成了含有分子內(nèi)三氫鍵及不含分子內(nèi)氫鍵的穩(wěn)定結(jié)構(gòu),由于氫原子的引入,分子內(nèi)的氫鍵無(wú)論從鍵數(shù)和鍵能上均被削弱.結(jié)果表明,在兩性離子、陽(yáng)離子、離子液體的穩(wěn)定構(gòu)型中,同一體系不同穩(wěn)定構(gòu)型,氫鍵數(shù)目越多,能量越大,構(gòu)型越穩(wěn)定.

    3)離子液體最穩(wěn)定構(gòu)型是陰離子與陽(yáng)離子功能基團(tuán)發(fā)生相互作用時(shí)形成的,而當(dāng)它們之間形成氫鍵時(shí),HSO3中的氫更易與HSO4中的氧結(jié)合,最后以兩性離子與H2SO4相互作用的形式存在.

    4)通過電荷分析,離子液體中陰離子及陽(yáng)離子功能基團(tuán)上的氫體現(xiàn)酸性,吡啶環(huán)中的11H也有一定的酸性.對(duì)于磺酸類功能化離子液體,當(dāng)離子液體中氫鍵被削弱時(shí),離子液體的酸性將會(huì)增強(qiáng).為對(duì)其進(jìn)行酸性調(diào)變并找到結(jié)構(gòu)對(duì)催化活性和目標(biāo)產(chǎn)物選擇性的影響規(guī)律提供理論基礎(chǔ).

    表4 離子液體 [C4SO3HPy]HSO4的Mulliken電荷布局分析Tab.4 Mulliken atomic chargesand population analysisof[C4SO3HPy]HSO4

    參考文獻(xiàn):

    [1]Seddon KR.Ionic Liquids for Clean Technology[J].Journalof Chem ical Technology&Biotechnology,1997,68(4):351-356.

    [2]Welton T.Room-temperature ionic liquids.Solvents for synthesisand catalysis[J].Chemical reviews,1999,99(8):2071-2084.

    [3]張鎖江,劉曉敏,姚曉倩,等.離子液體的前沿、進(jìn)展及應(yīng)用 [J].中國(guó)科學(xué)(B輯:化學(xué)),2009(10):1134-1144.

    [4]ColeAC,Jensen JL,Ntai I,etal.NovelBr?nsted acidic ionic liquidsand theiruseasdualsolvent-catalysts[J].Journalof theAmericanChem ical Society,2002,124(21):5962-5963.

    [5]Welton T.Ionic liquids in catalysis[J].Coordination Chemistry Reviews,2004,248(21):2459-2477.

    [6]Now icki J,MuszynskiM,Gryglew icz S.Novelbasic ionic liquids from cyclic guanidinesand amidines-new catalysts for transesterification of oleochem icals[J].Journalof Chem ical Technology and Biotechnology,2014,89(1):48-55.

    [7]Olivier-Bourbigou H,Magna L,Morvan D.Ionic liquidsand catalysis:recentprogress from know ledge toapplications[J].Applied CatalysisA:General,2010,373(1):1-56.

    [8]Wasserscheid P,Welton T.Ionic liquids in synthesis[M/OL].Wiley Online Library,2008.

    [9]Plechkova NV,Seddon KR.Applicationsof ionic liquids in the chem ical industry[J].Chem icalSociety Reviews,2008,37(1):123-150.

    [10]ZhangWL,Yue Y,SuWY,etal.Metalchloridesor sulfuric acid in ionic liquid solvents convertcatechol to p-tert-Butylcatechol[J].Catalysis Communications,2015,65:113-116.

    [11]TurnerEA,Pye CC,Singer RD.Useofab initio calculations toward the rationaldesign of room temperature ionic liquids[J].Journalof Physical Chem istry A,2003,107(13):2277-2288.

    [12]Wang Y,LiH,Han S.A theoreticalinvestigationof the interactionsbetweenwatermoleculesand ionic liquids[J].Journalof PhysicalChemistry B,2006,110(48):24646-24651.

    [13]Ding K,Zhang S,WangD,etal.Hydrogenbondsin im idazolium ionic liquids[J].JournalofPhysicalChemistry A,2006,110(31):9775-9782.

    [14]UmebayashiY,M itsugiT,F(xiàn)ujimoriT,eral.Lithium ion solvationn in room-temperature ionic liquids involving bis(trifluoromethanesulfonyl)im ideanion studied by Raman spectroscopy and DFT calculation[J].The Journalof PhysicalChem istry B,2007,111(45):13028-13032.

    [15]FredlakeCP,Crosthwaite JM,HertDG,etal.Thermophysicalpropertiesof imidazolium-based ionic liquids[J].JournalofChemical&Engineering Data,2004,49(4):954-964.

    [16]Urahata SM,Ribeiro MCC.Structure of ionic liquidsof 1-alkyl-3-methylim idazolium cations:A systematic computer simulation study[J].The Journalof Chem icalPhysics,2004,120(4):1855-1863.

    [17]Weingartner H.Understanding ionic liquids at themolecular level:facts,problems,and controversies[J].Angewandte Chem ie International Edition,2008,47(4):654-670.

    [18]蘇培峰,譚凱,吳安安,等.理論與計(jì)算化學(xué)研究進(jìn)展 [J].廈門大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,50(2):311-318.

    [19]魯萍,譚振江.對(duì)量子化學(xué)計(jì)算軟件性能的分析及對(duì)比 [J].吉林師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,32(2):107-109.

    [20]韓曉剛,廖奕.Gaussian和Multiw fn軟件在結(jié)構(gòu)化學(xué)教學(xué)中的應(yīng)用 [J].首都師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,33(5):23-26.

    [21]ZhangWL,DuW,Wang XW,etal.Progress in adjustmentof ionic liquidsacidity based on thealkylation[J].Advanced MaterialsResearch,2013,634-638:581-586.

    [責(zé)任編輯 田 豐 夏紅梅]

    Theoreticalstudy of N-(4-sulfonic acid butyl)pyridine bisulfate ionic liquids

    ZHANGWenlin,QIAN Haozhi,CHEN Yao,LIJunfei,LIChunli

    (Schoolof Chem ical Engineering,HebeiUniversity of Technology,Tianjin 300130,China)

    The precursor compound[C4SO3Py]±,cation[C4SO3HPy]+and the ion pairsof SO3H functionalized acidic ionic liquids[C4SO3HPy]HSO4werestudied systematically by theDensity FunctionalTheory(DFT)at6-31+G*level. Thestructureoptim ization,frequency analysisand energy calculationwereobtained,then interactionenergywascalculated and atom ic chargewasanalyzed.The resultsindicated thattheintramolecularhydrogen bondsand intermolecularhydrogen bonds single existed or co-existed in the stable configuration of the zw itterion,cation and ionic liquid.The strongest interaction in the ionic liquidwas found between theanion and the functionalgroup.When the strong intramolecularhydrogen bondswere formed,the hydrogen of HSO3was prone to combined w ith oxygen atom of HSO4.The final existing form was the interaction between zw itterion and H2SO4.Charge analysis showed that the hydrogen in anion and cation functioned groupswere acidic,besides,the hydrogen in the pyridine ring also show a certain acidic charge.

    ionic liquids;structuralanalysis;DFT;hydrogen bond;algorithm

    TE667

    A

    1007-2373(2016)04-0045-08

    10.14081/j.cnki.hgdxb.2016.04.008

    2016-03-10

    河北省高等學(xué)??茖W(xué)技術(shù)研究項(xiàng)目(ZD2015118)

    張文林(1968-),男(漢族),教授,博士.

    猜你喜歡
    兩性離子吡啶氫鍵
    教材和高考中的氫鍵
    兩性離子聚合物/多元醇復(fù)合井壁強(qiáng)化劑的研制與作用機(jī)理
    吡啶-2-羧酸鉻的制備研究
    云南化工(2021年10期)2021-12-21 07:33:28
    勘 誤
    今日農(nóng)業(yè)(2019年11期)2019-08-13 00:49:02
    甜菜堿類兩性離子共聚物涂層的構(gòu)建及釋藥行為研究
    二水合丙氨酸復(fù)合體內(nèi)的質(zhì)子遷移和氫鍵遷移
    含吡啶的大環(huán)席夫堿錳(Ⅱ)配合物:合成、表征及抗菌性質(zhì)
    銥(Ⅲ)卟啉β-羥乙與基醛的碳?xì)滏I活化
    兩性離子催化劑的研究進(jìn)展
    岛国毛片在线播放| 久久狼人影院| 亚洲三区欧美一区| 欧美日韩黄片免| 精品人妻1区二区| 欧美av亚洲av综合av国产av| 丁香欧美五月| 男女床上黄色一级片免费看| 在线天堂中文资源库| 男女之事视频高清在线观看| 色精品久久人妻99蜜桃| 18禁黄网站禁片午夜丰满| 两人在一起打扑克的视频| 久9热在线精品视频| 久久天躁狠狠躁夜夜2o2o| 免费观看a级毛片全部| 无人区码免费观看不卡 | 人人妻人人添人人爽欧美一区卜| 久久久久精品国产欧美久久久| 两个人看的免费小视频| 国产一区二区 视频在线| 色综合婷婷激情| 五月天丁香电影| 午夜两性在线视频| 午夜视频精品福利| 女人被躁到高潮嗷嗷叫费观| 亚洲专区中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人精品巨大| 十八禁网站网址无遮挡| 国产成人精品无人区| 99久久精品国产亚洲精品| 国产精品 欧美亚洲| 丝袜美足系列| 曰老女人黄片| 欧美日韩av久久| 他把我摸到了高潮在线观看 | 成年女人毛片免费观看观看9 | 人人妻人人澡人人看| 欧美日韩亚洲国产一区二区在线观看 | 国产精品国产av在线观看| 免费看a级黄色片| 熟女少妇亚洲综合色aaa.| 久久午夜亚洲精品久久| a在线观看视频网站| 在线观看免费日韩欧美大片| 精品国产一区二区久久| 国产亚洲欧美精品永久| 亚洲欧洲精品一区二区精品久久久| 看免费av毛片| 免费在线观看视频国产中文字幕亚洲| 老司机在亚洲福利影院| 每晚都被弄得嗷嗷叫到高潮| 色在线成人网| 久久毛片免费看一区二区三区| 在线天堂中文资源库| 香蕉丝袜av| 狠狠狠狠99中文字幕| 18禁裸乳无遮挡动漫免费视频| 午夜福利影视在线免费观看| 国产精品香港三级国产av潘金莲| 18在线观看网站| 国产精品一区二区在线不卡| 三上悠亚av全集在线观看| 亚洲熟妇熟女久久| 99久久国产精品久久久| 国产精品久久久久成人av| 波多野结衣av一区二区av| 香蕉丝袜av| 亚洲第一欧美日韩一区二区三区 | 99re在线观看精品视频| 亚洲欧洲日产国产| 亚洲情色 制服丝袜| 成年人午夜在线观看视频| 亚洲成人免费电影在线观看| 精品亚洲乱码少妇综合久久| 久久精品国产亚洲av香蕉五月 | 日韩视频一区二区在线观看| 母亲3免费完整高清在线观看| 免费在线观看黄色视频的| 人妻久久中文字幕网| 精品国产乱子伦一区二区三区| 丰满少妇做爰视频| 国产视频一区二区在线看| videos熟女内射| 99国产精品免费福利视频| 中亚洲国语对白在线视频| 国产片内射在线| 欧美成人午夜精品| 五月开心婷婷网| 热99re8久久精品国产| 大片免费播放器 马上看| 国产男靠女视频免费网站| 黄频高清免费视频| 亚洲情色 制服丝袜| 12—13女人毛片做爰片一| 国产欧美亚洲国产| 久久这里只有精品19| 99九九在线精品视频| 国产精品久久久av美女十八| 十八禁人妻一区二区| 黑人欧美特级aaaaaa片| 国产深夜福利视频在线观看| 男女免费视频国产| 欧美久久黑人一区二区| 亚洲专区中文字幕在线| 肉色欧美久久久久久久蜜桃| 免费av中文字幕在线| 国产日韩欧美视频二区| 国产男女超爽视频在线观看| svipshipincom国产片| 亚洲精品成人av观看孕妇| 日韩欧美三级三区| 免费不卡黄色视频| 女性被躁到高潮视频| 国产精品久久久人人做人人爽| 91国产中文字幕| 啦啦啦在线免费观看视频4| 自线自在国产av| 精品国产乱码久久久久久男人| 2018国产大陆天天弄谢| 亚洲精品中文字幕在线视频| 性色av乱码一区二区三区2| 国产男女内射视频| 女性被躁到高潮视频| 欧美 日韩 精品 国产| 国产精品免费大片| 欧美大码av| 美女国产高潮福利片在线看| 久久久国产欧美日韩av| 久久天堂一区二区三区四区| 高清黄色对白视频在线免费看| 可以免费在线观看a视频的电影网站| 精品国产乱子伦一区二区三区| av免费在线观看网站| 999久久久国产精品视频| 黄色视频不卡| 成人黄色视频免费在线看| 人人妻人人添人人爽欧美一区卜| 久久久久网色| 婷婷成人精品国产| 国产又爽黄色视频| 俄罗斯特黄特色一大片| 在线观看免费日韩欧美大片| av片东京热男人的天堂| 日日摸夜夜添夜夜添小说| 亚洲欧美激情在线| 国产精品久久久av美女十八| 国产一区二区三区综合在线观看| 一个人免费在线观看的高清视频| 久久久久久久大尺度免费视频| 亚洲成人国产一区在线观看| 国产成人精品无人区| 首页视频小说图片口味搜索| 中文字幕高清在线视频| 黄网站色视频无遮挡免费观看| 人妻 亚洲 视频| 国内毛片毛片毛片毛片毛片| 欧美日韩视频精品一区| 黑人欧美特级aaaaaa片| 亚洲五月色婷婷综合| 国产精品国产av在线观看| 亚洲国产欧美一区二区综合| 狂野欧美激情性xxxx| 亚洲欧美一区二区三区久久| 欧美+亚洲+日韩+国产| 国产精品香港三级国产av潘金莲| 欧美激情久久久久久爽电影 | 最新在线观看一区二区三区| 久久国产亚洲av麻豆专区| 蜜桃国产av成人99| 他把我摸到了高潮在线观看 | 国产精品一区二区精品视频观看| 欧美av亚洲av综合av国产av| 老司机影院毛片| 一进一出抽搐动态| 色视频在线一区二区三区| 精品国产亚洲在线| 久久中文看片网| 男女边摸边吃奶| 午夜成年电影在线免费观看| 久久精品国产a三级三级三级| 欧美黑人欧美精品刺激| 狂野欧美激情性xxxx| 黄网站色视频无遮挡免费观看| 午夜福利免费观看在线| 欧美大码av| 久久精品国产a三级三级三级| 亚洲黑人精品在线| 日日夜夜操网爽| 午夜福利视频在线观看免费| 国产熟女午夜一区二区三区| av免费在线观看网站| 亚洲精品国产精品久久久不卡| 欧美一级毛片孕妇| 国产欧美日韩一区二区三区在线| 久久99一区二区三区| 黄色怎么调成土黄色| 午夜福利视频在线观看免费| 波多野结衣一区麻豆| 不卡av一区二区三区| 深夜精品福利| 国产精品 欧美亚洲| 男人舔女人的私密视频| 黄片大片在线免费观看| 亚洲自偷自拍图片 自拍| 欧美黄色片欧美黄色片| 亚洲伊人久久精品综合| 欧美日韩精品网址| 99国产综合亚洲精品| 国产高清国产精品国产三级| 亚洲国产看品久久| 精品国产乱码久久久久久小说| 高清欧美精品videossex| 美女高潮喷水抽搐中文字幕| 亚洲欧洲日产国产| 欧美精品高潮呻吟av久久| 中文欧美无线码| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| 久久久精品94久久精品| 搡老岳熟女国产| 中文欧美无线码| 亚洲 欧美一区二区三区| 99国产极品粉嫩在线观看| 黑丝袜美女国产一区| 女人高潮潮喷娇喘18禁视频| 亚洲精品乱久久久久久| 亚洲五月色婷婷综合| 黄片小视频在线播放| 久久精品亚洲av国产电影网| 成人国语在线视频| 老司机在亚洲福利影院| 1024香蕉在线观看| 国产国语露脸激情在线看| 女人爽到高潮嗷嗷叫在线视频| 午夜日韩欧美国产| 窝窝影院91人妻| 国产又爽黄色视频| 亚洲中文av在线| 亚洲伊人久久精品综合| 一边摸一边抽搐一进一出视频| 欧美老熟妇乱子伦牲交| 久久九九热精品免费| 久久久久久久久免费视频了| 久久久久久久国产电影| 精品人妻熟女毛片av久久网站| 久久午夜亚洲精品久久| 免费在线观看日本一区| 亚洲国产成人一精品久久久| 视频区欧美日本亚洲| 成人国产av品久久久| 韩国精品一区二区三区| 免费高清在线观看日韩| 免费不卡黄色视频| 美女国产高潮福利片在线看| 欧美老熟妇乱子伦牲交| 嫁个100分男人电影在线观看| 国产日韩一区二区三区精品不卡| 亚洲国产成人一精品久久久| 亚洲av电影在线进入| 亚洲成国产人片在线观看| 777久久人妻少妇嫩草av网站| 欧美黑人欧美精品刺激| 一本色道久久久久久精品综合| 精品国产超薄肉色丝袜足j| 中文字幕人妻丝袜制服| 黑人猛操日本美女一级片| 久久人人97超碰香蕉20202| 色精品久久人妻99蜜桃| 大片免费播放器 马上看| 国产欧美日韩一区二区三| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看| 日本精品一区二区三区蜜桃| 国产欧美日韩综合在线一区二区| 欧美日韩国产mv在线观看视频| 老司机亚洲免费影院| 亚洲免费av在线视频| 黄网站色视频无遮挡免费观看| 成人免费观看视频高清| 免费不卡黄色视频| 变态另类成人亚洲欧美熟女 | 中文字幕人妻熟女乱码| 亚洲av片天天在线观看| av在线播放免费不卡| 国产精品成人在线| 日本av手机在线免费观看| 在线亚洲精品国产二区图片欧美| 久久久久精品人妻al黑| av又黄又爽大尺度在线免费看| 国产精品久久久久久精品古装| 俄罗斯特黄特色一大片| 2018国产大陆天天弄谢| 日韩中文字幕欧美一区二区| 19禁男女啪啪无遮挡网站| 色婷婷av一区二区三区视频| 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 午夜老司机福利片| 欧美精品啪啪一区二区三区| 91麻豆精品激情在线观看国产 | 人妻 亚洲 视频| 免费观看a级毛片全部| 午夜福利影视在线免费观看| 欧美午夜高清在线| 精品国产一区二区久久| 国产精品久久久久成人av| 亚洲视频免费观看视频| 国产精品美女特级片免费视频播放器 | 青青草视频在线视频观看| 精品人妻在线不人妻| 国产精品亚洲av一区麻豆| 国产深夜福利视频在线观看| 99香蕉大伊视频| av不卡在线播放| 国产淫语在线视频| 嫁个100分男人电影在线观看| 99久久99久久久精品蜜桃| 99re在线观看精品视频| 久久久久国产一级毛片高清牌| www.自偷自拍.com| 热99re8久久精品国产| 亚洲国产欧美日韩在线播放| 99精品久久久久人妻精品| 在线观看www视频免费| 1024视频免费在线观看| 亚洲欧美一区二区三区久久| 精品乱码久久久久久99久播| e午夜精品久久久久久久| 一级a爱视频在线免费观看| 男女床上黄色一级片免费看| 性少妇av在线| 黄色视频在线播放观看不卡| 精品人妻熟女毛片av久久网站| 夜夜夜夜夜久久久久| 精品少妇一区二区三区视频日本电影| 亚洲男人天堂网一区| 免费不卡黄色视频| 国产视频一区二区在线看| 免费高清在线观看日韩| 久久久精品94久久精品| www.自偷自拍.com| 亚洲精华国产精华精| 91精品国产国语对白视频| 久久精品熟女亚洲av麻豆精品| 日本av免费视频播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 国产aⅴ精品一区二区三区波| 亚洲国产看品久久| 国产精品 欧美亚洲| 一区福利在线观看| 亚洲欧洲精品一区二区精品久久久| 精品一品国产午夜福利视频| 黑人猛操日本美女一级片| 亚洲精品久久成人aⅴ小说| 高清av免费在线| 99re在线观看精品视频| 色综合欧美亚洲国产小说| 免费黄频网站在线观看国产| 亚洲第一欧美日韩一区二区三区 | 国产av又大| 不卡一级毛片| 久久午夜亚洲精品久久| 在线亚洲精品国产二区图片欧美| 日本wwww免费看| 久久精品国产a三级三级三级| 久久这里只有精品19| 黑人巨大精品欧美一区二区蜜桃| 99久久精品国产亚洲精品| 97在线人人人人妻| 一本色道久久久久久精品综合| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 国产精品久久久av美女十八| 亚洲欧洲日产国产| 国产精品免费一区二区三区在线 | 亚洲午夜理论影院| 精品国产乱码久久久久久男人| 香蕉久久夜色| 男女下面插进去视频免费观看| 亚洲午夜精品一区,二区,三区| 人成视频在线观看免费观看| 国产精品久久久久久精品电影小说| 国产精品久久久久久人妻精品电影 | 成年人免费黄色播放视频| 日日爽夜夜爽网站| 精品人妻熟女毛片av久久网站| 欧美大码av| 久久中文看片网| 国产精品一区二区免费欧美| 国产在线免费精品| 在线观看舔阴道视频| 香蕉丝袜av| 欧美乱码精品一区二区三区| 满18在线观看网站| 国产精品99久久99久久久不卡| 亚洲精品自拍成人| 999久久久精品免费观看国产| 男人操女人黄网站| 国产精品影院久久| 99热网站在线观看| 午夜福利欧美成人| 黑人猛操日本美女一级片| 18禁国产床啪视频网站| 两个人看的免费小视频| 80岁老熟妇乱子伦牲交| 熟女少妇亚洲综合色aaa.| 两个人免费观看高清视频| av又黄又爽大尺度在线免费看| 国产欧美日韩一区二区三区在线| bbb黄色大片| 97在线人人人人妻| 欧美日韩亚洲国产一区二区在线观看 | 久久精品国产a三级三级三级| 国产成人免费观看mmmm| 欧美精品高潮呻吟av久久| 熟女少妇亚洲综合色aaa.| 欧美午夜高清在线| 国产精品电影一区二区三区 | 日韩一卡2卡3卡4卡2021年| 亚洲精品在线观看二区| 自拍欧美九色日韩亚洲蝌蚪91| 999久久久国产精品视频| 日本vs欧美在线观看视频| 国产成人系列免费观看| 黑人猛操日本美女一级片| 国产伦理片在线播放av一区| 啦啦啦中文免费视频观看日本| 久久精品人人爽人人爽视色| 免费在线观看视频国产中文字幕亚洲| 人成视频在线观看免费观看| 日韩中文字幕欧美一区二区| 操出白浆在线播放| 波多野结衣av一区二区av| 国产男女内射视频| 国产伦理片在线播放av一区| 黄色毛片三级朝国网站| 91av网站免费观看| 黄色成人免费大全| 一区在线观看完整版| 国产在线一区二区三区精| 夜夜夜夜夜久久久久| 亚洲国产看品久久| 9色porny在线观看| 露出奶头的视频| 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 欧美日韩黄片免| 亚洲国产av新网站| 国产淫语在线视频| 精品一区二区三区av网在线观看 | 91老司机精品| 国产男女超爽视频在线观看| a级毛片黄视频| 亚洲人成电影观看| 嫩草影视91久久| av又黄又爽大尺度在线免费看| av视频免费观看在线观看| 国产精品一区二区在线不卡| 亚洲精品粉嫩美女一区| 老司机深夜福利视频在线观看| 久久天躁狠狠躁夜夜2o2o| 老司机靠b影院| 亚洲成a人片在线一区二区| 男女无遮挡免费网站观看| 两个人看的免费小视频| 免费高清在线观看日韩| 777久久人妻少妇嫩草av网站| 香蕉国产在线看| www.自偷自拍.com| 2018国产大陆天天弄谢| 精品少妇黑人巨大在线播放| 一本综合久久免费| 18在线观看网站| 18禁观看日本| 高清av免费在线| 国产高清激情床上av| 黄色a级毛片大全视频| 日韩视频一区二区在线观看| av免费在线观看网站| 午夜精品久久久久久毛片777| 色在线成人网| 黑人猛操日本美女一级片| 国产一区二区三区综合在线观看| 欧美人与性动交α欧美软件| 一区二区三区精品91| 老熟女久久久| 一级毛片电影观看| 亚洲精品美女久久av网站| 亚洲精品国产一区二区精华液| 亚洲第一青青草原| 王馨瑶露胸无遮挡在线观看| 一级,二级,三级黄色视频| 老熟妇仑乱视频hdxx| 亚洲精品乱久久久久久| 美女高潮喷水抽搐中文字幕| 国产亚洲精品一区二区www | 国产aⅴ精品一区二区三区波| 欧美日韩中文字幕国产精品一区二区三区 | 免费黄频网站在线观看国产| 国产精品一区二区在线不卡| 国产男女超爽视频在线观看| 少妇 在线观看| 99久久国产精品久久久| 国产又色又爽无遮挡免费看| 精品国内亚洲2022精品成人 | 成年人黄色毛片网站| 黄色 视频免费看| 狂野欧美激情性xxxx| 亚洲,欧美精品.| 午夜91福利影院| 麻豆av在线久日| 自线自在国产av| 亚洲欧洲日产国产| 亚洲av日韩精品久久久久久密| 曰老女人黄片| 日韩有码中文字幕| 搡老岳熟女国产| 午夜免费鲁丝| 黄频高清免费视频| 久久婷婷成人综合色麻豆| 可以免费在线观看a视频的电影网站| a级片在线免费高清观看视频| 大片免费播放器 马上看| 日韩人妻精品一区2区三区| 中文字幕制服av| 老司机午夜福利在线观看视频 | 久久久久国产一级毛片高清牌| 欧美黑人欧美精品刺激| 亚洲va日本ⅴa欧美va伊人久久| 又大又爽又粗| 国产精品影院久久| www.熟女人妻精品国产| 国产一区二区三区在线臀色熟女 | 午夜91福利影院| 电影成人av| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久| 五月开心婷婷网| 九色亚洲精品在线播放| 欧美日韩一级在线毛片| 久久精品aⅴ一区二区三区四区| 真人做人爱边吃奶动态| 久久久久久久精品吃奶| 亚洲国产成人一精品久久久| 国产精品久久电影中文字幕 | 99国产精品一区二区蜜桃av | 在线观看一区二区三区激情| 亚洲国产中文字幕在线视频| 麻豆乱淫一区二区| 黄片小视频在线播放| 精品一区二区三区av网在线观看 | 狠狠婷婷综合久久久久久88av| 老司机深夜福利视频在线观看| 天堂中文最新版在线下载| 超碰成人久久| 亚洲黑人精品在线| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 欧美 日韩 精品 国产| 黄色片一级片一级黄色片| 涩涩av久久男人的天堂| 一本—道久久a久久精品蜜桃钙片| 久久亚洲精品不卡| 青青草视频在线视频观看| av有码第一页| 亚洲自偷自拍图片 自拍| 国产精品 国内视频| 妹子高潮喷水视频| 18禁黄网站禁片午夜丰满| 午夜福利视频在线观看免费| 国产男女内射视频| 亚洲av美国av| 久久青草综合色| 丝袜美腿诱惑在线| 97人妻天天添夜夜摸| 欧美精品一区二区免费开放| 亚洲 国产 在线| av又黄又爽大尺度在线免费看| 他把我摸到了高潮在线观看 | 国产一区二区激情短视频| 咕卡用的链子| 电影成人av| 久久婷婷成人综合色麻豆| 国产区一区二久久| 人人妻人人爽人人添夜夜欢视频| av不卡在线播放| 国产国语露脸激情在线看| 亚洲成人国产一区在线观看| 国产精品熟女久久久久浪| 国产国语露脸激情在线看| 黄片小视频在线播放| 欧美日韩国产mv在线观看视频| 男女之事视频高清在线观看| av免费在线观看网站| 欧美成狂野欧美在线观看| 久久精品国产a三级三级三级| 亚洲欧洲日产国产| 丝袜美腿诱惑在线| 欧美 日韩 精品 国产| 黄色丝袜av网址大全| 黄频高清免费视频| 亚洲全国av大片| 国产一区二区 视频在线| 欧美在线一区亚洲| 人人妻,人人澡人人爽秒播| 人人妻人人添人人爽欧美一区卜| 电影成人av| 波多野结衣av一区二区av| 黄色毛片三级朝国网站| 大香蕉久久成人网| 在线观看免费视频网站a站| 一区二区三区乱码不卡18|