• <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

    久久精品久久久久久久性| a级片在线免费高清观看视频| 我的女老师完整版在线观看| 中文字幕亚洲精品专区| 欧美xxⅹ黑人| 亚洲精品视频女| 日本wwww免费看| 一本色道久久久久久精品综合| 欧美精品国产亚洲| 色视频在线一区二区三区| 免费日韩欧美在线观看| 麻豆乱淫一区二区| 日本欧美视频一区| 久久精品国产自在天天线| 内地一区二区视频在线| 一区二区日韩欧美中文字幕 | 日韩一本色道免费dvd| 中文乱码字字幕精品一区二区三区| 亚洲欧美精品自产自拍| a级毛片黄视频| 伊人亚洲综合成人网| 欧美xxxx性猛交bbbb| 成人毛片60女人毛片免费| 高清欧美精品videossex| 久久99一区二区三区| a 毛片基地| 蜜臀久久99精品久久宅男| 最近手机中文字幕大全| 精品熟女少妇av免费看| 精品一区二区三区视频在线| 母亲3免费完整高清在线观看 | 91国产中文字幕| 久久久久久人妻| 成人国产av品久久久| 又粗又硬又长又爽又黄的视频| 精品国产一区二区三区四区第35| 国产精品三级大全| 久久av网站| 少妇猛男粗大的猛烈进出视频| 国产高清不卡午夜福利| 国产精品久久久久成人av| 中文天堂在线官网| 欧美成人午夜精品| av免费观看日本| 亚洲国产成人一精品久久久| 国产成人免费观看mmmm| 又黄又爽又刺激的免费视频.| 有码 亚洲区| 国产成人一区二区在线| 久久精品久久精品一区二区三区| 男的添女的下面高潮视频| 国产麻豆69| 伦理电影免费视频| 久久99一区二区三区| 麻豆乱淫一区二区| 日本wwww免费看| 日韩欧美一区视频在线观看| 天堂8中文在线网| 一级毛片黄色毛片免费观看视频| 水蜜桃什么品种好| 18在线观看网站| 亚洲美女搞黄在线观看| 亚洲综合色网址| 人妻一区二区av| 久久午夜福利片| 欧美精品一区二区大全| 啦啦啦啦在线视频资源| 色吧在线观看| 久久ye,这里只有精品| 国产成人91sexporn| 超色免费av| 亚洲国产精品999| 另类精品久久| 欧美人与性动交α欧美精品济南到 | 97在线人人人人妻| 免费大片18禁| 成年美女黄网站色视频大全免费| 超碰97精品在线观看| 黑人高潮一二区| 欧美激情国产日韩精品一区| 99久久精品国产国产毛片| 在线观看三级黄色| 欧美日韩视频精品一区| 国产免费福利视频在线观看| 精品人妻熟女毛片av久久网站| 亚洲丝袜综合中文字幕| 18禁动态无遮挡网站| 99热6这里只有精品| 少妇人妻久久综合中文| 欧美+日韩+精品| 好男人视频免费观看在线| 亚洲av国产av综合av卡| 免费黄频网站在线观看国产| 桃花免费在线播放| 久久鲁丝午夜福利片| 成年人午夜在线观看视频| 搡女人真爽免费视频火全软件| 久热这里只有精品99| 在线观看免费高清a一片| 天天操日日干夜夜撸| 成人二区视频| 黑人猛操日本美女一级片| 九草在线视频观看| 亚洲欧美一区二区三区国产| 日本wwww免费看| 国产探花极品一区二区| 久久久久国产精品人妻一区二区| 一级毛片 在线播放| 九草在线视频观看| 亚洲人成网站在线观看播放| 久久 成人 亚洲| 日韩制服骚丝袜av| 免费人成在线观看视频色| 在线观看美女被高潮喷水网站| 大陆偷拍与自拍| 亚洲精品成人av观看孕妇| 成人国语在线视频| 大香蕉久久网| 欧美性感艳星| 亚洲成国产人片在线观看| av国产久精品久网站免费入址| 天天躁夜夜躁狠狠久久av| 韩国精品一区二区三区 | 国产亚洲精品久久久com| 国产 一区精品| 精品一品国产午夜福利视频| 午夜91福利影院| 日韩中字成人| 久久久欧美国产精品| 国产老妇伦熟女老妇高清| 国产成人91sexporn| 波多野结衣一区麻豆| 99久久中文字幕三级久久日本| 99香蕉大伊视频| a级毛片在线看网站| 欧美精品一区二区大全| 99视频精品全部免费 在线| 大香蕉久久网| 国产成人午夜福利电影在线观看| 黄色视频在线播放观看不卡| 在线观看三级黄色| 成人毛片60女人毛片免费| 国产色婷婷99| 中文乱码字字幕精品一区二区三区| 五月玫瑰六月丁香| 亚洲国产精品一区二区三区在线| 内地一区二区视频在线| 国产成人a∨麻豆精品| av免费观看日本| 又粗又硬又长又爽又黄的视频| av国产精品久久久久影院| 久久综合国产亚洲精品| 亚洲av福利一区| 久久久国产一区二区| 久久精品国产自在天天线| 欧美日韩视频精品一区| 菩萨蛮人人尽说江南好唐韦庄| 一级a做视频免费观看| 亚洲激情五月婷婷啪啪| 日韩成人伦理影院| 久久人人爽人人片av| 一级毛片电影观看| 日韩av在线免费看完整版不卡| 18禁国产床啪视频网站| 天堂俺去俺来也www色官网| 国产片内射在线| 婷婷色麻豆天堂久久| 亚洲,欧美精品.| 韩国高清视频一区二区三区| 成人漫画全彩无遮挡| 大香蕉久久成人网| 免费看光身美女| 久久久久久久大尺度免费视频| 久久韩国三级中文字幕| 99香蕉大伊视频| 美女脱内裤让男人舔精品视频| 老女人水多毛片| 国产亚洲精品久久久com| 自拍欧美九色日韩亚洲蝌蚪91| 国产熟女欧美一区二区| 国产在线一区二区三区精| 性色avwww在线观看| 在线观看三级黄色| 日韩av在线免费看完整版不卡| 国产毛片在线视频| 国产精品一二三区在线看| 91在线精品国自产拍蜜月| 少妇人妻 视频| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三卡| 久久人人爽人人片av| 啦啦啦中文免费视频观看日本| 日韩中字成人| 日本免费在线观看一区| 日韩欧美一区视频在线观看| 性色av一级| 男男h啪啪无遮挡| 国产精品国产三级国产av玫瑰| 成人毛片a级毛片在线播放| 国产精品不卡视频一区二区| 国产成人欧美| 最后的刺客免费高清国语| 欧美日韩视频精品一区| 久久女婷五月综合色啪小说| 91成人精品电影| 亚洲欧美中文字幕日韩二区| 免费大片18禁| av片东京热男人的天堂| 成人毛片a级毛片在线播放| 伦精品一区二区三区| 精品亚洲乱码少妇综合久久| 色婷婷久久久亚洲欧美| 精品第一国产精品| 国产国语露脸激情在线看| 日本免费在线观看一区| 免费看不卡的av| 国产爽快片一区二区三区| 欧美人与性动交α欧美软件 | 在线亚洲精品国产二区图片欧美| 咕卡用的链子| 国产亚洲最大av| 久久精品国产自在天天线| 成年美女黄网站色视频大全免费| 18+在线观看网站| 黑人巨大精品欧美一区二区蜜桃 | 日本黄色日本黄色录像| 91成人精品电影| 久久99热6这里只有精品| 99久久综合免费| 精品一区在线观看国产| 亚洲av成人精品一二三区| 交换朋友夫妻互换小说| 午夜福利,免费看| 一个人免费看片子| 国产精品国产三级国产专区5o| 美女内射精品一级片tv| 欧美另类一区| 亚洲国产精品专区欧美| 免费高清在线观看视频在线观看| 久久国产精品大桥未久av| 国产成人午夜福利电影在线观看| av线在线观看网站| 久久精品夜色国产| 美女内射精品一级片tv| 狂野欧美激情性xxxx在线观看| 亚洲国产精品专区欧美| 2022亚洲国产成人精品| 最近中文字幕2019免费版| 97超碰精品成人国产| 欧美精品一区二区大全| 最新中文字幕久久久久| 国产精品免费大片| 啦啦啦中文免费视频观看日本| 我要看黄色一级片免费的| 国产精品一二三区在线看| 久久精品国产a三级三级三级| 国产一级毛片在线| 亚洲国产色片| 99九九在线精品视频| 天美传媒精品一区二区| 多毛熟女@视频| 亚洲精品色激情综合| 精品国产一区二区久久| 国产极品天堂在线| 天天操日日干夜夜撸| 美女脱内裤让男人舔精品视频| 国产在线视频一区二区| 午夜免费观看性视频| 日韩av不卡免费在线播放| 国产综合精华液| videos熟女内射| 亚洲一码二码三码区别大吗| 搡老乐熟女国产| 婷婷成人精品国产| 在线天堂最新版资源| 国国产精品蜜臀av免费| 激情视频va一区二区三区| 久久国产精品大桥未久av| 久久人人爽人人爽人人片va| 久久99热6这里只有精品| 久久鲁丝午夜福利片| 我的女老师完整版在线观看| 日本欧美国产在线视频| 日韩电影二区| 亚洲熟女精品中文字幕| 卡戴珊不雅视频在线播放| av在线老鸭窝| 一级,二级,三级黄色视频| 国产精品一区二区在线不卡| 高清毛片免费看| 国产福利在线免费观看视频| 国产69精品久久久久777片| 黑人巨大精品欧美一区二区蜜桃 | 欧美国产精品va在线观看不卡| 高清视频免费观看一区二区| 日本wwww免费看| 97人妻天天添夜夜摸| 欧美日本中文国产一区发布| 最新中文字幕久久久久| 看非洲黑人一级黄片| 久久青草综合色| 日韩制服丝袜自拍偷拍| 视频在线观看一区二区三区| 午夜老司机福利剧场| 国产不卡av网站在线观看| 亚洲国产av影院在线观看| 丰满少妇做爰视频| 欧美老熟妇乱子伦牲交| 七月丁香在线播放| 97精品久久久久久久久久精品| 日韩熟女老妇一区二区性免费视频| 久久久久久伊人网av| 黄网站色视频无遮挡免费观看| 香蕉丝袜av| 黑人欧美特级aaaaaa片| 久久精品国产a三级三级三级| 男的添女的下面高潮视频| 亚洲,一卡二卡三卡| 在线观看人妻少妇| 国产精品偷伦视频观看了| 亚洲综合精品二区| 亚洲精品美女久久久久99蜜臀 | av电影中文网址| 免费看不卡的av| 亚洲国产精品专区欧美| 最近的中文字幕免费完整| 五月开心婷婷网| 中国三级夫妇交换| 麻豆精品久久久久久蜜桃| 精品酒店卫生间| 精品久久蜜臀av无| 国产黄色视频一区二区在线观看| 日本黄色日本黄色录像| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 精品国产一区二区三区久久久樱花| 狂野欧美激情性xxxx在线观看| 内地一区二区视频在线| 久久久久久久久久久久大奶| 黄色怎么调成土黄色| 国产高清三级在线| av在线观看视频网站免费| 一本久久精品| 免费看不卡的av| 国产一区二区在线观看av| 成年人午夜在线观看视频| 18+在线观看网站| 亚洲av.av天堂| 青春草国产在线视频| 全区人妻精品视频| 青青草视频在线视频观看| av视频免费观看在线观看| 国产激情久久老熟女| 精品视频人人做人人爽| 久久人人爽人人爽人人片va| 日韩欧美一区视频在线观看| av天堂久久9| 亚洲精品一二三| 搡老乐熟女国产| 蜜桃在线观看..| 国产精品一区二区在线不卡| 亚洲久久久国产精品| 黄片播放在线免费| 成年人免费黄色播放视频| 午夜免费鲁丝| 免费看光身美女| 看免费av毛片| 久久久久久久久久久免费av| 午夜老司机福利剧场| 99视频精品全部免费 在线| 人妻人人澡人人爽人人| 男女免费视频国产| 一二三四在线观看免费中文在 | 高清欧美精品videossex| 欧美日韩一区二区视频在线观看视频在线| 老司机亚洲免费影院| 亚洲欧美一区二区三区黑人 | 精品视频人人做人人爽| 国产片内射在线| 国产精品人妻久久久久久| 久久影院123| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美日韩视频高清一区二区三区二| 深夜精品福利| 在线观看免费高清a一片| 人人妻人人澡人人爽人人夜夜| 91精品国产国语对白视频| 国产欧美亚洲国产| 9191精品国产免费久久| 插逼视频在线观看| 热99国产精品久久久久久7| 亚洲一码二码三码区别大吗| 亚洲欧美日韩另类电影网站| 国产精品偷伦视频观看了| 亚洲成人一二三区av| 亚洲熟女精品中文字幕| 乱人伦中国视频| 最近最新中文字幕免费大全7| 国产一区有黄有色的免费视频| 五月玫瑰六月丁香| 国产午夜精品一二区理论片| 美女xxoo啪啪120秒动态图| 人妻一区二区av| 最近中文字幕高清免费大全6| av有码第一页| 夜夜骑夜夜射夜夜干| 中文字幕另类日韩欧美亚洲嫩草| a 毛片基地| 婷婷色麻豆天堂久久| 啦啦啦啦在线视频资源| 精品国产一区二区久久| 亚洲第一av免费看| 亚洲欧洲日产国产| a级片在线免费高清观看视频| 青青草视频在线视频观看| 久久久久国产网址| 日韩欧美精品免费久久| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 亚洲色图 男人天堂 中文字幕 | 国产精品偷伦视频观看了| 飞空精品影院首页| 美女国产视频在线观看| 夜夜爽夜夜爽视频| 成人18禁高潮啪啪吃奶动态图| 91精品三级在线观看| 少妇猛男粗大的猛烈进出视频| 99视频精品全部免费 在线| 五月天丁香电影| 亚洲激情五月婷婷啪啪| 欧美成人精品欧美一级黄| 久久这里只有精品19| 国产亚洲精品第一综合不卡 | 免费看光身美女| 美女内射精品一级片tv| 国产女主播在线喷水免费视频网站| 狂野欧美激情性bbbbbb| 哪个播放器可以免费观看大片| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕免费大全7| 久久狼人影院| 免费大片18禁| 色网站视频免费| 成年人午夜在线观看视频| 丝袜人妻中文字幕| 亚洲精品国产av蜜桃| 少妇人妻久久综合中文| h视频一区二区三区| 国产高清国产精品国产三级| 国产成人精品福利久久| 热99国产精品久久久久久7| 日韩av不卡免费在线播放| 国产无遮挡羞羞视频在线观看| 99热国产这里只有精品6| 久久久a久久爽久久v久久| 亚洲国产欧美在线一区| 欧美激情 高清一区二区三区| 日日摸夜夜添夜夜爱| 国产有黄有色有爽视频| 国产成人aa在线观看| 国产精品一国产av| 日日啪夜夜爽| 免费av中文字幕在线| 亚洲国产精品专区欧美| tube8黄色片| 精品一区在线观看国产| 五月开心婷婷网| 99久国产av精品国产电影| 久久午夜福利片| 亚洲一码二码三码区别大吗| 亚洲成人一二三区av| 国产成人91sexporn| h视频一区二区三区| 午夜免费鲁丝| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美日韩另类电影网站| 成人国产麻豆网| 欧美97在线视频| 成人18禁高潮啪啪吃奶动态图| 国产精品.久久久| 2021少妇久久久久久久久久久| 国产精品 国内视频| 丰满饥渴人妻一区二区三| 久久精品久久久久久噜噜老黄| 又黄又爽又刺激的免费视频.| 亚洲精品国产色婷婷电影| www.av在线官网国产| 国产在视频线精品| 超色免费av| 乱人伦中国视频| 亚洲欧洲日产国产| 国产一区二区激情短视频 | 久久精品熟女亚洲av麻豆精品| 免费大片18禁| 中文字幕另类日韩欧美亚洲嫩草| 一边摸一边做爽爽视频免费| 搡女人真爽免费视频火全软件| 美女中出高潮动态图| 日韩成人av中文字幕在线观看| 色网站视频免费| 51国产日韩欧美| 五月伊人婷婷丁香| 亚洲第一av免费看| 校园人妻丝袜中文字幕| 欧美xxⅹ黑人| 美女视频免费永久观看网站| 最新中文字幕久久久久| 亚洲三级黄色毛片| 久久青草综合色| 亚洲欧美精品自产自拍| 蜜桃国产av成人99| 精品午夜福利在线看| 欧美+日韩+精品| 高清欧美精品videossex| 亚洲精品,欧美精品| 国产日韩欧美亚洲二区| 久久影院123| 99精国产麻豆久久婷婷| 人妻系列 视频| 久久久久久伊人网av| 亚洲精华国产精华液的使用体验| 成人无遮挡网站| 亚洲精品自拍成人| 国产乱人偷精品视频| 亚洲精品国产色婷婷电影| 免费av不卡在线播放| 少妇 在线观看| 精品国产一区二区久久| 国产又色又爽无遮挡免| 一级毛片我不卡| 亚洲精品久久久久久婷婷小说| 99久久人妻综合| 精品一区在线观看国产| 欧美日韩综合久久久久久| 欧美成人午夜精品| 成人国产av品久久久| 亚洲精品一区蜜桃| 国产亚洲最大av| 亚洲成人av在线免费| 一边摸一边做爽爽视频免费| 桃花免费在线播放| 伦精品一区二区三区| 2021少妇久久久久久久久久久| 欧美日韩视频精品一区| 亚洲国产看品久久| 日日摸夜夜添夜夜爱| 一级片免费观看大全| 黑人巨大精品欧美一区二区蜜桃 | 久久久国产精品麻豆| 一二三四在线观看免费中文在 | 国产精品蜜桃在线观看| 午夜久久久在线观看| 久久久久久久久久久免费av| 免费黄色在线免费观看| 精品一区二区三区四区五区乱码 | 国产亚洲欧美精品永久| 欧美日本中文国产一区发布| 丝袜脚勾引网站| 国产精品秋霞免费鲁丝片| kizo精华| 欧美 亚洲 国产 日韩一| 69精品国产乱码久久久| 久热这里只有精品99| 9191精品国产免费久久| av播播在线观看一区| 久久国产精品男人的天堂亚洲 | 亚洲一码二码三码区别大吗| 18禁观看日本| 男女高潮啪啪啪动态图| 在线天堂最新版资源| 亚洲精品456在线播放app| 亚洲内射少妇av| videosex国产| 99热这里只有是精品在线观看| 日韩中文字幕视频在线看片| www.av在线官网国产| 少妇 在线观看| 成人国语在线视频| 看免费成人av毛片| 亚洲欧美一区二区三区国产| 女人精品久久久久毛片| 国产在线视频一区二区| 狂野欧美激情性bbbbbb| 最近最新中文字幕免费大全7| 不卡视频在线观看欧美| 午夜激情av网站| 女人精品久久久久毛片| 欧美日本中文国产一区发布| 最近2019中文字幕mv第一页| 菩萨蛮人人尽说江南好唐韦庄| 国产成人一区二区在线| 曰老女人黄片| av片东京热男人的天堂| 免费看光身美女| 51国产日韩欧美| 各种免费的搞黄视频| 少妇高潮的动态图| 秋霞在线观看毛片| xxx大片免费视频| 色94色欧美一区二区| 在线观看www视频免费| 免费久久久久久久精品成人欧美视频 | 日韩av在线免费看完整版不卡| 精品亚洲成国产av| 久久久久久久精品精品| 天天躁夜夜躁狠狠久久av| 丝袜喷水一区| 亚洲 欧美一区二区三区| 婷婷色麻豆天堂久久| 天天影视国产精品| 大片免费播放器 马上看| 国产精品无大码| 两个人免费观看高清视频| 高清欧美精品videossex| 亚洲欧美成人综合另类久久久|