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

    Electrochemical catalytic oxidationof p-aminophenol at p-phenol modifiedgraphite paste electrode by cyclicvoltammetry and PM6 semi-empiricalmolecular orbital theory

    2016-08-07 11:53:57ZHUYongchunJIANGQingxiaoSUGuitianCHENQingyang
    關(guān)鍵詞:化學(xué)修飾石墨粉氨基

    ZHU Yongchun, JIANG Qingxiao, SU Guitian, CHEN Qingyang

    (1.Institute of Catalysis for Energy and Environment, Shenyang Normal University, Shenyang 110034, China;2. College of Chemistry and Chemical Engineering, Shenyang Normal University, Shenyang 110034, China)

    ?

    Electrochemical catalytic oxidationof p-aminophenol at p-phenol modifiedgraphite paste electrode by cyclicvoltammetry and PM6 semi-empiricalmolecular orbital theory

    ZHU Yongchun1,2, JIANG Qingxiao2, SU Guitian2, CHEN Qingyang2

    (1.Institute of Catalysis for Energy and Environment, Shenyang Normal University, Shenyang 110034, China;2. College of Chemistry and Chemical Engineering, Shenyang Normal University, Shenyang 110034, China)

    The graphite powder surface is chemically modified by p-aminophenol with a diazo-reaction, and used to compose a carbon paste electrode. At the modified electrode, p-aminophenol shows an electrochemical catalytic oxidation with an oxidation potential negatively shifted for about 73 mV and an increase of the oxidation peak current for about 60% compared with that at a carbon paste electrode without the modification. In order to deeply understand the catalytic oxidation from molecular view, PM6 of a semi-empirical quantum chemical method in MOPAC2012 software was used in the study of electrochemical catalytic activity of the system based on the molecular cluster models of a piece of graphene with 38 carbon atoms as the surface of graphite powder, p-aminophenol and its diazo-reaction product in vacuum conditions. The calculation results show that all the molecular cluster models are thermodynamic stable with negative Gibbs free energies of the reactions. The analysis of frontier orbital gives out an electron transfer induced electrochemical catalysis. The electrochemical redox mechanism of p-aminophenol at modified and unmodified carbon paste electrodes were qualitatively explained. The exchangeable electron of p-aminophenol at p-aminophenol modified carbon paste electrode is about 64.8% larger than that at unmodified carbon paste electrode, which is accord well with the increase of oxidation peak current, and the anodic and cathodic potentials have a linear relation with the energies of HOMO and the LUMO of the molecular clusters of p-aminophenol.

    semi-empirical molecular orbital theory; p-aminophenol ; chemical modified graphite paste electrode; electrochemical catalysis; graphite powder

    0 Introduction

    Carbon paste electrode (CPE) being composed of graphite powder as solid phase and organic binder as liquid phase[1-3]has been used for many years in the electroanalysis and electrochemical studies. In order to increase its functionality, CPEs were modified by different methods such as adsorptions, chemically modifying or intercalating molecules with specific function groups into carbon material surfaces[4-8], as well as the encapsulations of modifiers or solvents into the organic binders[9-12]. However the theoretical or quantum chemical studies on the CPEs have very few reports in literature due to the difficulty in modeling the graphite powder surface[13-14]and its chemical modifications. Nowadays, more and more people pay more attention to the successful applications of quantum chemistry in studies of electrochemical systems, so it also motivate us to apply quantum chemical methods in studies on the CPE system.

    In the present paper,graphite powder surfaces were chemically modified with p-phenol by means of diazo-reaction of p-aminophenol, and used them to make a chemically modified CPE. The modified CPE shows an electrochemical catalytic behavior for the electrochemical oxidation of p-aminophenol with negatively peak potential shift and 60% increase of oxidation peak current compared with the unmodified CPE. In order to obtain an insight into the electron transfer mechanics from quantum chemical view, the graphite powder surface was modeled as a piece of graphene constructed with 38 carbons, and whole working electrode was modeled as molecular cluster composed of the piece of graphene, modifiers and the target molecules. The systems were calculated with semi-empirical quantum chemical method PM6 from MOPAC2012 software[22-23], and the the calculated results were analyzed in thermodynamic and frontier orbital theoretical concepts. Some interesting results are reported here.

    2 Experimental and computational details

    2.1 Materials and instruments

    All electrochemical experiments were carried out on a CHI620b electrochemical instrument (CHI Co., Shanghai, China) with three-electrode system, including a home-made platinum wire electrode as the counter electrode, a p-aminophenol modified carbon paste electrode as the working electrode and a KCl saturated calomel electrode (Model 232, Shanghai, China) as the reference electrode, and all potential reported here was respect to this reference electrode.

    Polyamide resin (Tianjin, China), epoxy resin (Wuxi, China) methyl silicone oil (Shenyang, China) and Graphite powder (Spectrographic pure) were all used for the preparation of basic carbon electrode. P-aminophenol (Analytical pure) was prepared into 0.05 M, and hydrogen chloride (Analytical pure) was prepared into 1.0 M, and sodium nitrite (Analytical pure) was prepared into 0.1 M. All of them were used for the modification of graphite powder. All solutions were prepared with ultra-pure water (18.2 M Ω, USA) for all the experiments. All experiments were carried out at room temperature.

    The basic carbon electrode was prepared as previous work[24]. Graphite powder, epoxy resin and polyamide resin was mixed into paste, filled into a clean glass tube about 6~7 cm in length and 5.0 mm in diameter. A piece of copper wire was inserted into the paste from other end of the tube as electronic lead. The electrode was solidified in air for about 72 hours. After the solidification, the top of the electrode was polished on sand paper, and tip out some paste from the top of the electrode, and left a cavity about 1 mm depth to hold carbon paste.

    The graphite powder surface was modified by diazo-reaction of p-aminophenol with sodium nitrite in acidic solution at lower temperature. 20.0 ml of 0.05 M p-aminophenol solution was mixed with 3 ml of 1.0 M HCl, 10 ml of 0.1 M sodium nitrite in 100 ml beaker. The beaker was put into another 200 mL beaker surround with smashed ice to keep the temperature inside the 100 mL beaker about zero. The mixture solution was stirred with a glassy bar, and 1.0 g of graphite powder was added into the mixture in the beaker while the stirring. The beaker was stored at about 0~4 ℃ for one day waiting for the reaction finished and cooled. The mixture solution was heated and stirred with glassy bar at temperature below 100 ℃ to evaporate the solvent until the graphite powder in the mixture solution become dry powder again. The resulted graphite powder was called p-phenol modified graphite powder. The modified graphite powder was mixed with methyl silicone oil in weight ratio of 2∶1 on a clean glass plate into a paste, the paste was filled into the basic carbon electrode cavity, polished on a piece of glassy paper into mirror surface as the working electrode. If the paste was prepared with unmodified graphite powder and methyl silicone oil, the prepared CPE was called the unmodified carbon paste electrode.

    2.2 Molecular cluster models and computational details

    All quantum chemical calculation in this work was performed with PM6 semi-empirical quantum chemical method in MOPAC2012 software[22-23]on personal computer. All molecules and molecular clusters were designed by PCMOD 7.0[25]and further optimized with MOPAC2012 software before the energy and molecular orbital calculations. The heat and entropy of formation and frontier orbitals were calculated by the semi-empirical calculations in vacuum condition. There were not considered about the influence of solvent effects in the calculation, partly due to the calculation complexity and the methyl silicone oil of organic phase mixed with the graphite powder and greatly reduced the influences of water molecules on the molecule and molecular clusters. All of the obtained data were analyzed and evaluated for the stability and the electron transfer properties of each molecular cluster with thermodynamic and quantum chemical concepts.

    Graphite powderis commonly considered as a stratification of multiple graphene sheets. The surface model of graphite powder in quantum chemical studies can be described as pieces of graphene sheet[26-27]. In our case, a piece of graphene sheet was constructed with 38 carbon atoms in match-join form of twelve six-rings, which was described as C38 in figure 1,C38. A molecule of p-aminophenol was also given in Fig.1 and shorted as pAP. The graphene sheet interact with p-minophenol was built up and optimized in configuration and described as the molecular cluster of C38+pAp in Fig.1. In the case of graphene surface modified with one to three p-phenol molecules through diazo-reaction from p-minophenol, the first p-phenol group was stood on one side of the graphene sheet, and it was continually modified in the same side, then we obtained molecular clusters, and described as C38-p, C38-p2and C38-p3, respectively. These molecular clusters may interact with one p-minophenol, and give out molecular clusters as C38-p+pAP, C38-p2+pAP, and C38-p3+pAP at modified CPE in Fig.1, respectively. After electrochemical oxidation of p-aminophenol, a semi-quinone was produced as an intermediate from the hydroxyl group, it may be tautomerized into semiquinone from amino group, which interacts with C38 modified and forms molecular clusters as C38+pAQ, C38+pNQ at unmodified CPE, and C38-p3+PAQ and C38-p3+pAQ at modified CPE, respectively.

    C38:graphene with 38 carbon atoms; pAP:p-aminophenol; C38+pAP graphene C38 with p-aminophenol; C38-p+pAP:C38-p,graphene C38 with one p-phenol group and one of p-aminophenol; C38-p2+pAP:C38 graphene with two p-phenol groups and one p-aminophenol; C38-p3+pAP:C38 graphen with three p-phenol groups and one p-aminophenol; C38-p3+pNQ: C38 graphene with three p-phenol groups and one p-hydroxylquinoneimine tautomerized from p-aminosemiquinone.

    Fig.1 Molecule or molecular clusters

    3 Results and discussions

    3.1 Electrocatalytic behavior of p-aminophenol at the modified carbon paste electrodes

    The cyclic voltammogram of p-aminophenol on the modified carbon paste electrode shows one oxidation peak located at 0.575 V(58.81 μA), and two reduction peaks at 0.356 V(41.671 μA), and 0.141 V (28.24 μA), respectively as shown in Fig.2A curve 1, while at the unmodified carbon paste electrode, there are also one oxidation peak at 0.648 V(36.71 μA), and two reduction peaks at 0.352 V(13.88 μA), and 0.151 V(4.96 μA) as shown in Fig.2 A curve 2. The modified electrode in electrolyte solution without p-aminophenol has no any peak at all as shown in Fig.2A curve 3.

    (a)—Cyclic voltammograms; (b)— oxidation peak current with square root of scan rate

    Fig.2 p-amino phenol at different carbon paste electrodes

    Compared the curve 1 and curve 2 in Fig.2(a), it is an obvious electrochemical catalytic oxidation of p-aminophenol at p-phenol modified carbon paste electrode with oxidation peak potential negatively shifted about 0.073 V and an increase in oxidation peak current for about 60%. The oxidation peak current was plotted against squared root of scan rate as shown in Fig.2(b), a linear line was obtained with a regression equation of,

    Thisrealationshipindicatestheoxidationisadiffusioncontrolledprocess.Thatmeanstheelectrochemicaloxidationrateisfasterthanthediffusionrate.Inordertounderstandtheelectrocatalyticbehavior,quantumchemicalcomputationalmethodswasperformed,andsomesignificantresultsweregivenasthefollowings.

    3.2 Quantum chemical calculation results and analysis

    3.2.1 Thermodynamic stability of molecule and molecular clusters

    Heat of formation of molecule or molecular clusters can be easily obtained from the PM6 calculation in MOPAC2012, and used to evaluate the stability of molecule or molecular clusters. According to thermodynamic principle[28], the Gibbs free energy change(ΔGr) of a reaction can be calculated with enthalpy change(ΔHr) and entropy change(ΔSr) of the reaction, and described as,

    The enthalpy change and entropy change of a reaction can be calculated fromheat and entropy of formations of products and reactants as,

    Where n and m are the numbers of products and reactants, respectively included in the reactions.

    The heats of formation of all the molecules and molecular clusters used in this work were calculated by PM6 semi-empirical quantum chemical method. The reaction heat of each cluster formation was also calculated according to equation (2~4) and listed in Table 1. All reaction heats of cluster formations are exothermal ones with negative values. If the influence of reaction entropy is too small compared with the reaction heat to be omitted, then the reaction Gibbs free energy is a negative too, so the reactions must be a spontaneous process. The heats of formation of clusters of C38 graphene with p-phenol groups are much higher, but after it interacts with another p-aminophenol with hydrogen bonds, and forms a new molecular cluster, the heat of formation is decreased rapidly. The new molecular cluster formation stabilizes the system, and makes the reaction occur spontaneously. The other molecular clusters follow the similar way with the similar tendency.

    Tab.1 The calculated thermodynamic parameters of molecule or molecular cluster

    All terms were the same as described in Fig.1.

    3.2.3 Frontier orbital energy analysis

    According to quantum chemical theory[29], there are n number of molecular orbitals linearly combined of m number of valence electron orbitals of the atoms in the molecule or molecular cluster as the following,

    The Hamiltonian,Hi, of the system can be described as,

    Wherehjis the kinetic energy part of the system;Vjis the potential energy part of the system. According to Sch?dinger equation,

    Whereφiis the molecular orbital,Eiis the energy of the molecular orbital.so by performing quantum chemical calculation, we can obtained the energy of system, energy and function of each molecular orbital. According to the Kohn-Sham Theory, the frontier molecular orbital (FMO) composed of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) among the orbitals are very important in description of the chemical reactive tendency of a molecule with approximations may be useful for evaluation of electrochemical activity. Some descriptors were theoretical studied and reported such as the frontier orbital gap (ΔE), chemical potential (μ)[30], molecular hardness (η)[31], molecular electron affinity ()[32-33], and charge exchange (ΔN) fromAtoBmolecules[34-36]. All of these terms can be calculated from the frontier orbital energy of the HOMO (εHOMO) and LUMO (εLUMO). All these descriptors are defined as,

    Here ΔNwas originally used to describe the electron transfer quantity from molecule A to molecule B, but in our case, it is used to describe the electron transfer from molecule A to another electrode or output electric circuit withμB, andηB=0, we called molecule B as environment, the new quantity may be called exchangeable electron. All energies of the frontier orbital and all the corresponding descriptors of the molecules and molecular clusters used in this work were calculated by PM6 in MOPAC2012 software, and listed in Table 2. In order to compare the electron transfer capacity of graphene C38 and modified C38 and other molecular clusters, we used the difference of the exchangeable electrons between any molecular cluster (ΔNMC) and graphene C38 (ΔNC38) to describes relative amount of exchangeable electrons included in the molecules or molecular clusters, we call it as difference of exchangeable electron (ΔΔN, as another descriptor) and may be defined as,

    All terms were the same as decribed in Fig.1. The terms with”+” indicates the molecular cluster was calculated with a positive charge.

    The difference of the exchangeable electrons values (ΔN) can be positive or negative. The positive value means the exchangeable electron of a molecular cluster is more than that of graphene, and may be likely to transfer the electron from the molecular cluster to graphene, and the molecular cluster may likely to be oxidized. In the negative case, the molecular cluster has less exchangeable electron than that of graphene, electron may be transferred from environment to the molecular cluster, and the molecular cluster through the graphene is likely to be reduced. So the difference of exchangeable electron may be used to evaluate the electron transfer possibility of a molecular cluster in the CPE system combined with the experimental results. From the calculation results in the Table 2, we can seen that all values from pAP to C38-p3+pAP are negative, that means the modifications and interaction of graphene C38 with the tested molecule have a tendency to reduce the exchangeable electron of the molecular cluster, and result in the reduction processes of the molecular clusters. This conclusion is not the case of the experiments result, where the oxidation process does really occur. In order to found the reasons of the unacceptable theoretical results, we further analyze the components of the frontier orbital of each molecule and molecular cluster in the next section.

    Tab.2 The calculated descriptors of molecule and molecular clusters

    3.2.4 Component analysis of frontier orbital

    The analysis of a molecular orbital components can help us to understand deeply the molecular activities in a chemical reactions, as well as electron transfer property. In our case, we can localize which group of the molecule taking part in the electron transfer process and which molecular orbital in the molecular cluster is mainly responsible for electrochemical oxidation reaction. All the frontier orbital including the LOMO and HOMO as well as the next HOMO were analyzed until we found out the group of our interested one, and all results were plotted in Fig.3.

    All of the LUMO and the first HOMOs of the molecule and molecular clusters are composed of carbon pz atomic orbital from C38 graphene except C38+pAP cluster, which is the composed of p orbitals of atoms of N, O an C in pAp molecule. The all carbon atoms in the graphene, C38, are sp2hybridized, and left pz atomic orbital with one p electron for the composite ofπmolecular orbital. These results indicate the frontier orbitals of graphene related clusters areπorπ*molecular orbitals with higher energies. The modifier of p-phenol group or interacting molecule of p-aminophenol only have small influences on the frontier orbital but not change the orbital components and the properties. The p-phenol group and p-phenol molecular orbials only insert into the groups of the C38 molecular orbitals with energy less than the frontier orbitals of C38. The next HOMO molecular orbitals of clusters including p-aminophenol molecular orbitals have smaller energy gap from the first HOMOs, so if the molecular cluster loss one electron and left one position for the electron transfer of p-aminophenol, then the situation will be changed. This is the electron transfer induced electrochemical catalysis. The out electric circuit draw one electron of graphene from the HOMO orbital composed of pz atomic orbital of carbon in graphene, and then electron in the next HOMO molecular orbital mainly composed of atomic orbitals of p-aminophenol becomes the first HOMO of the molecular cluster, and the next time to draw out one electron from the molecular cluster must be from the molecular orbital mainly composed of atomic orbitals in p-aminophenol. So we further perform the calculations of the systems with one positive charge. The all calculation results are listed in the last part of Table 2. In these cases, all the ΔEterms are much smaller than the corresponding systems with zero charge, which means the electron transfer from HOMO to LUMO is much easier than the corresponding zero charged molecular cluster. All ΔNterms of molecular clusters related to the oxidation of p-aminophenol are positive.

    The plot of energies and components of frontier orbital for the one positively charged molecular clusters was shown in Fig.4. All of the LUMO orbital of the clusters are composed of carbon pz atomic orbital except that of C38+pAP, of which the first and second LUMOs are composed of the p orbitals of C, O and N atoms in pAP molecule, and the HOMO is composed of pz atomic orbital of C atom in C38 due to the first HOMO at zero charge state of the molecular cluster. All of other HOMO of the rest molecular clusters are composed of p orbital of O, N, and C in pAP molecule, so the nest step of electron transfer just from the molecule of pAP molecule, which results in the oxidation of this molecule, and responsible for the oxidation process.

    From the comparison of the ΔNvalues of C38-p3+pAP with that of C38+pAP, It is very interested to found the value of ΔNis increased for about 64.8%, which is accord well with the oxidation peak current in CV curves in Fig.2 with 60.2% increased in oxidation peak current of p-amino-phenol at the p-phenol modified carbon paste electrode compared with that at the unmodified carbon paste electrode.

    3.2.5 Relationship between redox potential and energies of frontier orbitals

    The frontier orbitals of HOMO and LUMO are important in chemical reactions, especially in electrochemical reactions of receiving or losing electron from them.So as early as 1949, Maccoll showed the good relation between reduction potential of a conjugated system and the frontier orbital energy[37]. The oxidation potential,EOX, is proportional to the energy of HOMO,εHOMO, and the reduction potentialEREDis proportional to the energy of LUMO,εLUMO, and all together the expression can be written as[38],

    WhereEOX/REDis the experimental potentials of oxidation and reductions, anda,bare constants. In our system, there are two molecular clusters involuing in the electrochemical redox reactions, one is the C38+pAP for the unmodified CPE and the other is the C38-p3+pAP for the modified CPE. In the C38+pAP molecular cluster shows one pair of well defined redox peaks related to the first electron transfer process, so the frontier orbitals of C38+pAP+were used. In the C38-p3-pAP cluster shows one pair of redox peaks for the first electron transfer process. The other reduction peak belongs to the second electron transfer process. So the quantum chemical calculation was perform again on the molecular cluster of C38-p3+pAP with two positive charge, so the calculation result were also liested in table 2 in term of C38-p3+pAP2+. The energy of LUMO (εLUMO=-6.388 eV) is responsible for the first reduction peak in CV curve. The plot of redox peak potentials against the energies of frontier orbitals is a linear (as shown in Fig.5) with a regression equation of,

    Fig.5 The relationship between redox potential and energy of frontier orbital of p-aminophenol

    The good linear relationship of experimental peak potential and calculated frontier orbital energies further indicates that the quantum chemical calculations are accord well with the experimental results, and offers theoretical evidence for the electrocatalytic behavior of p-aminophenol at the p-aminophenol modified CPE. The second oxidation peak potential can be predicted from the energy of HOMO orbital in C38-p3+pAP2+(-11.814 eV) according to equation (11) as 0.824 V.

    3.2.6 Qualitative explanation of the redox mechanism of p-aminophenol

    In cyclic voltammogram in Fig.2, event the potential window was already up to 1.0 V, there is no second oxidation peak appear, however we obtained the second reduction peak with a relatively larger reduction peak current compared with the second one. The reason for this behavior may be that the first oxidation peak already includes the second oxidation step.

    From organic chemistry view we known, the oxidation product of p-aminophenol is p-amino semiquinone (pAQ), which may tautomerize into p-hydroxylquinoneimine (pNQ), and corresponds to the second reduction peak. From the thermodynamic calculation results in Tab. 1 we can see that the reaction free energy of pNQ is -33.26 kJ/mol, -3 449.8 kJ/mol with C38 and C38+p3, respectively, while the reaction free energy for pAQ is -0.947 8 kJ/mol and -3 158 kJ/mol, respectively. That means pNQ in the molecular clusters are more stable than pAQ. In the tautomerization process, pNQ is preferable to pAQ. In molecular orbital calculation results in Table 2, the energy of first LUMO of pAQ in the molecular clusters is so the -1.099 eV and -3.679 eV with C38 and C38+p3, respectively, while the ones of pNQ is -0.997 eV and -3.461 eV, respectively. That means in the two cases the pAQ is easier to be reduced than pNQ. So the first reduction peak belongs to pAQ, and second reduction peak belongs to pNQ. Because in the two cases the pNQ is more stable than pAQ, so after pAQ produced, it tautomerizes to pNQ mainly, that is responsible for the larger reduction peak current of second one than the first one. The electrochemical redox mechanism may be described as,

    In this way, we can explain qualitatively the all the redox peaks in the experimental CVs from thermodynamic points.

    4 Conclusions

    In the present paper, graphite powder surface was chemically modified with p-amino-phenol by diazo reaction, and shows electrochemical catalytic behavior for p-amino-phenol. Semi-empirical quantum chemical method of PM6 in MOPAC2012 software was used in quantum chemical calculations based on the molecular cluster models including graphene with 38 carbons, p-phenol groups and p-aminophenol. Some conclusions were summarized as the followings. The designed molecular cluster models are stable from the thermodynamic view. The catalytic behavior comes from electron transfer induced catalytic electrochemical oxidation by analysis of frontier orbital energy and composited components. The difference of exchangeable electron can be used to describe the electron transfer possibility which is accord well with the increase of oxidation peak current. The quantum chemical calculations combined with the molecular descriptors will help us to deeply understand the electrochemical catalytic behavior of the systems. The combination of redox potentials shows a linear relation with the calculated energies of frontier orbital, which can be used to predict the unknown redox potential of a molecule, and help us to understand the electrochemical mechanism of the molecule.

    Acknowledgements

    The authors would like to acknowledge the financial supports of the Chinese National Science Foundation (2004DFA02700), Liaoning Education Minister(2004-c022)and national key, and National Key Laboratory on Electroanalytical Chemistry (2006-06).

    [ 1 ]CORDERO-RANDO M M, HIDALGO-HIDALGO de C J L, BLANCO E,et al. The sonogel-carbon electrode as a Sol-Gel graphite-based electrode[J]. Anal Chem, 2002,74(10):2423-2427.

    [ 2 ]BELLIDO-MILLA D, CUBILLANA-AGUILERA L M, KAOUTIT M E, et al. Recent advances in graphite powder-based electrodes[J]. Anal Bioanal Chem, 2013,405(11):3525-3539.

    [ 3 ]HAO J, ZHU Yongchun, LI G. Electrochemical liquid membrane micro-extraction of bromide on ethyl benzoate coated graphite-epoxy composite electrode surface by cyclic voltammetry[J]. Microchim Acta,2009, 167(3):267-272.

    [ 4 ]VYTRAS K, ?倲VANCARA I, METELKA R. Carbon paste electrodes in electroanalytical chemistry[J]. J Serb Chem Soc, 2009,74(10):1021-1033.

    [ 5 ]SISWANA M, OZOEMENA K I, NYOKONG T. Electrocatalytic behaviour of carbon paste electrode modified with iron (Ⅱ) phthalocyanine (FePc) nanoparticles towards the detection of amitrole[J].Talanta, 2006, 69(5):1136-1142.

    [ 6 ]ZHAI Z Q, WU J, SUN W, et al. Direct electrochemistry of hemoglobin and its electrocatalysis based on a carbon nanotube paste electrode[J]. J Chinese Chem Soc, 2009(2):56:561-567.

    [ 7 ]PUNBUSAYAKULYAKUL N. carbon nanotubes architectures in electroanalysis[J]. Procedia Engin, 2012,32(1):683-689.

    [ 8 ]LI M, WU Z. A review of intercalation composite phase change material: Preparation, structure and properties[J]. Renew Sustain Energ Rev, 2012,6(4):2094-2101.

    [ 9 ]APETREI C, APETREI I M, DE SAJA J A, et al. Carbon paste electrodes made from different carbonaceous materials: application in the study of antioxidants[J]. Sens, 2011,11(2):1328-1344.

    [10]RAHBAR N, PARHAM H. Carbon baste electrode modified with Cuo-nanoparticles as a probe for square wave voltammetric determination of atrazine[J]. Jundishapur J Nat Pharm Prod, 2013,8(3):118-124.

    [11]AGLANL R F, MOHAMED G G, MOHAMED H A. Chemically modified carbon paste electrode for determination of cesium ion by potentiometric method[J]. Am J Anal Chem, 2012,3(8):576-586.

    [12]SIM C M, HONG Y J, KIM M H, et al. Voltammetry behavior of modified carbon paste electrode with cytochrome C and Mn2O3nanoparticles for hydrogen peroxide sensing[J]. Int J Electrochem Sci, 2012,7(12):12098-12109.

    [13]STRELKO V V, DUDARENKO V V, LAVRIENKO-OMETSINSKAYA E D. Quantum-chemical interpretation of some electrochemical properties of activated carbons[J]. Ukr Khim Zh, 1987,53(6):820-824.

    [14]LAVRINENKO-OMETSINSKAYA E D, REPYAKH I V, KARTEL N T, et al. States of a heteroatom in graphite-like structures of nitrogen-containing carbons[J].Teor kspo Khim, 1982,18(2):224-227.

    [15]JAMAL M, SARAC A S, MAGNER E. Conductive copolymer-modified carbon fibre microelectrodes: electrode characterisation and electrochemical detection of p-aminophenol[J]. Sens Actuat B, 2004,97(1):59-66.

    [16]SUN W, JIAO K, ZHANG S, et al. Electrochemical detection for horseradish peroxidase-based enzyme immunoassay using p-aminophenol as substrate and its application in detection of plant virus[J]. Anal Chim Acta, 2001,434(1):43-50.

    [17]HALSALL H B, HEINEMAN W R, JENKINS S H. Capillary immunoassay with electrochemical detection[J]. Clin Chem, 1988,34(9):1701-1702.

    [18]WIJAYAWARDHANA C A, HALSALL H B, HEINEMAN W R. Micro volume rotating disk electrode (RDE) amperometric detection for a bead-based immunoassay[J]. Anal Chim Acta, 1999,399(1):3-11.

    [19]GIL EP, Tang H, Halsall B,et al. Competitive heterogeneous enzyme immunoassay for theophylline by flow-injection analysis with electrochemical detection of p-aminophenol[J]. Clin Chem, 1990,36(4):662-665.

    [20]ENSAFI A A, HEYDARI-BAFROOEI E, REZAEI B. Simultaneous detection of hydroxylamine and phenol using p-aminophenol-modified carbonnanotubepaste electrode[J]. Chin J Catal, 2013,34(9):1768-1775.

    [21]AHMED T M F, SANKARAPAPAVINASAM S. Poly(para-aminophenol): a new soluble, electroactive conducting polymer[J]. J Electroanal Chem, 1992, 338(1/2):347-352.

    [22]KWON S J, YANG H, KWARK K J J. An electrochemical immunosensor using p-aminophenol redox cycling by NADH on a self-assembled monolayer and ferrocene-modified Au electrodes[J]. Analyst , 2008,133(11):1599-1604.

    [23]AGUILAR Z P, FRITSCH I. Immobilized enzyme-linked DNA-hybridization assay with electrochemical detection for cryptosporidium parvum hsp70 mRNA[J].Anal Chem, 2003,75:3890-3897.

    [24]SREWART J J P. Application of the PM6 method to modeling the solid state [J]. J Mol Model, 2008,14(6):499-535.

    [25]CLARK T. A handbook of computational chemistry: a practical guide to chemical structure and energy calculations[M]. New York:Wiley, John & Sons, Incorporated, 1985.

    [26]ZHU Yongchun, GUAN J, CAO L, et al. Electrochemical solid phase nano-extraction and determination of trace iodide in iodised table salt on silver sulfate-modified carbon paste electrode by differential pulse voltammetry[J]. Talanta, 2010,80(3):1234-1238.

    [27]STEWART J J P. Optimization of parameters for semiempirical methods V:modification of NDDO approximations and applications to 70 elements[J]. J Mol Model, 2007,13(12):1173-1213.

    [28]CALLEN H B. Thermodynamics and introduction to thermostatistics[M]. New York:Wiley, 1985.

    [29]MAGNASCO V. Methods of molecular quantum mechanics: an introduction to electronic molecular structure[M]. United Kingdom:John Wiley and sons, Ltd, Chichester, 2009.

    [30]CHOI Y, MEBANE DS, WANG J, et al. Continuum and quantum-chemical modeling of oxygen reduction on the cathode in a solid oxide fuel cell[J]. Top Catal, 2007,46(3):386-401.

    [31]CALLEN H B. Thermodynamics and introduction to thermostatistics[M]. New York: Wiley, 1985.

    [32]PARR RG, DONNELLY RA, LEVY M, et al. Electronegativity: the density functional viewpoint[J]. J Chem Phys, 1978,68(8):3801-3807.

    [33]PEARSON R G. Hard and soft acids and bases[J]. J Am Chem Soc, 1963,85(26):3533-3539.

    [34]PARR R G, SZENTPALY LV, LIU S. Electrophilicity index[J]. J Am Chem Soc, 1999,121(9):1922-1924.

    [35]MENDEZ F, GARCIA-GARIBAY MA. A hard soft acid base and DFT analysis of singlet triplet gaps and the addition of singlet carbenes to alkenes[J]. J Org Chem, 1999,64(19):7061-7066.

    [36]SURENDRA B N. Computational studies foroxidation reduction reactions of cinnoline-4(1H)-one, in aqueous phase by density functional theory[J]. Brit J Appl Sci Techn, 2014,4(3):465-467.

    [37]MACCOLL A. Reduction potentials of conjugated systems[J]. Nature, 1949,163:178-179.

    1673-5862(2016)03-0343-11

    對(duì)氨基酚在修飾碳糊電極上電氧化的PM6法研究

    朱永春1,2, 蔣清笑2, 蘇桂田2, 陳慶陽2

    (1. 沈陽師范大學(xué) 能源與環(huán)境催化研究所, 沈陽 110034;2. 沈陽師范大學(xué) 化學(xué)化工學(xué)院, 沈陽 110034)

    以重氮化反應(yīng)將對(duì)氨基化學(xué)修飾在石墨粉表面上并構(gòu)建成碳糊電極。在此修飾電極上對(duì)-氨基酚表現(xiàn)處電化學(xué)催化氧化行為,與裸電極相比,氧化峰電位負(fù)移了73 mV,氧化峰電流升高了60%。為了從分子水平上深入了解這種催化氧化行為,以38個(gè)碳構(gòu)成的石墨烯片段-對(duì)氨基-對(duì)氨基酚組成的分子簇模型,MOPAC2012軟件包內(nèi)的PM6半經(jīng)驗(yàn)分子軌道法研究了真空條件下體系的電化學(xué)催化活性。計(jì)算結(jié)果表明,所有建立的分子簇模型都是熱力學(xué)穩(wěn)定的,都具有負(fù)的反應(yīng)吉布斯自由能。前線分子軌道的分析給出了電化學(xué)催化誘導(dǎo)下的電子轉(zhuǎn)移過程,并定性地解釋了對(duì)氨基酚在修飾和未修飾電極上的電化學(xué)氧化機(jī)理。對(duì)氨基酚在修飾電極上的可交換電子量比在未修飾電極上的大了64.8%。此結(jié)果恰好與試驗(yàn)結(jié)果中的氧化峰電流的升高結(jié)果一致;試驗(yàn)結(jié)果中的氧化還原峰電位對(duì)計(jì)算結(jié)果中的相應(yīng)最高占有軌道和最低空軌道能呈現(xiàn)很好的線性關(guān)系。

    半經(jīng)驗(yàn)分子軌道理論; 對(duì)氨基酚; 化學(xué)修飾碳糊電極; 電化學(xué)催化; 石墨粉

    date: 2015-05-17.

    O646 Document code: A

    10.3969/ j.issn.1673-5862.2016.03.019

    Supported: Project supported by National Natural Science Foundation of China (20807563).

    Biography: ZHU Yongchun(1955-), male, was born in Huaide city of Jilin province, professor of Shenyang Normal University, doctor.

    猜你喜歡
    化學(xué)修飾石墨粉氨基
    石墨粉對(duì)磷石膏基自流平砂漿性能的影響
    高溫氣冷堆核燃料元件用石墨粉粒度表征方法的研究
    豬δ冠狀病毒氨基肽酶N的研究進(jìn)展
    可再生石墨粉染色圓領(lǐng)T恤
    超疏水功能材料用于食品非法添加的超靈敏SERS檢測(cè)
    聚乙二醇修飾重組人生長(zhǎng)激素的研究
    基于納米材料化學(xué)修飾電極在農(nóng)藥檢測(cè)中的應(yīng)用進(jìn)展
    相轉(zhuǎn)移催化合成2-氨基異煙酸
    鍍鎳碳纖維/鍍鎳石墨粉填充PC/ABS復(fù)合材料導(dǎo)電性能研究
    中國塑料(2015年2期)2015-10-14 05:34:16
    多糖化學(xué)修飾方法的研究進(jìn)展
    日本欧美视频一区| 丰满迷人的少妇在线观看| 99re6热这里在线精品视频| 亚洲欧美精品综合一区二区三区| 麻豆国产av国片精品| 精品少妇黑人巨大在线播放| 精品国产国语对白av| 亚洲五月婷婷丁香| 在线天堂中文资源库| av天堂在线播放| 久久久国产一区二区| 亚洲熟妇熟女久久| 精品国产超薄肉色丝袜足j| 国产精品98久久久久久宅男小说| 日韩有码中文字幕| 老熟女久久久| 中文欧美无线码| www日本在线高清视频| 日韩大码丰满熟妇| 亚洲三区欧美一区| 一级毛片电影观看| 国产av精品麻豆| 久久狼人影院| 丰满少妇做爰视频| 亚洲欧洲日产国产| 成人黄色视频免费在线看| 大陆偷拍与自拍| avwww免费| 国产精品一区二区免费欧美| 久久性视频一级片| 精品国产一区二区三区久久久樱花| 国产真人三级小视频在线观看| 亚洲色图av天堂| 精品国产乱子伦一区二区三区| 一本—道久久a久久精品蜜桃钙片| 香蕉国产在线看| 五月开心婷婷网| 男女午夜视频在线观看| 久久精品亚洲精品国产色婷小说| 色综合欧美亚洲国产小说| 国产黄频视频在线观看| 午夜福利乱码中文字幕| 成年人午夜在线观看视频| 高清在线国产一区| aaaaa片日本免费| 黄片大片在线免费观看| 99精国产麻豆久久婷婷| 日韩视频在线欧美| 久热爱精品视频在线9| 热re99久久国产66热| 性少妇av在线| 国产伦人伦偷精品视频| 另类亚洲欧美激情| 97在线人人人人妻| 国产真人三级小视频在线观看| 久久婷婷成人综合色麻豆| 国产精品99久久99久久久不卡| 熟女少妇亚洲综合色aaa.| 国产国语露脸激情在线看| 91字幕亚洲| 亚洲国产av新网站| 国产精品久久久人人做人人爽| 亚洲avbb在线观看| 国产精品国产av在线观看| 一进一出好大好爽视频| 两性午夜刺激爽爽歪歪视频在线观看 | 黄片播放在线免费| 亚洲国产欧美一区二区综合| 99re在线观看精品视频| 国产在线一区二区三区精| 超碰97精品在线观看| 国产精品国产高清国产av | 欧美黄色淫秽网站| 一级,二级,三级黄色视频| 国产精品av久久久久免费| 亚洲成人国产一区在线观看| 亚洲 国产 在线| 热99国产精品久久久久久7| av一本久久久久| 亚洲国产av影院在线观看| 亚洲一区二区三区欧美精品| 欧美激情极品国产一区二区三区| 精品视频人人做人人爽| 侵犯人妻中文字幕一二三四区| 脱女人内裤的视频| 99热国产这里只有精品6| 一级片免费观看大全| 国产免费现黄频在线看| 国产在线观看jvid| 老熟女久久久| 在线十欧美十亚洲十日本专区| 丰满迷人的少妇在线观看| 一级毛片电影观看| 久久天躁狠狠躁夜夜2o2o| 亚洲视频免费观看视频| 午夜精品国产一区二区电影| 精品国产国语对白av| 久久午夜综合久久蜜桃| 少妇的丰满在线观看| 精品福利永久在线观看| 十八禁网站免费在线| 巨乳人妻的诱惑在线观看| 午夜视频精品福利| 久久久久网色| 国产人伦9x9x在线观看| tocl精华| 国产精品久久久久久人妻精品电影 | 中文字幕人妻丝袜一区二区| 超碰97精品在线观看| 色老头精品视频在线观看| 一区二区三区乱码不卡18| 欧美日韩亚洲国产一区二区在线观看 | 老司机亚洲免费影院| 久久午夜综合久久蜜桃| 精品欧美一区二区三区在线| 伦理电影免费视频| 黄片大片在线免费观看| 这个男人来自地球电影免费观看| 中亚洲国语对白在线视频| 久久av网站| avwww免费| 日韩精品免费视频一区二区三区| 一区二区三区激情视频| 黑人操中国人逼视频| 午夜久久久在线观看| 国产成人欧美在线观看 | 黄频高清免费视频| 中文字幕人妻丝袜一区二区| 丝瓜视频免费看黄片| 在线观看www视频免费| 国产无遮挡羞羞视频在线观看| 欧美日韩国产mv在线观看视频| 高潮久久久久久久久久久不卡| 亚洲欧洲日产国产| 国产成人一区二区三区免费视频网站| 久9热在线精品视频| 超色免费av| 18禁美女被吸乳视频| 1024香蕉在线观看| 欧美午夜高清在线| 国产精品国产av在线观看| 国产高清激情床上av| h视频一区二区三区| 精品一区二区三区视频在线观看免费 | 多毛熟女@视频| 成人国语在线视频| 婷婷丁香在线五月| 久久中文看片网| 老司机影院毛片| 国产一区二区激情短视频| 中文字幕人妻丝袜一区二区| 国产高清激情床上av| 欧美变态另类bdsm刘玥| 91成年电影在线观看| 精品亚洲成国产av| 动漫黄色视频在线观看| 亚洲专区国产一区二区| 黄色 视频免费看| 天天躁夜夜躁狠狠躁躁| 亚洲欧美精品综合一区二区三区| av不卡在线播放| 黄色成人免费大全| 日韩中文字幕欧美一区二区| 一区二区三区激情视频| 变态另类成人亚洲欧美熟女 | 男女高潮啪啪啪动态图| 91麻豆av在线| 最新在线观看一区二区三区| 久久久久精品国产欧美久久久| 无限看片的www在线观看| 美国免费a级毛片| 亚洲中文日韩欧美视频| 一本综合久久免费| a级毛片黄视频| 国产精品98久久久久久宅男小说| www日本在线高清视频| 成人永久免费在线观看视频 | 美女高潮喷水抽搐中文字幕| 午夜福利一区二区在线看| 亚洲视频免费观看视频| 亚洲精品国产精品久久久不卡| 少妇裸体淫交视频免费看高清 | 亚洲熟妇熟女久久| 午夜福利在线免费观看网站| 19禁男女啪啪无遮挡网站| 99精国产麻豆久久婷婷| 精品福利永久在线观看| 一区二区三区精品91| 国产精品秋霞免费鲁丝片| 夜夜夜夜夜久久久久| 老司机福利观看| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品电影小说| 侵犯人妻中文字幕一二三四区| 一级黄色大片毛片| a级毛片黄视频| 亚洲精品成人av观看孕妇| 午夜两性在线视频| 亚洲精品国产区一区二| 热re99久久精品国产66热6| 久久久国产欧美日韩av| h视频一区二区三区| 五月天丁香电影| 天堂中文最新版在线下载| 欧美日韩中文字幕国产精品一区二区三区 | 一级黄色大片毛片| 午夜视频精品福利| 国产日韩欧美亚洲二区| 亚洲国产毛片av蜜桃av| 午夜福利,免费看| 国产亚洲欧美在线一区二区| 国产亚洲精品久久久久5区| 亚洲国产欧美一区二区综合| 99riav亚洲国产免费| 桃红色精品国产亚洲av| 久久精品亚洲av国产电影网| 丝袜人妻中文字幕| 精品熟女少妇八av免费久了| 视频区图区小说| 可以免费在线观看a视频的电影网站| 色综合欧美亚洲国产小说| 国产精品久久电影中文字幕 | 极品教师在线免费播放| 一级片免费观看大全| 成年人黄色毛片网站| 三级毛片av免费| 在线av久久热| 美女扒开内裤让男人捅视频| 国产aⅴ精品一区二区三区波| 精品国产乱子伦一区二区三区| 亚洲 欧美一区二区三区| 一区二区三区国产精品乱码| 亚洲国产欧美网| 欧美 亚洲 国产 日韩一| 国产精品免费视频内射| 在线亚洲精品国产二区图片欧美| 欧美成人午夜精品| 天堂俺去俺来也www色官网| 黑人巨大精品欧美一区二区蜜桃| 亚洲第一av免费看| tube8黄色片| 国产淫语在线视频| kizo精华| 五月开心婷婷网| 国产伦人伦偷精品视频| 最近最新中文字幕大全免费视频| 交换朋友夫妻互换小说| 一级,二级,三级黄色视频| 国产精品美女特级片免费视频播放器 | 淫妇啪啪啪对白视频| 国产成+人综合+亚洲专区| 欧美 亚洲 国产 日韩一| netflix在线观看网站| 精品国产超薄肉色丝袜足j| 国产日韩欧美在线精品| 性高湖久久久久久久久免费观看| 俄罗斯特黄特色一大片| 久久久久久久精品吃奶| 国产成人av教育| 女人精品久久久久毛片| 国产日韩一区二区三区精品不卡| 精品国产国语对白av| 妹子高潮喷水视频| 黑人巨大精品欧美一区二区mp4| 国产成人欧美| 国产xxxxx性猛交| 欧美黑人欧美精品刺激| 国产精品.久久久| 最近最新中文字幕大全免费视频| 99精品久久久久人妻精品| 亚洲熟女精品中文字幕| 国产精品秋霞免费鲁丝片| 三级毛片av免费| 欧美成人免费av一区二区三区 | 免费在线观看完整版高清| 国产伦理片在线播放av一区| 在线观看免费视频网站a站| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品 国内视频| 久久性视频一级片| 后天国语完整版免费观看| 亚洲人成电影免费在线| 老鸭窝网址在线观看| 国产精品香港三级国产av潘金莲| 精品一区二区三区视频在线观看免费 | 国产亚洲欧美精品永久| 啦啦啦在线免费观看视频4| 精品熟女少妇八av免费久了| 国产麻豆69| 精品国产乱子伦一区二区三区| 午夜两性在线视频| 国产精品国产av在线观看| 国产精品 国内视频| 可以免费在线观看a视频的电影网站| 午夜日韩欧美国产| www.自偷自拍.com| 午夜福利乱码中文字幕| 操美女的视频在线观看| 91麻豆av在线| 久久99热这里只频精品6学生| 最近最新中文字幕大全电影3 | 日韩中文字幕视频在线看片| 天天操日日干夜夜撸| 高清黄色对白视频在线免费看| 在线观看免费午夜福利视频| 一级片'在线观看视频| 国产免费视频播放在线视频| 亚洲午夜精品一区,二区,三区| 不卡一级毛片| 丁香欧美五月| 99re6热这里在线精品视频| 91av网站免费观看| 美女主播在线视频| 黄色丝袜av网址大全| 香蕉久久夜色| 免费在线观看影片大全网站| 国产淫语在线视频| 国产欧美亚洲国产| 精品福利永久在线观看| 桃花免费在线播放| xxxhd国产人妻xxx| 777久久人妻少妇嫩草av网站| 热99国产精品久久久久久7| 妹子高潮喷水视频| 蜜桃在线观看..| av网站免费在线观看视频| 国产av又大| 亚洲欧洲日产国产| 精品久久久精品久久久| 精品国产亚洲在线| 精品亚洲成a人片在线观看| 人妻一区二区av| 亚洲七黄色美女视频| 免费观看a级毛片全部| 青草久久国产| 日本wwww免费看| 欧美日韩视频精品一区| 国产男女内射视频| 香蕉久久夜色| 欧美日韩黄片免| 在线观看一区二区三区激情| 高清视频免费观看一区二区| 日本a在线网址| 亚洲天堂av无毛| 中文字幕最新亚洲高清| 亚洲va日本ⅴa欧美va伊人久久| 欧美精品亚洲一区二区| 在线观看一区二区三区激情| 天天操日日干夜夜撸| 久久久久精品国产欧美久久久| 亚洲成人免费电影在线观看| 国产精品二区激情视频| 免费观看a级毛片全部| 热re99久久国产66热| 99国产极品粉嫩在线观看| 免费观看av网站的网址| 亚洲色图 男人天堂 中文字幕| 一级毛片电影观看| 美国免费a级毛片| 狠狠精品人妻久久久久久综合| 欧美日本中文国产一区发布| 黄频高清免费视频| 99热网站在线观看| 一级毛片精品| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 超色免费av| 亚洲欧美日韩高清在线视频 | 亚洲全国av大片| 动漫黄色视频在线观看| 国产伦理片在线播放av一区| 欧美乱妇无乱码| 黑人猛操日本美女一级片| 亚洲欧美精品综合一区二区三区| 国产日韩欧美亚洲二区| 亚洲精品国产精品久久久不卡| 中文字幕人妻丝袜制服| 美女主播在线视频| 免费在线观看日本一区| 超碰97精品在线观看| av有码第一页| 高清黄色对白视频在线免费看| 交换朋友夫妻互换小说| 在线观看免费午夜福利视频| 日韩有码中文字幕| 国产精品 欧美亚洲| av电影中文网址| 男人舔女人的私密视频| 国产欧美亚洲国产| 黄片大片在线免费观看| 欧美中文综合在线视频| 欧美 亚洲 国产 日韩一| 亚洲国产看品久久| 这个男人来自地球电影免费观看| www.精华液| 日韩中文字幕欧美一区二区| 99re6热这里在线精品视频| 日本wwww免费看| 亚洲自偷自拍图片 自拍| 日本欧美视频一区| 欧美乱妇无乱码| 成人18禁高潮啪啪吃奶动态图| 每晚都被弄得嗷嗷叫到高潮| 黑人巨大精品欧美一区二区蜜桃| 日本a在线网址| 久久久久精品国产欧美久久久| 一区二区三区激情视频| 深夜精品福利| 欧美日韩视频精品一区| 久久 成人 亚洲| 精品亚洲成国产av| 久久久久久久精品吃奶| 欧美久久黑人一区二区| 亚洲 国产 在线| 国产免费福利视频在线观看| 久久久精品94久久精品| 999精品在线视频| 国产精品自产拍在线观看55亚洲 | 国产男女超爽视频在线观看| 国产单亲对白刺激| 国产精品1区2区在线观看. | 国产精品一区二区精品视频观看| 91麻豆精品激情在线观看国产 | 亚洲精品中文字幕在线视频| videosex国产| 久久精品国产99精品国产亚洲性色 | 国产精品免费一区二区三区在线 | 无遮挡黄片免费观看| 黄频高清免费视频| 男女高潮啪啪啪动态图| 蜜桃国产av成人99| 免费高清在线观看日韩| 成人三级做爰电影| 亚洲欧美精品综合一区二区三区| 另类亚洲欧美激情| 精品久久蜜臀av无| 精品福利观看| 欧美黄色片欧美黄色片| 别揉我奶头~嗯~啊~动态视频| 精品国产乱子伦一区二区三区| 99久久人妻综合| 纵有疾风起免费观看全集完整版| 成人影院久久| 伊人久久大香线蕉亚洲五| 这个男人来自地球电影免费观看| 亚洲精品av麻豆狂野| 99国产精品免费福利视频| 久久香蕉激情| 免费在线观看视频国产中文字幕亚洲| www.999成人在线观看| 菩萨蛮人人尽说江南好唐韦庄| 美女国产高潮福利片在线看| 少妇被粗大的猛进出69影院| 男人操女人黄网站| 黄色成人免费大全| 国产在线免费精品| 亚洲第一欧美日韩一区二区三区 | 中文字幕制服av| 高清黄色对白视频在线免费看| 在线播放国产精品三级| 欧美日韩亚洲高清精品| 欧美激情高清一区二区三区| 久久国产精品人妻蜜桃| 精品亚洲成a人片在线观看| 亚洲精品自拍成人| 国产免费视频播放在线视频| 国产淫语在线视频| 免费在线观看黄色视频的| 国产精品一区二区精品视频观看| 国产欧美亚洲国产| 亚洲综合色网址| 成人av一区二区三区在线看| 国产成人精品无人区| 国产成+人综合+亚洲专区| 久久久水蜜桃国产精品网| 国产成人啪精品午夜网站| 精品亚洲成国产av| 欧美亚洲 丝袜 人妻 在线| 免费在线观看视频国产中文字幕亚洲| 99re6热这里在线精品视频| 亚洲avbb在线观看| 91成人精品电影| 国产成人一区二区三区免费视频网站| 妹子高潮喷水视频| 色精品久久人妻99蜜桃| 日韩免费高清中文字幕av| 亚洲成人国产一区在线观看| 狂野欧美激情性xxxx| 在线亚洲精品国产二区图片欧美| 一级片免费观看大全| 欧美乱码精品一区二区三区| 欧美午夜高清在线| 日韩中文字幕视频在线看片| 久久久精品区二区三区| cao死你这个sao货| 免费在线观看视频国产中文字幕亚洲| 成在线人永久免费视频| 午夜91福利影院| 国产一区二区三区在线臀色熟女 | 精品视频人人做人人爽| 男人操女人黄网站| 少妇 在线观看| 啪啪无遮挡十八禁网站| 亚洲自偷自拍图片 自拍| 欧美激情极品国产一区二区三区| 精品第一国产精品| 搡老岳熟女国产| 久9热在线精品视频| 欧美精品亚洲一区二区| 在线看a的网站| 午夜91福利影院| 露出奶头的视频| 国产99久久九九免费精品| 80岁老熟妇乱子伦牲交| 欧美黄色淫秽网站| 成人av一区二区三区在线看| 亚洲av欧美aⅴ国产| 一区二区三区国产精品乱码| 亚洲精品美女久久av网站| 自拍欧美九色日韩亚洲蝌蚪91| 日本五十路高清| 午夜老司机福利片| 飞空精品影院首页| 18在线观看网站| 免费在线观看视频国产中文字幕亚洲| 91精品三级在线观看| 久久人妻av系列| 久久毛片免费看一区二区三区| 两个人看的免费小视频| 在线观看一区二区三区激情| 国产精品亚洲av一区麻豆| 欧美 日韩 精品 国产| 一边摸一边抽搐一进一小说 | netflix在线观看网站| 日韩三级视频一区二区三区| 正在播放国产对白刺激| 国产真人三级小视频在线观看| 老司机午夜十八禁免费视频| 国产精品久久久久久精品古装| 一边摸一边抽搐一进一小说 | 日日爽夜夜爽网站| 1024香蕉在线观看| 国产成人啪精品午夜网站| 日韩成人在线观看一区二区三区| 免费在线观看黄色视频的| 午夜免费成人在线视频| 精品久久蜜臀av无| 在线观看人妻少妇| 精品国产一区二区久久| 中文亚洲av片在线观看爽 | 天堂中文最新版在线下载| 丁香六月欧美| 99re6热这里在线精品视频| 亚洲欧洲精品一区二区精品久久久| av在线播放免费不卡| 久久精品国产99精品国产亚洲性色 | 嫁个100分男人电影在线观看| www日本在线高清视频| 狠狠婷婷综合久久久久久88av| 久久精品国产综合久久久| 久久精品人人爽人人爽视色| 少妇猛男粗大的猛烈进出视频| 99国产综合亚洲精品| 99国产精品一区二区三区| cao死你这个sao货| 日韩中文字幕欧美一区二区| 久久香蕉激情| 在线观看人妻少妇| 熟女少妇亚洲综合色aaa.| 国产不卡av网站在线观看| 国产精品一区二区免费欧美| 91精品三级在线观看| 啦啦啦中文免费视频观看日本| 成人免费观看视频高清| 99国产精品99久久久久| 99国产综合亚洲精品| 夜夜骑夜夜射夜夜干| 午夜福利免费观看在线| 久久国产亚洲av麻豆专区| av有码第一页| 午夜日韩欧美国产| 大陆偷拍与自拍| 超碰成人久久| 少妇裸体淫交视频免费看高清 | 久久久水蜜桃国产精品网| 精品第一国产精品| 日韩精品免费视频一区二区三区| 99热网站在线观看| 纯流量卡能插随身wifi吗| 久久久国产成人免费| 黄色 视频免费看| 99精品久久久久人妻精品| 国产高清视频在线播放一区| 久久香蕉激情| 亚洲精品中文字幕在线视频| 久9热在线精品视频| 999久久久国产精品视频| 久久久精品国产亚洲av高清涩受| 亚洲成人免费av在线播放| 桃花免费在线播放| 久久国产精品男人的天堂亚洲| svipshipincom国产片| 蜜桃在线观看..| 看免费av毛片| 性高湖久久久久久久久免费观看| 9191精品国产免费久久| 日韩欧美免费精品| 国产精品国产高清国产av | 狠狠婷婷综合久久久久久88av| 国产精品久久久久成人av| 亚洲一区二区三区欧美精品| 超碰97精品在线观看| 国产一区二区在线观看av| 涩涩av久久男人的天堂| 大片电影免费在线观看免费|