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

    Simulation of space heavy-ion induced primary knock-on atoms in bipolar devices

    2024-01-25 07:29:14BinZhang張彬HaoJiang姜昊XiaoDongXu徐曉東TaoYing應濤ZhongLiLiu劉中利WeiQiLi李偉奇JianQunYang楊劍群andXingJiLi李興冀
    Chinese Physics B 2024年1期
    關(guān)鍵詞:張彬

    Bin Zhang(張彬), Hao Jiang(姜昊), Xiao-Dong Xu(徐曉東), Tao Ying(應濤), Zhong-Li Liu(劉中利),Wei-Qi Li(李偉奇),3, Jian-Qun Yang(楊劍群), and Xing-Ji Li(李興冀),?

    1School of Physics,Harbin Institute of Technology,Harbin 150001,China

    2School of Material Science and Engineering,Harbin Institute of Technology,Harbin 150001,China

    3State Key Laboratory of Intense Pulsed Radiation Simulation and Effect,Xi’an 710024,China

    Keywords: Monte Carlo simulation,primary knock-on atom(PKA),space-heavy ion,radiation damage

    1.Introduction

    The influence of space radiation on electronic components is mainly due to the total ionizing dose, displacement damage and single-event effect.[1]The radiation effects of the space radiation environment on the electronic components inside a spacecraft have attracted extensive attention.[1,2]The above-mentioned radiation effects can lead to abnormality or failure of electronic devices and eventually the failure of the mission of the orbiting spacecraft.As silicon-based devices with good current drive capability,linearity,and other operating characteristics, bipolar junction transistors (BJTs) are often used in spacecraft.However, BJTs are sensitive to cosmic rays and may suffer displacement damage during on-orbit missions.[2–4]

    Displacement damage is caused by energetic particles colliding with target material atoms, resulting in the dislodging of atoms from their lattice position.[5]Displacement damage is hard to recover, which will seriously affect the electrical performance and lead to the degradation of semiconductor devices.[6,7]Primary knock-on atoms (PKAs), which are the first atoms to be displaced, usually have a higher energy level that can result in further cascade damage.Much work has been done on cascading damage in past studies, mainly on the microscopic mechanisms of displacement damage, to replace empirical methods such as non-ionizing energy loss(NIEL).[8–10]For example, using molecular dynamics (MD)simulation, Heet al.[6]studied the relationship of temperature with primary damage and found that temperature could increase the fraction of isolated vacancies in vacancy clusters because vacancy clusters may be thermally unstable.Hamedaniet al.[11]performed a comparative study of different materials to examine the damage production mechanisms for a relatively wide range of materials and recoil energies.The abovementioned MD simulation works were generally to select PKA incident materials with specific energy artificially to study the damage process of semiconductors.Researchers are trying to link the macroscopic electrical effects observed in the experiments with the microscopic damage structures generated at the atomic level,[12]and are attempting to combine different simulation techniques(e.g.,Monte Carlo,MD,ab initio).Obtaining a wealth of information on PKAs is crucial in this process.Researchers have extensively studied the PKAs generated by irradiated particles in materials.For example,Liuet al.[13]investigated radiation damage in a silicon drift detector (SDD)and found that the distribution of displacement damage caused by neutrons in SDDs is uniform.Shenet al.[14]studied the non-ionization energy loss of protons in silicon carbide and found that with an increase in incident proton energy the proportion of Si and C particles in the primary recoil atoms gradually decreased,from primary products to secondary products.

    Previous studies have primarily focused on radiation damage to materials by neutrons and protons.However,there are few reports on the PKAs induced by space heavy ions in semiconductors, despite heavy ions being common radiation particles in outer space and ground radiation sources.The simulation of PKAs induced by space heavy ions in semiconductors has important application scenarios and practical significance.In addition, the rich PKA information obtained by simulations can be used as important MD reference input parameters.Therefore, this study uses Monte Carlo simulation to study the PKAs induced by heavy ions in BJTs and qualitatively analyzes the radiation damage of typical space heavy ions in semiconductor materials.This research is of great significance as it provides a better understanding of the process of displacement damage caused by charged heavy ions in semiconductors.

    2.Models and physical processes

    2.1.Model setup

    In our previous research,[2]we conducted an experimental study on the NPN BJT 2N2222.The thickness of the base(p+) and epitaxial layer (n?) of the 2N2222 BJT is about 1.5 μm and 10 μm, respectively.In this simulation, we reduced the 2N2222 BJT to bulk Si cubes with dimensions of 10μm×10μm×10μm,focusing on the simulation of intrinsic Si,which is similar to a previous study.[12]In addition,to studying the distribution of PKAs in the BJT as a function of depth, the Si bulk material was divided into 100 parts vertically to count the number of PKAs in each layer.The model diagram is shown in Fig.1.

    Fig.1.The direction of ion incidence and the model of BJTs in this study.2N2222 BJTs were simplified as Si bulk material with 100 layers.The length and width of each layer is 10μm and the thickness is 0.1μm.

    2.2.Beam source energy setting

    Galactic cosmic rays are mainly composed of light ions(protons and He) and heavy ions.[7,15]Heavy ions are found to decrease considerably in relative abundance after Fe ions in space.Hence,we conducted our simulations with ions that have an atomic number lower than Fe.The ions used in the simulation were12C,14N,16O and56Fe.Energetic ions will dissipate energy due to electronic excitation and atomic displacement in host materials.Therefore,ion energy selection in simulations is mainly based on the ion range calculated using SRIM (stopping and range of ions in matter).[16]The calculated results for the ion ranges are shown in Fig.2.Previous studies have revealed that an increase in the energy of incident ions tends to weaken the displacement damage dose at the corresponding location within a device.[2]Therefore,in the simulation we only consider ion incidence within a certain energy range.In other words,we take the total thickness of the model BJTs (10 μm) as the maximum range of the ion and choose an appropriate energy as the ion energy used in our simulation.The selected ion kinetic energy parameters are shown in Table 1.

    Fig.2.The range of different types of ions (C, N, O and Fe) in Si bulk material.The horizontal axis represents the incident energies of ions.The red dotted line represents the total thickness of the BJTs.

    Table 1.The types of ions used in the simulation and the ion energy chosen for this simulation

    2.3.Simulation parameters

    Our research work uses one module of the software Extreme-environment Radiation Effect Technology Computer Aided Design(ERETCAD)developed by the present authors,which is mainly based on GEANT4 (GEometry ANd Tracking).GEANT4[17]is an open-source simulation toolkit that offers a wide set of physics models;it has a good performance in particle transport and is also used in the field of semiconductor irradiation and medical and space applications.[18,19]

    The calculation results of TRIM(transport of ions in matter) were used as the benchmark in this study.[16]TRIM is a computational subroutine of SRIM that employs Monte Carlo simulations to track the ion trajectory and determine the energy transferred to each target atom collision.In the following discussions we will refer to both SRIM and TRIM interchangeably but only use the name SRIM.After considering the convergence of the simulations and the storage of the data,we conducted simulations up to 100000 events in our code and up to 300 events using SRIM.The number of events remained unchanged for different types of ions.

    A proper physical list is crucial to obtaining accurate results of the simulation.We use a physical list called standardNR to simulate the transport of heavy ions in BJTs and the production of PKAs.A single Coulomb scattering process is included in standardNR.This physical list also explicitly integrates the classical equations of motion for scattering events,resulting in precise tracking of both the projectile and the recoil target nucleus.[20]In addition,accurate description of the transport process of heavy ions requires not only electromagnetic processes but also physical processes such as elastic and inelastic heavy ions.The physical process of the developed test application is based on the GEANT4 TestEm7 example.The NeutronHPphysics list is used to complete some neutron collision simulation studies.For the simulation of neutrons in materials, we refer to Raineet al.’s research work.[12]An option called‘detailed calculation with full damage cascades’was used in this simulation to achieve the accuracy of SRIM calculation results.This option follows every recoil until its energy drops below the lowest displacement energy of any target atom.Hence all collisional damage to the target is analyzed.For the extraction of PKA information in SRIM in this study we mainly refer to the research of Liuet al.[21]

    3.Results and analysis

    3.1.Energy spectrum of PKAs

    The energy distribution of PKAs can provide an important parameter for semiconductor defect evolution.The first parameter extracted from simulations was the energy spectrum of the PKAs.The results of calculation from SRIM are used as a reference because of the wide acceptance of SRIM in many fields.As shown in Fig.3,all the results are normalized to the total number of events used in this simulation.

    Fig.3.Energy spectra of PKAs induced by different heavy ions in BJTs: (a)C,(b)N,(c)O,(d)Fe.The simulation results of SRIM are shown with purple diamonds.The results for 0.7 MeV protons are shown in panel (a).In the inset of each figure, we present the details of the PKA energy spectrum around the peak position.

    As shown in Fig.3, it is found that PKAs induced by space ions of different types and energies have a similar energy spectral distribution trend.The energy spectrum of PKAs reaches its peak at 17 eV and drops rapidly thereafter.The variation range of the energy spectrum is mainly 20 eV–100 eV and the trend of each energy spectrum is consistent after 1000 eV.PKAs have different cut-off energies according to different incident energies and types of incident ions.The process of producing PKAs is an elastic collision process and the effects of electron excitation can be basically ignored.This process can be described by a binary collision approximation based on the conservation of momentum and energy[22]

    whereMnandMare the masses of the incident ion and the target atom,Enis the energy of the incident ion andθis the scattering angle of the recoil atom.As shown in Eq.(1), the energy attained by PKAs is affected by many factors,which is the reason why there is a maximum cut-off energy.

    The amplitude of the energy spectrum is greatly affected by the incident energy and ion type.As shown in Fig.3, we noted that they=axbfunction can better describe the change in the energy spectrum in the low-energy region (17 eV–1000 eV).In this paper,we use fitting parameters to intuitively analyze the impact of different types and energies of incident ions on the energy spectrum distribution of PKAs.To obtain a complete calculation and results analysis,ions with higher energy have also been considered in this study;the range of these ions exceeds the thickness of the BJTs(10μm)studied in this paper.The functional relationship between fitting parameters and incident ion energy is shown in Fig.4.

    Fig.4.The relationship between the fitting parameters obtained by fitting the PKA energy spectrum with functions and the incident ions energy (in the PKA energy range of 17 eV–1000 eV).(a) The fitting parameter a.(b)The fitting parameter b as a function of the initial energy of the incident ions.

    Coefficientadetermines the amplitude height of the function and coefficientbdetermines the steepness of the function.As shown in Fig.4(a),at the same incident energy coefficientais positively correlated with the atomic number(the mass of the incident ion),which means that the amplitude height of the PKA energy spectrum increases.The increase in atomic number will directly lead to an increase in the number of PKAs.This is mainly because the elastic scattering cross section of ions in a material is related to the ion mass.[7]When the range of ions is close to the size of the BJT,and reaches its maximum value,the peak height of the PKA energy spectrum reaches its maximum.However,when the ion range is larger than the size of the BJT, the number of PKAs begins to reduce.The main reason for this phenomenon is that the energy obtained by the target atom is negatively correlated with the velocity or energy of the incident ions.[23]The increase in the incident velocity of ions means that the reaction time between ions and the target decreases.Thus,the number of PKAs induced by high-energy ion incidence decreases.As shown in Fig.4(b),the value ofbis less affected by the type of incident ion and its energy; the value ofbis maintained at a value of around?1.5.The magnitude ofbalso determines the profile and trend of the energy spectrum, therefore the PKAs produced in BJTs by different types of ions have a similar energy spectrum distribution.

    The energy spectrum of space ion-induced PKAs is mainly in the low-energy region, and PKAs with low energy will not produce cascading damage.Therefore, we believe that the characteristics of PKAs can be used to quickly and qualitatively analyze the displacement damage of space ions in BJTs.Previous MD studies have shown that larger atomicscale structural damage is often caused by PKAs with greater energy.[6,24]In our simulation,information about high-energy PKAs was also recorded, but the probability of their generation is relatively small.Therefore, we believe that through the local simulation of these PKAs with fewer occurrences but higher energy we can complete a more complete theoretical analysis of displacement damage to semiconductor materials or devices.Finally, as shown in Fig.3, the result of the simulation is close to those of SRIM,so we believe that the code of the simulation in this study is effective; this allows us to further investigate the nature of PKAs induced by space ions in BJTs.We attribute the variation seen between our code and SRIM in the incident scenarios of 13 MeV N ions in Fig.3(b)and 37 MeV Fe ions in Fig.3(c)to the distinct principles of the two calculation methods: SRIM employs the binary collision approximation method while GEANT4 is more reliant on the collision cross section.

    3.2.PKA angle distribution

    Another important characteristic of PKAs, the angle distribution, is also worthy of attention.In previous MD studies[6,11]the PKAs were mostly chosen to fall vertically.However, previous reported studies have also shown that the incidence direction and distance of PKA may have a certain impact on the damage yield.[22,25]Therefore,the angle distribution of PKAs induced by heavy ions in BJTs was also studied in this paper.The scattering angle of PKAs was defined to be the angle between the velocity direction of the PKA and the initial direction of the incident ion.The results of the simulation are shown in Fig.5.

    The results of the simulation demonstrated that the velocity distribution of PKAs is from 0°to 180°.This is mainly because incident ions are kinds of charged particles that will be deflected by the Coulomb interaction of the target atoms during the ion transport process.Therefore,the PKAs induced by ions have the entire spatial distribution.The PKAs induced by typical heavy ions and protons in BJTs (Fig.5(a))had a distinct single-peak feature with a peak position around 88°and were not affected by ion type and energy.This also means that within a certain range of incident energies, the PKAs generated by the collision of heavy ions have a larger forward scattering angle.In addition,our proton irradiation results are very close to previous research conclusions.Whitlow and Nakagawa[26]calculated the angle distribution of PKAs produced by 1 MeV protons incident in Si, and their results showed that the peak position of the PKAs was around 89°.According to the binary collision theory(Eq.(1)),in the case of a larger scattering angle the energy obtained by PKAs is smaller,which is consistent with the energy spectrum distribution of PKAs obtained in our simulation.In addition,for this single-peak feature of PKA angle distribution, we have also done further research to explore whether other particle types can produce similar phenomena in BJTs.Previous studies[22]have shown that neutrons,as a typical neutral particle,can produce large displacement damage.Therefore, we investigated the neutron-induced angular distribution of PKAs in BJT.In our neutron simulation, the energy of incident neutrons can cover all physical reactions between neutrons and matter,and the simulation results are shown in Fig.6.

    Fig.5.Scattering angle distribution of PKAs induced by different heavy ions in BJTs: (a) C, (b) N, (c) O, (d) Fe.In (a) the purple symbols represent the PKA angle distribution induced by proton incidence.The red dotted line indicates that the scattering angle is 90°.

    It is important to note that we did not apply the same normalization to the neutron calculations when comparing the properties of neutron and space ion-induced PKAs.This is mainly because the charged particles and neutrons use different physical models in the simulation so they can only be qualitatively compared.As shown in Fig.6, the angular distribution of PKAs induced by neutrons can also be distributed in the entire spatial angle range.When low-energy neutrons are incident(1 MeV),the simulation results show that the angular distribution of PKAs has a wide range of peaks,and the peak position is about 64°.When neutrons with an energy of 15 MeV are incident,three angular peaks appear.The number of peaks is affected by the neutron energy.A similar phenomenon also appeared in the simulation by Caiet al.[22]when simulating neutron-induced PKAs in Zr-based alloys.The main reason for this phenomenon is that the interaction between neutrons and matter is accompanied by a nuclear reaction, resulting in a variety of nuclear fragments.For charged ions, the transport process of 1 MeV protons in SiC materials predicted by Shenet al.[14]showed that PKAs generated from nuclear reactions accounted for a small proportion of all PKAs.The fragmentation cross section of heavy ions is affected by the energy of the incident ions.In this study, the simulated energy of incident ions is small, and the PKAs generated come mostly from the elastic scattering process.The contribution of nuclear reactions can be ignored.Therefore,there is no multipeak phenomenon in the angular distribution of PKAs induced by heavy ions.

    Fig.6.Scattering angle distribution of PKAs induced by neutron in BJTs.A large neutron energy range is used for incidence on the BJTs,covering most of the reaction processes of neutrons and materials.

    From the simulation results, it can be seen that both charged ions and neutral particles have similar scattering angle distributions in Si-based BJTs,and the PKAs have a large forward scattering angle.Therefore,we believe that when using MD to study the evolution of defects induced by heavy ions in BJTs we can choose a certain angle for the PKA incidence direction instead of just vertical incidence.In addition,it can be seen from Fig.6 that neutron-induced PKAs in BJTs also have a large scattering angle (61°–81°).The main reason for this may be that the probability of head to head collision of particles with target atoms in the process of transmission of particles in the material is small.Charged ions are also affected by the Coulomb force, so the PKA angle induced by ions is larger.Neutral particles can produce more forward scattering.

    3.3.The depth distribution of PKAs

    The ranges of space ions with different energies vary in materials,resulting in displacement damage in varied areas of BJTs by charged ions.The localized defect evolution of BJTs can help us gain invaluable insights into the defect distribution in real irradiation environments.If the relationship between ion energy and PKA distribution in the device is established it can be used to guide ground experiments to study the effect of displacement damage in specific regions of the BJT on its electrical performance.The simulation results of PKA depth–count curves in a BJT chip are plotted in Fig.7.

    Fig.7.PKA depth–count curves induced by different heavy ions in BJTs: (a)C,(b)N,(c)O,(d)Fe.The pink diamonds represent the PKAs distribution of 0.7 MeV protons in BJTs.

    From the results of the simulation,it can be seen that the distribution of PKAs along the direction of the incident ions exhibits an obvious peak characteristic.The peak initially increases with increase in depth, and then decreases rapidly after reaching the peak.The displacement damage area is not evenly distributed.More displacement damage occurs near the peak position.This is because the energy loss of charged particles is characterized primarily by the square of the nuclear charge and the inverse square of the projectile velocity.Space ions are slower at the end of their range and therefore produce more PKAs.Similarly,this phenomenon also occurs at the incident position, and the number of PKAs decreases with increase in the incident particle energy.At the end of the ion range,the ion energy is not sufficient to induce PKAs,so the number of PKAs decreases rapidly.

    The relationship between the peak position of PKAs and the incident ion energy is shown in Fig.8.It can be seen from the calculation results that the depth position of the PKA peak is a linear function of the incident ion energy for the same incident ion.However,it should be noted that the peak position should occur within a certain energy range, that is, the range of ions is smaller than the thickness of the BJT.The peak position is affected by many factors, mainly the atomic number and the energy of the incident ion.At the same ion incident energy,the atomic number has a negative correlation with the peak of the incident depth.

    Fig.8.PKA depth–count curves induced by different heavy ions in BJTs:(a)C,(b)N,(c)O,(d)Fe.

    4.Conclusion

    Monte Carlo simulation was used to compute the displacement damage induced by space ions, focusing on the properties of primary knock-on atoms (PKAs).Our results demonstrate that the properties of PKAs are mainly affected by the type and energy of incident ions.The energy spectrum of PKAs is mainly concentrated in the low-energy region,which is insufficient to cause further cascade damage.Hence the PKAs generated by ions can effectively qualitatively describe the displacement damage in a BJT.Within a certain range of ion incident energies, the distribution of ion-induced PKAs in BJT is inhomogeneous,with an objective peak distribution trend.Our calculation results show that the peak position is a function of the incident ion and energy, which allows us to study the influence of displacement damage in specific areas of BJTs on their electrical performance by adjusting the type of incident ion and incident energy.Additionally,the research results reveal that PKA chooses a certain angle of incidence and localized defect evolution to confirm the displacement damage in each area of BJT can be more consistent with the real radiation damage scenario.In summary,our work supports the investigation of displacement damage induced by space ions in BJTs and the simulation of defect evolution in semiconductor materials.

    Acknowledgments

    Project supported by the National Natural Science Foundation of China(Grant Nos.11974091,51973046,U22B2044,and 21673025) and the Open Projects of State Key Laboratory of Intense Pulsed Radiation Simulation and Effect(Grant No.SKLIPR2020).

    猜你喜歡
    張彬
    Compared discharge characteristics and film modifications of atmospheric pressure plasma jets with two different electrode geometries
    復雜光照環(huán)境的車道線檢測方法
    EXISTENCE RESULTS FOR SINGULAR FRACTIONAL p-KIRCHHOFF PROBLEMS*
    女醫(yī)生喚醒沉睡愛人:攀山渡河筑一座愛城
    電力拖動實訓接線板的設(shè)計與研究
    一種相控陣天線波束指向角計算方法
    GEANT4 simulation study of over-response phenomenon of fiber x-ray sensor?
    酒鬼報警
    民間文學(2020年5期)2020-07-04 02:37:19
    患癌一瞞到底:我們笑著我們含淚炮制深情
    孤獨舞者情殤7年,“天邊的愛”就在盈盈眉眼處
    免费在线观看亚洲国产| 男人舔女人下体高潮全视频| 免费大片18禁| 久久99热6这里只有精品| 脱女人内裤的视频| 在线观看免费视频日本深夜| 我的女老师完整版在线观看| 久久午夜福利片| 婷婷丁香在线五月| 少妇熟女aⅴ在线视频| 午夜福利高清视频| 欧美性猛交黑人性爽| 日韩亚洲欧美综合| 国产成人a区在线观看| 国产精品98久久久久久宅男小说| 久久婷婷人人爽人人干人人爱| 最近中文字幕高清免费大全6 | 老司机午夜十八禁免费视频| 日日摸夜夜添夜夜添av毛片 | 欧美精品国产亚洲| 国模一区二区三区四区视频| 真实男女啪啪啪动态图| 国产一级毛片七仙女欲春2| 久久久成人免费电影| 欧美zozozo另类| 久久久国产成人精品二区| 国产伦精品一区二区三区四那| 国产精品一区二区性色av| 级片在线观看| 日本成人三级电影网站| 美女高潮喷水抽搐中文字幕| 欧美成人一区二区免费高清观看| 欧美国产日韩亚洲一区| 亚洲经典国产精华液单 | 久久人妻av系列| av福利片在线观看| 国产精品一区二区三区四区久久| 久99久视频精品免费| 国产色婷婷99| 亚洲精品亚洲一区二区| 久久久久久大精品| 欧美丝袜亚洲另类 | av天堂在线播放| 国产一区二区激情短视频| 亚洲国产色片| 18+在线观看网站| 嫩草影视91久久| 国产精品三级大全| 无人区码免费观看不卡| 麻豆一二三区av精品| 亚洲专区中文字幕在线| 在线观看美女被高潮喷水网站 | 成人av一区二区三区在线看| 亚洲av电影在线进入| 无人区码免费观看不卡| 日本熟妇午夜| 色综合亚洲欧美另类图片| 中亚洲国语对白在线视频| 日本三级黄在线观看| 国语自产精品视频在线第100页| 美女免费视频网站| 一个人观看的视频www高清免费观看| 精品日产1卡2卡| 在线十欧美十亚洲十日本专区| 久久国产乱子伦精品免费另类| 国内毛片毛片毛片毛片毛片| 国产精品久久久久久久久免 | 国产精品自产拍在线观看55亚洲| 男人狂女人下面高潮的视频| 国语自产精品视频在线第100页| 黄色女人牲交| 麻豆成人午夜福利视频| 色5月婷婷丁香| 国产男靠女视频免费网站| av在线天堂中文字幕| 又爽又黄无遮挡网站| 国产伦一二天堂av在线观看| 少妇熟女aⅴ在线视频| 91字幕亚洲| 精品人妻偷拍中文字幕| 欧美色欧美亚洲另类二区| 真人做人爱边吃奶动态| 色哟哟哟哟哟哟| 国产私拍福利视频在线观看| а√天堂www在线а√下载| 精品人妻熟女av久视频| 国产大屁股一区二区在线视频| 久久久久免费精品人妻一区二区| 日韩人妻高清精品专区| 亚洲国产欧洲综合997久久,| 午夜老司机福利剧场| 精品人妻一区二区三区麻豆 | 九色成人免费人妻av| 亚洲五月天丁香| 国产美女午夜福利| 国产精品伦人一区二区| av黄色大香蕉| 小蜜桃在线观看免费完整版高清| 国产精品乱码一区二三区的特点| 成年女人看的毛片在线观看| 久久久久亚洲av毛片大全| 热99在线观看视频| 欧美精品国产亚洲| 亚洲av熟女| 亚洲欧美日韩高清专用| 中出人妻视频一区二区| av专区在线播放| 欧美区成人在线视频| 村上凉子中文字幕在线| 99国产综合亚洲精品| 婷婷丁香在线五月| 悠悠久久av| 成人特级av手机在线观看| 麻豆av噜噜一区二区三区| 欧美精品啪啪一区二区三区| 色吧在线观看| 免费在线观看亚洲国产| 麻豆成人av在线观看| 亚洲av美国av| 日本黄色片子视频| 亚洲av五月六月丁香网| 亚洲av美国av| 免费看日本二区| 亚洲av五月六月丁香网| 国产一区二区亚洲精品在线观看| 如何舔出高潮| 悠悠久久av| 极品教师在线免费播放| 最近最新中文字幕大全电影3| 九九在线视频观看精品| 国产一区二区亚洲精品在线观看| 亚洲不卡免费看| 精品国产亚洲在线| 搡老妇女老女人老熟妇| 99热这里只有是精品50| 深夜a级毛片| 夜夜爽天天搞| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av香蕉五月| 成年人黄色毛片网站| 午夜日韩欧美国产| 亚洲综合色惰| 欧美xxxx性猛交bbbb| 成人永久免费在线观看视频| 美女高潮喷水抽搐中文字幕| 极品教师在线视频| 亚洲中文字幕日韩| 少妇裸体淫交视频免费看高清| 欧美国产日韩亚洲一区| 麻豆久久精品国产亚洲av| 黄色视频,在线免费观看| 99热这里只有精品一区| 波多野结衣高清无吗| 51午夜福利影视在线观看| 国产精品永久免费网站| 午夜福利在线在线| 色播亚洲综合网| 男女那种视频在线观看| 午夜影院日韩av| 国产精品伦人一区二区| 999久久久精品免费观看国产| 国产亚洲欧美在线一区二区| 国产精品永久免费网站| 日韩人妻高清精品专区| 97热精品久久久久久| 免费看日本二区| 国产av麻豆久久久久久久| 深夜a级毛片| a级毛片免费高清观看在线播放| 国产人妻一区二区三区在| 国产精品亚洲美女久久久| 婷婷六月久久综合丁香| 中文字幕精品亚洲无线码一区| 少妇高潮的动态图| 欧美黄色片欧美黄色片| 欧美中文日本在线观看视频| 变态另类成人亚洲欧美熟女| 国产精品女同一区二区软件 | 小说图片视频综合网站| 乱人视频在线观看| 九色国产91popny在线| 婷婷精品国产亚洲av在线| 亚洲aⅴ乱码一区二区在线播放| 又黄又爽又刺激的免费视频.| 亚洲精品亚洲一区二区| 夜夜爽天天搞| 搡老熟女国产l中国老女人| 国产欧美日韩一区二区三| 无人区码免费观看不卡| 一区二区三区激情视频| 免费观看人在逋| 又爽又黄a免费视频| 久久精品人妻少妇| 一边摸一边抽搐一进一小说| 欧美+日韩+精品| 我的女老师完整版在线观看| 特大巨黑吊av在线直播| 国产精品一区二区性色av| 精品免费久久久久久久清纯| 啦啦啦韩国在线观看视频| 少妇人妻一区二区三区视频| 国产精品美女特级片免费视频播放器| 九九在线视频观看精品| 超碰av人人做人人爽久久| 亚洲天堂国产精品一区在线| 五月伊人婷婷丁香| 99热6这里只有精品| 三级男女做爰猛烈吃奶摸视频| 欧美日韩中文字幕国产精品一区二区三区| 中文在线观看免费www的网站| 国产精品久久久久久亚洲av鲁大| 国产精华一区二区三区| 伦理电影大哥的女人| 有码 亚洲区| 亚洲av电影不卡..在线观看| 午夜免费成人在线视频| 国产欧美日韩一区二区精品| 亚洲精品一卡2卡三卡4卡5卡| 亚洲无线在线观看| 亚洲av成人精品一区久久| 少妇人妻一区二区三区视频| 非洲黑人性xxxx精品又粗又长| 一级黄色大片毛片| 日本免费a在线| 久久久久久久亚洲中文字幕 | 亚洲人成网站高清观看| 婷婷六月久久综合丁香| 99久久精品国产亚洲精品| 欧美日韩综合久久久久久 | 18美女黄网站色大片免费观看| 免费av观看视频| 国产精品久久久久久精品电影| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 一个人免费在线观看电影| 日韩大尺度精品在线看网址| 国产免费av片在线观看野外av| 国产日本99.免费观看| 国产精品98久久久久久宅男小说| 国产精品三级大全| 国产精品永久免费网站| 色噜噜av男人的天堂激情| 欧美xxxx性猛交bbbb| 国产老妇女一区| 免费大片18禁| 成人精品一区二区免费| 一二三四社区在线视频社区8| 精品99又大又爽又粗少妇毛片 | 一个人免费在线观看的高清视频| 免费观看精品视频网站| 高清日韩中文字幕在线| 亚洲avbb在线观看| 无遮挡黄片免费观看| 久久久色成人| 动漫黄色视频在线观看| 国模一区二区三区四区视频| 伦理电影大哥的女人| 一个人免费在线观看的高清视频| 丁香六月欧美| 亚洲熟妇熟女久久| 国产精品一区二区免费欧美| 午夜a级毛片| 精品无人区乱码1区二区| 动漫黄色视频在线观看| 18禁裸乳无遮挡免费网站照片| 国产伦在线观看视频一区| 少妇的逼好多水| 99久久成人亚洲精品观看| 久久国产乱子免费精品| 亚洲一区高清亚洲精品| 亚洲av成人不卡在线观看播放网| 成人av在线播放网站| 老熟妇乱子伦视频在线观看| 亚洲最大成人av| 淫妇啪啪啪对白视频| avwww免费| 欧美色欧美亚洲另类二区| 悠悠久久av| 日韩欧美精品免费久久 | 久久久久国内视频| 网址你懂的国产日韩在线| 偷拍熟女少妇极品色| av专区在线播放| 日韩欧美三级三区| 国产成人啪精品午夜网站| av天堂在线播放| 一级作爱视频免费观看| 身体一侧抽搐| 日韩欧美精品免费久久 | 久久久色成人| 不卡一级毛片| 一本一本综合久久| 欧美最新免费一区二区三区 | 18+在线观看网站| 免费高清视频大片| 精品久久久久久久久久免费视频| 久久婷婷人人爽人人干人人爱| 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| av女优亚洲男人天堂| 亚洲中文字幕日韩| 亚洲欧美清纯卡通| 午夜两性在线视频| 美女cb高潮喷水在线观看| 免费无遮挡裸体视频| 深爱激情五月婷婷| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 久久国产精品人妻蜜桃| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 搡女人真爽免费视频火全软件 | 狠狠狠狠99中文字幕| 18+在线观看网站| 国产一区二区在线av高清观看| 国产在视频线在精品| 男女下面进入的视频免费午夜| 夜夜夜夜夜久久久久| 欧美一区二区精品小视频在线| 欧美乱妇无乱码| 日本成人三级电影网站| 99riav亚洲国产免费| 欧美激情久久久久久爽电影| 69av精品久久久久久| 欧美激情国产日韩精品一区| 男女之事视频高清在线观看| 欧美绝顶高潮抽搐喷水| 美女xxoo啪啪120秒动态图 | 91久久精品国产一区二区成人| 国产熟女xx| 国产欧美日韩精品亚洲av| 国产 一区 欧美 日韩| 热99在线观看视频| 99久久成人亚洲精品观看| 国内揄拍国产精品人妻在线| 别揉我奶头 嗯啊视频| avwww免费| 亚洲av电影不卡..在线观看| 久久精品夜夜夜夜夜久久蜜豆| 又爽又黄无遮挡网站| 免费看a级黄色片| 非洲黑人性xxxx精品又粗又长| 日本免费一区二区三区高清不卡| 精品熟女少妇八av免费久了| 国产爱豆传媒在线观看| 午夜激情福利司机影院| 好男人电影高清在线观看| 精品一区二区三区视频在线| 国产高潮美女av| 欧美成狂野欧美在线观看| 亚洲人成网站在线播放欧美日韩| 嫁个100分男人电影在线观看| 女同久久另类99精品国产91| 国产精华一区二区三区| 免费在线观看影片大全网站| 好男人在线观看高清免费视频| 亚洲精品在线美女| 国产真实乱freesex| 每晚都被弄得嗷嗷叫到高潮| 亚洲人与动物交配视频| 亚洲真实伦在线观看| 99久久精品国产亚洲精品| 精品人妻1区二区| 免费在线观看亚洲国产| 亚洲欧美日韩无卡精品| 欧美日韩亚洲国产一区二区在线观看| 国产精品av视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 中文字幕人成人乱码亚洲影| 精品99又大又爽又粗少妇毛片 | 欧美xxxx黑人xx丫x性爽| 免费av不卡在线播放| 亚洲精品影视一区二区三区av| 国产精品98久久久久久宅男小说| 俺也久久电影网| 国产精品久久久久久久电影| 麻豆国产av国片精品| 日韩欧美 国产精品| 欧美日韩福利视频一区二区| 嫩草影视91久久| 亚洲人成伊人成综合网2020| 国产真实伦视频高清在线观看 | 在线十欧美十亚洲十日本专区| 两个人视频免费观看高清| 搞女人的毛片| 毛片一级片免费看久久久久 | 又粗又爽又猛毛片免费看| 国产精品美女特级片免费视频播放器| 日本五十路高清| 97超级碰碰碰精品色视频在线观看| 欧美日本亚洲视频在线播放| 精品日产1卡2卡| 日本 av在线| 亚洲精品乱码久久久v下载方式| 亚洲美女视频黄频| 九色国产91popny在线| 88av欧美| 赤兔流量卡办理| 亚洲性夜色夜夜综合| 久久久久久久久中文| 99国产综合亚洲精品| 欧美绝顶高潮抽搐喷水| 日本 欧美在线| 国产视频一区二区在线看| 免费观看的影片在线观看| 日本黄大片高清| 亚洲国产色片| 人人妻人人澡欧美一区二区| 色av中文字幕| 国产精品,欧美在线| 国产成人a区在线观看| 哪里可以看免费的av片| 色综合婷婷激情| 在线观看美女被高潮喷水网站 | 俄罗斯特黄特色一大片| 人妻制服诱惑在线中文字幕| 亚洲精品久久国产高清桃花| 日本黄大片高清| 日韩 亚洲 欧美在线| 人人妻人人澡欧美一区二区| 99热这里只有精品一区| 久久久久久久久久黄片| 少妇的逼水好多| 制服丝袜大香蕉在线| 久久久久久久久久成人| 中亚洲国语对白在线视频| 久久精品91蜜桃| 国产淫片久久久久久久久 | 久久久色成人| 国产av在哪里看| 精品午夜福利视频在线观看一区| 九九久久精品国产亚洲av麻豆| 尤物成人国产欧美一区二区三区| 美女免费视频网站| 国产成人av教育| 国产一级毛片七仙女欲春2| 婷婷六月久久综合丁香| 亚洲人成电影免费在线| 国产精品98久久久久久宅男小说| av在线天堂中文字幕| 精品国内亚洲2022精品成人| 少妇丰满av| 热99re8久久精品国产| 国产黄片美女视频| 一级毛片久久久久久久久女| 99久久精品一区二区三区| 免费av毛片视频| 午夜激情福利司机影院| 色哟哟哟哟哟哟| 一边摸一边抽搐一进一小说| 亚洲成人免费电影在线观看| 无人区码免费观看不卡| 午夜福利欧美成人| 精品久久久久久久久久免费视频| 久久久精品欧美日韩精品| 国内精品一区二区在线观看| 男人的好看免费观看在线视频| 久久久久久九九精品二区国产| 日本黄色片子视频| 亚洲国产精品sss在线观看| 99久久精品国产亚洲精品| 中文字幕av成人在线电影| 久久久色成人| 久久精品91蜜桃| 十八禁人妻一区二区| 婷婷精品国产亚洲av| 亚洲黑人精品在线| 欧美不卡视频在线免费观看| 性色av乱码一区二区三区2| 国产黄片美女视频| 国产真实乱freesex| 国产午夜精品久久久久久一区二区三区 | 日韩欧美国产一区二区入口| 中文字幕免费在线视频6| 国产成年人精品一区二区| 在线国产一区二区在线| 男女视频在线观看网站免费| 日韩欧美 国产精品| 又爽又黄无遮挡网站| 一进一出好大好爽视频| 中文字幕高清在线视频| 久久6这里有精品| 日本黄色视频三级网站网址| 美女 人体艺术 gogo| 色噜噜av男人的天堂激情| 久久久久久久亚洲中文字幕 | 国产精品国产高清国产av| 国产亚洲精品久久久com| 91字幕亚洲| 欧美另类亚洲清纯唯美| 五月玫瑰六月丁香| 特级一级黄色大片| 国产精华一区二区三区| 欧美日本亚洲视频在线播放| 亚洲天堂国产精品一区在线| or卡值多少钱| 变态另类成人亚洲欧美熟女| 无人区码免费观看不卡| 国产精品av视频在线免费观看| 99精品在免费线老司机午夜| а√天堂www在线а√下载| 亚洲国产高清在线一区二区三| 免费在线观看亚洲国产| 日韩欧美一区二区三区在线观看| 赤兔流量卡办理| 一个人观看的视频www高清免费观看| 精品一区二区三区人妻视频| 亚洲第一欧美日韩一区二区三区| 五月玫瑰六月丁香| 国产探花极品一区二区| 久久久色成人| 五月伊人婷婷丁香| 欧美最黄视频在线播放免费| 亚洲av熟女| 免费在线观看成人毛片| 黄色女人牲交| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 免费在线观看影片大全网站| 色综合欧美亚洲国产小说| 十八禁网站免费在线| 欧美日韩国产亚洲二区| 禁无遮挡网站| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 在线观看美女被高潮喷水网站 | 免费人成视频x8x8入口观看| 两个人的视频大全免费| 色综合欧美亚洲国产小说| 国产av不卡久久| 国产熟女xx| 嫁个100分男人电影在线观看| 免费看a级黄色片| 国产又黄又爽又无遮挡在线| 超碰av人人做人人爽久久| 亚洲第一区二区三区不卡| 国产单亲对白刺激| 久久人人爽人人爽人人片va | 18禁在线播放成人免费| 久久香蕉精品热| 日本黄色片子视频| 深夜精品福利| 两性午夜刺激爽爽歪歪视频在线观看| 青草久久国产| 日韩成人在线观看一区二区三区| 欧美一级a爱片免费观看看| 一个人免费在线观看电影| 身体一侧抽搐| 欧美另类亚洲清纯唯美| 天堂动漫精品| 成人国产综合亚洲| 三级毛片av免费| 国产欧美日韩一区二区三| 国产熟女xx| 欧美成人一区二区免费高清观看| 最好的美女福利视频网| 成人美女网站在线观看视频| 成人午夜高清在线视频| 中文字幕av成人在线电影| 高清在线国产一区| 深爱激情五月婷婷| 成年女人看的毛片在线观看| 精品午夜福利视频在线观看一区| 有码 亚洲区| 久久国产精品影院| 国产亚洲欧美98| 国产精品一区二区三区四区久久| 国产毛片a区久久久久| 成人av在线播放网站| 九九久久精品国产亚洲av麻豆| 91在线观看av| 亚洲av免费高清在线观看| www日本黄色视频网| 亚洲av中文字字幕乱码综合| 午夜老司机福利剧场| 国产真实伦视频高清在线观看 | 91午夜精品亚洲一区二区三区 | 亚洲 欧美 日韩 在线 免费| 麻豆av噜噜一区二区三区| 亚洲国产欧洲综合997久久,| 99久久九九国产精品国产免费| 午夜激情福利司机影院| 99热只有精品国产| 欧美激情久久久久久爽电影| 美女 人体艺术 gogo| 在线免费观看的www视频| 久久久久久久久久成人| 国产人妻一区二区三区在| 精品一区二区免费观看| 熟妇人妻久久中文字幕3abv| 亚洲人成电影免费在线| 久久久久久久久大av| 欧美绝顶高潮抽搐喷水| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 国产成+人综合+亚洲专区| 亚洲国产色片| 丰满的人妻完整版| 国产成+人综合+亚洲专区| 成人性生交大片免费视频hd| 九九久久精品国产亚洲av麻豆| 国产成+人综合+亚洲专区| 精品熟女少妇八av免费久了| 亚洲av.av天堂| 人妻制服诱惑在线中文字幕| 在线免费观看不下载黄p国产 | 国产人妻一区二区三区在| 国产欧美日韩一区二区三| 少妇熟女aⅴ在线视频| 99热6这里只有精品| 国产高清视频在线观看网站| 亚洲经典国产精华液单 | 亚洲,欧美精品.| 国产精品一区二区三区四区免费观看 | 色哟哟·www|