• <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年,“天邊的愛”就在盈盈眉眼處
    久久久久久九九精品二区国产| 日本一本二区三区精品| 97超视频在线观看视频| 久久鲁丝午夜福利片| 亚洲成a人片在线一区二区| a级毛片免费高清观看在线播放| 亚洲成人久久性| 日韩欧美精品v在线| av专区在线播放| 亚洲欧洲国产日韩| 免费看光身美女| 中文字幕熟女人妻在线| 九九在线视频观看精品| 免费人成在线观看视频色| 最近2019中文字幕mv第一页| 午夜亚洲福利在线播放| 精品熟女少妇av免费看| 男女那种视频在线观看| 久久韩国三级中文字幕| 九色成人免费人妻av| 91久久精品电影网| 看片在线看免费视频| 狂野欧美白嫩少妇大欣赏| 日日摸夜夜添夜夜爱| 亚洲精品成人久久久久久| 日本熟妇午夜| 精品人妻一区二区三区麻豆| 一级毛片电影观看 | 欧美高清成人免费视频www| 日韩 亚洲 欧美在线| 99九九线精品视频在线观看视频| 午夜a级毛片| 国产伦在线观看视频一区| 久久久久九九精品影院| 精品熟女少妇av免费看| 国产在线男女| 久久九九热精品免费| 成人三级黄色视频| 女同久久另类99精品国产91| 中文字幕av在线有码专区| 成年免费大片在线观看| 爱豆传媒免费全集在线观看| 少妇人妻一区二区三区视频| 最好的美女福利视频网| 色哟哟·www| 波野结衣二区三区在线| av黄色大香蕉| 免费搜索国产男女视频| 少妇人妻一区二区三区视频| 啦啦啦韩国在线观看视频| 最新中文字幕久久久久| 日韩欧美精品免费久久| 中文字幕久久专区| 人妻少妇偷人精品九色| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲91精品色在线| 你懂的网址亚洲精品在线观看 | 欧美日本视频| 狂野欧美白嫩少妇大欣赏| 最新中文字幕久久久久| 国产亚洲精品av在线| 亚洲av电影不卡..在线观看| 日本一二三区视频观看| 欧美一区二区国产精品久久精品| 此物有八面人人有两片| 美女 人体艺术 gogo| 91在线精品国自产拍蜜月| 免费av毛片视频| 欧美人与善性xxx| 精品一区二区三区视频在线| www.色视频.com| 国产伦一二天堂av在线观看| 精品久久久久久成人av| 少妇人妻精品综合一区二区 | 久久鲁丝午夜福利片| 尤物成人国产欧美一区二区三区| a级毛色黄片| 天堂√8在线中文| 男人舔女人下体高潮全视频| 又粗又爽又猛毛片免费看| 黄色欧美视频在线观看| 99久久成人亚洲精品观看| 亚洲va在线va天堂va国产| 国产午夜福利久久久久久| 午夜福利视频1000在线观看| 国内揄拍国产精品人妻在线| 午夜福利在线观看免费完整高清在 | 少妇熟女欧美另类| 国产黄色小视频在线观看| 男的添女的下面高潮视频| 内射极品少妇av片p| 精品人妻一区二区三区麻豆| 日韩一区二区三区影片| 好男人视频免费观看在线| 亚洲内射少妇av| 哪个播放器可以免费观看大片| 国产成人一区二区在线| 欧美一级a爱片免费观看看| 99热这里只有是精品在线观看| 国产精品伦人一区二区| 国产亚洲av嫩草精品影院| 免费看美女性在线毛片视频| 男女啪啪激烈高潮av片| 91狼人影院| 精品久久久久久久久久久久久| 亚洲精品国产av成人精品| 搡老妇女老女人老熟妇| 国产极品精品免费视频能看的| 国产成人一区二区在线| 午夜视频国产福利| 国产精品麻豆人妻色哟哟久久 | 男女那种视频在线观看| 国产黄色小视频在线观看| 日本黄大片高清| 天堂网av新在线| 黑人高潮一二区| 欧美zozozo另类| 精品久久久久久久久亚洲| 一进一出抽搐gif免费好疼| 亚洲在久久综合| 看黄色毛片网站| 国产在线男女| 亚洲欧美成人精品一区二区| 人人妻人人看人人澡| 成人美女网站在线观看视频| 久久久精品大字幕| 精品久久久久久久久久免费视频| 国产欧美日韩精品一区二区| 精品一区二区三区视频在线| 黄色日韩在线| 少妇丰满av| 久久精品国产清高在天天线| 日韩欧美三级三区| 日韩在线高清观看一区二区三区| 青春草亚洲视频在线观看| 丰满人妻一区二区三区视频av| 成人国产麻豆网| 黄片无遮挡物在线观看| av天堂中文字幕网| 久久久久国产网址| 日韩在线高清观看一区二区三区| 日韩大尺度精品在线看网址| 好男人视频免费观看在线| 欧美一区二区国产精品久久精品| 精品一区二区免费观看| www日本黄色视频网| 波多野结衣高清无吗| 亚洲一区二区三区色噜噜| 亚洲最大成人中文| 在线国产一区二区在线| 级片在线观看| 国产成人aa在线观看| 亚洲高清免费不卡视频| 能在线免费看毛片的网站| 我要搜黄色片| a级毛片a级免费在线| 亚州av有码| 一区福利在线观看| 国产极品天堂在线| 成人午夜精彩视频在线观看| 中文字幕av在线有码专区| av专区在线播放| 日韩大尺度精品在线看网址| 免费在线观看成人毛片| 亚洲成人久久爱视频| 免费看av在线观看网站| 国产亚洲精品久久久com| 亚洲国产欧美人成| 欧美+日韩+精品| 黄色欧美视频在线观看| 国产综合懂色| 六月丁香七月| 精品免费久久久久久久清纯| 深爱激情五月婷婷| 久久久精品94久久精品| 中文精品一卡2卡3卡4更新| 亚洲自偷自拍三级| 久久精品国产亚洲av香蕉五月| 乱人视频在线观看| 一个人免费在线观看电影| 国模一区二区三区四区视频| 亚洲av中文av极速乱| 老司机福利观看| 12—13女人毛片做爰片一| 可以在线观看毛片的网站| 日本撒尿小便嘘嘘汇集6| 青春草国产在线视频 | 永久网站在线| 亚洲无线观看免费| 成人午夜精彩视频在线观看| 国产欧美日韩精品一区二区| 啦啦啦观看免费观看视频高清| 日韩大尺度精品在线看网址| 99久久成人亚洲精品观看| 日韩欧美精品免费久久| 美女内射精品一级片tv| 两个人视频免费观看高清| 国产精品爽爽va在线观看网站| 色哟哟·www| 又粗又硬又长又爽又黄的视频 | 亚洲精华国产精华液的使用体验 | 2022亚洲国产成人精品| 精品不卡国产一区二区三区| .国产精品久久| 国内精品一区二区在线观看| 97超碰精品成人国产| 成人欧美大片| 99久久成人亚洲精品观看| 亚洲第一电影网av| 级片在线观看| 午夜久久久久精精品| 级片在线观看| 简卡轻食公司| 深爱激情五月婷婷| or卡值多少钱| 欧美一区二区亚洲| 亚洲av第一区精品v没综合| 搡老妇女老女人老熟妇| 欧美性猛交╳xxx乱大交人| 九九热线精品视视频播放| 一边亲一边摸免费视频| 在线免费观看的www视频| 精品一区二区免费观看| 日韩三级伦理在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲丝袜综合中文字幕| 一进一出抽搐动态| 在线免费十八禁| 精品一区二区三区视频在线| 性插视频无遮挡在线免费观看| 99在线人妻在线中文字幕| 久久精品综合一区二区三区| 久久久久久九九精品二区国产| 免费观看在线日韩| 22中文网久久字幕| av天堂在线播放| 九九在线视频观看精品| 99久久久亚洲精品蜜臀av| 日本黄色视频三级网站网址| 国产精品国产高清国产av| 亚洲av二区三区四区| 波多野结衣高清无吗| 最好的美女福利视频网| www日本黄色视频网| 久久精品久久久久久噜噜老黄 | 成人特级av手机在线观看| 狂野欧美激情性xxxx在线观看| 人体艺术视频欧美日本| 亚洲人成网站在线播| 日本一本二区三区精品| 国产精品蜜桃在线观看 | 亚洲av.av天堂| 精品人妻偷拍中文字幕| 欧美不卡视频在线免费观看| 久久久久久久久久成人| 亚洲欧美成人精品一区二区| 一级毛片aaaaaa免费看小| 99久久精品一区二区三区| 精品人妻一区二区三区麻豆| 级片在线观看| 欧美另类亚洲清纯唯美| 国产精品一二三区在线看| 日韩av不卡免费在线播放| 你懂的网址亚洲精品在线观看 | 免费观看在线日韩| 欧美最新免费一区二区三区| 99视频精品全部免费 在线| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久久久免| 免费看日本二区| 亚洲美女搞黄在线观看| 亚洲av不卡在线观看| 免费黄网站久久成人精品| 国产高清三级在线| 久久久久免费精品人妻一区二区| 久久久久久久久久久丰满| 欧美日韩国产亚洲二区| 国产成人freesex在线| 亚洲成人av在线免费| 91午夜精品亚洲一区二区三区| 国产精品福利在线免费观看| 男女下面进入的视频免费午夜| 久久久久久久午夜电影| 狂野欧美白嫩少妇大欣赏| 黄色欧美视频在线观看| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站| 99热这里只有是精品在线观看| 国产成人午夜福利电影在线观看| 久久欧美精品欧美久久欧美| 男女边吃奶边做爰视频| 亚洲欧美中文字幕日韩二区| 久久久久久久亚洲中文字幕| 色哟哟哟哟哟哟| 色吧在线观看| 18禁在线无遮挡免费观看视频| 亚洲人成网站在线播| 亚洲经典国产精华液单| 久久婷婷人人爽人人干人人爱| 国产成人福利小说| 国产不卡一卡二| 国产成人福利小说| 精品免费久久久久久久清纯| 日韩av在线大香蕉| 国产精品女同一区二区软件| 一级二级三级毛片免费看| 一本久久中文字幕| 久久久国产成人精品二区| 亚洲三级黄色毛片| 国产伦在线观看视频一区| 欧美激情国产日韩精品一区| 女同久久另类99精品国产91| 国产极品天堂在线| 少妇丰满av| 男人舔女人下体高潮全视频| 亚洲国产欧美人成| 99国产精品一区二区蜜桃av| 男人和女人高潮做爰伦理| 国产单亲对白刺激| 最近最新中文字幕大全电影3| 成人综合一区亚洲| av在线亚洲专区| 婷婷色综合大香蕉| 亚洲国产欧美在线一区| 国产精品女同一区二区软件| 午夜免费激情av| 简卡轻食公司| 成人亚洲精品av一区二区| 免费人成在线观看视频色| 亚洲在久久综合| 一个人看视频在线观看www免费| 成人特级黄色片久久久久久久| 2021天堂中文幕一二区在线观| 免费av毛片视频| 日本爱情动作片www.在线观看| 成人无遮挡网站| 少妇熟女欧美另类| 久久精品国产亚洲av天美| 久久这里有精品视频免费| 久久久久久久久久成人| 边亲边吃奶的免费视频| 免费观看a级毛片全部| 日韩欧美三级三区| 激情 狠狠 欧美| 黄色视频,在线免费观看| 亚洲五月天丁香| 一本一本综合久久| 日韩高清综合在线| 麻豆av噜噜一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 日韩视频在线欧美| 天天一区二区日本电影三级| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品合色在线| 午夜福利高清视频| 久久这里只有精品中国| 丰满乱子伦码专区| 色视频www国产| 青春草国产在线视频 | 青春草亚洲视频在线观看| 久久久国产成人精品二区| 国产精品一区二区在线观看99 | 久久久久网色| 亚洲av中文字字幕乱码综合| 国产精品综合久久久久久久免费| 国产精品精品国产色婷婷| 极品教师在线视频| 三级男女做爰猛烈吃奶摸视频| 中文精品一卡2卡3卡4更新| 韩国av在线不卡| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩精品成人综合77777| 一区二区三区免费毛片| 美女cb高潮喷水在线观看| 非洲黑人性xxxx精品又粗又长| 乱码一卡2卡4卡精品| 黄色日韩在线| 一区福利在线观看| 丰满乱子伦码专区| 国产成人影院久久av| 韩国av在线不卡| 中文字幕久久专区| 国产黄片美女视频| 国产一区二区三区av在线 | 深夜精品福利| 久久人人爽人人爽人人片va| 久久午夜福利片| 国产单亲对白刺激| kizo精华| 黄色日韩在线| 天天躁夜夜躁狠狠久久av| 亚洲综合色惰| 一夜夜www| 亚洲成a人片在线一区二区| 久久久a久久爽久久v久久| 欧美最新免费一区二区三区| 中文欧美无线码| 国产精品免费一区二区三区在线| 国产高清有码在线观看视频| 在现免费观看毛片| 99热这里只有精品一区| 亚洲欧美日韩东京热| 国产三级中文精品| 99热只有精品国产| 美女内射精品一级片tv| 亚洲va在线va天堂va国产| 色视频www国产| 国产激情偷乱视频一区二区| 日本-黄色视频高清免费观看| 12—13女人毛片做爰片一| 高清日韩中文字幕在线| 乱人视频在线观看| 亚洲国产色片| 直男gayav资源| 两性午夜刺激爽爽歪歪视频在线观看| 日产精品乱码卡一卡2卡三| 波多野结衣高清作品| 国产成人精品婷婷| 嫩草影院入口| 欧洲精品卡2卡3卡4卡5卡区| 搞女人的毛片| 波多野结衣高清无吗| 中文精品一卡2卡3卡4更新| 男女边吃奶边做爰视频| 亚洲精品国产av成人精品| 别揉我奶头 嗯啊视频| 国产av一区在线观看免费| 国产探花极品一区二区| 国产精品久久久久久亚洲av鲁大| 99在线人妻在线中文字幕| 丝袜美腿在线中文| 桃色一区二区三区在线观看| 青春草视频在线免费观看| 网址你懂的国产日韩在线| 久久久成人免费电影| 久久韩国三级中文字幕| 男人和女人高潮做爰伦理| 一级毛片aaaaaa免费看小| 久久久久久九九精品二区国产| a级毛片免费高清观看在线播放| 亚洲真实伦在线观看| 有码 亚洲区| 国产精品一区www在线观看| 一本精品99久久精品77| h日本视频在线播放| 久久这里只有精品中国| 一级毛片我不卡| 青春草国产在线视频 | 一进一出抽搐gif免费好疼| 中国美白少妇内射xxxbb| 三级经典国产精品| 亚洲av中文av极速乱| 久久人人爽人人爽人人片va| 欧美成人一区二区免费高清观看| 亚洲不卡免费看| 免费无遮挡裸体视频| 国产真实伦视频高清在线观看| 国产黄色小视频在线观看| 老女人水多毛片| 成人三级黄色视频| 一进一出抽搐动态| 久久99精品国语久久久| 99国产极品粉嫩在线观看| 亚洲精品国产av成人精品| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 简卡轻食公司| 人体艺术视频欧美日本| 校园人妻丝袜中文字幕| 中文在线观看免费www的网站| 最好的美女福利视频网| 午夜精品在线福利| 国产毛片a区久久久久| 久久精品影院6| 日韩欧美国产在线观看| 免费看a级黄色片| 亚洲精品久久国产高清桃花| 91狼人影院| 我的女老师完整版在线观看| 蜜桃亚洲精品一区二区三区| 精品人妻偷拍中文字幕| 日本欧美国产在线视频| 岛国毛片在线播放| 免费观看人在逋| 亚洲精华国产精华液的使用体验 | 一区二区三区四区激情视频 | 亚洲内射少妇av| 在线播放国产精品三级| 国产不卡一卡二| 国产亚洲欧美98| 日韩av不卡免费在线播放| videossex国产| 精品少妇黑人巨大在线播放 | 欧美精品国产亚洲| 一本精品99久久精品77| 男人舔女人下体高潮全视频| 99热这里只有是精品在线观看| ponron亚洲| 亚洲最大成人手机在线| 波野结衣二区三区在线| or卡值多少钱| 最近最新中文字幕大全电影3| 看非洲黑人一级黄片| 一本一本综合久久| 国产高清不卡午夜福利| 校园人妻丝袜中文字幕| 看片在线看免费视频| 黑人高潮一二区| 中文字幕制服av| 自拍偷自拍亚洲精品老妇| 国产片特级美女逼逼视频| 国产精品一区二区在线观看99 | av在线播放精品| 人人妻人人看人人澡| av又黄又爽大尺度在线免费看 | 狂野欧美激情性xxxx在线观看| 成人美女网站在线观看视频| 男女做爰动态图高潮gif福利片| 蜜桃亚洲精品一区二区三区| 久久久午夜欧美精品| 亚洲四区av| 伦理电影大哥的女人| 一夜夜www| 欧美+日韩+精品| 色综合站精品国产| 中文字幕人妻熟人妻熟丝袜美| 日韩精品有码人妻一区| 久久精品久久久久久久性| 亚洲精品日韩在线中文字幕 | 97超碰精品成人国产| 日韩三级伦理在线观看| 又爽又黄a免费视频| 搡女人真爽免费视频火全软件| 久久精品国产99精品国产亚洲性色| 欧美+亚洲+日韩+国产| 可以在线观看的亚洲视频| 高清午夜精品一区二区三区 | 国产成人午夜福利电影在线观看| 欧美又色又爽又黄视频| 中国国产av一级| 午夜免费男女啪啪视频观看| 成人国产麻豆网| 亚洲欧美日韩高清专用| av天堂中文字幕网| 听说在线观看完整版免费高清| 国产激情偷乱视频一区二区| 黑人高潮一二区| 日韩欧美国产在线观看| 99久久精品热视频| 极品教师在线视频| 亚洲精品久久久久久婷婷小说 | 中文字幕av成人在线电影| 亚洲精品影视一区二区三区av| 日韩视频在线欧美| 热99在线观看视频| 欧美成人免费av一区二区三区| 亚洲av中文av极速乱| 99热6这里只有精品| 欧美+亚洲+日韩+国产| 春色校园在线视频观看| 欧美变态另类bdsm刘玥| 国产探花极品一区二区| or卡值多少钱| 国产精品电影一区二区三区| 亚洲内射少妇av| 精品国内亚洲2022精品成人| 黑人高潮一二区| 丝袜喷水一区| 波多野结衣高清作品| 黄色欧美视频在线观看| 一边亲一边摸免费视频| 黄色欧美视频在线观看| 床上黄色一级片| 亚洲丝袜综合中文字幕| АⅤ资源中文在线天堂| 青春草视频在线免费观看| 精品久久国产蜜桃| 国产伦精品一区二区三区四那| 综合色av麻豆| 人人妻人人看人人澡| 男插女下体视频免费在线播放| 夫妻性生交免费视频一级片| 高清毛片免费观看视频网站| 国产成人精品久久久久久| 熟女人妻精品中文字幕| 欧美日本亚洲视频在线播放| 性欧美人与动物交配| 成人国产麻豆网| 男女那种视频在线观看| 久久精品综合一区二区三区| 亚洲,欧美,日韩| 欧美在线一区亚洲| 变态另类成人亚洲欧美熟女| 老司机影院成人| av在线观看视频网站免费| 欧美极品一区二区三区四区| 亚洲一级一片aⅴ在线观看| 国产av麻豆久久久久久久| 久久久精品欧美日韩精品| 人妻制服诱惑在线中文字幕| 一级毛片我不卡| 久久热精品热| 免费看日本二区| 久久久久久伊人网av| 日韩视频在线欧美| 成人漫画全彩无遮挡| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 91久久精品电影网| 国产精品一二三区在线看| av免费观看日本| 国产大屁股一区二区在线视频| 国产成人福利小说| 特级一级黄色大片| 一级毛片我不卡|