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

    Anharmonic Properties of Aluminum from Direct Free Energy Interpolation Method?

    2018-01-24 06:23:17ZhiGuoZhao趙志國JunShengSun孫俊生XiuLuZhang張修路HaiFengYang楊海峰andZhongLiLiu劉中利
    Communications in Theoretical Physics 2017年5期
    關(guān)鍵詞:海峰修路

    Zhi-Guo Zhao(趙志國),Jun-Sheng Sun(孫俊生),Xiu-Lu Zhang(張修路),Hai-Feng Yang(楊海峰), and Zhong-Li Liu(劉中利),*

    1College of Physics and Electric Information,Luoyang Normal University,Luoyang 471022,China

    2The Unit 63615 of People’s Liberation Army,Kuerle 841001,China

    3Laboratory for Extreme Conditions Matter Properties,Southwest University of Science and Technology,Mianyang 621010,China

    1 Introduction

    Like platinum and tungsten,aluminum is another frequently used pressure standard material in high-pressure experiments,such as the shock wave experiment.The pressure-volume-temperature equation of states(EOS)provides key information for it as a pressure standard.[1]Al is a typicalsp-bonded simple metal and crystallized in face-centered-cubic structure under ambient conditions.With the simple electronic and lattice structures,it is often taken as a typical prototype for theoretical calculations.[2?3]For the theoretical calculation of free energy,the accurate anharmonic effects are very hard to be taken into account in the theoretical EOS in previous methods.

    Usually,the quasi-harmonic approximation(QHA)was used to calculate the temperature effects and it only takes account into the volume dependence of phonon frequencies of lattice vibrations.However,in real materials the phonon frequencies depend largely on temperature,especially in the highly anharmonic solids or at the temperature close to melting.This is because that the QHA omits part of anharmonicity and is only accurate when temperature below Debye temperature.[4?6]In order to overcome these diffculties of QHA,several anharmonic free energy calculation techniques have been developed,such as the thermodynamic integration technique[7?8]and the self-consistent ab initio lattice dynamics(SCAILD)method.[9?10]

    Molecular dynamics(MD)simulations include anharmonicity to any order beyond Debye temperature,[11?12]especially applicable to high-temperature free energy calculations.Thus it can correct the free energy of the QHA considerably in the free energy calculations.In this paper,we apply the MD simulations to extract anharmonic free energy of Al and then construct more accurate EOS beyond the QHA.

    The rest of the paper is organized as follows.In Sec.2,we describe the details about the computation method.We illustrate the results in detail and make deep discussions in Sec.3.The conclusions are presented in Sec.4.

    2 Computational Method

    2.1 Quasi-Harmonic Approximation

    The free energy of Al was first calculated by the quasi-harmonic approximation and then corrected by taking into account the phonon-phonon interactions.The quasi-harmonic approximation(QHA)describes volumedependent thermal effects based on lattice vibration models known as phonon,neglecting the phonon-phonon interactions.In the framework of the QHA,the Helmholtz free energy of a crystal system is,

    whereEstatic(V)is the zero-temperature energy of a staticlattice with volumeV.The termFzpis the zero-point motion energy written as,

    whereg(ω,0K)is the 0 K phonon density of states(PDOS)of the frequencyω.

    The termFph(V,T)in Eq.(1)is the lattice vibrational free energy,and is calculated from,

    2.2 Anharmonicity Corrections

    In the harmonic approximation,the PDOSg(ω,0K)is calculated from the force constants using either the direct method[13?14]or the density functional perturbation theory.[5,15]In the harmonic approximation,the phononphonon interaction contribution to free energy are neglected,resulting in large errors of thermodynamic properties at high temperature.Therefore,the phonons interactions must be fully taken into account to include the rest of the anharmonic effects.However,the rest anharmonic free energy calculation is very challenging.As we know that the explicit anharmonicity are naturally implied in molecular dynamics(MD)simulations.We can take into account the rest anharmonic effects from MD simulations.The PDOS atT,g(ω,T),can be extracted from the Fourier transform of autocorrelation of atomic position function r(t)in a constant temperature MD simulation.[16]

    According to the Wiener–Khintchine theorem,[17?18]the autocorrelation of position function r(t)is given by

    The PDOS atTis,

    The Helmholtz free energy of a crystal system is then rewritten as,

    whereEstatic(V)is the total energy of a static lattice with volumeV.

    In Eq.(6),the termFzpis the zero-point motion energy of the lattice written as,

    The termFph(V)in Eq.(6)is calculated from,

    After obtaining the finite number of temperature dependent PDOS(TD-PDOS)data,we use a much straitforward and simple method named the direct free energy interpolation(DFEI)method,to construct full anharmonic effects at any temperature.

    The DFEI method extracts anharmonic free energy at any temperature based on the limited several TD-PDOS data from MD simulations.The free energiesF(V,T)can be directly deduced from the TD-PDOS data at the specific temperaturesTs with the temperature interval 50 K,according to Eq.(6).The free energies at any temperatures were obtained by interpolating these resulting free energy data using the numerical analysis technique of the basis spline(B-spline).Finally,the anharmonic effects of other thermodynamic properties were corrected.

    2.3 The Interatomic Potential of Al

    We perform classical MD simulations to simulate Al at different pressure and temperature.Then,we gather the trajectories of all the atoms in crystal to calculate phonon density of states.The embedded-atom-method(EAM)[19?20]potential developed by Mishin et al.[21?22]was used to describe the interatomic interactions of Al atoms.For a metal system containingNatoms,the total potential energyEtotis the sum of the embedding energyFand a pair potential?,

    whererijis the distance between the atomsiandj.The functionFi(ρi)is the energy to embed the atomiinto the background electron densityρi,which is the superposition of the atomic densities,

    TheMD simulationswereconducted with the large-scale atomic/molecular massively parallel simulator(LAMMPS)[23?24]package.The simulation box constructed from the multiplication 14×14×14 of the facecentered-cubic(fcc)conventional unit cell including 10976 atoms.The simulations were conducted for all the supercells with different volumes in the canonical ensemble(NVT).The periodic boundary condition was used for all the atoms in the simulation box.The time step was 1 fs and the total number of time steps were 10000.From Fig.1,we see that 10000 MD steps are able to converge the PDOS well.

    Fig.1 (Color online)The PDOS obtained using different numbers of MD steps.

    3 Results and Discussions

    3.1 High-Temperature Phonon Dispersion and DOS

    The phonon dispersion curve of Al were calculated using the direct method from molecular dynamics simulations.[25]The calculated phonon dispersion curve of Al are compared well with experimental data at 298 K[26]in Fig.2(a).The frequency of lattice vibration decreases with temperature until very small negative frequencies occur at 900 K as shown in the circle of Fig.2(a).Such very small negative frequencies often result from the numerical noises in the direct phonon calculation method with molecular dynamics simulations.[26]The good agreement of the calculated phonon dispersion curve with experimental data show the validity of the EAM potential for Al.The TD-PDOSs of Al with atomic volumeV0=16.61?A3are obtained by Fourier transform of autocorrelation of atomic position function and plotted in Fig.2(b).The PDOSs are from 0 to 900 K with the interval of 100 K.The density of states of higher frequencies increase with temperature,attributed to the increasing collective excitations of higher frequencies from the phonon-phonon interactions.On the contrary,the lower frequencies of vibrations are inhibited with increasing temperature.The 900 K PDOS does not show any negative frequency and thus does not affect the free energy calculation.

    The volumes of supercells used for PDOS calculations had volumes varying from 1.08V0to 0.64V0with the interval of 0.1V0,whereV0is the equilibrium volume at ambient pressure.The simulation temperature for each volume ranged from 0 K to 950 K with the interval of 50 K.The volume range and interval made a good sample for the description of the equation of states of Al(Fig.3(a))compared with the experimental data,[27]indicating that the anharmonic free energy interpolation based on the PDOS with the temperature interval of 50 K is suffcient.

    Fig.2 (a)The comparison of the calculated 0 GPa phonon dispersion curves of Al with experimental data at 298 K.[26]The circle indicates the negative frequencies at Γ point.(b)The phonon DOS of Al at atomic volume of 16.61?A3from 0 to 900 K with the interval of 100 K along the arrow direction.

    3.2 Thermal Expansion Coeffcients

    We calculate the Helmholtz free energy from the PDOS according to Eq.(1).The third-order Birch–Murnaghan equation of states[28]was fitted to the calculated free energyF–Vdata at each temperature. The third-order Birch–Murnaghan equation of states is given by,

    whereη=(V0/V)1/3,andF0,V0,B0andB′are fitting parameters.The thermodynamic properties and the anharmonic effects are automatically analyzed using the our PhaseGO package.[29?30]

    The accuracy of the calculated thermal properties was checked by the thermal expansion coeffcient defined by,The calculated thermal expansion coeffcients as a function of temperature are shown in Fig.3(b).The QHA results are in good agreement with experimental data[31]below 400 K,while when the temperature is above 400 K the zero-pressure thermal expansion coeffcient gradually deviates from experiment.This attributes to the neglect of anharmonicity caused by phonon-phonon interactions.This deviation is also observed in the calculations for other transition metals,such as Ta,[32]Pt,[33]and Pd.[34]However,the DFEI results show good agreement with the experimental data[31]in Fig.3(b).

    Fig.3 (a)The comparison of the calculated isothermal compression curve with experimental data.[27](b)The calculated thermal expansion coeffcients of Al compared with experimental data.[31]

    3.3 Heat Capacity

    Fig.4 The specific heat of Al at constant volume and at constant pressure compared with the experimental data.[35]

    whereFis the Helmholtz free energy calculated from Eq.(6).The thermal expansion caused by anharmonic effects results in a difference betweenCPandCV.The difference between the two can be written as

    whereαVis the volume thermal expansion coeffcient andB0is the bulk modulus.Because the heat capacities are

    The specific heat at constant volume was calculated by obtained from a second derivative of the Helmholtz free energy,they are more sensitive to the tiny errors in the Helmholtz free energy.

    Figure 4 showsCVandCPas a function of temperature at 0 GPa.The calculatedCVaccords reasonably well with experimental data[35]below 600 K,and first increases dramatically as pressure increases and then finally approaches to 3R.The DFEI and QHA results are very similar.The calculatedCPis also in good comparison with experimental data.[35]The better agreement with experimental data is also found after anharmonic correction(DFEI),especially whenTis beyond 800 K.

    3.4 Grüneisen Parameter and Bulk Modulus

    The thermodynamic Grüneisen parameter was achieved by

    whereBTis isothermal bulk modulus.The isothermal bulk modulus was calculated from

    The calculated isothermal bulk modulus versus temperature curves at different pressures are shown in Fig.5(a).The isothermal bulk modulus curve at 0 GPa is very good agreement with the experimental data.[36]The calculated thermodynamic Grüneisen parameter is compared with the ab initio data from Ref.[37]in Fig.5(b).The calculated thermodynamic Grüneisen parameter is also in very good agreement with the reference data from Ref.[38].

    Fig.5(a)The calculated isothermal bulk modulus of Al at different pressures,in comparison with zero-pressure data[36](open circles).(b)The calculated thermodynamic Grüneisen parameter compared with ab initio data from Ref.[37].

    4 Conclusion

    In conclusion,we calculate and correct the thermodynamic properties of Al using the direct calculation method of full free energy of lattice vibrations,i.e.,the DFEI method which makes accurate anharmonic corrections for the thermodynamic properties of materials beyond the QHA.The calculated thermal expansion coeffcients of Al using DFEI method agree better with experiment than the QHA results.Furthermore,the calculated isotherms,the constant volume and constant pressure heat capacities,the bulk moduli using DFEI are all in good agreement with experimental data.This indicates that the DFEI method is indeed an effcient and accurate method to determining the thermodynamic properties of materials from the high temperature phonon density of states.

    [1]C.J.Pickard and R.J.Needs,Nat.Mater.9(2010)624.

    [2]R.Martin,Nature(London)400(1999)117.

    [3]J.C.Boettger and S.B.Trickey,Phys.Rev.B 53(1996)3007.

    [4]U.Hansen,P.Vogl,and V.Fiorentini,Phys.Rev.B 60(1999)5055.

    [5]S.Baroni,S.D.Gironcoli,A.D.Corso,and P.Giannozzi,Rev.Mod.Phys.73(2001)515.

    [6]S.Baroni,P.Giannozzi,and E.Isaev,Rev.Miner.Geochem.71(2010)39.

    [7]D.Alf`e,G.D.Price,and M.J.Gillan,Phys.Rev.B 64(2001)45123.

    [8]D.Alf`e,G.D.Price,and M.J.Gillan,Phys.Rev.B 65(2002)165118.

    [9]P.Souvatzis,O.Eriksson,M.I.Katsnelson,and S.P.Rudin,Phys.Rev.Lett.100(2008)95901.

    [10]P.Souvatzis,O.Eriksson,M.I.Katsnelson,and S.P.Rudin,Comput.Mater.Sci.44(2009)888.

    [11]I.Errea,M.Calandra,and F.Mauri,Phys.Rev.B 89(2014)64302.

    [12]R.Car and M.Parrinello,Phys.Rev.Lett.55(1985)2471.

    [13]D.Alf`e,Comput.Phys.Commun.180(2009)2622.

    [14]A.Togo,F.Oba,and I.Tanaka,Phys.Rev.B 78(2008)134106.

    [15]X.Gonze,Phys.Rev.B 55(1997)10337.

    [16]M.Thomas,M.Brehm,R.Fligg,et al.,Phys.Chem.Chem.Phys.15(2013)6608.

    [17]N.Wiener,Acta Math.55(1930)117.

    [18]A.Khintchine,Math.Ann.109(1934)604.

    [19]M.S.Daw and M.I.Baskes,Phys.Rev.B 29(1984)6443.

    [20]S.M.Foiles,M.I.Baskes,and M.S.Daw,Phys.Rev.B 33(1986)7983.

    [21]J.M.Winey,A.Kubota,and Y.M.Gupta,Modell.Simul.Mater.Sci.Eng.17(2009)55004.

    [22]J.M.Winey,A.Kubota,and Y.M.Gupta,Modell.Simul.Mater.Sci.Eng.18(2010)29801.

    [23]S.J.Plimpton,J.Comput.Phys.117(1995)1.

    [24]http://lammps.sandia.gov.

    [25]L.T.Kong,Comput.Phys.Commun.182(2011)2201.

    [26]R.Stedman and G.Nilsson,Phys.Rev.145(1966)492.

    [27]A.Hnstrm and P.Lazor,J.Alloy Compd.305(2000)209.

    [28]F.Birch,Phys.Rev.71(1947)809.

    [29]Z.L.Liu,Comput.Phys.Commun.191(2015)150.

    [30]Z.L.Liu,Comput.Phys.Commun.197(2015)341.

    [31]F.C.Nix and D.MacNair,Phys.Rev.60(1941)597.

    [32]Z.L.Liu,L.C.Cai,X.R.Chen,et al.,J.Phys.:Condens.Matter 21(2009)95408.

    [33]T.Sun,K.Umemoto,Z.Wu,J.C.Zheng,and R.M.Wentzcovitch,Phys.Rev.B 78(2008)24304.

    [34]Z.L.Liu,J.H.Yang,L.C.Cai,F.Q.Jing,and D.Alf`e,Phys.Rev.B 83(2011)144113.

    [35]M.Forsblom,N.Sandberg,and G.Grimvall,Phys.Rev.B 69(2004)165106.

    [36]A.Dewaele,P.Loubeyre,and M.Mezouar,Phys.Rev.B 70(2004)94112.

    [37]J.L.Nie,L.Ao,F.A.Zhao,M.Jiang,and X.T.Zu,Can.J.Phys.93(2015)55004.

    [38]W.B.Holzapfel,M.Hartwig,and W.Sievers,J.Phys.Chem.Ref.Data 30(2001)515.

    猜你喜歡
    海峰修路
    以牙還牙
    科教新報(2024年51期)2024-12-11 00:00:00
    Progress and challenges in magnetic skyrmionics
    活著
    歌海(2022年1期)2022-03-29 21:39:55
    小豬琪琪修路
    倪海峰
    兒童大世界(2019年3期)2019-04-11 03:33:38
    開常古道上有個“修路節(jié)”
    中國公路(2017年8期)2017-07-21 14:26:20
    My School
    Nonlinear symbolic LFT model for UAV
    修路
    小說月刊(2015年8期)2015-04-19 02:41:24
    修路
    小說月刊(2014年8期)2014-04-19 02:39:18
    色吧在线观看| 国产高清激情床上av| 特级一级黄色大片| av黄色大香蕉| 在线观看免费午夜福利视频| 国产成人福利小说| 美女免费视频网站| 国产亚洲精品久久久久久毛片| 免费一级毛片在线播放高清视频| 亚洲第一欧美日韩一区二区三区| 欧美成狂野欧美在线观看| 亚洲av第一区精品v没综合| 国产精品免费一区二区三区在线| 国产三级在线视频| 亚洲一区二区三区色噜噜| 亚洲av中文字字幕乱码综合| 别揉我奶头~嗯~啊~动态视频| 69人妻影院| 热99在线观看视频| 久久久久亚洲av毛片大全| 一个人看的www免费观看视频| 免费人成在线观看视频色| 国产单亲对白刺激| 色精品久久人妻99蜜桃| 丰满的人妻完整版| 嫩草影视91久久| 小说图片视频综合网站| 色尼玛亚洲综合影院| 一级作爱视频免费观看| 精品免费久久久久久久清纯| 97超视频在线观看视频| 亚洲乱码一区二区免费版| 熟妇人妻久久中文字幕3abv| 亚洲av成人精品一区久久| 午夜激情欧美在线| 国产精华一区二区三区| 日本a在线网址| 精品乱码久久久久久99久播| 国产野战对白在线观看| 黄色日韩在线| 欧美绝顶高潮抽搐喷水| 亚洲国产精品久久男人天堂| 老司机午夜十八禁免费视频| 中国美女看黄片| 国产三级黄色录像| 99久久精品热视频| 欧美+亚洲+日韩+国产| 小说图片视频综合网站| 国产精品99久久久久久久久| 国产精品一区二区免费欧美| 国产v大片淫在线免费观看| 精品无人区乱码1区二区| 久久久国产成人免费| 69av精品久久久久久| 午夜福利免费观看在线| 男女下面进入的视频免费午夜| 国产三级中文精品| 麻豆久久精品国产亚洲av| 国产极品精品免费视频能看的| av天堂中文字幕网| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕av在线有码专区| 国产高清videossex| 欧美极品一区二区三区四区| 非洲黑人性xxxx精品又粗又长| 十八禁人妻一区二区| 久久久久久久精品吃奶| 欧美zozozo另类| 俄罗斯特黄特色一大片| 久久香蕉国产精品| 欧美性感艳星| 美女高潮喷水抽搐中文字幕| 国产野战对白在线观看| ponron亚洲| 变态另类成人亚洲欧美熟女| 一个人看视频在线观看www免费 | 午夜免费成人在线视频| 国产三级中文精品| 成人永久免费在线观看视频| 亚洲第一欧美日韩一区二区三区| av片东京热男人的天堂| 欧美日韩一级在线毛片| 1024手机看黄色片| 欧美乱妇无乱码| 免费人成视频x8x8入口观看| 俄罗斯特黄特色一大片| 男女床上黄色一级片免费看| 动漫黄色视频在线观看| 精品不卡国产一区二区三区| 小蜜桃在线观看免费完整版高清| 在线看三级毛片| 成人特级av手机在线观看| 天堂动漫精品| 亚洲,欧美精品.| 村上凉子中文字幕在线| 亚洲专区中文字幕在线| 三级男女做爰猛烈吃奶摸视频| 手机成人av网站| 亚洲欧美日韩东京热| 啦啦啦观看免费观看视频高清| 亚洲第一欧美日韩一区二区三区| 欧美成狂野欧美在线观看| 亚洲国产精品999在线| 午夜福利在线观看免费完整高清在 | svipshipincom国产片| 麻豆国产av国片精品| 香蕉丝袜av| 88av欧美| 亚洲精品色激情综合| 美女黄网站色视频| 国产免费av片在线观看野外av| av国产免费在线观看| 此物有八面人人有两片| 好男人电影高清在线观看| 免费在线观看影片大全网站| 午夜两性在线视频| 亚洲成人免费电影在线观看| 一个人看的www免费观看视频| 久久久精品大字幕| 国产色爽女视频免费观看| 国产精品永久免费网站| 中文字幕人妻熟人妻熟丝袜美 | 国产高潮美女av| 国产亚洲精品久久久久久毛片| 亚洲成人免费电影在线观看| netflix在线观看网站| 黑人欧美特级aaaaaa片| svipshipincom国产片| 少妇的丰满在线观看| 熟女少妇亚洲综合色aaa.| 桃色一区二区三区在线观看| 最新在线观看一区二区三区| 欧美黄色淫秽网站| av女优亚洲男人天堂| 一本久久中文字幕| 国产成人影院久久av| 欧美成狂野欧美在线观看| 两个人视频免费观看高清| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美一区二区三区黑人| 亚洲av日韩精品久久久久久密| 90打野战视频偷拍视频| 深爱激情五月婷婷| 亚洲性夜色夜夜综合| 久久久精品欧美日韩精品| 日本免费a在线| 免费人成视频x8x8入口观看| 啪啪无遮挡十八禁网站| 欧美一区二区亚洲| 亚洲性夜色夜夜综合| 在线免费观看的www视频| 国产高清视频在线观看网站| 免费av毛片视频| 午夜免费观看网址| 国产男靠女视频免费网站| 97人妻精品一区二区三区麻豆| 波多野结衣高清无吗| 精品国产超薄肉色丝袜足j| 国产精品1区2区在线观看.| 亚洲av五月六月丁香网| 亚洲人成网站在线播放欧美日韩| 一区二区三区激情视频| 一区二区三区激情视频| 精品一区二区三区视频在线观看免费| 亚洲真实伦在线观看| 精品国内亚洲2022精品成人| 级片在线观看| 国产精品 欧美亚洲| 国产在视频线在精品| 怎么达到女性高潮| 欧美高清成人免费视频www| 少妇的逼水好多| 精品久久久久久久人妻蜜臀av| 国产av麻豆久久久久久久| 国产亚洲精品综合一区在线观看| 国产成+人综合+亚洲专区| 欧美日韩精品网址| 手机成人av网站| 在线观看66精品国产| 中文亚洲av片在线观看爽| 亚洲中文字幕一区二区三区有码在线看| 精品人妻偷拍中文字幕| 搡老岳熟女国产| 日韩欧美一区二区三区在线观看| 日本精品一区二区三区蜜桃| 国产一区二区三区在线臀色熟女| 宅男免费午夜| 51国产日韩欧美| 国产69精品久久久久777片| www.熟女人妻精品国产| 亚洲精品国产精品久久久不卡| 精品国产美女av久久久久小说| 国产一区二区亚洲精品在线观看| 午夜亚洲福利在线播放| 一边摸一边抽搐一进一小说| 精品一区二区三区视频在线 | 欧美日韩精品网址| 每晚都被弄得嗷嗷叫到高潮| 一进一出抽搐gif免费好疼| 成人国产一区最新在线观看| 国产一区在线观看成人免费| 国产精品日韩av在线免费观看| 91在线精品国自产拍蜜月 | 在线观看舔阴道视频| 国内精品一区二区在线观看| 欧美区成人在线视频| 国产三级中文精品| 国产精品香港三级国产av潘金莲| 久久亚洲精品不卡| 97超视频在线观看视频| 亚洲国产欧美人成| 欧美在线一区亚洲| 男人和女人高潮做爰伦理| www国产在线视频色| 级片在线观看| 亚洲中文日韩欧美视频| 国产亚洲精品一区二区www| 国产亚洲精品久久久久久毛片| 日韩人妻高清精品专区| av欧美777| 国产av不卡久久| 天天躁日日操中文字幕| 可以在线观看的亚洲视频| 亚洲国产精品久久男人天堂| 观看免费一级毛片| 亚洲无线观看免费| 国产老妇女一区| 欧美一区二区精品小视频在线| 午夜福利在线在线| 精品人妻偷拍中文字幕| 亚洲美女黄片视频| 黄片小视频在线播放| 国产黄色小视频在线观看| 一本一本综合久久| 一卡2卡三卡四卡精品乱码亚洲| netflix在线观看网站| 午夜福利在线观看免费完整高清在 | 波多野结衣高清无吗| 日韩欧美 国产精品| 五月伊人婷婷丁香| 黄色视频,在线免费观看| 亚洲人成电影免费在线| 亚洲av电影不卡..在线观看| 99热这里只有是精品50| 人人妻人人看人人澡| 麻豆成人av在线观看| av中文乱码字幕在线| 哪里可以看免费的av片| 欧美激情在线99| 99久久精品国产亚洲精品| 亚洲人成网站高清观看| 无限看片的www在线观看| 动漫黄色视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲av成人av| a在线观看视频网站| 精品日产1卡2卡| 色av中文字幕| 一级毛片高清免费大全| 亚洲成人免费电影在线观看| 国产精品久久久久久亚洲av鲁大| 国产免费av片在线观看野外av| 欧美日韩中文字幕国产精品一区二区三区| 丝袜美腿在线中文| 非洲黑人性xxxx精品又粗又长| 丰满人妻熟妇乱又伦精品不卡| 九色成人免费人妻av| 18禁国产床啪视频网站| 国产一区二区三区在线臀色熟女| 亚洲av免费高清在线观看| 久99久视频精品免费| 国产aⅴ精品一区二区三区波| 夜夜看夜夜爽夜夜摸| 亚洲成人中文字幕在线播放| 精品久久久久久久毛片微露脸| 一本久久中文字幕| 日本黄色片子视频| 看黄色毛片网站| 国产一区二区三区视频了| 国产一区二区在线观看日韩 | 免费大片18禁| av女优亚洲男人天堂| 免费观看精品视频网站| 不卡一级毛片| 国产精品一区二区三区四区免费观看 | 亚洲精品国产精品久久久不卡| 露出奶头的视频| 亚洲性夜色夜夜综合| www.www免费av| 欧美黄色淫秽网站| 高潮久久久久久久久久久不卡| 国产毛片a区久久久久| 亚洲中文字幕一区二区三区有码在线看| 国产午夜精品久久久久久一区二区三区 | 欧美成人免费av一区二区三区| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| 首页视频小说图片口味搜索| 99热只有精品国产| 香蕉久久夜色| 老司机福利观看| 99国产综合亚洲精品| 日本黄大片高清| 99在线视频只有这里精品首页| 久久精品人妻少妇| 国产伦精品一区二区三区四那| 男插女下体视频免费在线播放| 国产精品,欧美在线| 国产视频内射| 精品不卡国产一区二区三区| 欧美性猛交黑人性爽| 99久久成人亚洲精品观看| 每晚都被弄得嗷嗷叫到高潮| 久久精品91无色码中文字幕| 久久精品国产综合久久久| 香蕉av资源在线| 69av精品久久久久久| 51国产日韩欧美| 国产欧美日韩精品亚洲av| 国产精品98久久久久久宅男小说| av黄色大香蕉| 午夜激情福利司机影院| 国模一区二区三区四区视频| netflix在线观看网站| 美女黄网站色视频| 欧美另类亚洲清纯唯美| 日本免费一区二区三区高清不卡| 一边摸一边抽搐一进一小说| 国产精品乱码一区二三区的特点| 综合色av麻豆| 人妻丰满熟妇av一区二区三区| 在线观看一区二区三区| 免费在线观看成人毛片| 999久久久精品免费观看国产| 99国产精品一区二区蜜桃av| 在线观看舔阴道视频| www.999成人在线观看| 午夜福利成人在线免费观看| 亚洲真实伦在线观看| 真人一进一出gif抽搐免费| h日本视频在线播放| 18禁在线播放成人免费| 精品熟女少妇八av免费久了| 国产精品自产拍在线观看55亚洲| 人人妻人人看人人澡| 国产精品日韩av在线免费观看| 日韩欧美 国产精品| 99热这里只有精品一区| 国产伦精品一区二区三区视频9 | 午夜福利成人在线免费观看| 欧美日韩福利视频一区二区| 十八禁人妻一区二区| 国产探花极品一区二区| 舔av片在线| 免费在线观看日本一区| 18禁国产床啪视频网站| 亚洲国产精品久久男人天堂| www日本黄色视频网| 在线免费观看不下载黄p国产 | 老司机在亚洲福利影院| 性色av乱码一区二区三区2| 一区二区三区国产精品乱码| АⅤ资源中文在线天堂| 真实男女啪啪啪动态图| 他把我摸到了高潮在线观看| 精品99又大又爽又粗少妇毛片 | 人妻夜夜爽99麻豆av| 国模一区二区三区四区视频| 婷婷精品国产亚洲av| 一夜夜www| 国产亚洲精品av在线| 有码 亚洲区| 久久久久久久亚洲中文字幕 | 中亚洲国语对白在线视频| 伊人久久精品亚洲午夜| 国产一区二区三区在线臀色熟女| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 精品久久久久久久久久久久久| 午夜福利欧美成人| 欧美一级毛片孕妇| 日韩成人在线观看一区二区三区| 国产亚洲欧美98| 成人国产一区最新在线观看| 精品熟女少妇八av免费久了| 12—13女人毛片做爰片一| 高清毛片免费观看视频网站| 在线观看66精品国产| 亚洲aⅴ乱码一区二区在线播放| 国产精品99久久久久久久久| 叶爱在线成人免费视频播放| 日韩欧美在线乱码| 成人特级av手机在线观看| 99热精品在线国产| 国产精品香港三级国产av潘金莲| 国产v大片淫在线免费观看| 男女午夜视频在线观看| 成年免费大片在线观看| 精品人妻1区二区| 午夜激情欧美在线| 亚洲精品亚洲一区二区| 黄色片一级片一级黄色片| 人妻夜夜爽99麻豆av| 午夜福利视频1000在线观看| 欧美不卡视频在线免费观看| 有码 亚洲区| 制服人妻中文乱码| 高清日韩中文字幕在线| 国产色爽女视频免费观看| 香蕉av资源在线| 日韩欧美国产在线观看| 欧美午夜高清在线| xxx96com| www.999成人在线观看| 欧美一级a爱片免费观看看| 十八禁网站免费在线| 美女被艹到高潮喷水动态| 久久精品综合一区二区三区| 韩国av一区二区三区四区| 日韩欧美免费精品| 婷婷亚洲欧美| a在线观看视频网站| 99视频精品全部免费 在线| 青草久久国产| 欧美一区二区精品小视频在线| 国产美女午夜福利| 日韩欧美精品v在线| 午夜福利成人在线免费观看| 国产私拍福利视频在线观看| 国产精品电影一区二区三区| 乱人视频在线观看| 99国产精品一区二区三区| 亚洲国产精品成人综合色| 日韩成人在线观看一区二区三区| 久久久久性生活片| 国产午夜精品久久久久久一区二区三区 | svipshipincom国产片| 午夜免费观看网址| 一级毛片女人18水好多| 90打野战视频偷拍视频| 国产免费男女视频| 天天躁日日操中文字幕| 夜夜看夜夜爽夜夜摸| 麻豆国产av国片精品| 大型黄色视频在线免费观看| 亚洲一区高清亚洲精品| 亚洲美女视频黄频| 亚洲av成人精品一区久久| 19禁男女啪啪无遮挡网站| 成人精品一区二区免费| 日本免费a在线| 婷婷亚洲欧美| 久久九九热精品免费| 免费在线观看影片大全网站| 禁无遮挡网站| 一区二区三区高清视频在线| 成人鲁丝片一二三区免费| 国产一区二区在线观看日韩 | 搡老岳熟女国产| 国产精品 欧美亚洲| 国产高清激情床上av| 国产 一区 欧美 日韩| 亚洲精品美女久久久久99蜜臀| 日日摸夜夜添夜夜添小说| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| 天天躁日日操中文字幕| 亚洲成人久久性| 99久久精品国产亚洲精品| 中文字幕人成人乱码亚洲影| 国产aⅴ精品一区二区三区波| www.999成人在线观看| 欧美高清成人免费视频www| 18禁美女被吸乳视频| 欧美最新免费一区二区三区 | 熟女人妻精品中文字幕| 搡老岳熟女国产| av黄色大香蕉| 欧美日韩中文字幕国产精品一区二区三区| 搞女人的毛片| 最后的刺客免费高清国语| 12—13女人毛片做爰片一| 日本黄色片子视频| 国产免费av片在线观看野外av| 国内精品久久久久精免费| 精品久久久久久久毛片微露脸| 欧美日韩瑟瑟在线播放| 亚洲五月婷婷丁香| 日韩欧美 国产精品| 真人一进一出gif抽搐免费| 亚洲电影在线观看av| 中文字幕精品亚洲无线码一区| 少妇人妻一区二区三区视频| 亚洲国产中文字幕在线视频| 亚洲熟妇中文字幕五十中出| 少妇的逼水好多| 身体一侧抽搐| 精品久久久久久久毛片微露脸| 一进一出抽搐动态| 国产精品永久免费网站| 午夜福利高清视频| 51午夜福利影视在线观看| 欧美+日韩+精品| 精品久久久久久久末码| 大型黄色视频在线免费观看| 精品久久久久久成人av| 亚洲精品美女久久久久99蜜臀| 久久久久免费精品人妻一区二区| 欧美成人a在线观看| 日本黄大片高清| 成年免费大片在线观看| 久久精品国产清高在天天线| 成人欧美大片| 99国产精品一区二区蜜桃av| 国产高潮美女av| 又爽又黄无遮挡网站| 在线免费观看的www视频| 日韩欧美在线乱码| 中国美女看黄片| 在线观看av片永久免费下载| 老司机福利观看| 久久精品影院6| 在线十欧美十亚洲十日本专区| 男插女下体视频免费在线播放| av天堂在线播放| 日本三级黄在线观看| 亚洲精品粉嫩美女一区| 久久久久国产精品人妻aⅴ院| 午夜精品在线福利| 搡老岳熟女国产| 国产色婷婷99| 亚洲七黄色美女视频| 成年女人看的毛片在线观看| 欧美黄色淫秽网站| 久久99热这里只有精品18| 丁香六月欧美| 夜夜爽天天搞| 老汉色∧v一级毛片| 欧美日韩一级在线毛片| 中文字幕av成人在线电影| 亚洲成人久久爱视频| 久久精品人妻少妇| 国产中年淑女户外野战色| 女人高潮潮喷娇喘18禁视频| 国产69精品久久久久777片| 美女免费视频网站| 欧美中文综合在线视频| 九九久久精品国产亚洲av麻豆| 亚洲国产中文字幕在线视频| 亚洲片人在线观看| 国产精品永久免费网站| 色综合亚洲欧美另类图片| 色尼玛亚洲综合影院| 国产蜜桃级精品一区二区三区| 在线播放国产精品三级| 他把我摸到了高潮在线观看| 国产综合懂色| 国产精品免费一区二区三区在线| 免费一级毛片在线播放高清视频| 99久久无色码亚洲精品果冻| 亚洲精品亚洲一区二区| 亚洲天堂国产精品一区在线| av在线天堂中文字幕| 国产精品久久久久久久电影 | 色噜噜av男人的天堂激情| 国产av麻豆久久久久久久| 男女视频在线观看网站免费| 久久精品人妻少妇| 午夜久久久久精精品| 乱人视频在线观看| 国产午夜福利久久久久久| 婷婷丁香在线五月| 变态另类丝袜制服| 亚洲精品一卡2卡三卡4卡5卡| 美女 人体艺术 gogo| 搡女人真爽免费视频火全软件 | 国产黄a三级三级三级人| 欧美不卡视频在线免费观看| 美女 人体艺术 gogo| 久久中文看片网| 亚洲欧美日韩东京热| 亚洲国产色片| 搡女人真爽免费视频火全软件 | 国产一区二区在线av高清观看| 欧美不卡视频在线免费观看| 51国产日韩欧美| 国产精品一及| 亚洲欧美日韩东京热| bbb黄色大片| 一进一出好大好爽视频| 一本综合久久免费| 久9热在线精品视频| 久久亚洲精品不卡| netflix在线观看网站| 久久九九热精品免费| 国产精品香港三级国产av潘金莲| 亚洲国产欧美网| 黄片大片在线免费观看| 91av网一区二区| 国产av麻豆久久久久久久| 最后的刺客免费高清国语| 97人妻精品一区二区三区麻豆| 国产精品电影一区二区三区| 18美女黄网站色大片免费观看| 国产野战对白在线观看| 啦啦啦观看免费观看视频高清| 国产精品久久久久久人妻精品电影| 12—13女人毛片做爰片一| 国产av一区在线观看免费| 观看美女的网站| 91久久精品电影网| 国产高清videossex| 亚洲欧美一区二区三区黑人| 大型黄色视频在线免费观看| 两个人看的免费小视频|