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

    HSH-C10: A new quasi-2D carbon allotrope with a honeycomb-star-honeycomb lattice

    2022-09-16 05:24:56QinGoLifuZhngCiynZhengShuliLeiShujunLiZhenpengHu
    Chinese Chemical Letters 2022年8期

    Qin Go, Lifu Zhng, Ciyn Zheng, Shuli Lei, Shujun Li, Zhenpeng Hu,?

    a School of Physics, Nankai University, Tianjin 300071, China

    b Hubei Key Laboratory of Low Dimensional Optoelectronic Materials and Devices, Hubei University of Arts and Science, Xiangyang 441053, China

    c Institute of Mathematics, Free University of Berlin, Arnimallee 6, Berlin D-14195, Germany

    ABSTRACT Two-dimensional (2D) materials with honeycomb, kagome or star lattice have been intensively studied because electrons in such lattices could give rise to exotic quantum effects.In order to improve structural diversity of 2D materials to achieve unique properties, here we propose a new quasi-2D honeycombstar-honeycomb (HSH) lattice based on first-principles calculations.A carbon allotrope named HSH-C10 is designed with the HSH lattice, and its mechanical properties have been intensively investigated through total energy, phonon dispersion, ab initio molecular dynamic simulations, as well as elastic constants calculations.Besides the classical covalent bonds, there is an interesting charge-shift bond in this material from the chemical bonding analysis.Additionally, through the analysis of electronic structure, HSH-C10 is predicted to be a semiconductor with a direct band gap of 2.89 eV, which could combine the desirable characteristics of honeycomb and star lattice.Importantly, by modulating coupling strength, a flat band near the Fermi level can be obtained in compounds HSH-C6Si4 and HSH-C6Ge4, which have potential applications in superconductivity.Insight into such mixed lattice would inspire new materials with properties we have yet to imagine.

    Keywords:Carbon allotrope Two-dimensional lattice Charge-shift bond Flat band First-principles

    The unique two-dimensional (2D) lattices, like honeycomb lattice [1], kagome lattice [2] and star lattice [3], have been intensively studied because electrons in such lattices could give rise to exotic quantum effects.The honeycomb lattice could induce Dirac bands, while the kagome lattice and star lattice could induce both Dirac bands and flat bands [1–5].They are interesting subjects to study various properties in quantum fields, such as quantum spin liquid, frustrated spin, and fractional quantum Hall effect [6].Especially, the kinetic energy of electrons is quenched in the flat bands, which is helpful to realize strongly correlated electronic states, such as superconductivity and Wigner crystal [7,8].These lattices not only provide ideal research platforms for understanding novel physics, but also have potential applications in electronics and spintronics [9–14].

    Many studies have focused on the practical application of materials with unique lattices based on the special properties[15–22] mentioned above.With honeycomb lattice, graphene is one of the most widely investigated 2D carbon allotropes [23].It has also actively encouraged and promoted the research on other 2D topological materials with honeycomb lattice, like silicene, germanene, and stanine [24–26].Besides structures with honeycomb lattice, various structures with kagome lattice have also been proposed [27–29].Materials with kagome lattice usually exhibit interesting properties.For example, it is reported that the designed phosphorus carbide P2C3with kagome lattice contains double kagome bands, which is named from the property that the flat band in one kagome band accidently crossing the Dirac point in the other kagome band [30].And the experimental evidence of flat band in a kagome lattice has been observed in twisted bilayer silicene, making it a possible candidate for realizing exotic quantum Hall states [31].Recently, a star carbon monolayer was also proposed with a flat band near the fermi level [32,33].By the tight-binding models on the decorated star lattices, Mizoguchi and his co-authors proposed a class of carbon-based materials with exotic electronic structures such as Dirac cones and flat bands [34].

    Fig.1.The top and side views of (a) bilayer honeycomb lattice, (b) star lattice, and(c) honeycomb-star-honeycomb lattice.

    To design novel 2D materials, a new route which mixes different lattices in structures has been devised recently.A novel honeycomb–kagome (HK) lattice was proposed, which could be regarded as a mixed lattice composed of honeycomb and kagome lattices [35–39].Comparing the single lattice, the mixed lattice significantly increases the structural diversity of 2D materials and could provide more opportunities to achieve unique properties [40].With experimental synthesis methods making rapid progress, practical methods are being explored to establish well-defined manufacturing techniques to realize the mixed lattices.For example, a spin-superstructure with honeycomb-triangular lattice was realized experimentally, which the magnetic layers in triangular and honeycomb lattices exhibit sequential magnetic orderings [41].It is noted that the novel 2D materials with a mixed lattice have quickly attracted significant attention and would be a potential research focus.

    In this study, a newquasi-2D honeycomb-star-honeycomb(HSH) lattice was proposed.The idea of the lattice was sprung up by the peculiar features of honeycomb lattice and star lattice.Material with such a mixed lattice would combine the desirable characteristics of both honeycomb and star lattices.The carbon-based materials would induce various characteristic electronic structures,and thus become our primary consideration to probe the physical and chemical properties.The dynamical, thermal and mechanical stability as well as the electronic properties are evaluated.It is expected that more new materials with exotic properties could be designed based on this mixed lattice.

    In this study, all the studies were conducted based on density functional theory (DFT), as implemented in the Vienna ab initio simulation package (VASP) [42,43].The Perdew?Burke?Ernzerhof(PBE) [44] exchange-correlation functional and PAW [45] pseudopotential were used.The electronic wave function was expanded in plane waves to a cut-off energy of 520 eV.For the unit cell and thep(2 × 2 × 1) supercell, Monkhorst-Pack meshes with 6 × 6 × 1 and 3 × 3 × 1 in Brillouin zone sampling were used, respectively.To avoid interactions between the repeated images, the vacuum thickness larger than 16was used.All the structures were relaxed with no symmetry constraint, and the convergence criteria for total energy and residual force on each ion were set to 10–5eV and 0.01 eV/, rescpectively.The phonon calculations were performed using PHONOPY package [46] with the forces calculated from density functional perturbation theory (DFPT) [47] as implemented in VASP.Because the convergence of phonon frequencies needs higher energy accuracy than normal electronic structure calculations, the total energy was converged to less than 10?8eV in all the phonon calculations.

    The star lattice can be considered as an expanded honeycomb lattice which replacing the sites of the honeycomb lattice by triangles.The HSH lattice is a mixed lattice with a star lattice sandwiched by two honeycomb lattices, as shown in Fig.1.Such lattice would significantly increase the structural diversity of 2D materials and provide more opportunities to achieve unique properties.Carbon is not only containing various allotropes, but also has the richest variety of compounds among all elements.As a preliminary exploration of the application in real materials, we designed a carbon allotrope based on the HSH lattice.The honeycomb lattice containing four atoms and the star lattice containing six atoms, the new carbon allotrope is named HSH-C10.

    Fig.2.(a) The top and (b) side views of HSH-C10, (c) the phonon dispersion along high symmetry directions of Brillouin zone, and (d–f) final structures with ab initio molecular dynamic simulations under the temperature of 500, 800 and 1000 K.

    The top and side views of HSH-C10are illustrated in Figs.2a and b.From a linear scan on the lattice constant (Fig.S1 in Supporting information), the carbon allotrope has a hexagonal cell with a lattice constant of 6.60and two different carbons with three different bond lengths.Six carbon atoms are marked from C1 to C6, as shown in Fig.2a and b.The lengths of three bonds marked as C2-C6, C1-C2, and C1-C3 are 1.30, 1.51 and 1.67, respectively.The carbon atom in honeycomb lattice is named CHoneycomb(C1,C3), while the one in star lattice is named CStar(C2, C4, C5, C6).For comparison, the energy of all the systems is defined as DFT computed ground-state energy per atom.The energies of graphene[48], two carbon allotropes (named as C-phase 1 and 2 here) [49],E-graphene [50], as well as graphdiyne [51] are also calculated,and details of all the structures are shown in Table 1.Graphene is the most energetically stable allotrope, and the energy of ?9.22 eV/atom is the same as previous DFT simulations [52].The energy of HSH-C10is ?7.80 eV/atom, which is higher than graphene and graphdiyne (?8.46 eV/atom), but lower than C-phase 1(?7.45 eV/atom), and a little higher than C-phase 2 (?7.93 eV/atom) and E-graphene (?7.93 eV/atom).The bond lengths of 1.51 and 1.67in HSH-C10are a little higher than the two kinds of bond lengths for C-phase 2 (1.50 and 1.59), while the bond length of 1.30is much smaller than that in C-phase 2.Furthermore, the recently reported T-carbon has an energy of?7.92 eV/atom [53], which is closed to the energy of C-phase 2.Obviously, the relative energy of HSH-C10structure is lower than C-phase 1 and much close to C-phase 2 and T-carbon.Since Cphase 1 is predicted to be stable at a low temperature [49] and T-carbon is already realized in the experiments [54], it is believed there would be a high possibility to realize the HSH-C10structure experimentally based on the energy data.

    From a structural analysis, the relatively high energy of this material comes from the formed hexahedral structure.The four bonds of CHoneycombatom point to the same side, leading to the inverted bond CHoneycomb-CHoneycomb(C1-C3) embedded inside a cage of six wing CHoneycomb-CStarbonds, thus generating relatively high strain energy.Interestingly, the hexahedral structure is the same as [1.1.1]propellane, which is a hydrocarbon with the molecular formula C5H6and could be synthesized in large quantities[55,56].It is reported that this molecule has remarkable stability despite the supposed strain energy from the inverted bridgehead bond [57–61].For comparison, the structure of [1.1.1]propellane is also shown in Fig.S2 (Supporting information).As listed in Table 1, the PBE calculated inverted bond length for [1.1.1]propellane is 1.58, consistent with the experimental results of 1.59[62], and little shorter than the CHoneycomb-CHoneycombbond length of 1.67in HSH-C10.As reported [61], the inverted bond is a charge-shift bond, which is regarded as a third kind of chemicalbond in addition to the ionic and covalent bonds.The charge-shift bond shows large covalent-ionic resonance interaction energy and depleted charge densities, but the bond itself could well be strong[57–61].From the structural connection, it is promising that HSHC10could be synthesized by the polymerization of [1.1.1]propellane.

    Table 1 The number of carbon atoms in unit cell (N), lattice constants (a and b), total energy per atom (Etot) and bond length (dC-C) for HSH-C10.For comparison, the corresponding parameters for graphene, two carbon allotropes (named as C-phase1 and 2), E-graphene, graphdiyne, and [1.1.1]propellane molecule are also provided here.

    To further confirm the dynamical stability of HSH-C10structure, the calculations of phonon dispersion are also performed.Ap(2 × 2 × 1) supercell is used to calculate the vibrational modes of HSH-C10along the high symmetry k-points throughout the first Brillouin zone.As shown in Fig.2c, the phonon dispersion of HSHC10has little imaginary frequency less than 4 cm-1close toГpoint.As the little imaginary frequency in long wavelength flexural acoustic (ZA) mode close toГpoint is the feature ofquasi-2D materials [63–66], the phonon dispersion indicates HSH-C10is dynamically stable.Besides, the highest frequency of HSH-C10reaches up to 58 THz (1935 cm-1), a little higher than that of graphene (1600 cm-1) [67], confirming the robust C–C bonding in 2D HSH-C10sheet.There also exists a large direct bandgap of 20 THz (~667 cm-1) at theГpoint in the phonon band structures,which would be useful to identify the structure of HSH-C10in experiments.

    Ab initiomolecular dynamic simulations in a Nose-Hoover thermostat are also performed to check the thermal stability.Ap(3 × 3 × 1) supercell containing 90 carbon atoms is built to reduce the constraint of periodic boundary condition and explore the possible structure reconstruction.Molecular dynamic simulations under the temperature of 500, 800, 1000 K for 10 ps are performed, and the final structures are depicted in Figs.2d–f.With no bonds breaking in the whole simulation, the structure is stable under the temperature of 500 K.Under the temperature of 800 K,the charge-shift bond (C1-C3) would fluctuate relatively greatly.When the temperature is heated to 1000 K, obvious structure reconstruction is found.The molecular dynamic simulations suggest that HSH-C10is separated by high energy barriers from other local minima on the potential energy surface.

    The mechanical properties of HSH-C10are also studied in this work.The elastic constants of materials based on theoretical calculations are not only important for the potential applications,but also could be used to investigate the mechanical stabilities.The finite distortion method is used to calculate the linear elastic constants [68].The independent elastic constantsCijcan be then expressed asCij=?F/?εiεj, and the strain energy density functionFis expressed asF=F0+F2ε2/2 +O(ε3),whereF0andF2ε2/2 are the static energy and the vibrational contribution, respectively.There are three independent elastic constants for the hexagonal structure, which areC11(C22=C11),C12andC44.Based on Born-Huang criteria [69], the mechanical stability for a 2D material needs to satisfyC11>C12andC44>0.By fitting the strain energy to the parabolic function of the axial and biaxial strains, it is obtained thatC22=C11= 394.14 N/m,C12= 384.42 N/m, andC44= 22.86 N/m, which also listed in Table 2.The elastic constants of HSH-C10satisfy the criteria tobe a mechanically stable material.Based on the linear elastic constant, Young’s modulus (Y) can be calculated using the equationY= (C112-C122)/C11.By multiplyingYand the length of nonperiodic direction of the calculation model, the in-plane stiffness(Ys) of HSH-C10is 20 N/m.It is in the same weight with that of germanene (43 N/m) and silicene (62 N/m) computed at the same theoretical level, which have been reported in our previous study[70].From the above discussion, it could be concluded that HSHC10also exhibits good mechanical property.

    Table 2 Independent elastic constants (Cij, kbar) and in-plane stiffness(Ys, N/m) of HSH-C10, germanene and silicene.

    Fig.3.(a) The top and (b) side views for electron localization function isosurface of HSH-C10 with an isovalue of 0.85, and (c) projected crystal orbital Hamilton population (pCOHP) analysis for the bonds of C2-C6, C1-C2, C1-C3.The Fermi level is set at 0 eV.

    To get a clear insight of the structural feature displayed in HSHC10, the chemical bonding is analyzed by the electron localization function (ELF) and projected crystal orbital Hamilton population(pCOHP) [71,72].ELF can be illustrated as a contour picture in real space with values ranging from 0 to 1.The region 0 means a low electron density area, while 1 indicates the perfect localization of electrons.It plays an important role in understanding the special characteristics of the chemical bonding.The typical kind of bonding in HSH-C10can be seen from ELF isosurface as shown in Figs.3a and b.Clearly, the isosurface of ELF at 0.85 mainly comprises three localization areas.The first is located in the middle of the CHoneycomb-CHoneycomb(C2-C6) within the 2D plane, the second is located in the middle of the CHoneycomb-CStar(C1-C2), and the third is located around the CStaratoms (C1, C3).The C2-C6 bond length is 1.30, considerably shorter than the C1-C2 bond with 1.51.Correspondingly, the ELF distribution around the middle of C2-C6 is higher than that of C1-C2, indicating a stronger covalent bonding.Interestingly, there is no evident electron localization between CStar-CStar(C1-C3).The charge is depleted from the center of the hexahedral cage and delocalized onto the bridge carbon atoms,which is a typical feature of charge-shift bond [57–61,73].Comparing with [1.1.1]propellane, the little longer CHoneycomb-CHoneycomb(C1-C3) bond length of 1.67could be attributed to the change in the hybridization of the orbitals of the bridgehead carbon atoms[73].Importantly, the charge-shift bond leads to the trend that more electrons are distributed outside thequasi-2D structure, thus suggesting the potential for applications in molecule adsorption and as a catalyst [74,75].

    Fig.4.The band structures for (a) HSH-C10 and (b) star carbon monolayer, and(c, d) atom-decomposed band structures along high symmetry direction of the Brillouin zone for HSH-C10.The circle sizes indicate the contribution of atoms in honeycomb or star lattice to the band structure.The Fermi level is set at 0 eV.

    To well display pCOHPs, the usual way is used by drawing negative contributions to the right and positive to the left, so bonding interaction goes to the right and antibonding interaction goes to the left.From pCOHPs shown in Fig.3c, it is noticed that the C2-C6 and C1-C2 bonds only fill bonding levels, while the inverted C1-C3 bond exhibits antibonding levels.This is well agreed with the previous observation in the literature that the charge-shift bond exhibits anti-bonding character [57,73].As a comparison, the pCOHP curves of [1.1.1]propellane is also shown in Fig.S3 (Supporting information), which exhibits the bonding features similar to HSHC10.Previous studies have reported that the structure could be well stabilized by partial population of the C1-C3σ-antibonding/bridgeπ-bonding orbitals [57,58].Combining the ELF to analyze, the C2-C6 and C1-C2 bonds are the classical covalent bond, while the inverted C1-C3 bond is the charge-shift bond.With ELF and pCOHP curves, the chemical bonding nature of HSH-C10is well captured.As the charge-shift bonding in previous reports normally happens in unique molecules, HSH-C10may open the door to the study of charge-shift bonding in the field of 2D materials.

    Fig.5.The top and side views of (a) HSH-C6Si4 and (b) HSH-C6Ge4, as well as (c, d)the band structures along high symmetry direction of the Brillouin zone.The Fermi level is set at 0 eV.

    The structures with a mixed lattice may combine the desirable characteristics of both the honeycomb and star lattices.From electronic band structures shown in Fig.4a, the PBE results indicate HSH-C10structure is a semiconductor with a direct band gap of 2.89 eV.It is noted that the PBE band gap for T-carbon is 2.25 eV,while it is 4.16 eV for diamond [53].The honeycomb lattice could induce Dirac bands with the Dirac cone just appear at the Fermi level [4,76], while the star lattice could induce both Dirac bands and flat bands near the Fermi level [3,32], as shown in Fig.4b.Comparing graphene with Dirac cone at the Fermi level, the band gap of HSH-C10is opened up by the coupling between CHoneycomband CStaratoms, whereas the two groups of Dirac-like bands separately located on the two sides of Fermi level.To obtain the correlation between the lattice and electronic structure, the atomdecomposed band structures (Figs.4c and d) along high symmetry direction of the Brillouin zone for HSH-C10are calculated, and the contributions of honeycomb and star lattices thus could be clearly observed.The valence band maximum (VBM) is mainly contributed from carbon atoms in star lattice while the conduction band minimum (CBM) is contributed from atoms in both star and honeycomb lattices.Three groups of Dirac-like bands and two flat bands are numbered from 1 to 8.The Dirac-like bands which numbered as 3, 4 and flat bands numbered as 8 are contributed from the carbon atoms in star lattice, while the Dirac-like bands numbered as 1, 2, 5, 6 and flat band numbered as 7 are contributed from the atoms in honeycomb lattice and star lattice.From the above discussion, the coupling of honeycomb and star lattices leads to the result that HSH-C10loses the metallic character in graphene and carbon allotrope with star lattice, but it retains some characteristics like flat band and Dirac bands.

    The electronic band structure is affected by the strength of coupling between atoms in honeycomb and star lattice, so it could be modulated by substituting atoms in the HSH lattice.Since silicon and germanium atoms have larger covalent radius and higher energy levels than carbon atoms, they are used to replace carbon atoms in honeycomb lattice, and the new structures are named as HSH-C6Si4and HSH-C6Ge4.The optimized crystal structures with top and side views are shown in Figs.5a and b, and the equilibrium lattice constants are estimated to be 7.55 and 7.81.The bond lengths of C-Si and C-Ge are 1.94 and 2.08, respectively.Phonon dispersion indicates that the structures are dynamically stable, as shown in Figs.S4a and b (Supporting information).The band structures would be acquired by including the couplings between the silicon (or germanium) atoms in honeycomb lattice and carbon atoms in star lattice.Different from HSH-C10with a large band gap, the HSH-C6Si4system exhibits semiconductor with a band gap of 0.21 eV, while HSH-C6Ge4is metallic.This is because the coupling between the carbon and silicon (or germanium)atoms is weaker than the CHoneycomb-CStarcoupling in HSH-C10.As illustrated in Figs.5c and d, a flat band occurs near the Fermi energy for HSH-C6Si4and HSH-C6Ge4.Interestingly, the flat band is attached to the upper Dirac-like band for HSH-C6Si4, while it is attached to the lower one for HSH-C6Ge4.The unique electronic structure renders HSH-C6Si4and HSH-C6Ge4potential 2D materials for application in nanoelectronics.The coexistence of flat and steep bands near the Fermi level is regarded as an important condition for the superconductivity of materials [77,78].All the results indicate that HSH lattice could combine the electronic characteristics of both the honeycomb and star lattices and expand the diversity of 2D materials.

    In summary, we designed a new HSH lattice and reported a newquasi-2D carbon allotrope named HSH-C10.Total energy,phonon dispersion, molecular dynamic simulations, as well as elastic constants calculations indicate the structure has good stability.Through chemical bonding analysis by ELF and pCOHP, the C2-C6 and C1-C2 bonds of HSH-C10are classical covalent bonds, while the inverted C1-C3 bond is described as a charge-shift bond, which is strong and shows large covalent-ionic resonance interaction energy.What is more, the charge-shift bond leads to the trend that more electrons are distributed outside the 2D plane of HSH-C10,suggesting a potential for applications in molecule adsorption and as a catalyst.It is expected that the HSH-C10could be synthesized by the polymerization of [1.1.1]propellane, and would be a fascinating 2D carbon allotrope for its unique structure and intriguing reactivity.As the first reported 2D material with charge-shift bond, HSH-C10may open the door for understanding charge-shift bonding in fundamental chemistry and physics, as well as other fields in application.The results of the calculated electronic band structures show that HSH-C10carbon is a semiconductor with a direct band gap of 2.89 eV and could combine the desirable characteristics of honeycomb and star lattice.To modulate the electronic band structure to the potential application in nanoelectronics, by replacing carbon element in honeycomb lattice with silicon and germanium, new structures named HSH-C6Si4and HSH-C6Ge4are also obtained.Different from HSH-C10with a large band gap,the HSH-C6Si4system exhibits semiconductor with a band gap of 0.21 eV, while HSH-C6Ge4is metallic.With different elements in the sub-lattice, combination of the desirable characteristics of both the honeycomb and star lattices could be realized, and the physical and chemical properties would be modulated.All the results indicate that HSH lattice would be used to derive the charge-shift bond family, increase the structural diversity of 2D materials and has potential in probing or engineering chemical and physical properties.This study may provide a creative platform to study the exotic properties of HSH lattice and stimulate more efforts on the research field of mixed lattice.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work is partially supported by the National Natural Science Foundation of China (Nos.21773124, 12134019), by the Fundamental Research Funds for the Central Universities (Nankai University, No.63213042), by the Ph.D.Candidate Research Innovation Fund of Nankai University, Science and Technology Research Project of Hubei Provincial Department of Education (No.D20212603) and Hubei University of Arts and Science (Nos.2020kypytd001, 2020kypytd002, XK2021024).This work is also supported by the Supercomputing Center of Nankai University(NKSC).

    Supplementary materials

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2021.11.027.

    大香蕉久久网| 久久99蜜桃精品久久| 亚洲七黄色美女视频| 身体一侧抽搐| 最近的中文字幕免费完整| 国内少妇人妻偷人精品xxx网站| 亚洲成a人片在线一区二区| 听说在线观看完整版免费高清| 伦理电影大哥的女人| 老女人水多毛片| 亚洲七黄色美女视频| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看| 黄色欧美视频在线观看| 欧美性猛交黑人性爽| 99视频精品全部免费 在线| 啦啦啦啦在线视频资源| 国产精品久久久久久精品电影| 尾随美女入室| 国产亚洲精品久久久com| 国产激情偷乱视频一区二区| 中国国产av一级| 免费观看在线日韩| 色尼玛亚洲综合影院| 激情 狠狠 欧美| 97超碰精品成人国产| 日韩三级伦理在线观看| 亚洲欧美成人精品一区二区| 免费人成视频x8x8入口观看| 男女啪啪激烈高潮av片| 午夜激情欧美在线| 爱豆传媒免费全集在线观看| 看免费成人av毛片| 久久久欧美国产精品| 中文在线观看免费www的网站| 日日干狠狠操夜夜爽| 精品午夜福利在线看| 色哟哟哟哟哟哟| 亚洲在线自拍视频| 亚洲国产精品久久男人天堂| 国产成人aa在线观看| 日韩成人av中文字幕在线观看| 神马国产精品三级电影在线观看| 亚洲国产精品合色在线| 日本一本二区三区精品| 乱系列少妇在线播放| 国产成人福利小说| 亚洲精品自拍成人| 九九在线视频观看精品| 亚洲熟妇中文字幕五十中出| 国产精品一区二区三区四区久久| 国内少妇人妻偷人精品xxx网站| 成年版毛片免费区| 成人高潮视频无遮挡免费网站| 中文字幕免费在线视频6| 国产中年淑女户外野战色| 国产精品精品国产色婷婷| 国产伦精品一区二区三区视频9| 日日啪夜夜撸| 午夜福利在线观看免费完整高清在 | 悠悠久久av| 内射极品少妇av片p| 中出人妻视频一区二区| 春色校园在线视频观看| 国产视频首页在线观看| 亚洲美女视频黄频| 国产老妇女一区| 美女脱内裤让男人舔精品视频 | 国产精品福利在线免费观看| 深夜a级毛片| 久久99蜜桃精品久久| 日本与韩国留学比较| 欧美又色又爽又黄视频| 国产av在哪里看| 国产精品久久久久久亚洲av鲁大| 精品熟女少妇av免费看| 亚洲av二区三区四区| 国产精品电影一区二区三区| 欧美高清性xxxxhd video| 夜夜夜夜夜久久久久| 国产中年淑女户外野战色| 99热这里只有精品一区| 免费在线观看成人毛片| 99在线视频只有这里精品首页| 午夜福利视频1000在线观看| 亚洲精品久久国产高清桃花| 性色avwww在线观看| 欧美成人a在线观看| 狠狠狠狠99中文字幕| 91久久精品国产一区二区三区| 亚洲熟妇中文字幕五十中出| 在线免费观看的www视频| avwww免费| 日本免费一区二区三区高清不卡| 国产伦精品一区二区三区四那| 免费看av在线观看网站| 国产精品一区www在线观看| 熟女人妻精品中文字幕| 熟女人妻精品中文字幕| 国产淫片久久久久久久久| 日韩高清综合在线| av卡一久久| av卡一久久| 老司机影院成人| 三级毛片av免费| 三级毛片av免费| 性欧美人与动物交配| 美女高潮的动态| 亚洲一区高清亚洲精品| 又爽又黄无遮挡网站| 国产精品永久免费网站| 女的被弄到高潮叫床怎么办| 免费看av在线观看网站| 成年免费大片在线观看| 亚洲精品456在线播放app| 日本黄色片子视频| 一卡2卡三卡四卡精品乱码亚洲| 久久精品久久久久久久性| 免费看av在线观看网站| 中文字幕免费在线视频6| 亚洲国产精品国产精品| 亚洲天堂国产精品一区在线| 国产精品电影一区二区三区| 啦啦啦观看免费观看视频高清| 国产高清三级在线| 又爽又黄a免费视频| 午夜久久久久精精品| 久久久精品大字幕| 久久久精品大字幕| 十八禁国产超污无遮挡网站| 十八禁国产超污无遮挡网站| 国产色婷婷99| 别揉我奶头 嗯啊视频| 亚洲人成网站在线播| 中文字幕免费在线视频6| 亚洲久久久久久中文字幕| 校园春色视频在线观看| 国产在线精品亚洲第一网站| 国内精品久久久久精免费| 国产精品不卡视频一区二区| www.av在线官网国产| 神马国产精品三级电影在线观看| 可以在线观看毛片的网站| 亚洲成av人片在线播放无| 成人高潮视频无遮挡免费网站| 亚洲成人av在线免费| 国产黄片美女视频| 黄色日韩在线| 欧美一区二区精品小视频在线| 两个人的视频大全免费| 99国产极品粉嫩在线观看| 男女下面进入的视频免费午夜| 久久精品国产自在天天线| 国内揄拍国产精品人妻在线| 99热网站在线观看| 99久久中文字幕三级久久日本| 美女国产视频在线观看| 亚洲av成人av| 中文精品一卡2卡3卡4更新| 一夜夜www| 久久九九热精品免费| av女优亚洲男人天堂| 久久精品国产亚洲网站| 精品久久久久久成人av| 午夜福利成人在线免费观看| 高清在线视频一区二区三区 | a级毛片免费高清观看在线播放| 毛片一级片免费看久久久久| 国产真实乱freesex| 麻豆精品久久久久久蜜桃| 特大巨黑吊av在线直播| 欧美高清性xxxxhd video| 成人永久免费在线观看视频| 国产欧美日韩精品一区二区| 成人欧美大片| 岛国毛片在线播放| 国产精品久久久久久久电影| 国产乱人视频| 26uuu在线亚洲综合色| 精品一区二区免费观看| 久久久精品94久久精品| 在线观看美女被高潮喷水网站| 一本久久中文字幕| 亚洲人与动物交配视频| 女的被弄到高潮叫床怎么办| 精品人妻一区二区三区麻豆| 国产精品野战在线观看| 日本熟妇午夜| 女同久久另类99精品国产91| 欧美又色又爽又黄视频| 精品国产三级普通话版| 欧美精品国产亚洲| 亚洲国产欧美在线一区| 成人漫画全彩无遮挡| 日韩高清综合在线| 精品无人区乱码1区二区| 乱码一卡2卡4卡精品| 亚洲最大成人中文| 麻豆久久精品国产亚洲av| 深夜a级毛片| 日韩欧美三级三区| 搡老妇女老女人老熟妇| 午夜精品在线福利| 色尼玛亚洲综合影院| 十八禁国产超污无遮挡网站| 久久久久久久久久成人| 国产极品精品免费视频能看的| 国产精品福利在线免费观看| 丰满人妻一区二区三区视频av| 亚洲最大成人av| 好男人视频免费观看在线| 欧美色欧美亚洲另类二区| 亚洲美女搞黄在线观看| 色5月婷婷丁香| 午夜福利在线观看免费完整高清在 | 亚洲成人精品中文字幕电影| 在线观看一区二区三区| 日韩强制内射视频| 春色校园在线视频观看| 男人狂女人下面高潮的视频| 成人一区二区视频在线观看| 国产伦一二天堂av在线观看| 亚洲欧美精品自产自拍| 欧美日韩一区二区视频在线观看视频在线 | 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 亚洲第一电影网av| 成人永久免费在线观看视频| 成人毛片60女人毛片免费| 少妇丰满av| 亚洲最大成人av| 亚洲av免费高清在线观看| 成人鲁丝片一二三区免费| 内地一区二区视频在线| 91狼人影院| 永久网站在线| 国内揄拍国产精品人妻在线| 99在线人妻在线中文字幕| 国产精品久久久久久av不卡| 亚洲成人久久爱视频| 成年av动漫网址| 亚洲无线观看免费| 成人毛片a级毛片在线播放| 热99在线观看视频| 免费观看人在逋| 日本撒尿小便嘘嘘汇集6| 少妇被粗大猛烈的视频| 成人毛片60女人毛片免费| 国产亚洲5aaaaa淫片| 国产精品人妻久久久影院| 免费看美女性在线毛片视频| 免费观看的影片在线观看| 男人狂女人下面高潮的视频| 成人午夜精彩视频在线观看| 亚洲成人av在线免费| www日本黄色视频网| 有码 亚洲区| 亚洲国产精品成人久久小说 | 男的添女的下面高潮视频| 99久久精品热视频| 国产黄片视频在线免费观看| 国产亚洲av嫩草精品影院| 校园春色视频在线观看| 国产伦在线观看视频一区| 天天躁夜夜躁狠狠久久av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 寂寞人妻少妇视频99o| 亚洲综合色惰| 一级毛片久久久久久久久女| 可以在线观看的亚洲视频| 色尼玛亚洲综合影院| 看免费成人av毛片| www日本黄色视频网| 最近中文字幕高清免费大全6| 久久精品国产亚洲网站| 看黄色毛片网站| 亚洲国产高清在线一区二区三| 国产蜜桃级精品一区二区三区| 97在线视频观看| 久久午夜亚洲精品久久| 免费黄网站久久成人精品| 91精品一卡2卡3卡4卡| 国产高潮美女av| 国产高清视频在线观看网站| 色播亚洲综合网| 性欧美人与动物交配| 在线国产一区二区在线| 国产免费一级a男人的天堂| 在线免费观看的www视频| 亚洲婷婷狠狠爱综合网| 国产69精品久久久久777片| 97超碰精品成人国产| 欧美日韩一区二区视频在线观看视频在线 | 国产一级毛片在线| 国产精品电影一区二区三区| 国产女主播在线喷水免费视频网站 | 免费观看人在逋| 一级av片app| 欧美+日韩+精品| 午夜福利高清视频| 国产黄a三级三级三级人| 国产成人影院久久av| 国产在视频线在精品| 天天一区二区日本电影三级| 中文欧美无线码| 欧美日韩乱码在线| 午夜激情福利司机影院| 啦啦啦观看免费观看视频高清| 国产精品一区二区性色av| 波多野结衣高清作品| 国产精品久久视频播放| 99久国产av精品| 久久久久免费精品人妻一区二区| 日韩av不卡免费在线播放| 91久久精品国产一区二区三区| 99热这里只有是精品在线观看| 中文字幕av在线有码专区| 日韩av不卡免费在线播放| 综合色av麻豆| 国产69精品久久久久777片| 国产免费男女视频| 亚洲自拍偷在线| 九色成人免费人妻av| 18禁在线无遮挡免费观看视频| 亚洲一区高清亚洲精品| 精品少妇黑人巨大在线播放 | 人人妻人人澡欧美一区二区| 哪里可以看免费的av片| 久久久久网色| 亚洲在线自拍视频| 日本-黄色视频高清免费观看| 99热只有精品国产| 亚洲国产精品成人久久小说 | 色综合亚洲欧美另类图片| 日韩高清综合在线| 最近最新中文字幕大全电影3| 中文亚洲av片在线观看爽| 国产精品,欧美在线| 性插视频无遮挡在线免费观看| 男女那种视频在线观看| 你懂的网址亚洲精品在线观看 | 美女黄网站色视频| 给我免费播放毛片高清在线观看| 国产又黄又爽又无遮挡在线| 国产极品精品免费视频能看的| 国产一区二区在线av高清观看| 国产美女午夜福利| 国产午夜精品论理片| 欧美bdsm另类| 成人特级av手机在线观看| 美女cb高潮喷水在线观看| 午夜精品国产一区二区电影 | 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av涩爱 | 亚洲国产色片| 午夜视频国产福利| 久久99蜜桃精品久久| av国产免费在线观看| 69人妻影院| 真实男女啪啪啪动态图| 欧美一区二区精品小视频在线| a级毛片a级免费在线| 人妻系列 视频| 色5月婷婷丁香| 成人国产麻豆网| 国内精品宾馆在线| 久久久欧美国产精品| 哪个播放器可以免费观看大片| 成人午夜高清在线视频| 国产又黄又爽又无遮挡在线| 日本-黄色视频高清免费观看| 男女边吃奶边做爰视频| 少妇熟女欧美另类| 免费看光身美女| 国产精品一及| 成人性生交大片免费视频hd| 伦精品一区二区三区| 欧美日韩精品成人综合77777| 免费搜索国产男女视频| 色尼玛亚洲综合影院| 日本-黄色视频高清免费观看| 欧美一区二区国产精品久久精品| 老司机福利观看| 中国美白少妇内射xxxbb| 免费观看a级毛片全部| 久久午夜福利片| 联通29元200g的流量卡| 男女啪啪激烈高潮av片| 精品免费久久久久久久清纯| 亚洲精品乱码久久久v下载方式| 婷婷精品国产亚洲av| 精品久久久久久久久av| 啦啦啦啦在线视频资源| 国产精品1区2区在线观看.| 国产男人的电影天堂91| 午夜激情福利司机影院| 中文资源天堂在线| www日本黄色视频网| 国产老妇伦熟女老妇高清| 欧美人与善性xxx| 丰满的人妻完整版| 精品人妻熟女av久视频| 一个人看视频在线观看www免费| 日韩制服骚丝袜av| 午夜老司机福利剧场| 国产精品人妻久久久久久| 成人国产麻豆网| 可以在线观看的亚洲视频| 少妇熟女aⅴ在线视频| 九九久久精品国产亚洲av麻豆| 久久久国产成人免费| 午夜免费男女啪啪视频观看| 少妇人妻一区二区三区视频| 日韩高清综合在线| 日韩在线高清观看一区二区三区| 久久精品国产亚洲av香蕉五月| 亚洲四区av| 亚洲,欧美,日韩| 深夜精品福利| 少妇熟女欧美另类| 精品少妇黑人巨大在线播放 | 日韩一本色道免费dvd| 91在线精品国自产拍蜜月| 一本精品99久久精品77| 国产乱人视频| 成人漫画全彩无遮挡| 国产精品永久免费网站| 悠悠久久av| 午夜亚洲福利在线播放| 直男gayav资源| 精品人妻偷拍中文字幕| 狠狠狠狠99中文字幕| 丝袜喷水一区| 亚洲成人精品中文字幕电影| 免费搜索国产男女视频| 99热精品在线国产| 男人舔女人下体高潮全视频| 亚洲欧美日韩卡通动漫| 国产成人a区在线观看| 自拍偷自拍亚洲精品老妇| 日韩高清综合在线| 久久这里有精品视频免费| 日韩中字成人| 午夜免费激情av| 国产午夜精品一二区理论片| 97人妻精品一区二区三区麻豆| 国产高清有码在线观看视频| av在线老鸭窝| 哪个播放器可以免费观看大片| 好男人在线观看高清免费视频| 99热这里只有精品一区| 国产精品人妻久久久影院| 国产一区二区三区在线臀色熟女| 九九爱精品视频在线观看| 欧美+亚洲+日韩+国产| 欧美成人a在线观看| 免费一级毛片在线播放高清视频| avwww免费| 免费在线观看成人毛片| av天堂中文字幕网| 免费看美女性在线毛片视频| 日本一本二区三区精品| 国产高潮美女av| 久久人人爽人人爽人人片va| 免费看日本二区| 欧美不卡视频在线免费观看| 日本色播在线视频| 久99久视频精品免费| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 精品久久久噜噜| 日韩精品有码人妻一区| 非洲黑人性xxxx精品又粗又长| 亚洲人成网站在线播放欧美日韩| 激情 狠狠 欧美| 天堂网av新在线| 美女黄网站色视频| 麻豆av噜噜一区二区三区| 亚洲三级黄色毛片| 男的添女的下面高潮视频| 能在线免费看毛片的网站| 久久精品国产亚洲av涩爱 | 特大巨黑吊av在线直播| 欧美+亚洲+日韩+国产| 日韩欧美精品免费久久| 亚洲精华国产精华液的使用体验 | 久久久久久久久久久丰满| 亚洲自拍偷在线| 国产黄色视频一区二区在线观看 | 午夜免费男女啪啪视频观看| 国产精品国产高清国产av| 校园人妻丝袜中文字幕| 国产伦一二天堂av在线观看| 亚洲av成人av| 亚洲国产精品久久男人天堂| 在线观看66精品国产| 成人鲁丝片一二三区免费| 好男人在线观看高清免费视频| 欧美+日韩+精品| 亚洲熟妇中文字幕五十中出| 免费黄网站久久成人精品| 91av网一区二区| 久久99热6这里只有精品| 国产高清不卡午夜福利| 校园人妻丝袜中文字幕| 99热精品在线国产| 久久久久久久久久黄片| 国产精品久久电影中文字幕| 亚洲av中文av极速乱| 亚洲精品影视一区二区三区av| 两性午夜刺激爽爽歪歪视频在线观看| 色哟哟·www| 老女人水多毛片| 精品久久久久久久久亚洲| 免费在线观看成人毛片| 国产一级毛片七仙女欲春2| 亚洲欧美日韩卡通动漫| av专区在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av天堂中文字幕网| 免费观看人在逋| 波多野结衣巨乳人妻| 亚洲一级一片aⅴ在线观看| 亚洲美女搞黄在线观看| 亚洲精品日韩av片在线观看| 午夜激情福利司机影院| 国产高潮美女av| 午夜福利视频1000在线观看| 国产高清不卡午夜福利| 国产精品久久视频播放| 亚洲丝袜综合中文字幕| 国产精品国产三级国产av玫瑰| 中文字幕免费在线视频6| 国产 一区精品| 欧美3d第一页| 黄片无遮挡物在线观看| 欧美色欧美亚洲另类二区| 日韩欧美 国产精品| 男插女下体视频免费在线播放| 三级经典国产精品| 99久久无色码亚洲精品果冻| 亚洲欧美日韩东京热| 舔av片在线| 日日摸夜夜添夜夜添av毛片| 色综合站精品国产| 国产午夜精品一二区理论片| 国产精品三级大全| 成人欧美大片| 男女边吃奶边做爰视频| 欧美最黄视频在线播放免费| 亚洲综合色惰| 在线播放国产精品三级| 欧美bdsm另类| 晚上一个人看的免费电影| 亚洲精品成人久久久久久| 性欧美人与动物交配| 亚洲精品日韩在线中文字幕 | 黄色日韩在线| 卡戴珊不雅视频在线播放| 黄色一级大片看看| 日韩成人av中文字幕在线观看| 久久久欧美国产精品| 亚洲精品粉嫩美女一区| 国产成人精品一,二区 | 色尼玛亚洲综合影院| 亚洲欧美成人综合另类久久久 | 午夜福利成人在线免费观看| 免费看a级黄色片| 国产日本99.免费观看| 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产| 日韩av不卡免费在线播放| av女优亚洲男人天堂| 亚洲国产欧洲综合997久久,| 在线免费观看不下载黄p国产| 神马国产精品三级电影在线观看| 99热6这里只有精品| av在线观看视频网站免费| 亚洲一区高清亚洲精品| 精品人妻视频免费看| 女的被弄到高潮叫床怎么办| 美女大奶头视频| 国产精品伦人一区二区| 免费观看a级毛片全部| 女人被狂操c到高潮| 国产 一区精品| 精品久久久久久成人av| 欧美日韩精品成人综合77777| 99热这里只有精品一区| 亚洲美女视频黄频| 直男gayav资源| 久久久久久久久大av| 久久久精品大字幕| 精品久久久久久久久av| 天堂影院成人在线观看| 国内揄拍国产精品人妻在线| 国内久久婷婷六月综合欲色啪| av免费在线看不卡| 国产 一区 欧美 日韩| 黑人高潮一二区| 99热精品在线国产| 黄色视频,在线免费观看| 内射极品少妇av片p| 天天躁日日操中文字幕| 欧美区成人在线视频| 久久国内精品自在自线图片| ponron亚洲| 青春草国产在线视频 | 国产一区二区激情短视频| 午夜免费男女啪啪视频观看| 青青草视频在线视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 99国产精品一区二区蜜桃av| 日本成人三级电影网站| 亚洲av免费高清在线观看|