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

    Charge Transfer Properties of Organic Semiconductor Molecules of Perylene Derivatives①

    2015-10-13 14:21:05TANYingXiongLIUJianBoLIQuanZHAOKeQing
    結(jié)構(gòu)化學(xué) 2015年3期

    TAN Ying-Xiong LIU Jian-Bo LI Quan ZHAO Ke-Qing

    ?

    Charge Transfer Properties of Organic Semiconductor Molecules of Perylene Derivatives①

    TAN Ying-Xiong LIU Jian-Bo LI Quan②ZHAO Ke-Qing

    (610066)

    The charge transfer rates of perylene and its four derivatives were studied at the level of B3LYP/6-31G** by density functional theory. The results showed that the perylene and its four derivatives belonged to the semiconductor molecules, which released energy when electron was injected. Therefore, they were suitable to be used as the electron injection material. The introduction of OH group can improve the electron transfer rate significantly. The formations of intramolecular hydrogen bonds were unfavorable to the hole transfer, but conducive to the electron transfer. The perylene derivatives, 2,5-3,4,5-(trifluorophenyl)ethynyl-8,11-3,4,5-trihydroxyphenyl ethynyl, designed in this article had the hole transfer rate of 1.57 cm2/V·s-1. Therefore, this kind of material will be potential hole transfer material with high transfer efficiency.

    perylene derivatives, density functional theory, hole transfer,electron transfer, organic semiconductors;

    1 INTRODUCTION

    The organic-conjugated materials have been widely used in the field of optoelectronics and microelectronics, such as organic light-emitting dio- des (OLEDs), organic field-effect transistors (OFETs) and organic photovoltaic devices (OPVs)[1-3]. The organic semiconductor materials can be divided into organic crystal, conjugated polymer material and liquid crystal semiconductor. The liquid crystal semiconductor has the characteri- stics such as high degree of ordering, molecular self-assembly and self-repair, and can be made by solution processing. Compared with the crystal material, it has the advantages of high charge transfer rate and low cost of device processing. The organic semiconductor materials can be designed and synthesized by molecules. The performance can be further enhanced. The organic compounds con- taining conjugated-electrons are typical semicon- ductor materials, which have received extensive attention. Polyacene and polythiophene are polyme- rized organic semiconductors. The triphenylene, truxene and perylene derivatives are the discoid organic semiconductor materials with bright pro- spect[4-6]. Perylene derivatives are a group of organic compounds with larger-electron structure. They have been proved to have excellent photoelectric properties recently. For these properties, they have found wide applications in the fields of charge transfer, organic light emitting and solar cells, and there are extensive studies in synthesis, characteri- zation and computational of perylene derivatives[7-9].

    Charge transfer is an important factor influencing the behavior of organic semiconductor materials, and charge carrier mobility is at the center of organic electronic devices[10-12]. According to the energy band theory, the movement of charge carrier, known as the charge transfer or charge migration, is an important factor for the characterization of semi- conductor materials. The holes, electron injections and transfer properties of four perylene derivatives were discussed through theoretical calculations in this article. The four perylene derivatives were perylene, 3,9-difluorine perylene, 3,9-dihydroxy perylene, 3,9-dihydroxy-4 and 10-difluorine per- ylene. The relationship between charge transfer properties and the structures of four perylene derivatives is studied, which contributes to the design and synthesis of perylene derivatives with high performance.

    2 THEORETICAL METHODS

    2. 1 Charge injection

    The process of electron injection is described with the following formula[13].

    The electron affinity of semiconductor molecule is expressed as, where(M-) and(M) denote the energy associated with negative ion and molecular geometry, respectively.represents that a molecule gains the energy released by an electron. The smaller the(the more negative), the easier the electron injection will be. This kind of material is suitable as the candidate material for the electron transfer layer.

    The process of hole injection is generally repre- sented by electron stripping.

    The performance of hole injection is generally described with the ionization energy

    where(M+) and(M) denote the energy asso- ciated with cation and molecular geometry, respect- tively. The ionization energyrepresents the energy required by a molecule to ionize one electron. The magnitude of ionization energy has a negative correlation with hole injection properties. The smaller the ionization energy is, the easier the ionization will be, and more suitable to be used as the material for the hole transfer layer.

    2. 2 Charge transfer

    The mobility of charge carrier can be obtained by the Einstein equation:

    The mobility of charge carrier is expressed as

    =e/(kB) (1)

    where e is the electron charge, kBis the Boltzmann constant,is the absolute temperature, andis the average diffusion coefficient of the charge starting from a molecule and towards all directions in three-dimensional space, as given by Formula (2):

    The diffusion coefficient is expressed as

    =(2)

    The small-molecule discotic liquid crystal material, as a kind of organic semiconductor ma- terial, has the properties of one-dimensional charge carrier migration. The average diffusion coefficient can be simplified as=2, whereis the disk spacing between adjacent discotic molecules, andis the constant of charge transfer rate between adjacent molecules. The carrier mobility is obtained by substituting into formula (1).

    According to the semi-classical model of Marcus charge transfer, the constant of charge transfer rate between adjacent molecules is expressed as the following formula[14, 15]:

    =(42h)2(4πkB)-0.5exp[-/(4kB)](4)

    whereis the charge transfer matrix element,is the absolute temperature,is the charge reor- ganization energy, h is the Planck constant and kBis the Boltzmann constant.

    The reorganization energy is calculated by the adiabatic potential energy surface. That is, the reorganization energy of hole transfer+and reorganization energy of electron transfer-are calculated as follows[10-12, 16]:

    +=(PE+/PE) –(PE+/PE+) +

    (PE/PE+) –(PE/PE) (5)

    -=(PE ̄/PE) –(PE ̄/PE ̄) +

    (PE/PE ̄) –(PE/PE) (6)

    where(PE+/PE) represents the total energy of PE+ions in PE configuration optimization.(PE+/PE+) represents the total energy of PE+ion configuration optimization. The denotations of other symbols are similar.

    The charge transfer matrix element characterizes the coupling strength of electron-electron interaction, and several methods have been proposed to evaluate the transfer integral within a molecular dimer. The simplest way is the frontier orbit energy level splitting method[10-12,16]. That is, the closed shell system was formed by adding one electron to the molecular/molecular cation system. The energy level splitting of HOMO and HOMO-1 at the transition state was calculated. One half is the hole transfer matrix element t+. The closed shell system is formed by reducing one electron from the molecular/molecular anion system. The energy level splitting of LUMO and LUMO+1 at the transition state is calculated. One half is the negative charge transfer matrix element t-.

    3 RESULTS AND DISCUSSION

    3. 1 Molecular structure and charge injection

    The geometrical structure optimization and energy calculation for the five molecules and molecular ions of perylene derivatives as shown in Fig. 1 were performed by the B3LYP/6-31G** the- oretical level in the Gaussian 09[17]. The results are shown in Table 1. D molecule was the perylene molecule not substituted; e molecule was 3,9- difluorine perylene; f molecule was 3,9-dihydroxy perylene; and g1 and g2 molecules (g2 contains F···H–O intramolecular hydrogen bond) are the derivatives, in which the Nos. 3 and 9 positions are substituted by F, and Nos. 4 and 10 positions by hydroxyl group. The binding energy of hydrogen bond was calculated as 19.69 kJ/mol. The bond lengths of O–H in the g1 and g2 molecules and their vibration frequencies are shown in Table 2.

    Table 1. Energy Gap Eg of Perylene Derivatives Molecules. Ionization Energy E(IP), Electron Affinity E(EA) and Recombination Energy (kJ/mol)

    TH, F-TH and OH-TH are triphenylene, fluorine-substituted triphenylene, and hydroxyl-substituted triphenylene respectively, as shown in literature [16].

    Table 2. Hydroxyl Bond Length R (O-H) (nm) in the g1 and g2 Molecules and Stretching Vibrational Frequency V (O-H) (cm-1)

    d:R1=R2=H;e:R1=H,R2=F;f:R1=H,R2=OH

    g1:R1=F,R2=OH; g2:R1=F,R2=OH

    Fig. 1 . Molecular structures of perylene derivatives

    As seen from the data in Table 1, the energy gap of five perylene derivatives molecules is about 285 kJ/mol (2.95 eV). The derivative molecules be- longed to the organic semiconductor molecules. A hole injection for the molecules will absorb higher energy (600 kJ/mol). Therefore, the molecules of derivatives were not suitable to be used as the hole injection materials. The electron affinities were ne- gative, which indicated that the energy release was required for electron injection. The amount of released energy was related to the types of sub- stituent. The substitution of e molecules of perylene derivatives by the electron-withdrawing fluorine atoms released the most energies, while the sub- stitution of f molecules by electron-donating hydroxyl group released the least energies. It is easier to achieve the electron injection of g2 molecule than that of g1 molecule, due to the formation of intramolecular hydrogen bonds. Therefore, the perylene derivatives were suitable to be used as electron injection materials. In terms of recombination energy, the introductions of elec- tron-withdrawing F or electron-donating OH were not conducive to electronic transfer or hole transfer, which was also the case for triphenylene derivatives. For either the hole transfer or electron transfer, the recombination energies of perylene and its deriva- tive molecules were smaller than those of tripheny- lene and its derivatives. It was thus indicated that perylene and its derivatives were more suitable to be used as charge transfer material than tripheny- lene and its derivatives.

    As shown from Table 2, the formation of F···H-O hydrogen bond in the molecules leng- thened the O–H bond, and the stretching vibrational frequency decreases, which presented the typical characteristics of red-shifted hydrogen bonds.

    3. 2 Charge transfer

    When the molecules and ions with optimized structures overlapped, and the disk spacing was 0.45 nm, the molecules and ions were rotated by 0~180o, considering the symmetric structure of each molecule. The difference was 20o. The total energy E of molecular ion dimer and energy level splitting value were calculated in each rotation angle (ΔHOMO and ΔLUMO). Fig. 2 shows the relationship between the energy level splitting in charge transfer and rotation angle for each molecule.

    Fig. 2 . Relationship between the energy level splitting value and rotation angle of perylene derivatives

    The energy level splitting values (ΔHOMO or ΔLUMO) are determined by the couplings of molecules and ions (HOMO or LUMO). For ΔLUMO, when the molecules and anions were overlapped at certain angle and distance to form the transition state molecules, ΔLUMO was the dif- ference between LUMO+1 and LUMO. If LUMO of molecules and ions had more overlapping in one-dimensional space, the molecular ion dimer formed had a lower LUMO. As a result, the accumulation of molecules and ions showed more prominent bonding characteristics. The ΔLUMO was also larger. Correspondingly, the charge transfer matrix element and the charge transfer rate increased. As shown in Fig. 2, the energy level splitting values of the hole or electron transfer presented regular changes with the rotation angle. The rotation angle was in the range of 0~180o. The images of energy level splitting values of charge transfer were symmetrical. When the disk spacing of molecules was fixed, the maximum values of energy level splitting of charge transfer for five molecules emerged within the range of rotation angle of 0o or 180o. The hole and electron transfer matrix elements were the maximum, given the maximum overlapped area of molecules and ions. When the rotation angle was in the range of 20~60o, the hole transfer matrix element was greater than that of the electron transfer. Within the range of 80~100o, the electron transfer matrix element was greater than that of the hole transfer. When the rotation angle was 0o, the hole transfer matrix element decreased successively from d, f, g2, e to the g1 molecules (the energy level splitting values are 0.147, 0.133, 0.111, 0.105 and 0.101 eV, respectively). It was indicated that the hole transfer matrix element for the perylene substituted by electron-donating group was greater than that of perylene substituted by electron-withdrawing groups. By comparison between g2 and g1, it was found that the hole transfer matrix element of g2 molecule was greater than that of g1 molecule, due to the better electronic delocalization of g2 mole- cule. There was a great difference with the previous research result of triphenylene. When the rotation angle was 0o, the electron transfer matrix element reduced successively from d, e, g2, f to the g1 mo- lecule (the energy level splitting values are 0.199, 0.188, 0.171, 0.132 and 0.114 eV, respectively). It was indicated that the introduction of electron- withdrawing group was more conducive to electron transfer than that of the electron-donating group. The formation of intramolecular hydrogen bonds was conducive to the electronic delocalization of rigid nucleus, which would increase the electron transfer matrix element.

    According to the lowest energy principle, the overlapped disk spacing and rotation angle between molecules and molecular ions were selected. Then, the charge transfer matrix element and charge transfer rate were calculated. The calculation results are shown in Table 3.

    Table 3. Relationship between Disk Spacing dis, the Lowest Energy of Charge Transfer (nm) and Rotation Angle θ (o). Relationship between the Charge Transfer Matrix Element T (kJ/mol) and Charge Transfer Rate (cm2/V·S-1)

    For perylene molecule d, when the distance was 0.55 nm and the rotation angle was 100o, the molecule and molecular ions overlapped with the lowest energy. At this moment, the accumulation was in the optimal state. The electron transfer rate was 2.76 times greater than the hole transfer rate. Compared with the hole transfer rates of tripheny- lene and its derivatives[12], those of perylene and its hydroxyl or fluorine-substituted derivatives decrea- sed obviously. But the electron transfer rates were equivalent. The introduction of electron-withdraw- ing fluorine atoms or electron-donating hydroxyl group would result in the increase of hole transfer matrix element of perylene, thereby leading to the increase of hole transfer rate of perylene. The introduction of electron-donating OH increased the electron transfer matrix element significantly, which improved the electron transfer rate. Comparing g1 and g2 molecules, it was found that the formation of intramolecular hydrogen bonds was not conducive to the hole transfer. For the g2 molecule, the formation of intramolecular hydrogen bond was apparently conducive to the electron transfer.

    3. 3 Design of a new type of perylene derivatives

    According to the above research result, the h molecule of perylene derivative is shown in Fig. 3.

    Fig. 3 . Structure of the h molecule

    In h molecule, at Nos. 2 and 5 positions, the 3,4,5-(trifluorophenyl) ethynyl was introduced, and at Nos. 8 and 11 positions, the 3,4,5-trihydroxy- phenyl ethynyl was introduced. When the molecules and ions rotated by 180o, the large-interaction between the molecules and ions emerged due to the overlap of electron-deficient trifluorophenyl and electron-rich trihydroxyphenyl. Therefore, the elec- tron coupling increased.

    The energy gap of h molecule was calculated on the basis of the optimized structures of molecules and molecular ions, which was 270.24 kJ/mol. The ionization energy and electron affinity were 590.64 and –151.39 kJ/mol, respectively. The recombina- tion energies of hole and electron transfer were 12.64 and 14.21 kJ/mol, respectively. It was indica- ted that the h molecule had semiconductor pro- perties. The performance of hole injection was poor, while the electron injection was accompanied by energy release. The recombination energies of hole and electron transfer were lower than those of d (perylene) molecule, which was beneficial to the charge transfer.

    When the h molecules and molecular ions overlapped, and the disk spacing was set at 0.45 nm, the molecular ion dimer was rotated at a difference of 20o within the range of 0~180o. Then the energy and energy level splitting value at the transition state were calculated. The results are shown in Figs. 4 and 5, respectively.

    As shown in Figs. 4 and 5, when the h molecule transferred the holes and electrons, the molecules and cations as well as the molecules and anions accumulated. The lowest energy value emerged at 40° and 140°. At this moment, the matrix element was the largest (energy level splitting value was also the largest), which was conducive to both the hole transfer and electron transfer. When the rotation angle of molecular ion dimer was set at 140°, and the distance between molecules and ions was changed to 0.34~0.45 nm, the energy of molecular ion system was calculated. The distance corresponding to the lowest energy was obtained. Moreover, the hole transfer and electron transfer matrix elements at this state were calculated, with the values being 5.78 and 3.81 kJ/mol, respectively. The hole transfer and electron transfer rates were 1.57 and 0.506 cm2/V·s-1, respectively. For either hole transfer or electron transfer, the transfer rate of h molecule was much greater than that of d, e, f and g molecules. Especially, the hole transfer properties were greatly improved. The reason was that the coupling degree of the electron-electron interaction was enhanced, while the recombination energy of transfer decreased, leading to the obvious increase of transfer matrix element. The hole transfer rate of h molecule was 3.1 times larger than the electron transfer rate. Thus, the h molecule can be used as the hole transfer material with high efficiency.

    Fig. 4 . Relationship between energy and rotation angle of h molecules and molecular ion dimers

    Fig. 5 . Relationship between energy level splitting value in charge transfer and rotation angle for the h molecule

    4 CONCLUSION

    The charge transfer rates of perylene and hy- droxyl and/or fluorine-substituted perylene deriva- tives were calculated theoretically. The following conclusions were derived:

    (1) The molecules of perylene derivatives were organic semiconductor molecules and therefore unsuitable to be used as the hole injection material but suitable as the electron injection material. The substitution properties and the formations of intramolecular hydrogen bonds would all influence the electron injection ability.

    (2) Compared with the hole transfer rates of triphenylene and its derivatives, those of perylene and its hydroxyl or fluorine-substituted derivatives were much lower. But the electron transfer rates were equivalent. The introductions of electron- withdrawing fluorine atoms or electron-donating hydroxyl groups could increase the hole transfer rate. The introduction of electron-donating OH would significantly improve the electron transfer rate. The formation of intramolecular hydrogen bonds was unfavorable to the hole transfer, but conducive to the electron transfer.

    (3) The perylene derivative h molecules of 2,5-3,4,5-(trifluorophenyl)ethynyl-8,11-3,4,5- trihy- droxyphenyl ethynyl designed in this article had the hole transfer rate of 1.57 cm2/V·s-1, with the electron transfer rate to be 0.506 cm2/V·s-1.

    (1) Schmidt-Mende, L.; Fechtenk?tter, A.; Müllen, K.; Moons, E.; Friend, R. H.; MacKenzie, J. D. Self-organized discotic liquid crystals for high-efficiency organic photovoltaics.2001, 293, 1119-1122.

    (2) Sergeyev, S.; Pisula, W.; Geerts, Y. H. Discotic liquid crystals: a new generation of organic semiconductors.2007,36, 1902-1929.

    (3) Laschat, S.; Baro, A.; Steinke, N.; Giesselmann, F.; H?gele, C.; Scalia, G.; Judele, R.; Kapatsina, E.; Sauer, S.; Schreivogel A.; Tosoni, M. Discotic liquid crystals: from tailor-made synthesis to plastic electronics.. 2007, 46, 4832-4887.

    (4) Feng, X.; Marcon, V.; Pisula, W.; Hansen, M. R.; Kirkpatrick, J.; Grozema, F.; Andrienko, D.; Kremer, K.; Müllen, K. Towards high charge-carrier mobilities by rational design of the shape and periphery of discotics.2009, 8, 421-422.

    (5) Bai, Y. F.; Zhao, K. Q.; Hu, P.; Wang, B. Q.; Shimizu, Y. Synthesis of amide group containing triphenylene derivatives as discotic liquid crystals and organic gelator.2009, 509, 60-76.

    (6) Ji, H.; Zhao, K. Q.; Yu, W. H.; Wang, B. Q.; Hu, P. Synthesis and mesomorphism of diacetylene-bridged triphenylene discotic liquid crystal dimers.2009, 52, 975-985.

    (7) Peter, O.; Frank, W. Effects of bay substituents on the racemization barriers of perylene bisimides: resolution of atropo-enantiomers.2007, 129, 14319-14326.

    (8) Wei, W.; Andrew, D. S.; Alexander, D. Q. L. Twisted perylene stereodimers reveal chiral molecular assembly codes.2008, 130, 8271-8279.

    (9) Tanya, N.; Singh, R.; Felix, N. C. Triplet sensitized red-to-blue photon upconversion.2010, 1, 195-200.

    (10) Chen, X. K.; Guo, J. F.; Zou, L. Y.; Ren, A. M.; Fan, J. X. A promising approach to obtain excellent n-type organic field-effect transistors: introducing pyrazine ring.2011, 115, 21416-21428.

    (11) Chem, X. K.; Zou, L. Y.; Ren, A. M.; Fan, J. X.How dual bridging atoms tune structural and optoelectronic properties of ladder-type heterotetracenes-a theoretical study.2011, 13, 19490-19498.

    (12) Chen, X. K.; Zou, L. Y.; Huang, S.; Min, C. G.; Ren, A. M.; Fen, J. K.; Sun, C. C.Theoretical investigation of charge injection and transport properties of novel organic semiconductor materials-cyclic oligothiophenes.2011, 12, 1198-1210.

    (13) Feng, J. K.2010, 26, 224. Feng, J. K. The molecular design investigation on OLED materials.. 2010, 26, 224-245.

    (14) Lemaur, V.; Filho, D. A.; Coropceanu, V.; Lehmann, M.; Geerts, Y.; Piris, J.; Debije, M. G.; Graats, A. M.; Senthilkumar, K.; Siebbeles, L. D. A.; Warman, J. M.; Bredas, J. L.; Cornil, J. Charge transport properties in discotic liquid crystals: a quantum-chemical insight into structure-property relations.2004, 126, 3271-3279.

    (15) Cornil, J.; Lemaur, V.; Calbert, J. P.; Bredas, J. L. Charge transport in discotic liquid crystals: a molecular scale description.2002, 14, 726-729.

    (16) Hu, J. D.; Li, Q.; Zhao, K. Q. Study on structure and charge transport properties of the hexahydroxy-triphorylene.2007, 65, 1784-1788.

    (17) Frisch, M. J.; Trucks G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Zakrzewski, V. G.; Montgomery, J. A.; Stratmann, J. R. E.; Burant, J. C.; Dapprich, S.; Millam, J. M.; Daniels, A. D.; Kudin, K. N.; Strain, M. C.; Farkas, O.; Tomasi, J.; Barone, V.; Cossi, M.; Cammi, R.; Mennucci, B.; Pomelli, C.; Adamo, C.; Clifford, S.; Ochterski, J.; Petersson, G. A.; Ayala, P. Y.; Cui, Q.; Morokuma, K.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Cioslowski, J.; Ortiz, J. V.; Baboul, A. G.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Andres, J. L.; Gonzalez, C.; Head-Gordon, M.; Replogle, E. S.; Pople, G. A., Revision B. 05. Pittsburgh, PA: Gaussian, Inc. 2009.

    31 March 2014; accepted 5 January 2015

    ①Supported by the National Natural Science Foundation of China (No. 51273133), the Opening Project of National Key Laboratory of Theoretical Chemical Computation (No. K1202), and the Department of Education in Sichuan Province (No. 11ZB086)

    . Li Quan, born in 1966, professor, Tel: 13808003581, E-mail: liquan6688@163.com

    10.14102/j.cnki.0254-5861.2011-0312

    av天堂中文字幕网| 在现免费观看毛片| 国产精品福利在线免费观看| 亚洲经典国产精华液单| 直男gayav资源| 男插女下体视频免费在线播放| 欧美一区二区国产精品久久精品| 亚洲国产欧美人成| 美女免费视频网站| 精品人妻偷拍中文字幕| 此物有八面人人有两片| 国产高清不卡午夜福利| 欧美激情国产日韩精品一区| 人妻制服诱惑在线中文字幕| .国产精品久久| 12—13女人毛片做爰片一| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区免费观看 | 夜夜夜夜夜久久久久| videossex国产| 岛国在线免费视频观看| 校园人妻丝袜中文字幕| 国产伦精品一区二区三区视频9| 婷婷精品国产亚洲av在线| 日本黄色视频三级网站网址| 国产精品野战在线观看| 麻豆精品久久久久久蜜桃| 久久久久国内视频| 99九九线精品视频在线观看视频| 日本精品一区二区三区蜜桃| 国产成人91sexporn| 91av网一区二区| 中文亚洲av片在线观看爽| 91狼人影院| 国产在视频线在精品| 熟女人妻精品中文字幕| 成人午夜高清在线视频| 国国产精品蜜臀av免费| 亚洲成人精品中文字幕电影| 深夜a级毛片| 久久亚洲精品不卡| www.色视频.com| 黑人高潮一二区| 亚州av有码| 国产亚洲91精品色在线| 亚洲精品在线观看二区| 免费看光身美女| 国内精品一区二区在线观看| 午夜精品国产一区二区电影 | 成熟少妇高潮喷水视频| 男人和女人高潮做爰伦理| 日韩av在线大香蕉| 嫩草影院精品99| 婷婷色综合大香蕉| www日本黄色视频网| 一进一出抽搐gif免费好疼| 搞女人的毛片| 波多野结衣巨乳人妻| 99精品在免费线老司机午夜| 中文字幕久久专区| 国产探花在线观看一区二区| 内射极品少妇av片p| 内射极品少妇av片p| 国产精品国产高清国产av| 亚洲色图av天堂| 午夜精品一区二区三区免费看| 国产单亲对白刺激| 精品一区二区三区视频在线| av视频在线观看入口| 欧美性猛交╳xxx乱大交人| 国产亚洲91精品色在线| 免费看a级黄色片| 国产毛片a区久久久久| 九九爱精品视频在线观看| 可以在线观看的亚洲视频| 51国产日韩欧美| 俺也久久电影网| 一级毛片久久久久久久久女| 国产av不卡久久| 色综合亚洲欧美另类图片| 日韩精品青青久久久久久| 国产精品野战在线观看| 成人亚洲精品av一区二区| 免费在线观看成人毛片| 精品无人区乱码1区二区| 亚洲欧美清纯卡通| 97超级碰碰碰精品色视频在线观看| 女的被弄到高潮叫床怎么办| 老司机午夜福利在线观看视频| 成人二区视频| 久久精品国产亚洲av天美| 国产 一区 欧美 日韩| 欧美另类亚洲清纯唯美| 日本撒尿小便嘘嘘汇集6| 久久99热6这里只有精品| 国产淫片久久久久久久久| 最后的刺客免费高清国语| 免费搜索国产男女视频| 男人和女人高潮做爰伦理| videossex国产| 一边摸一边抽搐一进一小说| 女人十人毛片免费观看3o分钟| 成人漫画全彩无遮挡| 国产精品日韩av在线免费观看| 欧美性猛交╳xxx乱大交人| 国产精品日韩av在线免费观看| 最近中文字幕高清免费大全6| 69av精品久久久久久| 亚洲国产精品成人综合色| 在线观看美女被高潮喷水网站| 听说在线观看完整版免费高清| 69人妻影院| 久久午夜福利片| 亚洲欧美日韩高清专用| 国产中年淑女户外野战色| 久久久国产成人精品二区| 亚洲久久久久久中文字幕| 国产在线精品亚洲第一网站| 国产精品电影一区二区三区| 男女那种视频在线观看| 国产色爽女视频免费观看| 亚洲国产精品久久男人天堂| 男女之事视频高清在线观看| 欧美性猛交╳xxx乱大交人| 国产国拍精品亚洲av在线观看| 国产一区二区三区在线臀色熟女| 久久久久久国产a免费观看| 少妇猛男粗大的猛烈进出视频 | 夜夜夜夜夜久久久久| 久久99热6这里只有精品| 中国美女看黄片| 校园人妻丝袜中文字幕| 丰满的人妻完整版| 亚洲av中文av极速乱| 2021天堂中文幕一二区在线观| 看非洲黑人一级黄片| 悠悠久久av| 亚洲中文字幕一区二区三区有码在线看| 久久人人精品亚洲av| 精品久久久久久久久av| a级毛片免费高清观看在线播放| 久久亚洲精品不卡| 一a级毛片在线观看| 国产精品免费一区二区三区在线| 久久人人爽人人爽人人片va| 性插视频无遮挡在线免费观看| 99riav亚洲国产免费| 老司机影院成人| 99热全是精品| 97超级碰碰碰精品色视频在线观看| 午夜久久久久精精品| 老师上课跳d突然被开到最大视频| 午夜精品国产一区二区电影 | 91狼人影院| 搞女人的毛片| 国产精品国产高清国产av| 国产精品人妻久久久久久| 亚洲精品在线观看二区| 国产精品亚洲一级av第二区| av.在线天堂| 又黄又爽又免费观看的视频| 大香蕉久久网| 九九在线视频观看精品| 日本成人三级电影网站| 99久久精品国产国产毛片| 能在线免费观看的黄片| 国产人妻一区二区三区在| 麻豆成人午夜福利视频| 国产精品嫩草影院av在线观看| 国产精品嫩草影院av在线观看| 成人av在线播放网站| 少妇人妻精品综合一区二区 | 国产一区二区在线观看日韩| 九九久久精品国产亚洲av麻豆| 久久精品国产清高在天天线| 精品午夜福利在线看| 色哟哟·www| 深爱激情五月婷婷| 春色校园在线视频观看| 国产精品一区www在线观看| 亚洲精品456在线播放app| 亚洲不卡免费看| 久久久a久久爽久久v久久| 国产av麻豆久久久久久久| 亚洲中文日韩欧美视频| 搡老熟女国产l中国老女人| 亚洲av五月六月丁香网| 最好的美女福利视频网| 国产真实伦视频高清在线观看| av免费在线看不卡| 欧美性猛交黑人性爽| 99国产精品一区二区蜜桃av| 亚洲五月天丁香| 能在线免费观看的黄片| 国产综合懂色| 免费av观看视频| 精品欧美国产一区二区三| 午夜福利视频1000在线观看| 日韩欧美精品免费久久| 色综合亚洲欧美另类图片| 嫩草影视91久久| 大型黄色视频在线免费观看| eeuss影院久久| 无遮挡黄片免费观看| 国产精品福利在线免费观看| 在线观看一区二区三区| 国内揄拍国产精品人妻在线| 香蕉av资源在线| 高清日韩中文字幕在线| 亚洲国产欧洲综合997久久,| 黑人高潮一二区| 长腿黑丝高跟| 亚洲性夜色夜夜综合| 男女视频在线观看网站免费| 欧美成人一区二区免费高清观看| 欧美成人a在线观看| 亚洲va在线va天堂va国产| 国产亚洲av嫩草精品影院| 人妻少妇偷人精品九色| 秋霞在线观看毛片| 成人二区视频| 内地一区二区视频在线| 国产精品一区二区免费欧美| 欧美+日韩+精品| .国产精品久久| 99久久久亚洲精品蜜臀av| 看非洲黑人一级黄片| 久久久久久伊人网av| 香蕉av资源在线| 九九热线精品视视频播放| 欧美激情在线99| 久久久色成人| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 久久精品91蜜桃| 色在线成人网| 我的女老师完整版在线观看| 中国国产av一级| 深爱激情五月婷婷| 男插女下体视频免费在线播放| 最新在线观看一区二区三区| 久久精品国产鲁丝片午夜精品| 97超碰精品成人国产| 久久精品久久久久久噜噜老黄 | 久久99热6这里只有精品| 一卡2卡三卡四卡精品乱码亚洲| 免费观看精品视频网站| 日韩,欧美,国产一区二区三区 | 别揉我奶头~嗯~啊~动态视频| 校园春色视频在线观看| 亚洲人成网站在线观看播放| 国产黄片美女视频| 国产亚洲精品久久久久久毛片| 日韩人妻高清精品专区| 国产精品,欧美在线| 国产精品精品国产色婷婷| videossex国产| 久久精品国产亚洲网站| 久久久久精品国产欧美久久久| 成人美女网站在线观看视频| 人妻制服诱惑在线中文字幕| ponron亚洲| 亚洲av中文av极速乱| 久久精品国产清高在天天线| АⅤ资源中文在线天堂| 男女之事视频高清在线观看| 性色avwww在线观看| 国产精品一区www在线观看| 精品久久久久久成人av| 真实男女啪啪啪动态图| 婷婷亚洲欧美| 日韩av不卡免费在线播放| 午夜a级毛片| АⅤ资源中文在线天堂| 99久国产av精品国产电影| 中文在线观看免费www的网站| 人妻制服诱惑在线中文字幕| 村上凉子中文字幕在线| 俄罗斯特黄特色一大片| a级毛色黄片| 国产在线男女| 91午夜精品亚洲一区二区三区| 综合色av麻豆| 99热全是精品| 少妇的逼水好多| 男人舔奶头视频| 国国产精品蜜臀av免费| 亚洲国产色片| 午夜久久久久精精品| 在线a可以看的网站| 国内揄拍国产精品人妻在线| 成人漫画全彩无遮挡| 韩国av在线不卡| 国产伦一二天堂av在线观看| 国产av在哪里看| 亚洲av美国av| 午夜精品在线福利| 中文字幕精品亚洲无线码一区| 亚洲图色成人| 久久精品国产亚洲av涩爱 | 色尼玛亚洲综合影院| 亚洲av五月六月丁香网| 国产真实伦视频高清在线观看| 乱人视频在线观看| 亚洲最大成人手机在线| 老师上课跳d突然被开到最大视频| 白带黄色成豆腐渣| 国产蜜桃级精品一区二区三区| 美女内射精品一级片tv| 高清毛片免费看| 国产白丝娇喘喷水9色精品| 五月伊人婷婷丁香| 精品午夜福利在线看| 久久天躁狠狠躁夜夜2o2o| 日本三级黄在线观看| 亚洲av五月六月丁香网| 久久99热这里只有精品18| 亚洲欧美清纯卡通| 亚洲av成人精品一区久久| 成人特级黄色片久久久久久久| 免费电影在线观看免费观看| 亚洲精品国产av成人精品 | 亚洲精品一区av在线观看| 我要搜黄色片| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 18禁黄网站禁片免费观看直播| 嫩草影院新地址| 内地一区二区视频在线| 成年版毛片免费区| 亚洲美女视频黄频| 18禁在线无遮挡免费观看视频 | 精品人妻视频免费看| а√天堂www在线а√下载| 精品久久久久久成人av| 欧美激情国产日韩精品一区| 久久人人精品亚洲av| 男女视频在线观看网站免费| 中文字幕免费在线视频6| 99久久精品一区二区三区| 男人舔奶头视频| 一个人免费在线观看电影| 日本在线视频免费播放| 网址你懂的国产日韩在线| 深夜a级毛片| 亚洲国产精品成人综合色| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av| 香蕉av资源在线| 俺也久久电影网| 啦啦啦观看免费观看视频高清| 12—13女人毛片做爰片一| АⅤ资源中文在线天堂| 成人二区视频| 尾随美女入室| 男人舔奶头视频| 淫秽高清视频在线观看| 色av中文字幕| 精品福利观看| 亚洲精品色激情综合| 毛片女人毛片| 成人美女网站在线观看视频| 精品熟女少妇av免费看| 热99在线观看视频| 18禁在线播放成人免费| 久久人人爽人人片av| 麻豆一二三区av精品| 亚洲精品亚洲一区二区| 久久久精品欧美日韩精品| a级毛色黄片| 91狼人影院| 一级a爱片免费观看的视频| 日本爱情动作片www.在线观看 | 国产白丝娇喘喷水9色精品| 黄片wwwwww| 嫩草影院新地址| 精品久久久久久成人av| 三级经典国产精品| 99久久精品热视频| 国产真实乱freesex| 免费看光身美女| 久久亚洲精品不卡| 国内精品宾馆在线| av专区在线播放| 老熟妇仑乱视频hdxx| 免费看美女性在线毛片视频| 精品久久久久久成人av| 久久久久国产精品人妻aⅴ院| 深夜a级毛片| 久久精品国产自在天天线| 五月玫瑰六月丁香| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 国产白丝娇喘喷水9色精品| 国产黄a三级三级三级人| 麻豆av噜噜一区二区三区| 国产午夜精品论理片| 国产 一区精品| 欧美激情在线99| 日韩欧美国产在线观看| 亚洲精品日韩在线中文字幕 | 欧美一区二区亚洲| 欧美另类亚洲清纯唯美| 黄色一级大片看看| 国产精品人妻久久久影院| 国产精品一区www在线观看| 午夜免费男女啪啪视频观看 | 亚洲一级一片aⅴ在线观看| 特级一级黄色大片| 草草在线视频免费看| 1024手机看黄色片| 日日摸夜夜添夜夜添av毛片| av国产免费在线观看| 99riav亚洲国产免费| 成人一区二区视频在线观看| 国产 一区精品| 亚洲四区av| 久久久久国产网址| 最新中文字幕久久久久| 日韩欧美精品免费久久| 亚洲三级黄色毛片| 两个人的视频大全免费| 成年女人永久免费观看视频| 欧美丝袜亚洲另类| 免费看光身美女| 菩萨蛮人人尽说江南好唐韦庄 | 蜜桃久久精品国产亚洲av| 一个人看的www免费观看视频| 一进一出好大好爽视频| 在现免费观看毛片| 国产成人一区二区在线| 亚州av有码| 晚上一个人看的免费电影| 日韩精品中文字幕看吧| 国产亚洲91精品色在线| 色av中文字幕| 日本黄大片高清| 亚洲国产精品成人综合色| 亚洲av成人精品一区久久| 97人妻精品一区二区三区麻豆| 97热精品久久久久久| 中文在线观看免费www的网站| 在线天堂最新版资源| 综合色丁香网| 日韩欧美国产在线观看| 人人妻人人澡欧美一区二区| 最后的刺客免费高清国语| 亚洲人成网站在线观看播放| 搡老岳熟女国产| 国产精华一区二区三区| 中文资源天堂在线| 我的女老师完整版在线观看| 精品无人区乱码1区二区| .国产精品久久| 综合色av麻豆| 亚洲美女搞黄在线观看 | 99久久无色码亚洲精品果冻| 国产色婷婷99| 18禁黄网站禁片免费观看直播| 久久九九热精品免费| 色5月婷婷丁香| av中文乱码字幕在线| а√天堂www在线а√下载| 欧美潮喷喷水| 日本色播在线视频| 国产视频内射| 欧美+亚洲+日韩+国产| 一本一本综合久久| 日韩精品中文字幕看吧| 天堂√8在线中文| 中文字幕免费在线视频6| 久久久久久九九精品二区国产| 亚洲四区av| 老熟妇乱子伦视频在线观看| 久久久久久久午夜电影| 亚洲国产精品国产精品| 国产精品福利在线免费观看| 亚洲欧美日韩高清专用| 黄色配什么色好看| 亚洲精品国产成人久久av| 国产在视频线在精品| 少妇人妻一区二区三区视频| 亚洲欧美日韩高清专用| 白带黄色成豆腐渣| 中国国产av一级| 天堂网av新在线| 晚上一个人看的免费电影| 日本黄大片高清| 色播亚洲综合网| 少妇的逼水好多| 床上黄色一级片| 精品久久久噜噜| 免费av毛片视频| 一级av片app| 神马国产精品三级电影在线观看| 俺也久久电影网| 网址你懂的国产日韩在线| 国产黄色小视频在线观看| 午夜免费男女啪啪视频观看 | 亚洲第一电影网av| 少妇被粗大猛烈的视频| 狂野欧美白嫩少妇大欣赏| 欧美三级亚洲精品| 欧美在线一区亚洲| 久久亚洲国产成人精品v| 日本精品一区二区三区蜜桃| 国产乱人偷精品视频| av国产免费在线观看| 亚洲成a人片在线一区二区| 22中文网久久字幕| 在现免费观看毛片| 女同久久另类99精品国产91| а√天堂www在线а√下载| 亚洲精华国产精华液的使用体验 | 久久精品国产亚洲av香蕉五月| 啦啦啦啦在线视频资源| 久久精品人妻少妇| 老司机影院成人| 中文亚洲av片在线观看爽| 免费人成视频x8x8入口观看| 搡老熟女国产l中国老女人| 国产一区二区在线av高清观看| 十八禁国产超污无遮挡网站| 我的女老师完整版在线观看| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 久久久久九九精品影院| 免费无遮挡裸体视频| 亚洲欧美日韩卡通动漫| 国产精品99久久久久久久久| 久久久精品94久久精品| 亚洲欧美日韩高清在线视频| www.色视频.com| 91精品国产九色| 亚洲真实伦在线观看| 国产片特级美女逼逼视频| 亚洲成人av在线免费| 狠狠狠狠99中文字幕| 人人妻人人看人人澡| 久久人人爽人人爽人人片va| 老师上课跳d突然被开到最大视频| av国产免费在线观看| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 一本久久中文字幕| 亚洲av二区三区四区| 一夜夜www| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 成人美女网站在线观看视频| 免费av毛片视频| 内地一区二区视频在线| 少妇人妻一区二区三区视频| 一级黄片播放器| 22中文网久久字幕| 最近在线观看免费完整版| 高清午夜精品一区二区三区 | 亚洲最大成人av| 亚洲丝袜综合中文字幕| 亚洲18禁久久av| 国产综合懂色| 午夜a级毛片| 精品日产1卡2卡| 91av网一区二区| 免费看av在线观看网站| 国产av在哪里看| 精品无人区乱码1区二区| 国产av在哪里看| 一级av片app| 久久99热这里只有精品18| 淫妇啪啪啪对白视频| 丰满乱子伦码专区| 成人一区二区视频在线观看| 国产欧美日韩精品亚洲av| 国产色婷婷99| 久久精品综合一区二区三区| 内射极品少妇av片p| or卡值多少钱| 成人永久免费在线观看视频| 成人漫画全彩无遮挡| 国产探花极品一区二区| 欧美性感艳星| 亚洲性久久影院| 综合色丁香网| 不卡视频在线观看欧美| 亚洲av一区综合| 搞女人的毛片| 99国产精品一区二区蜜桃av| 网址你懂的国产日韩在线| 18禁在线无遮挡免费观看视频 | 亚洲七黄色美女视频| 国产精品乱码一区二三区的特点| 亚洲国产精品成人综合色| 在线观看免费视频日本深夜| 麻豆成人午夜福利视频| 此物有八面人人有两片| 老司机午夜福利在线观看视频| 美女大奶头视频| 免费黄网站久久成人精品| 日韩欧美一区二区三区在线观看| 尤物成人国产欧美一区二区三区| 91麻豆精品激情在线观看国产| 久久久a久久爽久久v久久| 国产高清视频在线观看网站| 精品国内亚洲2022精品成人| 有码 亚洲区| 国产亚洲精品av在线| 国产黄色小视频在线观看| 亚洲国产精品合色在线| 夜夜爽天天搞| 国产精品久久久久久久电影| 久久午夜亚洲精品久久| 22中文网久久字幕| 韩国av在线不卡| 日本色播在线视频|