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

    基于密度泛函理論的尼克酸分子振動(dòng)光譜分析

    2016-06-30 03:01:01陶亞萍韓禮剛韓運(yùn)俠劉照軍
    光散射學(xué)報(bào) 2016年2期
    關(guān)鍵詞:密度泛函理論

    陶亞萍,韓禮剛,韓運(yùn)俠,劉照軍

    (洛陽(yáng)師范學(xué)院物理與電子信息學(xué)院,洛陽(yáng) 471022)

    基于密度泛函理論的尼克酸分子振動(dòng)光譜分析

    陶亞萍,韓禮剛,韓運(yùn)俠,劉照軍*

    (洛陽(yáng)師范學(xué)院物理與電子信息學(xué)院,洛陽(yáng)471022)

    摘要:實(shí)驗(yàn)測(cè)量了尼克酸分子的拉曼和紅外光譜,用B3LYP混合泛函和cc-PVDZ基組計(jì)算了尼克酸分子的平衡構(gòu)型、振動(dòng)頻率、拉曼和紅外強(qiáng)度。采用GAR2PED程序?qū)δ峥怂岱肿舆M(jìn)行了簡(jiǎn)正振動(dòng)分析,依據(jù)所得勢(shì)能分布對(duì)尼克酸分子的振動(dòng)頻率進(jìn)行了理論歸屬,彌補(bǔ)了以往文獻(xiàn)缺少對(duì)其振動(dòng)模式貢獻(xiàn)進(jìn)行定量研究的不足,提供了更多的振動(dòng)光譜信息。

    關(guān)鍵詞:尼克酸;密度泛函理論;振動(dòng)光譜;簡(jiǎn)正振動(dòng)分析;頻率歸屬

    1引言

    尼克酸(Niacin,C6H5NO2)為白色或淡黃色結(jié)晶性粉末,屬于水溶性B族維生素,與煙酰胺一起合稱為維生素PP。最初它是通過(guò)氧化尼古丁得到的,因此也稱為尼古丁酸。從結(jié)構(gòu)上看,它是由吡啶在3號(hào)位置連接一個(gè)羧基組成,亦稱為吡啶-3-甲酸?;瘜W(xué)上一般把具有生物活性的全部吡啶-3-甲酸及其衍生物統(tǒng)稱為煙酸。尼克酸在自然界中分布很廣,尤其在動(dòng)物肝臟、酵母、蛋黃、谷物、花生、魚(yú)類和綠色蔬菜中含量較高。它被廣泛用于食品、化學(xué)和醫(yī)藥等行業(yè)中。在食品工業(yè)中,可用做去味劑、防腐劑、保色劑和保鮮劑等;在飼料加工行業(yè)中,可提高雞鴨的產(chǎn)蛋率、奶牛產(chǎn)奶量及魚(yú)、牛、羊等禽畜的抗病能力和肉的質(zhì)量;在日用化學(xué)工業(yè)中,它可用于合成頭發(fā)生長(zhǎng)促進(jìn)劑、染發(fā)助劑、聚合物穩(wěn)定劑、感光材料的抗氧化劑與抗灰霧劑、電鍍光亮劑及PVC塑料的熱穩(wěn)定劑等;在醫(yī)藥行業(yè)中,它有促進(jìn)細(xì)胞新陳代謝、擴(kuò)張血管、緩解末稍血管痙攣的效果并且有保護(hù)皮膚中的膠原纖維的功效,可用于防治糙皮病和類似的維生素缺乏癥、動(dòng)脈硬化等癥,它作為醫(yī)藥中間體可用于合成治療中樞性呼吸及循環(huán)系統(tǒng)衰竭癥、偏頭痛、高血酯癥、冠心病等疾病的煙酸鉻、煙酰胺、尼克剎米、煙酸肌醇脂、癸煙酯等數(shù)十種藥物[1-3]。因此對(duì)尼克酸的研究具有重要的現(xiàn)實(shí)意義。

    目前有關(guān)尼克酸的研究主要集中在其含量的色譜法測(cè)定、晶體結(jié)構(gòu)[4-5]、二聚體結(jié)構(gòu)與性質(zhì)的分析方面[6],同時(shí)對(duì)尼克酸的表面增強(qiáng)拉曼光譜也有較多的研究[7-10]。振動(dòng)光譜是研究分子結(jié)構(gòu)的重要手段,而密度泛函理論可以從理論上對(duì)分子的振動(dòng)光譜進(jìn)行預(yù)測(cè),因此密度泛函理論在分子光譜的研究領(lǐng)域中應(yīng)用越來(lái)越普遍[11-12]。在尼克酸的振動(dòng)光譜分析方面,一直備受關(guān)注,已有一些相關(guān)的研究,比如:2001年,O.Sala測(cè)量了尼克酸水溶液1800~200 cm-1波段的拉曼光譜,并對(duì)此波段的振動(dòng)模式進(jìn)行了初步的指認(rèn)[13];2003年,P.Koczoń給出了尼克酸固體粉末3500~3000和1800~1300 cm-1波段的紅外、拉曼光譜的實(shí)驗(yàn)頻率及理論歸屬[14];2004年,莫衛(wèi)民等人對(duì)尼克酸的紅外光譜實(shí)驗(yàn)結(jié)果進(jìn)行了初步分析[15];2006年,Li-Ran Wang測(cè)量了尼克酸溶液1700~400 cm-1波段的拉曼光譜,并指出當(dāng)溶液呈中性或弱堿性時(shí),尼克酸以離子形式存在于溶液中[16];2009年,韓禮剛也曾結(jié)合Gauss View[17]軟件顯示的分子振動(dòng)的動(dòng)畫(huà)對(duì)尼克酸的拉曼光譜的振動(dòng)頻率進(jìn)行了初步歸屬[18];2011年,M.Kumar和R.A.Yadav對(duì)尼克酸及其N-氧化物的振動(dòng)光譜進(jìn)行了對(duì)比分析,也結(jié)合Gauss View軟件給出了經(jīng)驗(yàn)的歸屬[19]。雖然相關(guān)研究較多,但未見(jiàn)到結(jié)合密度泛函理論與簡(jiǎn)正振動(dòng)分析程序?qū)δ峥怂岬睦⒓t外光譜的全部波段進(jìn)行全面歸屬的文獻(xiàn),前人的討論均未能定量的給出每種基團(tuán)特征振動(dòng)模式的貢獻(xiàn)百分比?;诖?在實(shí)驗(yàn)測(cè)定尼克酸分子的紅外和拉曼光譜的基礎(chǔ)上,結(jié)合密度泛函理論和簡(jiǎn)正振動(dòng)模式分解的分析方法,來(lái)獲得各基團(tuán)特征振動(dòng)模式的勢(shì)能分布,從而對(duì)尼克酸分子的振動(dòng)光譜進(jìn)行了比較全面的歸屬。

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

    實(shí)驗(yàn)中用于測(cè)試的尼克酸樣品購(gòu)于北京化工廠,純度大于98%,未經(jīng)進(jìn)一步提純。紅外光譜采用的是KBr壓片制作,使用BIO-RAD 60 V 型傅立葉變換紅外光譜儀測(cè)量,光譜分辨率為1 cm-1。拉曼光譜使用LabRam HR800型共焦顯微拉曼光譜儀測(cè)量,激光光源為氬離子激光器的514.5 nm,到達(dá)樣品表面的功率為3 mW,50倍物鏡,1800 gr/mm的高分辨率閃耀衍射光柵,數(shù)據(jù)采集時(shí)間為30秒。

    3計(jì)算方法

    尼克酸是由吡啶在3號(hào)位被羧酸基團(tuán)取代而成,雖然只有1個(gè)取代基團(tuán),但由于羧酸基團(tuán)相對(duì)于吡啶環(huán)、羥基基團(tuán)相對(duì)于羰基基團(tuán)朝向的不同,尼克酸會(huì)出現(xiàn)4種構(gòu)型。為了找到尼克酸分子的最穩(wěn)定構(gòu)型,采用Gaussian 09[20]程序包、密度泛函理論方法(DFT)和B3LYP/6-311G(d,p)基組,計(jì)算了4種構(gòu)型的能量和頻率,發(fā)現(xiàn)除了(d)構(gòu)型有虛頻外,其它三種構(gòu)型均無(wú)虛頻。圖1是計(jì)算所得尼克酸分子4種構(gòu)型的示意圖。每種構(gòu)型對(duì)應(yīng)的能量分別為:(a)-436.865910 a.u.;(b)-436.865554 a.u.;(c)-436.853413 a.u.;(d)-436.855636 a.u.,通過(guò)比較,可以看到構(gòu)型(a)能量最低、最穩(wěn)定,李會(huì)學(xué)等人[6]通過(guò)過(guò)渡態(tài)的分析也確定了構(gòu)型(a)最穩(wěn)定,同時(shí)M.P.Gupta[4]、W.B.Wright[5]等人報(bào)道的晶體結(jié)構(gòu)也得出了同樣的結(jié)論,故此本文以下分析和討論是基于構(gòu)型(a)展開(kāi)的。

    Fig.1 The optimized equilibrium geometry of Niacin

    為了對(duì)簡(jiǎn)正振動(dòng)頻率做出正確的歸屬,采用更大的基組B3LYP/cc-PVDZ對(duì)構(gòu)型(a)進(jìn)行了進(jìn)一步的結(jié)構(gòu)優(yōu)化和頻率計(jì)算。尼克酸分子共有14個(gè)原子,42個(gè)笛卡爾坐標(biāo)和36個(gè)簡(jiǎn)正振動(dòng)頻率。根據(jù)分子的對(duì)稱性和P.Pulay提出的方法[21-22],構(gòu)造出一套共36個(gè)獨(dú)立、完備的局域?qū)ΨQ坐標(biāo),如表1所示。最后由GAR2PED 軟件[23]讀取Gaussian程序的輸出文件,得到每個(gè)簡(jiǎn)正振動(dòng)模式在各對(duì)稱坐標(biāo)上的勢(shì)能分布百分比,判斷出每種基團(tuán)特征振動(dòng)模式的貢獻(xiàn),從而完成各簡(jiǎn)正頻率的歸屬指認(rèn)。

    Tab.1 Local symmetry coordinates of Niacin

    4結(jié)果與討論

    4.1尼克酸分子的平衡構(gòu)型

    據(jù)我們所知,尼克酸分子的結(jié)構(gòu)在氣相中尚未確定,但固體狀態(tài)的晶體結(jié)構(gòu)已經(jīng)由M.P.Gupta等人[4]進(jìn)行了報(bào)道,理論和實(shí)驗(yàn)結(jié)構(gòu)參數(shù)的對(duì)比如表2所示。尼克酸是一個(gè)含氮原子的六元雜環(huán)化合物,是吡啶衍生物的一種。由于存在氮原子,使得分子中的六元環(huán)狀結(jié)構(gòu)略顯不規(guī)則,C1-N6、C5-N6鍵長(zhǎng)較C1-C2、 C2-C3、C3-C4、C4-C5略短;鍵角C4-C5-N6、N6-C1-C2略大于C1-C2-C3、C2-C3-C4、C3-C4-C5、C5-N6-C1。換言之,吡啶環(huán)上圍繞著氮原子的鍵長(zhǎng)和鍵角與苯環(huán)相比,均有較大變化。這一方面是因?yàn)榈幼陨淼挠绊?另一方面是因?yàn)榱u基與氮原子形成O-H…N分子間氫鍵的緣故。

    忽略羧基,吡啶環(huán)上C-C、C-N鍵長(zhǎng)的理論值與實(shí)驗(yàn)值的均方根偏差(RMSD)是0.019 ?,鍵角的偏差是0.72°??紤]羧基,尼克酸分子鍵長(zhǎng)的均方根偏差是0.028 ?,鍵角的偏差是1.82°。偏差最大的兩個(gè)角是C4-C7-O8和C4-C7-O9,這是因?yàn)橛捎诜肿娱g的相互作用,在固態(tài)中尼克酸的羧基會(huì)形成O-H…O氫鍵,會(huì)對(duì)羧基的鍵長(zhǎng)和鍵角產(chǎn)生較大的影響,這也是造成均方根偏差明顯增大的主要原因。同時(shí)C4-C7單鍵的鍵長(zhǎng)實(shí)驗(yàn)值比標(biāo)準(zhǔn)的C-C單鍵鍵長(zhǎng)(1.54 ?)明顯縮短,C7-O9鍵長(zhǎng)的實(shí)驗(yàn)值明顯比理論值小,亦可證明O-H…O氫鍵的存在。正是由于固體狀態(tài)中分子間會(huì)形成O-H…O和O-H…N兩種氫鍵,而理論計(jì)算的是單分子的氣態(tài)平衡構(gòu)型,未考慮到氫鍵的影響,所以理論計(jì)算的結(jié)果與實(shí)驗(yàn)數(shù)據(jù)會(huì)稍微有一些偏差。從表2的計(jì)算結(jié)果可以看到,所有的二面角都是0.0° 或180.0°,即羧酸基團(tuán)和吡啶環(huán)是共面的,表明尼克酸分子是平面結(jié)構(gòu)。

    Tab.2 Geometry parameters of Niacin

    4.2振動(dòng)頻率和歸屬

    為了對(duì)光譜圖進(jìn)行直觀的比較,我們根據(jù)文獻(xiàn)[24-25]把Gaussian09計(jì)算所得的拉曼活性Si轉(zhuǎn)換為譜峰的強(qiáng)度Ii,其理論和實(shí)驗(yàn)光譜的比較如圖2和3所示。計(jì)算結(jié)果表明尼克酸分子的36個(gè)簡(jiǎn)正振動(dòng)由25個(gè)面內(nèi)振動(dòng)和11個(gè)面外振動(dòng)組成,具有Cs對(duì)稱性,所有振動(dòng)模既有拉曼活性又有紅外活性。為了使密度泛函理論獲得的理論值與實(shí)驗(yàn)值更加吻合,J.P.Merrick 等人[26]在充分考慮電子結(jié)構(gòu)、非諧性等影響的基礎(chǔ)上,通過(guò)1066個(gè)理論頻率與相應(yīng)實(shí)驗(yàn)頻率的對(duì)比計(jì)算和分析,得出了B3LYP/cc-PVDZ 算法的頻率校正因子為0.9717。于是筆者把尼克酸分子的理論計(jì)算頻率與該校正因子相乘,經(jīng)過(guò)校正后,計(jì)算值與實(shí)驗(yàn)觀測(cè)值較好的吻合,其中存在的少許差別可能來(lái)自于理論模擬的是單分子理想化模型,而實(shí)驗(yàn)測(cè)量的樣品存在分子間相互作用,如表3所示。同時(shí)表3顯示了由簡(jiǎn)正振動(dòng)分析軟件得到的各種基團(tuán)特征振動(dòng)模式的貢獻(xiàn)百分比(勢(shì)能占比大于5%的振動(dòng))。

    通常自由的O-H鍵伸縮振動(dòng)應(yīng)該出現(xiàn)在3500 cm-1以上[27],但在尼克酸中O-H鍵與C=O雙鍵是相連的,在羰基和分子間形成的O-H…O氫鍵的雙重影響下,O-H鍵伸縮的振動(dòng)頻率會(huì)減小。 因此在本研究中,我們將位于3416 cm-1的紅外峰歸屬為O-H伸縮振動(dòng),是比較合理的。尼克酸分子屬于單取代雜環(huán)芳香烴分子,因此它應(yīng)該有4個(gè)C-H伸縮振動(dòng),且它們均屬于典型的含氮雜環(huán)芳香烴的C-H伸縮振動(dòng),應(yīng)該出現(xiàn)在3100~3010 cm-1區(qū)域[27],它們的振動(dòng)相對(duì)獨(dú)立,與其它鍵伸縮或彎曲振動(dòng)之間耦合較小,而本研究中位于3088、 3073、3042 cm-1的拉曼峰和位于3097、3071、3041 cm-1的紅外峰被指認(rèn)為C-H伸縮振動(dòng),與文獻(xiàn)報(bào)道的范圍基本一致;位于1708 cm-1的紅外強(qiáng)峰和1694 cm-1的拉曼強(qiáng)峰,從勢(shì)能分布的數(shù)值可以看到C=O伸縮振動(dòng)的比例高達(dá)92%,因此被歸屬為C=O伸縮振動(dòng),這與文獻(xiàn)[28]報(bào)道的羧基的特征基團(tuán)C=O伸縮振動(dòng)通常位于1715~1680 cm-1是一致的。事實(shí)上,M.Kumar和R.A.Yadav[19]根據(jù)Gauss View軟件的動(dòng)畫(huà)將位于3414 cm-1的紅外強(qiáng)峰歸屬為尼克酸分子的O-H伸縮振動(dòng),將位于3089、3073、3041 cm-1的拉曼峰與位于3099、3075、3040 cm-1的紅外峰歸屬為C-H伸縮振動(dòng),將位于1711 cm-1的紅外強(qiáng)峰和1695 cm-1的拉曼強(qiáng)峰歸屬為C=O伸縮振動(dòng),這與我們的歸屬都是一致的,可見(jiàn)對(duì)這種相對(duì)獨(dú)立、耦合較小的振動(dòng),歸屬并不存在爭(zhēng)議,差別在于我們是通過(guò)特征基團(tuán)貢獻(xiàn)百分比的數(shù)值進(jìn)行了定量化的歸屬。

    分子的振動(dòng)是分子中所有原子共同參與的簡(jiǎn)正振動(dòng),很多情況下表現(xiàn)為多個(gè)振動(dòng)耦合的形式。對(duì)于尼克酸分子的振動(dòng),在頻率較低的部分(<1700 cm-1),這種耦合現(xiàn)象就很普遍。M.Kumar和R.A.Yadav[19]指出對(duì)于這種強(qiáng)烈耦合振動(dòng)的歸屬是非常困難的,因此他們只給出了部分經(jīng)驗(yàn)歸屬,未展現(xiàn)耦合振動(dòng)部分。而本文依據(jù)簡(jiǎn)正振動(dòng)分析方法對(duì)尼克酸分子的振動(dòng)頻率能夠進(jìn)行全面的歸屬,展現(xiàn)了耦合振動(dòng)過(guò)程中各基團(tuán)振動(dòng)的貢獻(xiàn)百分比。如:從表3很容易看出來(lái)位于1417、1309、1187、1085 cm-1處的拉曼峰和位于1416、1297、1186、1088 cm-1處的紅外峰歸屬于C-H面內(nèi)彎曲振動(dòng),其中含有C-C伸縮振動(dòng)的成分;位于1002、975、954、812 cm-1處的拉曼峰和954、813 cm-1處的紅外峰歸屬于C-H面外彎曲振動(dòng);同時(shí)位于1600、1589 cm-1處的拉曼峰和1599、1583 cm-1處的紅外峰歸屬于C-C伸縮振動(dòng),這與以往文獻(xiàn)[29]報(bào)道環(huán)的C-C伸縮振動(dòng)一般應(yīng)出現(xiàn)在1625~1430 cm-1范圍內(nèi)是一致的;位于1321、1242 cm-1處的拉曼峰和1322、1243 cm-1處的紅外峰歸屬于C-N伸縮振動(dòng),亦耦合有C-C伸縮振動(dòng)成分;位于1115 cm-1處的拉曼和紅外峰歸屬于C-O伸縮振動(dòng);位于1031、751 cm-1處的拉曼峰和1039 cm-1處的紅外峰歸屬于環(huán)面內(nèi)變形,695、387 cm-1處的拉曼峰和695 cm-1處的紅外峰歸屬于環(huán)面外變形振動(dòng)。與P.Koczoń[16]的理論歸屬相比,差別主要在苯環(huán)上的C-C和C-N伸縮振動(dòng)的歸屬上。P.Koczoń[14]把環(huán)上各鍵的伸縮振動(dòng)統(tǒng)一稱為環(huán)伸縮振動(dòng),沒(méi)有具體給出環(huán)上哪個(gè)鍵的伸縮振動(dòng)所占比例更多一些,而我們給出的勢(shì)能分布結(jié)果很清晰的給出了各個(gè)振動(dòng)模式所占的比例。

    5結(jié)論

    用量子化學(xué)方法計(jì)算了尼克酸分子的4種可能構(gòu)型和相應(yīng)的能量,通過(guò)能量的對(duì)比、文獻(xiàn)調(diào)查和實(shí)驗(yàn)的佐證確定了尼克酸分子的最穩(wěn)定構(gòu)型,計(jì)算了最穩(wěn)定構(gòu)型的簡(jiǎn)正振動(dòng)頻率和相應(yīng)的紅外與拉曼強(qiáng)度,并對(duì)各簡(jiǎn)正振動(dòng)頻率進(jìn)行了全面歸屬。尼克酸分子有兩個(gè)氫鍵受體(-N、C=O)和一個(gè)供體(-OH),在固體狀態(tài)下存在O-H…O和O-H…N兩種氫鍵,會(huì)對(duì)其結(jié)構(gòu)和光譜產(chǎn)生一定的影響。M.Kumar和R.A.Yadav[19]曾經(jīng)指出,因?yàn)榄h(huán)上各鍵振動(dòng)之間存在耦合,與取代基的振動(dòng)也存在耦合,因此對(duì)尼克酸的吡啶環(huán)振動(dòng)的歸屬以及與環(huán)相連的C-H鍵振動(dòng)的歸屬都是非常困難的。而本文依據(jù)的簡(jiǎn)正振動(dòng)分析方法,通過(guò)比較勢(shì)能分布百分比的方式,可以容易地了解由哪種特征振動(dòng)模式所主導(dǎo),這種定量方法比以往的經(jīng)驗(yàn)歸屬更準(zhǔn)確。P.Koczoń[14]雖然也采用了簡(jiǎn)正振動(dòng)分析方法,但它只給出了部分波段的歸屬,而且沒(méi)有具體給出吡啶環(huán)上各振動(dòng)的貢獻(xiàn)比例。本文利用簡(jiǎn)正振動(dòng)分析方法給出了尼克酸分子全部譜帶的勢(shì)能分布,并對(duì)其進(jìn)行了全面的歸屬。尼克酸廣泛用于食品、化學(xué)和醫(yī)藥等行業(yè)中,對(duì)該物質(zhì)分子結(jié)構(gòu)的分析和振動(dòng)頻率的全面歸屬,為與其相關(guān)的藥學(xué)、生物學(xué)等方面的深入研究提供了重要的參考價(jià)值。

    Tab.3 Calculated and observed vibrational frequencies with PED of Niacin

    Fig.2 Infrared spectra of Niacin

    Fig.3 Raman spectra of Niacin

    參考文獻(xiàn)

    [1]趙東江,馬松艷.煙酸的合成方法與開(kāi)發(fā)前景[J].江蘇化工,2005,33(1):53-57.(Zhao Dongxiang,Ma Songyan.Synthesis Methods and Development Prospects of Nicotinic Acid [J].Jiangsu Chem Ind,2005,33(1):53-57.)

    [2]路百程.煙酸的合成與應(yīng)用[J].遼寧化工,1994,3:9-12.(Lu Baicheng,Synthesis and Application of Nicotinic Acid [J].Liaoning Chem Ind,1994,3:9-12.)

    [3]李艷云,尹振晏,宮彩紅.煙酸的研究進(jìn)展[J].北京石油化工學(xué)院學(xué)報(bào),2006,14(1):59-64.(Li Yanyun,Yin Zhenyan,Gong Caihong.Development of Nicotinic Acid [J].J Beijing Ins.Petrochem Tech,2006,14(1):59-64.)

    [4]Gupta M P,Kumar P.Nicotinic Acid:a Redetermination [J].Cryst Struct Commun,1975,4:365.

    [5]Wright W B,King G S.The crystal structure of nicotinic acid [J].Acta Crystal,1953,6(4):305-317.

    [6]李會(huì)學(xué),王曉峰,董小寧,等.煙酸二聚體的結(jié)構(gòu)與性質(zhì)[J].物理化學(xué)學(xué)報(bào),2009,25(1):161-168.(Li Huixue,Wang Xiaofeng,Dong Xiaoning.etal.Structures and Properties of Nicotinic Acid Dimer [J].Acta Phys-Chim Sin,2009,25(1):161-168.)

    [7]Park S M,Kim K.Adsorption of picolinic and niacins on a silver sol surface investigated by Raman spectroscopy [J].J Mol Stru,1995,344:195-203.

    [8]Baia M,Astilean S.Raman and Sers Investigations of Pharmaceuticals [M].Traian Iliescu Published by Springer,2008:137-142.

    [9]Zhang L S,Fang Y,Zhang P.Experimental and DFT theoretical studies of SERS effect on gold nanowires array [J].Chem Phy Lett,2008,451(1):102-105.

    [10]Wen R,Fang Y.Adsorption of pyridine carboxylic acids on silver surface investigated by potential-dependentSERS [J].Vib Spec,2005,39(1):106-113.

    [11]陶亞萍,吳義芳,韓禮剛,等.檸檬酸分子振動(dòng)光譜研究[J].光散射學(xué)報(bào),2015,27(3):256-262.(Tao Yaping,Wu Yifang,Han Ligang,etal.Investigation on the vibrational spectra of citric acid [J].J Light Scatter,2015,27(3):256-262)

    [12]陶亞萍,韓禮剛,梁會(huì)琴,等.2-巰基-5-甲基-1,3,4-噻二唑分子的拉曼、紅外光譜和簡(jiǎn)正振動(dòng)分析[J].光散射學(xué)報(bào),2012,24(4):381-387.(Tao Yaping,Han Ligang,Liang Huiqin,etal.Raman,FT-IR Spectra and Normal Mode Analysis of 2-Mercapto-5-methyl-1,3,4-thiadiazole [J].J Light Scatter,2012,24(4):381-387.)

    [13]Sala O,Goncalves N S,Noda L K.Vibrational analysis of Nicotinic species based on ab initio molecular orbital calculations [J].J Mol Stru,2001,566:411-416.

    [14]Koczoń P,Dobrowolski J C,Lewandowski W,etal.Experimental and theoretical IR and Raman spectra of picolinic,nicotinic and isoniacins [J].J Mol Stru,2003,655(1):89-95.

    [15]莫衛(wèi)民,薛建,胡寶祥,等.紅外光譜在煙酸系列化合物制備過(guò)程中的應(yīng)用[J].光譜實(shí)驗(yàn)室,2004,21(4):677-679.(Mo Weimin,Xue Jian,Hu Baoxiang,etal.The Application of FT-IR in the Synthesis of Nicotinic Derivatives [J].J Spec Lab,2004,21(4):677-679.)

    [16]Wang L R,Fang Y.Experimental(UV-Raman)and theoretical(DFT)study of pyridine carboxylic acid in aqueous solution [J].Chem Phy,2006,323(2):376-382.

    [17]Dennington R,Keith T,Millam J.GaussView 5.0.9[CP].Semichem Inc,Shawnee Mission,KS,2009.

    [18]韓禮剛.煙酸分子的簡(jiǎn)正振動(dòng)與拉曼光譜[J].洛陽(yáng)師范學(xué)院學(xué)報(bào),2009,28(2):44-46.(Han Ligang.Normal Vibration and Raman Spectrum of Nicotinic Acid [J].J Luoyang Norm Univ,2009,28(2):44-46.)

    [19]Kumar M,Yadav R A.Experimental IR and Raman spectra and quantum chemical studies of molecular structures,conformers and vibrational characteristics of nicotinic acid and its N-oxide [J].Spectrochim Acta Part A,2011,79(5):1316-1325.

    [20]Frisch M J,Trucks G W,Schlegel H B,etal.Gaussian 09 D.01[CP].Gaussian Inc,Wallingford CT,2013.

    [21]Pulay P,Fogarasi G,Pang F,etal.Systematic ab Initio Gradient Calculation of Molecular Geometries,Force Constants,and Dipole Moment Derivatives [J].J Amer Chem Soc,1979,101(10):2550-2560.

    [22]Pulay P,Paizs B.Newtonian molecular dynamics in general curvilinear internal coordinates [J].Chem Phy Lett,2002,353(5):400-406.

    [23]Martin J M,Alsenoy C V.Gar2ped [CP].University of Antwerp,1995.

    [24]Krishnakumar V,Surumbarkuzhali N,Muthunatesan S.Scaled quantum chemical studies on the vibrational spectra of 4-bromo benzonitrile [J].Spectrochim Acta,Part A,2009,71(5):1810-1813.

    [25]Chalmers J M,Griffiths P R.Handbook of Vibrational Spectroscopy [M].John Wiley & Sons Ltd,2002,71.

    [26]Merrick J P,Moran D,Radom L.An evaluation of harmonic vibrational frequency scale factors [J].J Phy Chem A,2007,111(45):11683-11700.

    [27]Socrates G.Infrared and Raman characteristic group frequencies,tables and charts [M].Wiley Chichester,2001:68-235.

    [28]Monicka J C,James C.FT-Raman and FTIR spectra,DFT investigation of the structure and vibrational assignment of mefenacet [J].J Mol Stru,2015,1095:1-7.

    [29]Tao Y,Han L,Han Y,etal.A combined experimental and theoretical analysis on molecular structure and vibrational spectra of 2,4-dihydroxybenzoic acid [J].Spectrochim Acta,Part A,2015,137:1078-1085.

    Analysis of Vibrational Spectra of Niacin Based on Density Functional Theory

    TAO Ya-ping,HAN Li-gang,HAN Yun-xia,LIU Zhao-jun*

    (CollegeofPhysicsandElectronicInformation,Luoyangnormaluniversity,Luoyang,471022)

    Abstract:The Raman and infrared spectra of niacin were measured experimentally.The molecular equilibrium geometries,vibrational frequencies,Raman and infrared intensity of niacin were calculated using the B3LYP level and cc-PVDZ basis set.Normal mode analysis was carried out using the program GAR2PED and the assignment of fundamental vibrations for niacin was obtained according to the potential energy distributions(PED).This study provides us with more quantitative vibrational spectral information which was not mentioned in previous literatures.

    Key words:niacin;DFT;vibrational spectrum;normal mode analysis;frequencies assignments

    收稿日期:2015-07-30; 修改稿日期:2015-12-21

    基金項(xiàng)目:NSFC-河南人才培養(yǎng)聯(lián)合基金(U1204109)河南省科技計(jì)劃項(xiàng)目(142102210109)

    作者簡(jiǎn)介:陶亞萍(1980-),女,內(nèi)蒙古通遼人,講師,主要從事分子光譜的理論和實(shí)驗(yàn)研究。E-mail:taoyaping2001@gmail.com 通訊作者:劉照軍(1966-),男,教授,從事拉曼光譜實(shí)驗(yàn)研究。E-mail:zhaojunliu@gmail.com

    文章編號(hào):1004-5929(2016)02-0175-07

    中圖分類號(hào):O657.37

    文獻(xiàn)標(biāo)志碼:A

    doi:10.13883/j.issn1004-5929.201602014

    猜你喜歡
    密度泛函理論
    理論計(jì)算Mn3Ga合金的電子結(jié)構(gòu)
    第一性原理研究鐵銅合金催化氨硼烷水解脫氫
    Cu/FePt/MgO多層體系的結(jié)構(gòu)和磁性研究
    氣相中Ni+活化環(huán)己烷C—H和C—C鍵的理論研究
    苯環(huán)間的扭轉(zhuǎn)角對(duì)分子器件整流特性的影響
    一種鐵配合物催化的烯烴氫化反應(yīng)機(jī)理的研究
    電子理論在材料科學(xué)中的應(yīng)用
    科技傳播(2016年8期)2016-07-13 23:46:24
    InnN (n=1~13) 團(tuán)簇的幾何結(jié)構(gòu)和電子性質(zhì)
    科技視界(2015年35期)2016-01-04 09:39:10
    重金屬離子對(duì)致癌性NDMA形成影響的實(shí)驗(yàn)與理論研究
    金屬離子的配位作用影響芳胺抗氧化性能的理論研究
    潤(rùn)滑油(2014年5期)2015-01-27 10:02:16
    国产男女超爽视频在线观看| 亚洲精品日韩在线中文字幕| 亚洲经典国产精华液单| 国产精品99久久久久久久久| 成人免费观看视频高清| 国产成人精品久久久久久| 国内少妇人妻偷人精品xxx网站| 在线观看免费日韩欧美大片 | 国产色爽女视频免费观看| 少妇精品久久久久久久| 久久这里有精品视频免费| 免费少妇av软件| 午夜免费男女啪啪视频观看| 在线观看av片永久免费下载| 成人影院久久| 干丝袜人妻中文字幕| 在线播放无遮挡| 91久久精品电影网| 久久久久人妻精品一区果冻| 亚洲精品456在线播放app| 色5月婷婷丁香| 中文字幕av成人在线电影| 2018国产大陆天天弄谢| 99久国产av精品国产电影| 亚洲国产精品国产精品| 国产成人精品久久久久久| 各种免费的搞黄视频| 精品久久久精品久久久| 在线观看人妻少妇| 精品国产露脸久久av麻豆| 精品久久久久久电影网| 久久久色成人| 日韩一区二区视频免费看| 日韩欧美 国产精品| 又粗又硬又长又爽又黄的视频| 国内精品宾馆在线| 80岁老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 夫妻午夜视频| 久久久久性生活片| 尾随美女入室| 亚洲精品,欧美精品| 偷拍熟女少妇极品色| 天堂俺去俺来也www色官网| 亚洲欧美成人精品一区二区| 欧美日韩精品成人综合77777| 亚洲精品久久午夜乱码| 欧美高清性xxxxhd video| 国产精品.久久久| 精品国产三级普通话版| 男人和女人高潮做爰伦理| 欧美97在线视频| 欧美xxⅹ黑人| 国产精品免费大片| 热re99久久精品国产66热6| 国产伦理片在线播放av一区| 精品人妻视频免费看| 青青草视频在线视频观看| 国产精品一二三区在线看| 国产成人一区二区在线| 亚洲成人中文字幕在线播放| 王馨瑶露胸无遮挡在线观看| 国产免费又黄又爽又色| 中文字幕av成人在线电影| 丰满乱子伦码专区| 99久国产av精品国产电影| 国产爽快片一区二区三区| 另类亚洲欧美激情| 人体艺术视频欧美日本| 亚洲欧美一区二区三区黑人 | 我要看黄色一级片免费的| 国产一区有黄有色的免费视频| 天堂中文最新版在线下载| 五月玫瑰六月丁香| 亚洲精品第二区| 欧美变态另类bdsm刘玥| 最黄视频免费看| 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 欧美最新免费一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 免费观看的影片在线观看| 精品久久久噜噜| 两个人的视频大全免费| 国产亚洲最大av| 国产精品一二三区在线看| 内射极品少妇av片p| 亚洲成人中文字幕在线播放| 欧美成人a在线观看| 精品人妻一区二区三区麻豆| 亚洲精品自拍成人| 亚洲aⅴ乱码一区二区在线播放| 丝瓜视频免费看黄片| 99视频精品全部免费 在线| 久久久久视频综合| 少妇 在线观看| 91狼人影院| 精品久久久久久久久亚洲| av在线播放精品| 国产爽快片一区二区三区| 色网站视频免费| 亚洲国产欧美在线一区| 欧美丝袜亚洲另类| 女性生殖器流出的白浆| 国产精品一及| 麻豆乱淫一区二区| 99精国产麻豆久久婷婷| 91久久精品国产一区二区三区| 久久99热6这里只有精品| 免费看光身美女| 国内揄拍国产精品人妻在线| 中文字幕精品免费在线观看视频 | 色视频在线一区二区三区| 中文字幕免费在线视频6| 黑丝袜美女国产一区| 国产一区亚洲一区在线观看| 热re99久久精品国产66热6| 伦理电影大哥的女人| 99久久精品热视频| 精品人妻熟女av久视频| 亚洲国产精品999| 精品国产乱码久久久久久小说| 亚洲国产最新在线播放| 人人妻人人看人人澡| 国产真实伦视频高清在线观看| 91精品国产国语对白视频| 乱码一卡2卡4卡精品| 边亲边吃奶的免费视频| 成人18禁高潮啪啪吃奶动态图 | 美女福利国产在线 | 国产一区二区三区av在线| 91久久精品国产一区二区成人| 国产大屁股一区二区在线视频| 边亲边吃奶的免费视频| 一二三四中文在线观看免费高清| 国产淫语在线视频| 少妇被粗大猛烈的视频| 精品人妻熟女av久视频| 三级经典国产精品| 大码成人一级视频| 18禁动态无遮挡网站| 精品国产露脸久久av麻豆| 有码 亚洲区| 婷婷色av中文字幕| 寂寞人妻少妇视频99o| 狂野欧美激情性xxxx在线观看| 特大巨黑吊av在线直播| 亚洲av欧美aⅴ国产| 男女国产视频网站| 日韩 亚洲 欧美在线| 亚洲av在线观看美女高潮| 亚洲av二区三区四区| 精品午夜福利在线看| 99热全是精品| 精品熟女少妇av免费看| 亚洲精品久久午夜乱码| 联通29元200g的流量卡| 婷婷色综合大香蕉| 精品人妻熟女av久视频| 亚洲欧美日韩东京热| 中国三级夫妇交换| 国产 精品1| 在线观看一区二区三区激情| 国产探花极品一区二区| 纵有疾风起免费观看全集完整版| 亚州av有码| 97热精品久久久久久| 秋霞伦理黄片| www.色视频.com| 熟女人妻精品中文字幕| 波野结衣二区三区在线| 精品久久国产蜜桃| 久久精品国产鲁丝片午夜精品| 91久久精品电影网| 久久韩国三级中文字幕| 夜夜爽夜夜爽视频| 成人国产麻豆网| 亚洲国产色片| 国产精品一区二区在线观看99| 日韩av免费高清视频| 国产人妻一区二区三区在| 青春草国产在线视频| 色婷婷av一区二区三区视频| 高清av免费在线| 最近中文字幕高清免费大全6| 国产黄色视频一区二区在线观看| 亚洲美女搞黄在线观看| 国产精品.久久久| 日韩中文字幕视频在线看片 | 国产高潮美女av| 亚洲精品久久午夜乱码| 国产av一区二区精品久久 | 麻豆乱淫一区二区| 国产成人精品久久久久久| 国产精品女同一区二区软件| 久久婷婷青草| 免费播放大片免费观看视频在线观看| 一级黄片播放器| 黄色一级大片看看| 国产亚洲最大av| 丝瓜视频免费看黄片| 中文字幕av成人在线电影| 久久韩国三级中文字幕| 国产精品爽爽va在线观看网站| 免费在线观看成人毛片| 身体一侧抽搐| 欧美成人一区二区免费高清观看| 不卡视频在线观看欧美| 麻豆国产97在线/欧美| 欧美日韩亚洲高清精品| 精品熟女少妇av免费看| 嫩草影院入口| 伊人久久精品亚洲午夜| 丰满乱子伦码专区| 成人美女网站在线观看视频| 亚洲欧美成人综合另类久久久| 熟女人妻精品中文字幕| 精品熟女少妇av免费看| 六月丁香七月| 亚洲最大成人中文| 久久久久人妻精品一区果冻| 久久精品熟女亚洲av麻豆精品| 看非洲黑人一级黄片| 乱系列少妇在线播放| 成人免费观看视频高清| 色视频在线一区二区三区| 黄色欧美视频在线观看| 国产成人免费观看mmmm| 亚洲人成网站在线播| 亚洲精品乱码久久久久久按摩| videossex国产| 中国美白少妇内射xxxbb| 一本—道久久a久久精品蜜桃钙片| 国产视频内射| 在线观看av片永久免费下载| av不卡在线播放| 一个人看视频在线观看www免费| 爱豆传媒免费全集在线观看| 老熟女久久久| 欧美日韩在线观看h| 精品一区二区三区视频在线| 你懂的网址亚洲精品在线观看| 日本av手机在线免费观看| 国产淫语在线视频| 各种免费的搞黄视频| 香蕉精品网在线| 精品人妻偷拍中文字幕| 久久久色成人| 国产精品久久久久久久电影| 18禁在线无遮挡免费观看视频| 亚洲精品一二三| 色吧在线观看| 免费黄网站久久成人精品| 日日啪夜夜撸| 水蜜桃什么品种好| 国产免费福利视频在线观看| 亚洲婷婷狠狠爱综合网| 精品酒店卫生间| 成人影院久久| 80岁老熟妇乱子伦牲交| 高清av免费在线| 高清在线视频一区二区三区| 亚洲欧美成人综合另类久久久| a级毛色黄片| 男女免费视频国产| 春色校园在线视频观看| 欧美日韩综合久久久久久| 人人妻人人添人人爽欧美一区卜 | 国产成人免费观看mmmm| 亚洲图色成人| 国内揄拍国产精品人妻在线| 18禁动态无遮挡网站| 国产黄色免费在线视频| 国产美女午夜福利| 永久免费av网站大全| 成人无遮挡网站| 欧美激情国产日韩精品一区| 成人漫画全彩无遮挡| 精品视频人人做人人爽| 国产欧美另类精品又又久久亚洲欧美| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人中文字幕在线播放| 亚洲精品中文字幕在线视频 | 欧美bdsm另类| 欧美成人精品欧美一级黄| 婷婷色麻豆天堂久久| 女性被躁到高潮视频| 色网站视频免费| 国产av国产精品国产| 日韩在线高清观看一区二区三区| 男女国产视频网站| av视频免费观看在线观看| av又黄又爽大尺度在线免费看| 高清毛片免费看| 街头女战士在线观看网站| 国产淫片久久久久久久久| 欧美区成人在线视频| 久久精品久久精品一区二区三区| 日韩国内少妇激情av| 激情 狠狠 欧美| 在线观看免费视频网站a站| 免费不卡的大黄色大毛片视频在线观看| 成年av动漫网址| 亚洲美女黄色视频免费看| 亚洲国产精品成人久久小说| 亚洲精品aⅴ在线观看| 日韩免费高清中文字幕av| 国产黄片美女视频| 欧美成人a在线观看| 国产成人freesex在线| 青春草亚洲视频在线观看| 久久精品熟女亚洲av麻豆精品| 国产精品国产三级国产av玫瑰| 国产欧美日韩一区二区三区在线 | 国内少妇人妻偷人精品xxx网站| 我要看黄色一级片免费的| 国产男女超爽视频在线观看| av国产久精品久网站免费入址| 少妇高潮的动态图| 国产免费视频播放在线视频| 国产白丝娇喘喷水9色精品| 国产黄片视频在线免费观看| 尾随美女入室| 天堂8中文在线网| 国产又色又爽无遮挡免| 噜噜噜噜噜久久久久久91| 秋霞伦理黄片| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 成人影院久久| 国产毛片在线视频| 色视频在线一区二区三区| 欧美老熟妇乱子伦牲交| 18禁在线播放成人免费| 在线精品无人区一区二区三 | 成人亚洲欧美一区二区av| 一二三四中文在线观看免费高清| 亚洲色图av天堂| 亚洲,一卡二卡三卡| 十八禁网站网址无遮挡 | 久热久热在线精品观看| 国产一区二区在线观看日韩| 熟妇人妻不卡中文字幕| 亚洲国产av新网站| 精品亚洲成a人片在线观看 | 久久亚洲国产成人精品v| 伦理电影免费视频| 欧美日本视频| 18禁在线播放成人免费| 精品少妇黑人巨大在线播放| 亚洲,一卡二卡三卡| 日韩,欧美,国产一区二区三区| 观看美女的网站| 国产伦在线观看视频一区| 亚洲国产精品一区三区| 亚洲在久久综合| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品电影小说 | 日本vs欧美在线观看视频 | 欧美日韩综合久久久久久| 六月丁香七月| 全区人妻精品视频| 一区二区三区乱码不卡18| 欧美区成人在线视频| 亚洲国产最新在线播放| 成年美女黄网站色视频大全免费 | 久久国产亚洲av麻豆专区| 不卡视频在线观看欧美| 久久韩国三级中文字幕| 欧美区成人在线视频| 性色avwww在线观看| 久久国产精品大桥未久av | 涩涩av久久男人的天堂| av又黄又爽大尺度在线免费看| 国产一区有黄有色的免费视频| 一本色道久久久久久精品综合| 亚洲国产毛片av蜜桃av| 有码 亚洲区| 国产伦精品一区二区三区四那| 一级a做视频免费观看| av专区在线播放| 夫妻午夜视频| 久久av网站| 午夜免费鲁丝| h日本视频在线播放| 好男人视频免费观看在线| 色网站视频免费| 性色avwww在线观看| 全区人妻精品视频| 简卡轻食公司| 国产亚洲5aaaaa淫片| 啦啦啦啦在线视频资源| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 99久久中文字幕三级久久日本| 一个人免费看片子| 国产亚洲午夜精品一区二区久久| 国产伦理片在线播放av一区| 亚洲欧美日韩卡通动漫| 国产男女超爽视频在线观看| 亚洲四区av| 三级国产精品欧美在线观看| 亚洲欧美日韩另类电影网站 | 国产 精品1| 国产黄色免费在线视频| 国产精品久久久久久久久免| 最黄视频免费看| av线在线观看网站| 精品午夜福利在线看| 成年女人在线观看亚洲视频| 成人18禁高潮啪啪吃奶动态图 | 欧美成人精品欧美一级黄| 久久久久精品久久久久真实原创| 国产亚洲av片在线观看秒播厂| 99热网站在线观看| 久久精品熟女亚洲av麻豆精品| a级毛片免费高清观看在线播放| 欧美日本视频| 91久久精品国产一区二区成人| av天堂中文字幕网| 亚洲精品日本国产第一区| 永久免费av网站大全| 97精品久久久久久久久久精品| 自拍偷自拍亚洲精品老妇| 国产精品久久久久久久久免| 99热网站在线观看| 精品一品国产午夜福利视频| 免费大片18禁| kizo精华| 2018国产大陆天天弄谢| 毛片一级片免费看久久久久| 国语对白做爰xxxⅹ性视频网站| 久久婷婷青草| 高清日韩中文字幕在线| 特大巨黑吊av在线直播| 国产黄频视频在线观看| 青春草国产在线视频| 亚洲激情五月婷婷啪啪| 亚洲av免费高清在线观看| 精品熟女少妇av免费看| 精品久久久精品久久久| 久久国产精品大桥未久av | 色视频在线一区二区三区| 日本欧美视频一区| 亚洲国产精品专区欧美| 美女视频免费永久观看网站| 亚洲精品视频女| 在线 av 中文字幕| 99精国产麻豆久久婷婷| 男男h啪啪无遮挡| 晚上一个人看的免费电影| 午夜福利在线在线| av女优亚洲男人天堂| 欧美日韩综合久久久久久| 欧美一区二区亚洲| 尾随美女入室| 成人特级av手机在线观看| 大又大粗又爽又黄少妇毛片口| 少妇人妻久久综合中文| 国产91av在线免费观看| 国产黄片视频在线免费观看| 亚洲人成网站在线播| 这个男人来自地球电影免费观看 | 亚洲精品乱码久久久v下载方式| 2022亚洲国产成人精品| 日韩免费高清中文字幕av| 国产黄片美女视频| 嫩草影院新地址| 18+在线观看网站| 久久久成人免费电影| 欧美丝袜亚洲另类| 亚洲av欧美aⅴ国产| av不卡在线播放| 人体艺术视频欧美日本| 91精品国产国语对白视频| 日韩国内少妇激情av| 一级毛片电影观看| 国产乱人偷精品视频| 最黄视频免费看| 晚上一个人看的免费电影| 国产 一区精品| 国产美女午夜福利| 久久久久久九九精品二区国产| 中文乱码字字幕精品一区二区三区| 老司机影院毛片| 欧美成人a在线观看| 久久国产精品大桥未久av | 黄色视频在线播放观看不卡| 欧美精品一区二区大全| 最近2019中文字幕mv第一页| 午夜福利网站1000一区二区三区| 99久久精品国产国产毛片| 精品国产一区二区三区久久久樱花 | 久久国产亚洲av麻豆专区| 成人免费观看视频高清| 亚洲aⅴ乱码一区二区在线播放| 久久久久国产网址| 一级毛片久久久久久久久女| 国产在线视频一区二区| 嫩草影院新地址| 免费看光身美女| 观看美女的网站| 18+在线观看网站| 亚洲精品国产色婷婷电影| 欧美人与善性xxx| 91在线精品国自产拍蜜月| 内射极品少妇av片p| 国产黄片视频在线免费观看| 少妇高潮的动态图| 日韩欧美精品免费久久| 丝瓜视频免费看黄片| 日韩人妻高清精品专区| 老司机影院成人| 欧美日韩一区二区视频在线观看视频在线| 夜夜骑夜夜射夜夜干| 精品国产露脸久久av麻豆| 91在线精品国自产拍蜜月| 在线精品无人区一区二区三 | 免费人成在线观看视频色| av一本久久久久| 久久av网站| 久久99热这里只频精品6学生| 欧美高清成人免费视频www| 午夜激情久久久久久久| 亚洲av成人精品一区久久| 啦啦啦在线观看免费高清www| 国产日韩欧美在线精品| 啦啦啦在线观看免费高清www| 欧美变态另类bdsm刘玥| 国产综合精华液| www.av在线官网国产| 王馨瑶露胸无遮挡在线观看| 成年美女黄网站色视频大全免费 | 啦啦啦在线观看免费高清www| 久热这里只有精品99| 精品一区在线观看国产| 在线精品无人区一区二区三 | 99精国产麻豆久久婷婷| 最近手机中文字幕大全| 午夜福利网站1000一区二区三区| 日韩 亚洲 欧美在线| 亚洲精品,欧美精品| 欧美日韩亚洲高清精品| 日韩制服骚丝袜av| 国产一区有黄有色的免费视频| 一级a做视频免费观看| a级一级毛片免费在线观看| 超碰av人人做人人爽久久| 亚洲国产日韩一区二区| 亚洲婷婷狠狠爱综合网| 精品熟女少妇av免费看| 乱系列少妇在线播放| 亚洲av欧美aⅴ国产| 亚洲精品日本国产第一区| 1000部很黄的大片| 黄色视频在线播放观看不卡| 街头女战士在线观看网站| 亚洲精品456在线播放app| 国产精品av视频在线免费观看| 成人美女网站在线观看视频| 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看 | 毛片一级片免费看久久久久| 国产免费福利视频在线观看| 最近2019中文字幕mv第一页| 只有这里有精品99| 亚洲精品456在线播放app| 最黄视频免费看| 免费黄频网站在线观看国产| 777米奇影视久久| 久久av网站| 国模一区二区三区四区视频| 国产成人精品久久久久久| 乱系列少妇在线播放| 日韩av在线免费看完整版不卡| 国产淫片久久久久久久久| 婷婷色麻豆天堂久久| 91精品国产九色| 亚洲欧洲国产日韩| tube8黄色片| 精品久久国产蜜桃| 内射极品少妇av片p| 能在线免费看毛片的网站| 男女免费视频国产| 三级国产精品片| 日本午夜av视频| 老师上课跳d突然被开到最大视频| 免费不卡的大黄色大毛片视频在线观看| 成人一区二区视频在线观看| 亚洲精品国产成人久久av| 岛国毛片在线播放| 亚洲综合精品二区| av网站免费在线观看视频| 韩国av在线不卡| 观看av在线不卡| 乱码一卡2卡4卡精品| 国产国拍精品亚洲av在线观看| 另类亚洲欧美激情| 久久99热这里只频精品6学生| 久久人妻熟女aⅴ| 亚洲精品国产av成人精品| 人人妻人人澡人人爽人人夜夜| 亚洲欧美日韩卡通动漫| 夜夜看夜夜爽夜夜摸| 日本欧美国产在线视频| 99视频精品全部免费 在线| 99久久精品国产国产毛片| 天天躁夜夜躁狠狠久久av| 久热久热在线精品观看| 少妇人妻 视频| 国产亚洲av片在线观看秒播厂| 亚洲内射少妇av| 国产在线视频一区二区| 一区二区三区乱码不卡18| 边亲边吃奶的免费视频|