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

    CeAuGa3的力學(xué)性質(zhì)及磁性的第一性原理計(jì)算

    2021-10-18 08:39:46程耀世劉本瓊
    關(guān)鍵詞:中國(guó)工程物理研究院核物理第一性

    程耀世, 劉本瓊

    (中國(guó)工程物理研究院核物理與化學(xué)研究所, 綿陽(yáng) 621999)

    1 Introduction

    It is very attractive to study the Cerium-based intermetallic compounds, since they exhibit lots of interesting behaviors resulting from the interplay between the single-impurity Kondo effect and the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions. Especially, the stoichiometric compounds with the general formula CeTX3or CeTyX4-y(Tis a transition metal,Xis a simple metal) have recently received great attention, because important new features have been observed in addition to their novel ground state properties which have been studied extensively. For examples, there is relatively low Cu-Al atomic site disorder[1-2]in CeCuAl3, which leads to the previous argument about its true crystal structure of either the ThCr2Si2-type[3]or the BaNiSn3-type[4-5], and indeed can influence the electronic properties. CeNiGe3undergoes antiferromagnetic (AFM) phase transition belowTN= 5.5 K at ambient pressure and exhibits a quantum critical point at about 6~7 GPa. The superconductivity in this material exists in a wide range of 1.7~9.3 GPa[6-7], and the spin-fluctuations may be not the hidden mechanism at the lower pressure region. Kimuraetal.[8]presented an unconventional pressure-induced superconductor CeRhSi3, where the superconductivity and the AFM state coexist, withTN=1.6 K at ambient pressure. Another similar example is the antiferromagnet CeIrSi3, where superconductivity appears in 1.8~3.5 GPa[9]. Although most of CeTX3compounds show AFM ordering at low temperature, CeAgAl3[10]is a scarce example of ferromagnetic compound withTC= 3.0 K. Further, selected members of this class offer the new insights that magnetoelastic coupling could result in the formation of a bound state between low-lying phonons and a crystal electric field (CEF) excitation. Adrojaetal. first reported the presence of an additional magnetic excitation in the neutron spectra of CeCuAl3measured in the paramagnetic phase[11], which cannot be explained by a pure CEF model. By generalizing the theoretical model proposed by Thalmeier and Fulde for the cubic system CeAl2[12-13]to the tetragonal point-group symmetry, a satisfactory interpretation for the magnetic inelastic neutron scattering spectra of CeCuAl3can be achieved. In recent studies of CeAuAl3[14-15], except for the pronounced vibronic bound state, a well-resolved anti-crossing of the CEF excitations with the acoustic phonon branch at zero magnetic field has also been observed by neutron scattering.

    In this work, we present detailed first-principles calculations of the structural and mechanical properties of CeAuGa3, which has been rarely studied and may behave qualitatively similarly to CeAuAl3. Total energies of different spin arrangements of CeAuGa3have been compared, and our results show that the total energy of ferromagnetic (FM) configuration is indeed the lowest,i.e., in the ground state a FM ordering state is favored. By further analyzing the paths of spin exchange coupling, it is found that all the spin exchange interactions are positive ferromagnetic type. The couplings between the Ce ions parallel to theabplane are comparable, while the third-nearest-neighbour coupling between the Ce ions within adjacent layers is relatively weak.

    2 Computational methods

    Several different spin configurations of CeAuGa3were examined by evaluating the exchange coupling constants between the Ce ions via mapping analysis based on first principles calculations. The Viennaabinitiosimulation package (VASP)[16-18]was used, with the standard frozen-core projector augmented wave (PAW) method, and the generalized gradient approximations (GGA) for the exchange and correlation corrections. The plane-wave cutoff energy was tested and set to 550 eV. According to the Monkhorst-Pack scheme, thek-point meshes in the Brillouin zone were sampled by 10×10×5 grids. For the elastic property calculations, a denserk-mesh sampling of 12×12×6 was adopted for accuracy. In order to obtain the phonon dispersion, Phonopy code[19]was used and a 2×2×2 supercell was constructed with thek-point meshes of 4×4×4. Other parameters included a convergence criterion of the Hellmann-Feynman force of 0.01 eV/?, and 10-8eV for the total energy tolerance.

    3 Discussion

    3.1 Spin exchange parameters

    According to Mocketal.[20], CeAuGa3crystallizes in the non-centrosymmetric BaNiSn3-type tetragonal structure (space group is I4mm, No.107). There are 10 atoms in the conventional unit cell, as shown in Fig. 1. The lattice parameters area= 4.337 ?,c= 10.66 ?, giving ac/aratio of 2.458. When doing geometry optimization, we start with the experimental geometries and calculate the dependence of the total energyEon the volumeV. For the tetragonal structure, we first calculate total energies for different volumes by fixing thec/aratio and changing the lattice parametera, and do a least-squares fit of theE-Vcurve to the third-order Birch-Murnaghan equation of state (EOS)[21],

    Fig.1 Ordered spin arrangements of Ce ions in the (a) FM, (b) AFM1 and (c) AFM2 states of CeAuGa3 (color online)The red dashed lines in (a) refer to the three spin exchange paths J1, J2 and J3 of Ce ions, respectively. For the ordered spin arrangement in AFM2 (c), only magnetic Ce sites are displayed for clarity in the 2 × 2 × 1 supercell

    (1)

    withE0being the equilibrium energy,V0the equilibrium volume, andB0′ the pressure derivative of the bulk modulusB0. Secondly, we fix the value of lattice parametera= 4.368 ? at the energy minimum, and change the volume by adjusting thec/aratio, as shown in Fig. 2. After the second EOS fit, the lattice parameters of CeAuGa3are finally determined asa= 4.384 ?,c= 10.655 ?,c/a= 2.43, and the equilibrium volumeV0= 204.75 ?3, which are in good agreement with the experimental data[20], as listed in Tab. 1. Among the different spin configurations of CeAuGa3(i.e., non-magnetic state NM, ferromagnetic state FM, and two different AFM states), the FM structure has the lowest total energy, indicating that the ground state prefers to be ferromagnetic ordering.

    Tab.1 Calculated lattice parameters (a,c), the c/a ratio, the equilibrium volume (V0) and relative energies (E0) of CeAuGa3, compared with previous experimental results

    Fig.2 The relationship between the total energy E and the unit cell volume V of CeAuGa3 (color online)

    The lines represent a least-squares fit to the third-order Birch-Murnaghan equation of state. The black stars denote that the data are obtained by fixingc/a= 2.458, and the red cycles correspond to the case of fixeda= 4.368 ?

    Based on the total energy calculations, one may analyze the magnetic state in terms of three spin-exchange parameters,i.e., the exchange couplingJ1between the nearest-neighbour Ce ions along thea(b) direction, the exchange couplingJ2between the next-nearest-neighbour Ce ions along the diagonal path within theab-plane, and the exchange couplingJ3between the third-nearest-neighbour Ce ions between the adjacent //ablayers, as denoted in Fig.1(a). To extract the values of the spin-exchange parameters, we express the total spin exchange interaction energies of the three ordered spin states in terms of the effective Heisenberg Hamiltonian,

    (2)

    whereJijdenotes the spin exchange parameter (fromJ1toJ3) for the spin exchange interaction between the spin sitesiandj,SiandSjare the spin angular momentum operators. By applying the energy expressions obtained for spin dimers[22-24], the total spin exchange energies of the three ordered spin states can be obtained as

    EFM=E0-50J1-25J2-50J3

    (3)

    EAFM1=E0-50J1-25J2+50J3

    (4)

    EAFM2=E0-175J1-62.5J2-200J3

    (5)

    By mapping the relative energies of the three-ordered spin states determined by first-principles calculations onto the corresponding relative energies determined from the above spin-exchange energies, one can obtainJ1= 0.572 meV,J2= 0.414 meV andJ3= 0.043 meV listed in Tab. 2. All the spin exchange parametersJ1,J2,J3are positive, indicating ferromagnetic interactions between the Ce ions, and further supporting the ferromagnetic ground state of CeAuGa3. The exchange coupling between the nearest-neighbour pairsJ1and that of the next-nearest-neighbour pairsJ2dominate, while the third-nearest-neighbour exchange couplingJ3may be negligible.

    ab.2 Ce-Ce distances associated with the spin exchange paths J1~J3 of CeAuGa3, and the values of J1~J3 determined by mapping analysis based on the present first-principles calculations

    3.2 Elastic properties

    For tetragonal structures, there are six independent elastic constants,i.e.,C11=C22,C12,C13=C23,C33,C44=C55, andC66. Under small deformations, the total energyEof a crystal associated with a strainecan be expressed as follows:

    (6)

    whereeiandejare components of the straine= (e1,e2,e3,e4,e5,e6). According to Eq. (6), all the elastic constantsCijcan be deduced by fitting the strain energy for particular deformations. We employ six different strain states for elastic constants calculations as listed in Tab. 3. For each strain type, 14 symmetric values of strains fromδ=-0.0175 to 0.0175 are taken to make the strain energy fit with a step of 0.0025. The corresponding energy-strain data is shown in Fig.3, and the calculated elastic constants are summarized in Tab. 4.

    Tab.3 Strain states used to determine the elastic constants of CeAuGa3

    Fig.3 Energy-strain data for CeAuGa3, lines connecting data points are shown to guide the eye only(color online)

    Tab.4 The calculated elastic constants Cij, bulk modulus B, shear modulus G, Young’s modulus Y and Poisson’s ratio ν of CeAuGa3

    One can find that all the predicted elastic constants fulfill the well-known Born-Huang stability criteria[25],

    C11-C12>0

    (7)

    C11+C33-2C13>0

    (8)

    C11>0,C33>0,C44>0,C66>0

    (9)

    2C11+C33+2C12+4C13>0

    (10)

    which confirms that CeAuGa3is mechanically stable. We have also checked the dynamical stability by the full phonon dispersion relation calculations[26]using the small displacement method[27], as shown in Fig. 4. There are in total 15 phonon branches, and all the calculated phonon frequencies along the high-symmetry directions Г-X-R-S0-Г and Г-M-S-N-P-G-M in the Brillouin zone are positive, suggesting a dynamically stable structure.

    Fig.4 Phonon dispersion curves of CeAuGa3 along high symmetric directions

    Based on the obtained six single-crystal elastic constants of CeAuGa3, one can further obtain its polycrystalline average isotropic bulk modulusBand shear modulusGaccording to the Voigt-Reuss-Hill approximation[28-30]:

    (11)

    (12)

    (13)

    (14)

    6/(C11-C12)+6/C44+3/C66]

    (15)

    The Young’s modulusYand Poisson’s ratioνcan be obtained as

    (16)

    The calculated bulk modulus, the shear modulus, Young’s modulus and Poisson’s ratio for CeAuGa3are also summarized in Tab. 4. For polycrystalline phases, a high (low) ratio of bulk to shear modulusB/Gis associated with ductility (brittleness) of the material[31]. The value of CeAuGa3is as high as 6.37, suggesting a high ductility of this compound.

    4 Conclusions

    Based on first-principles calculations, we have studied the structural properties for several different magnetic states of CeAuGa3. The total energy comparison shows that its ground state prefers to be ferromagnetic ordering. By mapping analysis of different spin states, one can find that all the exchanges between Ce ions are positive (ferromagnetic type). The exchange coupling between the nearest-neighbour pairs and that of the next-nearest-neighbour pairs of Ce ions dominate, while the third-nearest-neighbour interaction between the Ce spins within adjacent //ablayers is relatively weak. In additions, we have checked the mechanical and dynamical stabilities by computing the elastic properties and phonon dispersion, respectively. Our results imply that CeAuGa3is both mechanically and dynamically stable, since all the calculated elastic constants satisfy the well-known Born-Huang stability criteria, and without any imaginary phonon modes for all the wave vectors in the Brillouin zone.

    猜你喜歡
    中國(guó)工程物理研究院核物理第一性
    基于目標(biāo)航跡的引導(dǎo)誤差校正方法研究
    中國(guó)工程物理研究院
    軍工文化(2023年3期)2023-04-28 08:39:41
    DAAzF對(duì)DAAF熱性能的影響
    含能材料(2022年7期)2022-07-13 00:16:30
    科學(xué)文化評(píng)論(2021年3期)2021-11-22 00:51:14
    AuBe5型新相NdMgNi4-xCox的第一性原理研究
    SO2和NO2在γ-Al2O3(110)表面吸附的第一性原理計(jì)算
    宇宙中的“煉金術(shù)”
    《核物理介紹》教學(xué)中構(gòu)建學(xué)生自主學(xué)習(xí)激勵(lì)機(jī)制的探索性研究
    基于四傳感器的弱信號(hào)源定位方法
    傳感器世界(2019年9期)2019-03-17 18:52:46
    核物理研究在我國(guó)的主要發(fā)展動(dòng)向
    啦啦啦视频在线资源免费观看| 亚洲av美国av| av又黄又爽大尺度在线免费看| 男人舔女人的私密视频| 午夜av观看不卡| 久久人人爽人人片av| 欧美精品人与动牲交sv欧美| 热99国产精品久久久久久7| www.av在线官网国产| 国产精品欧美亚洲77777| 亚洲国产欧美日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美一区二区三区国产| 最近最新中文字幕大全免费视频 | 欧美成人午夜精品| 国产精品 国内视频| www日本在线高清视频| 国产有黄有色有爽视频| 一级毛片我不卡| av有码第一页| 丝袜美足系列| 欧美激情极品国产一区二区三区| 国产成人欧美| 国产无遮挡羞羞视频在线观看| 天堂8中文在线网| 成人手机av| 精品视频人人做人人爽| 中文字幕亚洲精品专区| 国产一区二区在线观看av| 亚洲国产毛片av蜜桃av| 国产精品免费大片| 熟女av电影| 久久影院123| 99国产精品免费福利视频| 国产99久久九九免费精品| 黄色一级大片看看| 丝袜美足系列| 人妻人人澡人人爽人人| 久久精品久久精品一区二区三区| 91精品三级在线观看| 色婷婷久久久亚洲欧美| 国产精品免费大片| 成年人黄色毛片网站| 午夜精品国产一区二区电影| 老司机午夜十八禁免费视频| 亚洲自偷自拍图片 自拍| 久久亚洲精品不卡| 国产高清国产精品国产三级| 精品少妇黑人巨大在线播放| 99国产精品99久久久久| 中文字幕另类日韩欧美亚洲嫩草| 一本综合久久免费| 亚洲黑人精品在线| 热re99久久精品国产66热6| 操出白浆在线播放| 91成人精品电影| 免费久久久久久久精品成人欧美视频| 午夜福利,免费看| 精品国产国语对白av| 国产成人精品无人区| 久热爱精品视频在线9| 最黄视频免费看| 精品国产一区二区三区四区第35| 中国国产av一级| 亚洲av日韩精品久久久久久密 | 国产黄色视频一区二区在线观看| 大香蕉久久成人网| 韩国精品一区二区三区| 啦啦啦在线观看免费高清www| a级片在线免费高清观看视频| av在线播放精品| 久久99一区二区三区| 麻豆乱淫一区二区| 亚洲av成人精品一二三区| 国产日韩欧美在线精品| 搡老乐熟女国产| videosex国产| 一边摸一边抽搐一进一出视频| 国产黄频视频在线观看| 十八禁高潮呻吟视频| 一级毛片电影观看| 日本黄色日本黄色录像| 又黄又粗又硬又大视频| 国产黄色免费在线视频| 欧美日韩一级在线毛片| 少妇的丰满在线观看| 免费观看a级毛片全部| 男女午夜视频在线观看| 亚洲情色 制服丝袜| 丝袜美腿诱惑在线| 麻豆乱淫一区二区| 日韩大片免费观看网站| 99re6热这里在线精品视频| 国产精品久久久久成人av| 亚洲人成电影观看| 极品少妇高潮喷水抽搐| 在线观看国产h片| 波野结衣二区三区在线| 亚洲欧美精品综合一区二区三区| av福利片在线| 赤兔流量卡办理| 久久精品国产综合久久久| 婷婷色麻豆天堂久久| 伊人久久大香线蕉亚洲五| 日本wwww免费看| 亚洲精品国产av成人精品| 91国产中文字幕| 男女无遮挡免费网站观看| 精品国产国语对白av| 免费在线观看完整版高清| 男男h啪啪无遮挡| 久久久久久久久久久久大奶| 日韩制服丝袜自拍偷拍| 男女国产视频网站| 色94色欧美一区二区| 搡老乐熟女国产| 久久国产精品男人的天堂亚洲| 啦啦啦视频在线资源免费观看| 亚洲精品日韩在线中文字幕| 日本91视频免费播放| 最近手机中文字幕大全| 美女视频免费永久观看网站| 欧美日韩av久久| 国产免费一区二区三区四区乱码| 久久久国产一区二区| 观看av在线不卡| 黄色片一级片一级黄色片| 99re6热这里在线精品视频| 一级黄片播放器| 国产亚洲午夜精品一区二区久久| 国产精品av久久久久免费| 成人国产一区最新在线观看 | 亚洲图色成人| 黄色怎么调成土黄色| 亚洲精品日韩在线中文字幕| av欧美777| 一级黄色大片毛片| 人人妻,人人澡人人爽秒播 | 极品少妇高潮喷水抽搐| 美女视频免费永久观看网站| 狂野欧美激情性bbbbbb| 婷婷色综合www| 一区二区三区精品91| av在线老鸭窝| 国产在视频线精品| 一级毛片 在线播放| 夫妻性生交免费视频一级片| 91字幕亚洲| 91精品国产国语对白视频| 如日韩欧美国产精品一区二区三区| 欧美老熟妇乱子伦牲交| 九色亚洲精品在线播放| 亚洲精品国产一区二区精华液| 国产免费福利视频在线观看| 新久久久久国产一级毛片| 91精品伊人久久大香线蕉| 亚洲一卡2卡3卡4卡5卡精品中文| 深夜精品福利| 老司机靠b影院| 中文字幕av电影在线播放| 欧美日韩一级在线毛片| 国产在线免费精品| 亚洲国产成人一精品久久久| 黄色毛片三级朝国网站| 欧美日韩一级在线毛片| 午夜福利免费观看在线| 99久久人妻综合| 人妻人人澡人人爽人人| 各种免费的搞黄视频| 51午夜福利影视在线观看| av在线老鸭窝| 久久久久网色| 精品久久蜜臀av无| 一级毛片黄色毛片免费观看视频| 久久久久久久国产电影| 国产成人免费无遮挡视频| 亚洲伊人久久精品综合| 狠狠婷婷综合久久久久久88av| 在线观看人妻少妇| 又黄又粗又硬又大视频| 又粗又硬又长又爽又黄的视频| 久久久久久久久久久久大奶| 赤兔流量卡办理| 热re99久久精品国产66热6| 亚洲精品久久午夜乱码| 久久精品久久久久久噜噜老黄| 国产精品一区二区在线观看99| 国产日韩欧美亚洲二区| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品999| 97精品久久久久久久久久精品| 国产精品一区二区免费欧美 | 国产亚洲欧美在线一区二区| 免费在线观看影片大全网站 | 国产精品熟女久久久久浪| 又粗又硬又长又爽又黄的视频| 色精品久久人妻99蜜桃| 久久亚洲国产成人精品v| 999精品在线视频| 精品少妇久久久久久888优播| 操出白浆在线播放| 国产高清不卡午夜福利| 一区二区三区激情视频| 国产精品三级大全| 视频在线观看一区二区三区| 人妻一区二区av| 欧美 日韩 精品 国产| 精品少妇久久久久久888优播| 国产欧美日韩综合在线一区二区| 老司机影院毛片| 国产xxxxx性猛交| 侵犯人妻中文字幕一二三四区| 一级毛片女人18水好多 | 国产精品99久久99久久久不卡| 91老司机精品| 菩萨蛮人人尽说江南好唐韦庄| 男女免费视频国产| 久久av网站| 国产野战对白在线观看| 狂野欧美激情性bbbbbb| 啦啦啦中文免费视频观看日本| 最新在线观看一区二区三区 | 日本av免费视频播放| 成在线人永久免费视频| 日韩欧美一区视频在线观看| 亚洲精品在线美女| 午夜91福利影院| 中文字幕色久视频| av国产久精品久网站免费入址| 丝袜美足系列| 欧美国产精品va在线观看不卡| 97人妻天天添夜夜摸| 蜜桃国产av成人99| 乱人伦中国视频| 秋霞在线观看毛片| 王馨瑶露胸无遮挡在线观看| 啦啦啦在线观看免费高清www| 日本色播在线视频| 精品免费久久久久久久清纯 | 久久精品久久精品一区二区三区| 91麻豆av在线| 女警被强在线播放| 一级片免费观看大全| 91精品三级在线观看| 丝袜人妻中文字幕| 国产片内射在线| 国产亚洲一区二区精品| 亚洲国产看品久久| av网站免费在线观看视频| 色综合欧美亚洲国产小说| 免费在线观看视频国产中文字幕亚洲 | 精品国产乱码久久久久久男人| 欧美日韩视频精品一区| 色综合欧美亚洲国产小说| 曰老女人黄片| 国产成人一区二区在线| 18在线观看网站| www.熟女人妻精品国产| 成年人午夜在线观看视频| 少妇精品久久久久久久| 性高湖久久久久久久久免费观看| 男女边摸边吃奶| 99国产精品一区二区蜜桃av | 纯流量卡能插随身wifi吗| 永久免费av网站大全| 亚洲美女黄色视频免费看| 国产真人三级小视频在线观看| 日本欧美视频一区| 精品人妻1区二区| 少妇精品久久久久久久| 啦啦啦 在线观看视频| 麻豆乱淫一区二区| 精品高清国产在线一区| 男人爽女人下面视频在线观看| 悠悠久久av| 日韩,欧美,国产一区二区三区| 色94色欧美一区二区| 中文字幕另类日韩欧美亚洲嫩草| 最黄视频免费看| 国产男人的电影天堂91| 亚洲精品成人av观看孕妇| 嫩草影视91久久| 中文字幕人妻丝袜一区二区| 亚洲av国产av综合av卡| 精品一区二区三区av网在线观看 | 国产欧美日韩一区二区三 | 国产免费现黄频在线看| 日韩一卡2卡3卡4卡2021年| 又黄又粗又硬又大视频| 人妻 亚洲 视频| 大香蕉久久成人网| 欧美精品av麻豆av| 国产欧美日韩一区二区三区在线| 婷婷色麻豆天堂久久| 精品久久久精品久久久| 无遮挡黄片免费观看| 国产伦理片在线播放av一区| 亚洲成人手机| 精品久久久精品久久久| 久久人人爽人人片av| 熟女av电影| 最新的欧美精品一区二区| 精品一区在线观看国产| a 毛片基地| 亚洲,欧美,日韩| 精品一区二区三区四区五区乱码 | 性色av一级| 久久毛片免费看一区二区三区| 亚洲av成人精品一二三区| 侵犯人妻中文字幕一二三四区| 肉色欧美久久久久久久蜜桃| 国精品久久久久久国模美| 国产欧美日韩一区二区三区在线| 成年人免费黄色播放视频| 在线观看人妻少妇| 国产黄频视频在线观看| 一区福利在线观看| 日本91视频免费播放| 宅男免费午夜| 少妇精品久久久久久久| 国产精品 国内视频| 中国国产av一级| 悠悠久久av| 90打野战视频偷拍视频| 18禁裸乳无遮挡动漫免费视频| 亚洲色图综合在线观看| 亚洲精品久久成人aⅴ小说| 日韩一卡2卡3卡4卡2021年| 三上悠亚av全集在线观看| 一区二区三区精品91| 一级黄色大片毛片| 最新的欧美精品一区二区| 国产精品一区二区精品视频观看| 日韩一本色道免费dvd| 不卡av一区二区三区| 亚洲九九香蕉| 我要看黄色一级片免费的| 国产精品.久久久| 亚洲 欧美一区二区三区| 亚洲第一青青草原| 久久影院123| 国产欧美日韩一区二区三 | 日本av免费视频播放| 国产片特级美女逼逼视频| 国产精品一区二区在线观看99| 亚洲精品久久午夜乱码| 人人妻人人爽人人添夜夜欢视频| 黄片播放在线免费| 亚洲成人国产一区在线观看 | 久久精品熟女亚洲av麻豆精品| 欧美日韩亚洲综合一区二区三区_| av有码第一页| 大陆偷拍与自拍| 大码成人一级视频| 欧美国产精品va在线观看不卡| 亚洲精品一二三| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩av久久| 午夜福利在线免费观看网站| 欧美中文综合在线视频| 日韩中文字幕欧美一区二区 | 一级毛片女人18水好多 | 天天躁日日躁夜夜躁夜夜| 熟女av电影| 欧美大码av| 国产精品久久久av美女十八| 欧美日韩福利视频一区二区| 亚洲图色成人| 少妇人妻久久综合中文| 亚洲精品一区蜜桃| 久久精品成人免费网站| 最近手机中文字幕大全| 亚洲九九香蕉| 五月天丁香电影| 香蕉丝袜av| 亚洲自偷自拍图片 自拍| 叶爱在线成人免费视频播放| 老司机亚洲免费影院| 美女中出高潮动态图| 日韩熟女老妇一区二区性免费视频| 黑人欧美特级aaaaaa片| 亚洲久久久国产精品| 丰满人妻熟妇乱又伦精品不卡| 赤兔流量卡办理| 国产精品欧美亚洲77777| 两个人免费观看高清视频| 亚洲欧洲日产国产| 午夜免费男女啪啪视频观看| 18禁裸乳无遮挡动漫免费视频| 涩涩av久久男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| 精品人妻1区二区| 日韩制服骚丝袜av| 亚洲欧美色中文字幕在线| 国产成人精品久久久久久| 人人澡人人妻人| 亚洲男人天堂网一区| 中国国产av一级| 欧美激情 高清一区二区三区| 久久久久久久久免费视频了| 免费在线观看完整版高清| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 国产成人一区二区在线| 精品一区二区三区四区五区乱码 | 国产真人三级小视频在线观看| 日韩伦理黄色片| 黄网站色视频无遮挡免费观看| 亚洲欧美精品自产自拍| 久久影院123| 人人妻,人人澡人人爽秒播 | 国产高清国产精品国产三级| 精品亚洲乱码少妇综合久久| 两性夫妻黄色片| 精品卡一卡二卡四卡免费| 久久久精品免费免费高清| 成人国语在线视频| 女人被躁到高潮嗷嗷叫费观| 婷婷色麻豆天堂久久| 又黄又粗又硬又大视频| 免费女性裸体啪啪无遮挡网站| 国产片特级美女逼逼视频| 欧美精品人与动牲交sv欧美| 国产亚洲一区二区精品| av福利片在线| 波多野结衣av一区二区av| 亚洲av男天堂| 亚洲av在线观看美女高潮| 这个男人来自地球电影免费观看| 咕卡用的链子| 色网站视频免费| 黑人欧美特级aaaaaa片| 亚洲成人免费av在线播放| 丝袜人妻中文字幕| 精品高清国产在线一区| 亚洲人成网站在线观看播放| 婷婷色综合大香蕉| 性色av一级| 国产成人精品久久久久久| 少妇猛男粗大的猛烈进出视频| 久久狼人影院| 精品久久久久久电影网| 亚洲欧美激情在线| 国产高清不卡午夜福利| 日日夜夜操网爽| 亚洲精品一二三| 一边亲一边摸免费视频| 久久人人爽人人片av| 伊人亚洲综合成人网| 又黄又粗又硬又大视频| 亚洲国产日韩一区二区| 成人午夜精彩视频在线观看| 日韩中文字幕视频在线看片| 蜜桃在线观看..| 黑人欧美特级aaaaaa片| 一本一本久久a久久精品综合妖精| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| 一本大道久久a久久精品| 天堂中文最新版在线下载| 午夜激情av网站| 99精国产麻豆久久婷婷| 精品少妇内射三级| 亚洲国产日韩一区二区| 美女主播在线视频| www.熟女人妻精品国产| 国产精品麻豆人妻色哟哟久久| 午夜老司机福利片| 国产黄色视频一区二区在线观看| 亚洲国产精品一区二区三区在线| 精品国产国语对白av| videosex国产| 国精品久久久久久国模美| 在线观看免费视频网站a站| 日韩一本色道免费dvd| 欧美日韩亚洲国产一区二区在线观看 | 国产免费现黄频在线看| 国产成人欧美在线观看 | 一级毛片黄色毛片免费观看视频| 国产免费视频播放在线视频| 国产一区二区三区综合在线观看| 亚洲成国产人片在线观看| videos熟女内射| 欧美日韩视频精品一区| 国产欧美日韩一区二区三区在线| 午夜av观看不卡| 丝袜在线中文字幕| 操出白浆在线播放| 亚洲成人免费电影在线观看 | 下体分泌物呈黄色| 日本a在线网址| 人妻人人澡人人爽人人| 国产免费现黄频在线看| 热re99久久国产66热| 黄色一级大片看看| 久久青草综合色| 久久鲁丝午夜福利片| 黄色a级毛片大全视频| 国产又色又爽无遮挡免| 午夜久久久在线观看| videos熟女内射| 视频区图区小说| 日韩制服丝袜自拍偷拍| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站 | 啦啦啦中文免费视频观看日本| 精品人妻熟女毛片av久久网站| 久久影院123| 国产主播在线观看一区二区 | 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 午夜福利乱码中文字幕| 各种免费的搞黄视频| 日本欧美视频一区| 后天国语完整版免费观看| 老鸭窝网址在线观看| 黄色怎么调成土黄色| 老司机靠b影院| 美女脱内裤让男人舔精品视频| 日韩免费高清中文字幕av| 色网站视频免费| 国产成人av激情在线播放| 一本一本久久a久久精品综合妖精| cao死你这个sao货| 国产男人的电影天堂91| 久久久国产欧美日韩av| 国产黄频视频在线观看| a级毛片黄视频| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲熟妇少妇任你| 午夜激情av网站| 久久99一区二区三区| 在现免费观看毛片| 在线天堂中文资源库| 色综合欧美亚洲国产小说| 精品亚洲成a人片在线观看| 国产成人一区二区三区免费视频网站 | 中国国产av一级| 深夜精品福利| 亚洲欧美日韩高清在线视频 | 国语对白做爰xxxⅹ性视频网站| 欧美日韩综合久久久久久| 黄色毛片三级朝国网站| 老司机影院毛片| 99久久99久久久精品蜜桃| 国产免费现黄频在线看| 咕卡用的链子| 久久久国产精品麻豆| 精品国产国语对白av| 免费女性裸体啪啪无遮挡网站| 亚洲欧美清纯卡通| 成在线人永久免费视频| 久久精品亚洲熟妇少妇任你| 国产女主播在线喷水免费视频网站| a级毛片黄视频| 最新的欧美精品一区二区| 丰满饥渴人妻一区二区三| 日本色播在线视频| 纯流量卡能插随身wifi吗| 女警被强在线播放| 精品亚洲成a人片在线观看| 狂野欧美激情性xxxx| 国产精品久久久av美女十八| 国产精品三级大全| kizo精华| 国产亚洲午夜精品一区二区久久| 亚洲,欧美精品.| 亚洲精品久久午夜乱码| 久久女婷五月综合色啪小说| 一级毛片我不卡| 亚洲七黄色美女视频| 亚洲精品第二区| 青草久久国产| 亚洲国产欧美日韩在线播放| 黄色片一级片一级黄色片| 男女午夜视频在线观看| 日韩大片免费观看网站| 男人舔女人的私密视频| 少妇被粗大的猛进出69影院| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三卡| 黑人猛操日本美女一级片| 国产91精品成人一区二区三区 | 午夜日韩欧美国产| 午夜福利在线免费观看网站| 美国免费a级毛片| 天天操日日干夜夜撸| 精品少妇一区二区三区视频日本电影| 91国产中文字幕| 老熟女久久久| av又黄又爽大尺度在线免费看| 久久久欧美国产精品| 两个人看的免费小视频| 岛国毛片在线播放| 精品久久久精品久久久| 亚洲七黄色美女视频| 亚洲成人免费电影在线观看 | 欧美日韩亚洲高清精品| 久久久久久亚洲精品国产蜜桃av| av在线老鸭窝| 免费av中文字幕在线| 久久鲁丝午夜福利片| 亚洲欧美一区二区三区久久| 国产精品久久久久久人妻精品电影 | 黑丝袜美女国产一区| 在线 av 中文字幕| 亚洲精品自拍成人| 中文字幕亚洲精品专区| 成人国语在线视频| 久久狼人影院| 色综合欧美亚洲国产小说| 国产精品国产av在线观看| 国产精品免费大片| 久久毛片免费看一区二区三区| 搡老乐熟女国产|