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

    A Benzimidazole-pyridine-2,3-dicarboxylic Acid Bridged Zinc(II) Coordination Complex–crystal Structure, Quantum Chemistry and Luminescence①

    2018-10-12 03:50:14WANGJiJunSUNHnYngLIChunBi
    結(jié)構(gòu)化學(xué) 2018年9期

    WANG Ji-Jun SUN Hn-Yng LI Chun-Bi, b

    ?

    A Benzimidazole-pyridine-2,3-dicarboxylic Acid Bridged Zinc(II) Coordination Complex–crystal Structure, Quantum Chemistry and Luminescence①

    WANG Jia-Juna②SUN Han-YangaLI Chuan-Bia, b②

    a(130103)b(136000)

    The structure of a zinc(II)coordinationcomplex (1),[C14H10N3O5Zn1.5]nor [Zn1.5(bzim)(pydc)(H2O)]n(H2pydc = pyridine-2,3-dicarboxylic acid, Hbzim = benzimidazole), has been determined by X-ray crystallography and characterized by elemental analysis, IR spectrum and luminescence.Chemical formula: C14H10N3O5Zn1.5.It crystallizes in the monoclinic system, space group21/with= 12.303(4),= 12.052(4),= 10.212(3) ?,= 104.147(4),= 1468.3(8) ?3,= 4,M= 398.30,D= 1.802 g/cm3,(000) = 800,= 2.501 mm-1and= 1.000.The 2-D network architecture of 1 is constructed from benzimidazole, zinc and pyridine- 2,3-dicarboxylic acid.The quantum-chemical calculations have been performed on ‘molecular fragments’ extracted from the crystal structure using the B3LYP method in Gaussian 09.The luminescence spectrum shows that complex 1 emits blue luminescence.

    pyridine-2,3-dicarboxylic acid, benzimidazolezinc(II) complex, quantum chemistry study, luminescence spectrum;

    1 INTRODUCTION

    Crystal engineering of supramolecular architec- tures based on metal organic frameworks (MOFs) have been fast expanding in the last few years owing to their potential applications such as absorp- tion[1], catalysis[2], luminescence[3],.The history ofMOFs used in the field of removing the highly toxic heavy metal ions is not long but has rapid development and broad prospects[4-9].MOFs exhi- biting luminescence emission in solid state is useful, especially, the fluorescence of the zinc(II) ion containing MOFs can be used for not only selective detection and differentiating Fe(III) and Cr(VI) ions[10], but also high iodine capture and nitro- explosive detection[11, 12, 13].

    Recently, many high-dimensional coordination complexes have been designed and prepared through molecular self-assembly process[14].Struc- turally, the weaker intermolecular forces, especially hydrogen bonds, play an important role in funda- mental biological processes in the field of bioche- mistry and supramolecular chemistry[15].The multicarboxylate ligands have been proved to be good candidates because they can be regarded as hydrogen-bonding accepters and hydrogen-bonding donors, dependent upon the number of deprotonated carboxylic groups.The imidazole-like ligand can coordinate to metal center, and it can even be deprotonated and act as an anion[16, 17].

    The hetero-ring carboxylic acid, pyridine-2,3- dicarboxylic acid (H2pydc) is the precursor of adenine-nicotinamide dinucleotide[18].It has great potential for coordinative interactions and hydrogen bonding, and hence may result in a large diversity of supramolecular architectures.Benzimidazole (Hbzim) is a fungicide[19].As a unique ligand it shows binding affinity to transition metals[20].Moreover, stable complexes are known to be formed as a neutral molecule or act as a benzimida- zole anion (bzim)[21].

    In this paper, in view of the weak water solubility of the ligand used, we adopt hydrothermal methods to carry out our experiment.We report a new coordination complex, [C14H10N3O5Zn1.50]nor [Zn1.5(bzim)(pydc)(H2O)]n(1, H2pydc = pyridine- 2,3-dicarboxylic acid, Hbzim = benzimidazole).Complex 1 exhibits a two-dimensional (2D)hydro- gen bonding network structure in theplane in the crystal lattice.

    2 EXPERIMENTAL

    2.1 Synthesis of complex 1

    Zn(Ac)2·2H2O (0.220 g, 1 mmol),benzimidazole(Hbzim, 0.118 g, 1 mmol), pyridine-2,3-dicar- boxylic acid(H2pydc, 0.167 g, 1 mmol) and 16 mL water were mixed with stirring followed by adjusting the pH value to 8 with an aqueous solution of NaOH.Then the mixture was sealed in a 25 mL Teflon-lined stainless-steel reactor and heated at 110 ℃ for 96 h to give colorless crystals of the title complex after cooling.The C, H and N contents were determined by elemental analysis: Calcd.(%) forC14H10N3O5Zn1.5: C, 42.18; H, 2.51; N, 10.54; Zn 24.63.Found (%): C, 41.95; H, 2.47; N, 10.39; Zn 24.98.IR(KBr, cm-1) 3395-3478br-vs, 3434vs, 2115m, 2101-1834w, 1826s, 1679s, 1625m, 1605s, 1596w, 1510m, 1365m, 1263m, 839m, 783scm-1.

    2.2 Structure determination and physical measurements

    A colorless block crystal for 1 with dimensions of 0.32mm × 0.27mm × 0.23mm was chosen for X-ray diffraction analysis.Crystal structure measurement was performed on a Bruker SMART APEX II CCD diffractometer equipped with a graphite-monochro- matic Mo(= 0.71073 ?) radiation by using anscan mode at 292(2) K.Absorption corrections were applied with a multi-scan mode[22].A total of 12414 reflections were obtained in the rangeof 2.40≤≤26.06o, of which 2907 were independent (int= 0.0722) and 2303 observed reflections with> 2()were employed for structuredetermination and refinement.The structure was solved by direct methods with SHELXS-97[23]and refined by full-matrix least-squares techniques using SHELXL-97 program[24]within WINGX[25].All non-hydrogen atoms were refined anisotropically.All H atoms on C atoms were positioned geometri- cally and refined as riding, with C–H = 0.93 andiso(H) = 1.2eq(C).The final= 0.0501,= 0.0743 (= 1/[2(F2) + (0.0265)2+ 1.5096], where= (F2+ 2F2)/3), (Δ/)max=0.000,= 1.000, (Δ)max= 0.362 and (Δ)min= –0.464 e/?3.The FT-IR spectrum was recorded from KBr pellets in the range of 4000~400 cm–1on a Perkin-Elmer 240C spectrometer.TG was performed using a Perkin-Elmer TG-7 analyzer in nitrogen with a heating rate of 10oC/min.The luminescence spectra for the powdered solid samples were measured at room temperature and the spectra were collected with a Perkin-Elmer LS-55 fluorescence spectro- meter.Other reagents were of analytical grade.

    2.3 Quantum chemistry calculation

    The quantum chemistry calculation of complex 1 was performed with Gaussian09 program[26]at the B3LYP/GenECP level, the 6-31+G(d) basis set for C, H, O, N and LANL2DZbasis set for Zn.For modeling the initial guess of the title complex was obtained from the X-ray refinement data (cif).In order to save computing time, we adopt the mini- mum structural unit of the complex to make an approximate structure model in the calculation, and use the neutral Hbzim instead of the deprotonatedbzim anion and the “one hydrogen atom contained” H2pydcanion ligand (that is, Hpydc) rather than the pydc dianionfor simplifying the polymeric complex 1 to a “monomer”, andkeep the charge balance in the model at the same time.

    3 RESULTS AND DISCUSSION

    3.1 Crystal structure of [Zn1.5(bzim)(pydc)(H2O)]n (1)

    The structure of complex 1 is described in Fig.1.It has two Zn centers, however, the coordination manners of them are different.Zn(1) is four-coor- dinated by two N and two O atoms, forming a tetrahedral coordination geometry.The two N atoms are from two different benzimidazole and the two O atoms are from two different pyridine-2,3-dicar- boxylic acids; while the Zn(2) is located at an inversion center and it is six-coordinated by two N and four O atoms, with two N and two O atoms from two different pyridine-2,3-dicarboxylic acids.The Npyridineand 2 site carboxylic O atoms from one ligand coordinated to Zn(2) in a chelated mode, while the other two O atoms are from two coor- dinated water molecules.For the Zn2-center, the two chelate O(4), N(1) atoms and two chelate O(4i), N(1i) atoms (symmetry code i:1–, –1–, 1–) are all located at the equatorial positions, while Owand Owioccupy two axial sites.This coordination manner formed an octahedral geometry (Fig.1).The bond lengths and bond angles are listed in Table 1.From Table 1, we can see that the Zn–N(bipy)distance is 2.067 ?, which agrees with the normal Zn–N(bipy)distance (2.064 ?)[27].The Zn–N(imidazole)distances are 1.986 and 1.990 ?, also in accord with the-imidazole type Zn–N(imidazole)distances (1.986 and 1.994 ?[28]).Zn···O(carboxyl)distances vary from 1.976 to 2.077 ?, which consist with the typical Zn–O(carboxyl)distance (2.010 ?[27]), and the Zn–Owbond (2.210 ?)is slightly longer than Zn–O(carboxyl), indicating that the former is stronger than the latter.With the connecting manner described above, every benzimidazole is deprotonized and connects two Zn(1) atoms, and every pyridine-2,3-dicarboxylic acid connects three Zn atoms.The two are Zn(1) atoms with each carboxyl in a single dentate coordination manner and another one is Zn(2) atom chelated by N and O.In this connecting mode, complex 1 formed a two-dimensional network[29]structure in theplane.In the network structure, there are O(w1)– H(W1A)???O(2iv) and O(w1)– H(W1B)???O(1v) (iv: 1–, –1–, –, v: 1–, –0.5+, 0.5–, in Table 2) hydrogen bonds which may stabilize[30]the structure of 1 extending into a2Dlayered supramolecular architecture (Fig.2).

    Fig.1. Coordination geometry of the coordination complex 1.

    (Symmetry codes for i: 1–, –1–, 1–;ii: 1–, –, 1–;iii:, 0.5–, 0.5+)

    Fig.2. 2D network supramolecular structure in theplane of complex 1 and hydrogen bonds (dotted line)

    Table 1. Selected Bonds Lengths (?) and Bond Angles (°)

    S ymmetry transformation:i:1–, –1–, 1–; ii: 1–, –, 1–; iii:, 0.5–, 0.5+

    Table 2. Hydrogen Bond for Complex 1 (? and °)

    Symmetry transformations used to generate the equivalent atoms: iv: 1–, –1–, –; v: 1–, –0.5+, 0.5–

    3.2 Quantum chemistry study

    The calculation covered 85 atoms, 1098 basis functions, 1978 primitive gaussians, 211electrons and 211electrons for the model of complex 1.As the model does not have a single electron, the spin multiplicity is 1.The total molecular energy is–2917.983 a.u., the energies of HOMO and LUMO are –0.227 and –0.112 a.u., respectively, with the Δ(ELUMO–EHOMO) value to be –0.115 a.u., which shows the complex is stable[31]in the ground state.The HOMO and LUMO are presented inFig.3, from which we can see the HOMO electron cloud is mainly located at bzim and the LUMO electron cloud at the pydc ligand.

    Selected atom net charges and electronic con- figuration of the title complex at the B3LYP/(6- 31+G(d) (for C, H, N and O) and Lanl2dz (for Zn) levels are listed in Table 3.The calculation results show that electronic configurations of the central Zn atoms are 40.25, 0.3139.9840.30, 0.34, and those of O and N atoms are 21.68~1.7425.02~5.3230.0130~0.01and 21.33~1.3724.18~4.3440.01, 0.02.The valence of zinc is +2.The zinc partially obtains electrons from the pydc and bzim ligands, illustrating that the charge of the central Zn metals for complex 1 has positive charge, and it is 1.41307 for Zn(1) and 1.41897 for Zn(2).In Table 3, the net charges of the coordinated O and N atoms are all negative.The charges of O(1), O(3ii), N(2), N(3iii) around Zn(1) are –0.75115, –0.74627, –0.62411, –0.72322, and those of N(1), N(1i), O(4), O(4i), O(w1), O(w1i) around Zn(2) are –0.54178, –0.55373, –0.75661, –0.85343, –0.99498, –1.05467, respectively.

    Fig.3. Molecular orbitals of complex 1

    Table 3. Selected Atom Net Charges and Electronic Configuration of Complex 1 at the B3LYP/(6-31+G(d) (for C, H, N and O) and Lanl2dz (for Zn) Levels

    Symmetry codes:i: 1–, –1–, 1–; ii: 1–, –, 1–; iii:, 0.5–, 0.5+

    3.4 Luminescence spectra

    Luminescence is an important property[32].To examine the luminescent properties of10metal complexes, the solid-state emission and excitation (inset) spectra of complex 1 at room temperature are depicted in Fig.4.The obvious emissions are observed at 466 nm (= 520 nm) for 1.To under- stand the nature of the emission spectra, the luminescence properties of the free Hbzim and H2pydc ligands under the same conditions were recorded for comparison.The Hbzim and H2pydc ligands exhibit emission bands located at 396 and 394 nm upon 570 and 560 nm excitation, respec- tively.These ligands are observed without emission bands above 460 cm-1.The emissions of complex 1 may be assigned to the metal-to-ligand charge- transfer (MLCT).

    The luminescence spectrum of complex 1 exhibits a red shift in solid state, compared with that of the two ligands (H2pydc and Hbzim).The reason may be that the10metal ion [zinc(II)] coordinates to the ligands, which probably forms the back-coupling-bond between the metal and ligands, and decreases the electron transition energy of intraligand charge transfer.Otherwise, the ligand coordinated with metal ions also forms additional five-membered rings, which also increases the* conjugation length and the conformational coplanarity of the ligand, accor- dingly reduces the energy gap between theand* molecular orbitals of the ligand[33, 34].

    Fig.4. Solid-state emission and excitation (inset) spectra of complex 1 and ligands at room temperature

    4 CONCLUSION

    The structure of a zinc(II)coordinationcomplex (1, [C14H10N3O5Zn1.50]n) has been determined by X-ray crystallography and characterized by the luminescence.The conclusions of this article are as follows: (a) It is a 2-D layered network architecture constructed from benzimidazole, zinc and pyridine- 2,3-dicarboxylic acid.(b) The calculations on ‘molecular fragments’ are extracted from the crystal structure using B3LYP method.The calculation reveals that the HOMO and LUMO orbitals are located at the bzim and pydc ligands, respectively.The atom net charges and electronic configuration are also shown.(c) The luminescence spectrum shows that complex 1 emits blue luminescence, and the luminescence mechanism of the emissions may be assigned to the metal-to-ligand charge-transfer (MLCT).

    (1) Chen, B.L.; Ockwig, N.W.; Millward, A.R.; Contreras, D.S.; Yaghi, O.M.High H2adsorption in a microporous metal-organic framework with open metal sites.2005, 44, 4745–4749.

    (2) Farrusseng, D.; Aguado, S.; Pinel, C.Metal-organic frameworks: opportunities for catalysis.2009, 48, 7502–7513.

    (3) Nagarkar, S.S.; Joarder, B.; Chaudhari, A.K.; Mukherjee, S.; Ghosh, S.K.Highly selective detection of nitro explosives by a luminescent metal-organic framework.2013, 52, 2881–2885.

    (4) Yu, C.X.; Shao, Z.C.; Hou, H.W.A functionalized metal-organic framework decorated with O?groups showing excellent performance for lead(II) removal from aqueous solution.2017, 8, 7611–7619.

    (5) Wang, H.R.; Meng, W.; Wu, J.; Ding, J.; Hou, H.W.; Fan, Y.T.Crystalline central-metal transformation in metal-organic frameworks.2016, 307, 130–146.

    (6) Meng, W.; Li, H.J.; Xu, Z.Q.; Du, S.S.Li, Y.X.; Zhu, Y.Y.; Han, Y.; Hou, H.W.; Fan, Y.T.; Tang, M.S.New mechanistic insight into stepwise metal-center exchange in a metal-organic framework based on asymmetric Zn4clusters.2014, 20, 2945–252.

    (7) Han, Y.; Chilton, N.F.; Li, M.; Huang, C.; Xu, H.; Hou, H.W.; Moubaraki, B.; Langley, S.K.; Batten, S.R.; Fan, Y.T.; Murray, K.S.Post-synthetic monovalent central-metal exchange, specific I2sensing, and polymerization of a catalytic [3×3] grid of [CuII5CuI4L6](I)2·13H2O.2013, 19, 6321–6328.

    (8) Mi, L.W.; Hou, H.W.; Song, Z.Y.; Han, H.Y.; Fan, Y.T.Polymeric zinc ferrocenyl sulfonate as a molecular aspirator for the removal of toxic metal ions.2008, 14, 1814–1821.

    (9) Zhao, J.A.; Mi, L.W.; Hu, J.Y.; Hou, H.W.; Fan, Y.T.Cation exchange induced tunable properties of a nanoporous octanuclear Cu(II) wheel with double-helical structure.2008, 130, 15222–15223.

    (10) Lv, R.; Li, H.; Su, J.; Fu, X.; Yang, B.; Gu, W.; Liu, X.Zinc metal-organic framework for selective detection and differentiation of Fe(III) and Cr(VI) ions in aqueous solution.2017, 56, 12348–12356.

    (11) Yao, R.X.; Cui, X.; Jia, X.X.; Zhang, F.Q.; Zhang, X.M.A luminescent zinc(II) metal-organic framework (MOF) with conjugated-electron ligand for high iodine capture and nitro-explosive detection.2016, 55, 9270–9275.

    (12) Wang, H.; Yang, W.T.; Sun, Z.M.Mixed-ligand Zn-MOFs for highly luminescent sensing of nitro compounds.2013, 8, 982–989.

    (13) Xing, S.; Bing, Q.; Qi, H.; Liu, J.; Bai, T.; Li, G.; Shi, Z.; Feng, S.; Xu, R.Rational design and functionalization of a zinc metal-organic framework for highly selective detection of 2,4,6-trinitrophenol.2017, 9, 23828–23835.

    (14) Avci, C.; Imaz, I.; Carné-Sánchez, A.; Pariente, J.A.; Tasios, N.; Pérez-Carvajal, J.; Alonso, M.I.; Blanco, A.; Dijkstra, M.; López, C.; Maspoch, D.Self-assembly of polyhedral metal-organic framework particles into three-dimensional ordered superstructures.2018, 10, 78–84.

    (15) Lee, H.R.; Helquist, S.A.; Kool, E.T.; Johnson, K.A.Importance of hydrogen bonding for efficiency and specificity of the human mitochondrial DNA polymerase.2008,283, 402–14410.

    (16) Huang, H.W.; Tsou, C.C.; Kuo, T.S.; Liaw, W.F.New members of a class of dinitrosyliron complexes (DNICs): interconversion and spectroscopic discrimination of the anionic {Fe(NO)2}9[(NO)2Fe(C3H3N2)2]?and [(NO)2Fe(C3H3N2)(SR)]?(C3H3N2= deprotonated imidazole, R =Bu, Et, Ph).2008, 47, 2196–2204.

    (17) Lai, S.W.; Chan, M.C.W.; Peng, S.M.; Che, C.M.Self-assembly of predesigned trimetallic macrocycles based on benzimidazole as nonlinear bridging motifs: crystal structure of a luminescent platinum(II) cyclic trimer.1999, 38, 669–671.

    (18) Heyes, M.P.Metabolism and neuropathologic significance of quinolinic acid and kynurenic acid.1993, 21, 83–89.

    (19) Davidse, L.C.Benzimidazole fungicides: mechanism of action and biological impact..1986, 24, 43–65.

    (20) Toth, A.; Floriani, C.; Chiesi-Villa, A.; Guastinil, C.Copper(I)-benzimidazole adducts: from mononuclear to polymeric complexes.1987, 26, 3897–3902.

    (21) Carina, R.F.; Verzegnassi, L.; Bernardinelli, G.; Williams, A.F.Modulation of iron reduction potential by deprotonation at a remote site.1998, 2681–2682.

    (22) Higashi, T..Rigaku Corporation, Tokyo, Japan 1995.

    (23) Sheldrick, G.M..University of G?ttingen, Germany 1997.

    (24) Sheldrick, G.M..University of G?ttingen, Germany 1997.

    (25) Farrugia, L.J.; Wing, X.A., University of Glasgow, Glasgow, UK 1988.

    (26) Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H.P.; Izmaylov, A.F.; Bloino, J.; Zheng, G.; Sonnenberg, J.L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, Jr.J.A.; Peralta, J.E.; Ogliaro, F.; Bearpark, M.; Heyd, J.J.; Brothers, E.; Kudin, K.N.; Staroverov, V.N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J.C.; Iyengar, S.S.; Tomasi, J.; Cossi, M.; Rega, Millam, N.J.; Klene, M.; Knox, J.E.; Cross, J.B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.E.; Stratmann, O.; Yazyev, A.J.; Austin, R.; Cammi, C.; Pomelli, J.W.; Ochterski, R.; Martin, R.L.; Morokuma, K.; Zakrzewski, V.G.; Voth, G.A.; Salvador, P.; Dannenberg, J.J.; Dapprich, S.; Daniels, A.D.; Farkas, O.; Foresman, J.B.; Ortiz, J.V.; Cioslowski, J.; Fox, D.J.Gaussian, Inc.,2009.Gaussian 09, Revision D.01.

    (27) Orpen, A.G.; Brammer, L.; Allen, F.H.; Kennard, O.; Watson, G.; Taylor, R.Tables of bond lengths determined by X-ray and neutron diffraction.Part 2.Organometallic compounds and Co-ordination complexes of the d- and f-block metals.1989, S1–S83.

    (28) Ashby, C.I.H.; Cheng, C.P.; Duesler, E.N.; Brown, T.L.Crystal structure and14N nuclear quadrupole resonance spectrum of catena--imidazolato-bis(imidazo1e)zinc nitrate.Donor characteristics of coordinated imidazolate.1978, 100, 6063–6067.

    (29) Kong, Z.G.; Liu, D.X.; Song, M.Y.; Zhang, L.Y.; Feng, S.Y.; Xu, Z.L.A two-dimensional supramolecular complex constructed by 1,10-phenanthroline derivative ligand: synthesis, structure and luminescence.2017, 36, 1859–1863.

    (30) Jurecka, P.; Hobza, P.True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine···cytosine, adenine···thymine, and their 9- and 1-methyl derivatives:? complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment.2003, 125, 15608–15613.

    (31) Fukui, K.Role of frontier orbitals in chemical reactions.1982, 218, 747–754.

    (32) Kong, Z.G.; Liu, D.X.; Li, R.; Jiang, Y.; Wang, L.J.; Hu, B.A new one-dimensional coordination polymer based on 3,5-dinitro-salicylic acid: synthesis, crystal structure, luminescent property and theoretical calculation.2017, 36, 841–847.

    (33) Perkovic, M.W.Allosteric manipulation of photoexcited state relaxation in (bpy)2RuII(binicotinic acid).2000, 39, 4962–4968.

    (34) Yuan, G.Z.; Huo, Y.P.; Nie, X.L.; Jiang, H.; Liu, B.; Fang, X.M.; Zhao, F.H.Controllable supramolecular structures and luminescent properties of unique trimeric Zn(II) 8-hydroxyquinolinates tuned by functional substituents.2013, 42, 2921–2929.

    18 January 2018;

    26 April 2018 (CCDC 689913)

    ① Supported by the Jilin Province Science and Technology Development Plan Item (No.20140204080GX), the Project of the Education Department of Jilin Province, China (No.JJKH20180777KJ) and the Science and Technology Development Projects of(No.2017057)

    Wang Jia-Jun, Tel: +86 0431 81765105.E-mail: jiajunwang@jlnu.edu.cn.Li Chuan-Bi, E-mail: li_c_b@163.com

    10.14102/j.cnki.0254-5861.2011-1955

    av在线天堂中文字幕| netflix在线观看网站| 午夜激情欧美在线| 久久九九热精品免费| 一边摸一边抽搐一进一小说| cao死你这个sao货| 成人午夜高清在线视频| 久久久久久久午夜电影| 国产三级黄色录像| 久久中文字幕人妻熟女| 综合色av麻豆| 国产精品爽爽va在线观看网站| 99在线人妻在线中文字幕| 日韩欧美三级三区| 97超视频在线观看视频| 久久久久久久久久黄片| 亚洲18禁久久av| 在线免费观看的www视频| 亚洲自拍偷在线| 亚洲乱码一区二区免费版| 午夜福利在线在线| 日本五十路高清| 亚洲av五月六月丁香网| 国产毛片a区久久久久| 国产毛片a区久久久久| 全区人妻精品视频| ponron亚洲| 男女午夜视频在线观看| 男女下面进入的视频免费午夜| 欧美黑人欧美精品刺激| 亚洲国产欧美人成| 国产精品精品国产色婷婷| 两个人视频免费观看高清| 最近最新免费中文字幕在线| www.www免费av| 精品熟女少妇八av免费久了| 国产美女午夜福利| 免费在线观看亚洲国产| 99久国产av精品| 色av中文字幕| 成熟少妇高潮喷水视频| 日韩欧美在线二视频| 国内久久婷婷六月综合欲色啪| 久久久精品大字幕| 国产成人精品久久二区二区免费| 丁香六月欧美| 国产av麻豆久久久久久久| 国产伦精品一区二区三区四那| 亚洲 欧美一区二区三区| 99精品在免费线老司机午夜| 高潮久久久久久久久久久不卡| 亚洲人成网站高清观看| 少妇的逼水好多| 啦啦啦韩国在线观看视频| 脱女人内裤的视频| 亚洲av片天天在线观看| 亚洲av片天天在线观看| 亚洲一区二区三区不卡视频| 99国产精品99久久久久| 看片在线看免费视频| 99精品在免费线老司机午夜| 亚洲中文日韩欧美视频| 亚洲国产日韩欧美精品在线观看 | 香蕉丝袜av| 亚洲九九香蕉| 日本a在线网址| 最近最新免费中文字幕在线| 一本久久中文字幕| 国产真实乱freesex| 国产美女午夜福利| 久久久久精品国产欧美久久久| 在线国产一区二区在线| 熟妇人妻久久中文字幕3abv| 黑人操中国人逼视频| 99久久99久久久精品蜜桃| 欧美精品啪啪一区二区三区| 国产精品99久久久久久久久| 网址你懂的国产日韩在线| 午夜福利在线观看免费完整高清在 | 我要搜黄色片| 国产欧美日韩精品亚洲av| 国产亚洲欧美在线一区二区| 久久久国产精品麻豆| а√天堂www在线а√下载| 母亲3免费完整高清在线观看| 国产精华一区二区三区| 国产成人一区二区三区免费视频网站| 国产69精品久久久久777片 | 日本五十路高清| 中文亚洲av片在线观看爽| 婷婷六月久久综合丁香| 观看免费一级毛片| 色av中文字幕| 精品电影一区二区在线| 看片在线看免费视频| 日韩欧美在线乱码| 精品福利观看| 国产午夜精品久久久久久| 国产爱豆传媒在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 哪里可以看免费的av片| 青草久久国产| 别揉我奶头~嗯~啊~动态视频| 欧美色欧美亚洲另类二区| 特大巨黑吊av在线直播| 国产人伦9x9x在线观看| 免费观看精品视频网站| 91在线观看av| 免费无遮挡裸体视频| 国产极品精品免费视频能看的| 亚洲精品在线美女| 日韩欧美国产一区二区入口| 国产一区二区在线观看日韩 | 成人特级黄色片久久久久久久| 免费在线观看影片大全网站| 婷婷精品国产亚洲av在线| 五月伊人婷婷丁香| 精品乱码久久久久久99久播| 麻豆国产97在线/欧美| 中文字幕熟女人妻在线| 亚洲av熟女| 欧美成狂野欧美在线观看| e午夜精品久久久久久久| 亚洲国产精品久久男人天堂| 最新美女视频免费是黄的| 亚洲午夜理论影院| 国产精品久久久久久久电影 | 9191精品国产免费久久| 中文字幕高清在线视频| 一级黄色大片毛片| 大型黄色视频在线免费观看| 国产一区二区在线观看日韩 | 日韩欧美 国产精品| 999久久久国产精品视频| 此物有八面人人有两片| 99re在线观看精品视频| 亚洲av免费在线观看| 观看美女的网站| 国产亚洲精品久久久com| 国产视频内射| 久久午夜综合久久蜜桃| 99久久99久久久精品蜜桃| 又黄又爽又免费观看的视频| 亚洲中文字幕日韩| 神马国产精品三级电影在线观看| 99re在线观看精品视频| 成年女人永久免费观看视频| 国产成人精品久久二区二区91| 午夜精品在线福利| 熟女人妻精品中文字幕| 欧美3d第一页| 丁香六月欧美| 国产精品98久久久久久宅男小说| 国产高清视频在线播放一区| 日韩 欧美 亚洲 中文字幕| 欧美黄色片欧美黄色片| 舔av片在线| 在线观看美女被高潮喷水网站 | 国产成人福利小说| 国模一区二区三区四区视频 | 老汉色∧v一级毛片| 最新中文字幕久久久久 | 欧美色欧美亚洲另类二区| 国产 一区 欧美 日韩| 中文资源天堂在线| 亚洲中文字幕一区二区三区有码在线看 | 色综合亚洲欧美另类图片| 琪琪午夜伦伦电影理论片6080| 午夜免费激情av| 欧美丝袜亚洲另类 | 国产av麻豆久久久久久久| 又紧又爽又黄一区二区| 日本五十路高清| 长腿黑丝高跟| 久久久久久久精品吃奶| 别揉我奶头~嗯~啊~动态视频| 精品午夜福利视频在线观看一区| 怎么达到女性高潮| 又黄又爽又免费观看的视频| 又粗又爽又猛毛片免费看| 国产一区二区在线观看日韩 | 12—13女人毛片做爰片一| 日韩 欧美 亚洲 中文字幕| 欧美性猛交黑人性爽| 久久午夜亚洲精品久久| 国产精品久久电影中文字幕| 国产亚洲精品av在线| av在线蜜桃| 午夜a级毛片| 国产一区二区三区在线臀色熟女| 久久天躁狠狠躁夜夜2o2o| 黄色片一级片一级黄色片| 日韩欧美在线二视频| 亚洲中文日韩欧美视频| 国产精品亚洲一级av第二区| 97碰自拍视频| 天堂√8在线中文| 国产精品久久视频播放| 在线看三级毛片| 美女被艹到高潮喷水动态| 亚洲人成网站在线播放欧美日韩| 一a级毛片在线观看| 久久久久久大精品| 亚洲成av人片免费观看| 国产v大片淫在线免费观看| 看片在线看免费视频| 夜夜看夜夜爽夜夜摸| 婷婷六月久久综合丁香| 国产伦人伦偷精品视频| 日韩有码中文字幕| 欧美一级a爱片免费观看看| 在线a可以看的网站| 丁香欧美五月| 国产一区二区在线av高清观看| 成人无遮挡网站| 色av中文字幕| 在线观看舔阴道视频| 黄色 视频免费看| 亚洲国产欧洲综合997久久,| a在线观看视频网站| 国产精品永久免费网站| 欧美日韩精品网址| 日本五十路高清| 精品久久久久久久毛片微露脸| 国产成人aa在线观看| 亚洲 欧美一区二区三区| 动漫黄色视频在线观看| 男女做爰动态图高潮gif福利片| 成在线人永久免费视频| 无人区码免费观看不卡| 国产综合懂色| 一个人看的www免费观看视频| 久久久国产欧美日韩av| 国产黄a三级三级三级人| 男人的好看免费观看在线视频| 亚洲第一电影网av| 婷婷亚洲欧美| 国产野战对白在线观看| 久久久精品欧美日韩精品| 国产淫片久久久久久久久 | 手机成人av网站| 人人妻人人澡欧美一区二区| 亚洲18禁久久av| 亚洲精品久久国产高清桃花| 校园春色视频在线观看| 成人永久免费在线观看视频| 日韩欧美国产一区二区入口| 九九热线精品视视频播放| 99精品久久久久人妻精品| 脱女人内裤的视频| 久久香蕉精品热| 九九久久精品国产亚洲av麻豆 | 国产成人欧美在线观看| 91在线观看av| 大型黄色视频在线免费观看| 哪里可以看免费的av片| 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 激情在线观看视频在线高清| 婷婷六月久久综合丁香| 搡老妇女老女人老熟妇| 少妇的丰满在线观看| 亚洲自偷自拍图片 自拍| 久久欧美精品欧美久久欧美| 美女高潮喷水抽搐中文字幕| 最近最新免费中文字幕在线| 国产视频内射| 天堂av国产一区二区熟女人妻| 精品国产亚洲在线| 亚洲国产欧洲综合997久久,| 亚洲五月婷婷丁香| xxx96com| 成人特级黄色片久久久久久久| 国产蜜桃级精品一区二区三区| 色精品久久人妻99蜜桃| 国产精品永久免费网站| 国语自产精品视频在线第100页| 18禁黄网站禁片免费观看直播| 欧美性猛交╳xxx乱大交人| 99久久国产精品久久久| 天堂av国产一区二区熟女人妻| 免费在线观看影片大全网站| 亚洲精品色激情综合| 一本久久中文字幕| 丰满人妻熟妇乱又伦精品不卡| 欧美三级亚洲精品| 99热只有精品国产| 国产亚洲精品一区二区www| 老汉色∧v一级毛片| 天堂√8在线中文| 中亚洲国语对白在线视频| 欧美日韩亚洲国产一区二区在线观看| 搡老熟女国产l中国老女人| 男女午夜视频在线观看| 少妇裸体淫交视频免费看高清| 欧美丝袜亚洲另类 | 亚洲欧美日韩东京热| 99re在线观看精品视频| 国产精品 国内视频| 国产久久久一区二区三区| 亚洲专区国产一区二区| 亚洲美女视频黄频| 亚洲五月天丁香| 色尼玛亚洲综合影院| 网址你懂的国产日韩在线| 午夜a级毛片| 嫩草影院入口| 欧美中文综合在线视频| 亚洲av五月六月丁香网| 一进一出好大好爽视频| 这个男人来自地球电影免费观看| 亚洲第一电影网av| 久久精品91蜜桃| 国产高清视频在线观看网站| 久久午夜亚洲精品久久| 99热精品在线国产| 国产乱人视频| 99国产综合亚洲精品| 免费看日本二区| 婷婷亚洲欧美| 久久精品国产清高在天天线| 香蕉国产在线看| 看黄色毛片网站| 波多野结衣高清无吗| 日本a在线网址| 国内久久婷婷六月综合欲色啪| 999久久久精品免费观看国产| 色综合站精品国产| 亚洲成a人片在线一区二区| 日韩精品中文字幕看吧| 精品乱码久久久久久99久播| 国产淫片久久久久久久久 | 欧美av亚洲av综合av国产av| 国产一区二区三区在线臀色熟女| 亚洲成a人片在线一区二区| 久久久久国产精品人妻aⅴ院| 久久99热这里只有精品18| bbb黄色大片| 在线观看舔阴道视频| 成人三级黄色视频| 久久精品91无色码中文字幕| 国产在线精品亚洲第一网站| 免费在线观看影片大全网站| 最近视频中文字幕2019在线8| 欧美一级毛片孕妇| av福利片在线观看| 亚洲18禁久久av| 桃色一区二区三区在线观看| 成人永久免费在线观看视频| 男女做爰动态图高潮gif福利片| 欧美乱妇无乱码| 搡老熟女国产l中国老女人| 一a级毛片在线观看| 日本一二三区视频观看| 精品一区二区三区av网在线观看| 久久久久久九九精品二区国产| 99国产综合亚洲精品| 欧美中文日本在线观看视频| 无限看片的www在线观看| 久久精品国产亚洲av香蕉五月| 午夜免费激情av| 国产精品亚洲美女久久久| 在线观看一区二区三区| 国产亚洲精品久久久com| 国产精品野战在线观看| 亚洲乱码一区二区免费版| 99久久无色码亚洲精品果冻| 中文在线观看免费www的网站| 亚洲国产高清在线一区二区三| 国产一区二区三区在线臀色熟女| 法律面前人人平等表现在哪些方面| 全区人妻精品视频| or卡值多少钱| 中出人妻视频一区二区| 动漫黄色视频在线观看| av在线天堂中文字幕| www.熟女人妻精品国产| 一区二区三区高清视频在线| 91av网站免费观看| 高清在线国产一区| 伊人久久大香线蕉亚洲五| 岛国在线免费视频观看| 久久午夜综合久久蜜桃| 日韩欧美国产一区二区入口| 老鸭窝网址在线观看| 99热精品在线国产| 亚洲精品在线观看二区| 国模一区二区三区四区视频 | 88av欧美| 18禁动态无遮挡网站| 免费观看在线日韩| 亚洲av一区综合| 天堂影院成人在线观看| 国产午夜福利久久久久久| 久久99热这里只频精品6学生 | .国产精品久久| 久久久精品欧美日韩精品| 中文在线观看免费www的网站| 只有这里有精品99| 亚洲欧美日韩卡通动漫| av在线老鸭窝| 国产一区二区三区av在线| 亚洲国产色片| 波野结衣二区三区在线| 99九九线精品视频在线观看视频| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版 | 人人妻人人澡欧美一区二区| 午夜福利视频1000在线观看| 久久精品国产亚洲av涩爱| 免费不卡的大黄色大毛片视频在线观看 | 国产三级中文精品| 六月丁香七月| 亚洲成av人片在线播放无| 国产精品人妻久久久影院| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| 国产精品久久视频播放| 欧美日韩国产亚洲二区| 亚洲av成人av| 18禁动态无遮挡网站| 中文字幕av成人在线电影| 卡戴珊不雅视频在线播放| av线在线观看网站| 蜜臀久久99精品久久宅男| 午夜免费激情av| 精品少妇黑人巨大在线播放 | 老女人水多毛片| 深夜a级毛片| 国产精品嫩草影院av在线观看| 日本熟妇午夜| 国产精品爽爽va在线观看网站| av黄色大香蕉| 日本与韩国留学比较| 亚洲四区av| 高清视频免费观看一区二区 | 麻豆国产97在线/欧美| 亚洲图色成人| 国产极品天堂在线| 亚洲四区av| 看十八女毛片水多多多| 丰满乱子伦码专区| 国产成人a区在线观看| 3wmmmm亚洲av在线观看| 午夜激情欧美在线| 夜夜看夜夜爽夜夜摸| 99国产精品一区二区蜜桃av| 在线播放无遮挡| 九色成人免费人妻av| 又爽又黄a免费视频| 国产精品福利在线免费观看| 国产淫片久久久久久久久| 小说图片视频综合网站| 午夜亚洲福利在线播放| 国产免费又黄又爽又色| 九九久久精品国产亚洲av麻豆| 观看美女的网站| 亚洲自偷自拍三级| 女的被弄到高潮叫床怎么办| 神马国产精品三级电影在线观看| 午夜福利视频1000在线观看| 国产黄片美女视频| 免费观看人在逋| a级一级毛片免费在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 91狼人影院| 亚洲av成人av| 尤物成人国产欧美一区二区三区| av国产免费在线观看| 国产精品国产三级国产专区5o | 日日干狠狠操夜夜爽| 免费av观看视频| 99在线人妻在线中文字幕| 99久久人妻综合| 亚洲成人中文字幕在线播放| 韩国av在线不卡| 亚洲中文字幕一区二区三区有码在线看| 亚洲婷婷狠狠爱综合网| 看免费成人av毛片| 欧美一区二区亚洲| 看黄色毛片网站| 国产亚洲午夜精品一区二区久久 | 美女xxoo啪啪120秒动态图| 久久这里有精品视频免费| 我的女老师完整版在线观看| 建设人人有责人人尽责人人享有的 | 精品欧美国产一区二区三| 99久久九九国产精品国产免费| 美女xxoo啪啪120秒动态图| 18+在线观看网站| 亚洲国产欧洲综合997久久,| 久久精品国产鲁丝片午夜精品| 日韩中字成人| 久久精品久久久久久久性| 久久欧美精品欧美久久欧美| 1024手机看黄色片| 亚洲电影在线观看av| 99国产精品一区二区蜜桃av| 亚洲,欧美,日韩| 最近视频中文字幕2019在线8| 一级毛片aaaaaa免费看小| 亚洲欧美成人综合另类久久久 | 91aial.com中文字幕在线观看| 三级国产精品片| 日日啪夜夜撸| 简卡轻食公司| 91久久精品电影网| 麻豆成人午夜福利视频| 大话2 男鬼变身卡| 久久午夜福利片| h日本视频在线播放| 久久韩国三级中文字幕| 亚洲性久久影院| 欧美性感艳星| 91狼人影院| 欧美精品国产亚洲| 亚洲人成网站高清观看| 免费电影在线观看免费观看| 亚洲久久久久久中文字幕| 色噜噜av男人的天堂激情| 老司机福利观看| 男人舔女人下体高潮全视频| 91精品国产九色| 自拍偷自拍亚洲精品老妇| 热99re8久久精品国产| 亚洲精品乱码久久久久久按摩| 亚洲欧美中文字幕日韩二区| 国产精品一区二区性色av| 久久久色成人| 国内精品美女久久久久久| 日本一二三区视频观看| 日本免费在线观看一区| 18禁裸乳无遮挡免费网站照片| 欧美区成人在线视频| 久久午夜福利片| 熟女电影av网| 国产亚洲av嫩草精品影院| 欧美区成人在线视频| 日本一二三区视频观看| 国产午夜福利久久久久久| 国产亚洲精品久久久com| 亚洲第一区二区三区不卡| 最后的刺客免费高清国语| 熟女电影av网| 日本欧美国产在线视频| 天堂中文最新版在线下载 | 麻豆成人av视频| 中文字幕熟女人妻在线| 中国美白少妇内射xxxbb| 一级毛片久久久久久久久女| 尾随美女入室| 在线观看美女被高潮喷水网站| kizo精华| 日本免费a在线| 日本熟妇午夜| 卡戴珊不雅视频在线播放| 日韩,欧美,国产一区二区三区 | 国产午夜精品论理片| 18禁在线播放成人免费| av在线亚洲专区| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 免费搜索国产男女视频| 成人亚洲欧美一区二区av| 久久午夜福利片| 在线天堂最新版资源| 男人和女人高潮做爰伦理| 久久久亚洲精品成人影院| 日本wwww免费看| 99久国产av精品| 国产 一区精品| 国产单亲对白刺激| 91av网一区二区| 女的被弄到高潮叫床怎么办| 热99在线观看视频| 精品久久久久久成人av| 亚洲av.av天堂| 午夜精品一区二区三区免费看| 麻豆国产97在线/欧美| 可以在线观看毛片的网站| 国产精品爽爽va在线观看网站| 色5月婷婷丁香| 丝袜美腿在线中文| 精品不卡国产一区二区三区| 免费观看的影片在线观看| 1000部很黄的大片| 少妇熟女欧美另类| 哪个播放器可以免费观看大片| 欧美极品一区二区三区四区| 草草在线视频免费看| 国产大屁股一区二区在线视频| 国产精品人妻久久久影院| 看片在线看免费视频| 舔av片在线| 精品少妇黑人巨大在线播放 | 国产成人freesex在线| 久久韩国三级中文字幕| 日本免费在线观看一区| 乱系列少妇在线播放| 国产精品福利在线免费观看| 亚洲人成网站在线播| 久久久久久伊人网av| 国产日韩欧美在线精品| 亚洲欧洲日产国产| 特级一级黄色大片| 午夜福利在线观看吧| 国产一级毛片七仙女欲春2| 免费黄色在线免费观看| 91av网一区二区| 天堂影院成人在线观看| 在线免费观看的www视频| 精品久久久久久久人妻蜜臀av| 少妇裸体淫交视频免费看高清| 国产精品1区2区在线观看.| 国产精品熟女久久久久浪|