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

    Investigation of decay modes of superheavy nuclei

    2021-12-15 10:31:52ManjunathaSowmyaDamodaraGuptaSridharNagarajaSeenappaAlfredCecilRaj
    Nuclear Science and Techniques 2021年11期

    H.C.Manjunatha ?N.Sowmya?P.S.Damodara Gupta?K.N.Sridhar ?A.M.Nagaraja,?L.Seenappa?S.Alfred Cecil Raj

    Abstract A detailed investigation of different decay modes,namely alpha decay,beta decay,cluster decay,including heavy particle emission(Zc>28),and spontaneous fission,was carried out,leading to the identification of new cluster and beta-plus emitters in superheavy nuclei with 104≤Z≤126.For the first time,we identified around 20 beta-plus emitters in superheavy nuclei.Heavy-particle radioactivity was observed in superheavy elements of atomic number in the range 116≤Z≤126.292-293Og were identified as 86Kr emitters,and 298122 and 300122 were identified as 94Zr emitters, whereas heavy-particle radioactivity from 91Y was also observed in 299123.Furthermore,the nuclei 300124 and 306126 exhibit 96Mo radioactivity.The reported regions of beta-plus and heavyparticle radioactivity for superheavy nuclei are stronger than those for alpha decay.The identified decay modes for superheavy nuclei are presented in a chart.This study is intended to serve as a reference for identifying possible decay modes in the superheavy region.

    Keywords Alpha decay.Beta decay.Heavy-particle radioactivity.Branching ratios

    1 Introduction

    The most important unanswered questions in Nuclear Physics are to determine the heaviest superheavy nucleithat can exist,and to investigate whether very-long-lived superheavy nuclei exist in nature.The past ten years have been marked by remarkable progress in the science of superheavy elements and nuclei.The existence of superheavy nuclei above Z=103 can be studied in terms of whether they can occur naturally or can be synthesized in the laboratory.There are no definitive conclusions regarding the existence of superheavy nuclei in nature.In contrast,such superheavy nuclei,with half-lives ranging between days toμs,can be synthesized using cold and hot fusion reactions.Cold fusion reactions involve either lead or bismuth as targets[1],whereas hot fusion reactions includE48Ca beamson various actinidetargets[2,3].Many theoretical predictions,such as microscopic–macroscopic[4](single-particle potential)and self-consistent approaches,including nucleus–nucleuspotential[5,6],relativistic field models[7,8],and multinucleon transfer reactions[9],provide information regarding the nucleus structure,shell closurelocation,and decay modesin heavy and superheavy nuclei.

    The discovery of superheavy elements[10,11]points to the island of stability.Boilley et al.[12]predicted the evaporation residue cross sections in superheavy elements and the influence of shell effects[13].The entrance channel dynamics were studied using48Ca as a projectile and208Pb as target[14].In 1966,two groups of researchers,namely Mayers and Swiatecki,and Viola and Seaborg[15],separately predicted the presence of heavy nucleinear the island of stability.Later,Sobiczewski et al.[16]predicted that the nucleus Z=114 will be the center of the island of stability,with neutron number N=184.In 1955,Nilsson[17]proposed a shell model which includes deformation property of the nuclei.Bender et al.[18]used a Skyrme energy density functional model and studied the deformation properties of closed proton and neutron shells.The nuclear mass,radius,and spectroscopy far away from the valley of stability were experimentally analyzed earlier[19].The investigation of isomers of the superheavy nucleus254No is a stepping stone toward the island of stability[20].Previous researchers[21]analyzed the nuclear shell structure and discovered additional stability near magic nuclei.The present scenario is almost near the center of the presumed island of stability,but the final landing is yet to be completed,and the intriguing question is how these superheavy nuclei are still accessible.

    The identification of superheavy nuclei is based on observations of decay chains.Superheavy element region 114≤Z≤118 were observed owing to their consistent decay chains,which end in the isotopes of rutherfordium(Rf)and dubnium(Db).Spontaneous fission andα-decay are the dominant decay modes in superheavy nuclei and limit their stability.Furthermore,newly synthesized superheavy elementsare primarily identified by their decay chains from unknown nuclei to known daughter nuclei by using the parent-and-daughter correlation.

    Thecompetition between different decay modes,such as ternary fission,spontaneous fission,cluster decay,proton decay,β-decay,and α-decay,in the heavy and superheavy region,has been extensively studied using various theoretical models,such as Coulomb and proximity potential models,modified generalized liquid drop models,effective liquid drop models,and temperature-dependent proximity potential models[22–33].The possible decay modes in the superheavy nuclei Z=119 and 120 are predicted in Ref.[34].From Ref.[35],it isclearly observed that theisotopes of the superheavy nuclei Z=104–112 have α-decay and spontaneous fission as dominant decay modes.However,onlyα-decay is dominant in the isotopes of superheavy nuclei Z=113,115–118.The isotopes of the superheavy nucleus Z=114 have spontaneous fission as the dominant decay mode in thenucleus284Fl,α-decay isdominant in the nuclei286-289Fl,andβ+is dominant in the nucleus290Fl.Furthermore,the concept of heavy-particle radioactivity[36]in thesuperheavy region hasimportant applicationsin the synthesis of superheavy nuclei.Despite the significant experimental and theoretical progress,there are many unanswered questions related to the decay modes of superheavy nuclei.Until now,onlyα-decay and spontaneous fission have been successfully observed in experiments.

    Experimental results suggest a considerable increase in the lifetime of nuclei as they approach closed proton and neutron shells[37].The lifetimes of most known superheavy nuclei are governed by the competition betweenαdecay and spontaneous fission.The existence of the island of stability has been confirmed experimentally in the previous decade[38].Some theoretical studies reveal that superheavy elements with 114 and 164 protons are stable against fission as well as alpha and beta decay[39].Various phenomenological and microscopic models,such as the fission model[40],the cluster model[41],generalized liquid drop model[42],and the unified model for alpha decay and alpha capture[43],are available to study the different decay modes of superheavy nuclei.In addition,many studies have been concerned with the alpha decay and spontaneous fission of superheavy nuclei[44–46].Simple empirical formulas are also available for determining the decay half-lives[47].The possible isotopes of new superheavy elements are identified by studying the competition between different probable decay modes,such as α-decay, β-decay,cluster decay,and spontaneous fission.This study focuses on the different decay modes of superheavy nuclei in the atomic number range 104≤Z≤126.After a detailed investigation of the competition between different decay modes,the possible isotopes and their decay modes with branching ratios are identified in the superheavy nuclei region.Hence,the contribution of this study is in the prediction of the most possible decay mode in superheavy nuclei,and in the identification of possible emitters in this superheavy region.The formalism is explained in Sect.2.The analysis of different decay modes and possible emitters in the superheavy region is explained in Sect.2.4.The paper is concluded in Sect.3.

    2 Theory

    2.1 Alpha decay and cluster decay

    In the effective liquid-drop model(ELDM),theα-decay half-life is computed using the relation

    where ν0is the assault frequency on the barrier,and ν0=1.8×1022s-1[48].Pαis the preformation factor,which is closely related to the shell structure[49].The empirical formula for Pαis expressed as

    where N,Z,and A are the neutron,charge,and mass number of the parent nucleus,respectively,Z1and Z2are the proton magic numbers around Z(Z1≤Z≤Z2),and N1and N2are the neutron magic numbers around N(N1≤N≤N2).p1,p2,and p3correspond to parameters in the region even(Z)-even(N),even(Z)-odd(N),odd(Z)-even(N),and odd(Z)-odd(N).They are presented in Table I of Ref.[50].P is the Gamow penetrability factor,given by the expression

    Table 1 Cluster-decay half-lives obtained from present study(PS)and available experiments(exp)

    Table 2 Comparison of logarithm half-lives(years)of spontaneous fission in the superheavy region 104≤Z≤114 from present study with those from available experiments

    Table 3 Comparison of alpha-decay half-lives from the present study(PS)and those from available experimental(Exp)values

    Table 4 Identified cluster emitters in the superheavy nuclei region

    Table 5 Identified alpha emitters in the superheavy nuclei region

    Table 5 continued

    Table 6 Identifiedβ+emitters in the superheavy nuclei region

    Table 7 Superheavy nuclei with dual decay mode and branching ratio

    Table 7 continued

    where μ is the inertial coefficient resulting from the Werner–Wheeler approximation[51].The limits of integration ζ0and ζcare the inner and outer turning points,expressed as ζ0=Rp-and.Rp is the radius of theparent nucleus,andis the final radius of the emitted cluster.In the ELDM,the total potential has been demonstrated to be the sum of Coulomb,proximity,and centrifugal potential[52,53].Hence,we can use the effective one-dimensional total potential energy as follows:

    To evaluate the Coulomb contribution in terms of the deformation parameter,we used V c as defined in Ref.[54]:

    with

    where S1and S2are the surface areas of the spherical fragments.σeffis the effective surface tension,which is defined as

    where R2is the final radius of the daughter fragment.The centrifugal potential energy is determined by

    where ? is the angular momentum of the emitted alpha/cluster and is calculated using the selection rules.In the case of alpha/cluster decay[55,56],the selection rules follow the condition

    In ELDM,a system with two intersecting spherical nuclei with different radii is considered[52].A schematic diagram for the representation of four independent coordinates,namely R1,R2,ζ,and ξ,is shown in Fig.1.Three constraints are used to reduce the four-dimensional spherical problem to an equivalent one-dimensional problem.The geometric constraint given below is introduced so that the spherical segments remain in contact:

    Fig.1 Schematic presentation of molecular phase of the dinuclear system(the daughter nucleus and the emitted(smaller)fragments).The distance between their geometrical centers and the distance between the center of the heavier fragment and the circular sharp neck of radius a are denoted byζandξ,respectively

    The variables ζ and ξ represent the distance between the geometrical centers and the distance between the center of the heavier fragment and the circular sharp neck of the radius,respectively[53,57].Assuming that nuclear matter is incompressible,the constraint for the conservation of the total volume of the system is

    where R=r0A1/3is the radius of the parent nucleus(r0=1.34 fm is an adjustable parameter),with A being the mass number of the parent.

    The radius of the α particle,R1,is assumed to be constant in the varying mass asymmetry shape description:

    2.2 Beta decay

    For all types ofβprocesses,the expression for the halflife Tβis given by[58]

    Here,EC isthe electron capture.For aparticular typeofβdecay,the half-life is expressed as follows:

    2.2.1 β±decay

    The Fermi function forβ-decay is expressed as

    Here,P(E)is the momentum of the particle,and F(E,Z)can be computed at the nuclear surface using themagnitude of the radial electron/positron wave function.The first approximation of F(E,Z)is

    At the surface of the nucleus(forβ+decay),the orbital electron screening effect has a significant impact on theβ electron/positron wave function.Thus,F(E,Z)becomes

    The expression for the energy released in β+decay is

    Similarly,for β-decay,

    2.2.2 Electron capture

    The value of Q for electron capture is given by

    Here,Be is the electron binding energy.Hence,even for the forbiddenβ+decay,electron capture is allowed.The captureof electrons of the K-shell for lower Z,and of the L-shell for higher Z is the major contributor to electron capture.The contributions of the electrons of higher shells are negligible.Thus,the Fermi function becomes

    In general,for any shell,

    Here,EXis the total energy of the electron:

    where ZXis the effective charge,which considers the screening of the Coulomb field of the nucleus by other electrons[61]:

    The nonzero components of the radial parts(gX&fX)of the wave function of the relativistic electron of orbit X are

    2.3 Spontaneous fission

    Spontaneous fission decay is studied by employing the quantum tunneling effect through the potential barrier.The decay constant of spontaneous fission is expressed as

    whereν,S,and P s are model-dependent quantities,namely assault frequency,preformation probability,and barrier penetrability,respectively.In the above equation,P=SP s and the spontaneous-fission half-lives are calculated as

    where h is the Planck constant,and Ev=hν/2 is the zeropoint vibration energy.The penetration probability is evaluated using the action integral K:

    and hence,the decimal logarithm of T(s)is given by

    If Eν=0.5 MeV,then the above equation becomes log10T=-log10P-20.5426.The action integral K is evaluated as follows:

    The term E(R)is the macroscopic energy in terms of the surface,volume,Coulomb,proximity energy,shell correction,and pairing energy term[62],and m is the rest mass of the neutron.A few superheavies are spherical,the rest are deformed,primarily prolate or oblate.To include this effect,deformations are also involved in the calculation of E(r),which is adopted from Ref.[62].In the above equation,R isthe separation distance between the center of the fission fragments,and Raand Rbare the turning points,which are evaluated using the boundary conditions E(Ra)and E(Rb)=Q.However,the term B(r)is the inertia with respect to r and is evaluated using the semi-empirical model for inertia[63]:

    whereμand k arethe reduced massof the fission fragments and a semi-empirical constant(k=14.8),respectively.R sph is the distance between the center of mass of the fission fragments,set as R sph/R0=0.75 in the symmetric case.The decay constant(λ)and the total fission decay constant are evaluated as described in Ref.[62].

    2.4 Results and discussion

    The mass excess values play a major role in the prediction of the decay mode and the corresponding half-lives.The predicted half-lives are sensitive to the Q-values,and small changes in the Q-values result in a notable change in the half-lives,with a magnitude of order 101to 102[36].Mass excess tables such as WS4[64],EBW[65],HFB28 and HFB29[66],DZ10[67],KTUY[68],finite-range droplet model(FRDM)[69],and AME16[70]areavailable in the literature.In the present study,we used the updated AME16[70]mass excess values up to Z=118,and above Z>118,themassexcessvaluesaretaken from the FRDM[69].The dominant decay mode is identified by studying the competition between different decay modes:α-decay,β-decay,cluster decay,and spontaneous fission in the superheavy nuclei region 104≤Z≤126.

    A detailed literature review indicates that there is no experimental evidence for cluster radioactivity in the superheavy region.Furthermore,experimental studies of cluster decay in the actinide region are available.To validatethepresent study,thecluster-decay half-livesobtained in the present study in the actinide region were compared with the experiments,and good agreement was observed.With this confidence,we studied cluster decay in the superheavy region,and theresultsarepresented in Table 1.Similarly,Table 2 shows a comparison of the studied logarithmic half-lives(in years)of spontaneous fission from the present study with those from available experiments.It can be seen that the cluster-decay and spontaneous-fission half-lives obtained in the present study are close to those of the experiments.

    As a part of this investigation,we studied theα-decay properties of superheavy nuclei using the formalism explained in the theory section.The predicted alpha-decay half-lives were validated by comparison with those from available experiments in the superheavy region.The results are given in Table 3.

    From the comparison,it is observed that the predicted half-lives are in good agreement with those of the experiments.With this confidence,we obtained the alpha-decay half-livesof superheavy nucleiin theregion 104≤Z≤126 Fig.2.

    Fig.2 (Color online)Map of nuclei reflecting the logarithmicαdecay half-lives for the isotopes of elements from Z=104 to 126.The Q-values were estimated using AME16 and FRDM95.The vertical line on the right side of the figure shows an increase in the log T1/2 values from the navy-blue region to the brown region

    shows a wide range ofα-decay half-lives.For a given superheavy nucleus,the alpha decay half-lives increase as the neutron number of its isotopes increases.For instance,theα-decay half-lives are of the order of nanoseconds at N/Z=1.307692 for Rutherfordium,whereas for the same superheavy element,theα-decay half-lives are of the order of 102s at N/Z=1.504762.Similarly,all neutron-rich superheavy nuclei have comparably longerα-decay halflives,which is in agreement with the report available in Ref.[73].The obtainedα-decay half-lives of all possible superheavy nuclei are presented in the heat map in Fig.2.The right vertical bar shows the magnitude of the log T1/2values.The color variation from navy blue to wine indicates values in the range 10-10–102s.The contrast in the blue region lies between 10-10s and 10-7s,in the green region,it lies in the range 10-6–10-4s,and the range 10-4-10-3s is presented in the yellow region.Finally,the red-to-wine region shows higher half-lives in the range 10-2-102s.The inset of Fig.2 on the top left side provides information on the magnified portion ofα-decay halflives in the superheavy region Z=104-114,whereas the bottom-right inset provides information on the magnified portion of the superheavy region Z=115-126.After the detailed investigation of theα-decay,a search wasmade to identify the cluster emitters in the superheavy region.Cluster radioactivity is energetically favorable if the Q-values are positive.We studied the possibility of cluster decay with 3≤Z c≤45 in the superheavy region 104≤Z≤126.For a given parent nucleus,the half-lives corresponding to various cluster emission were evaluated,and the cluster corresponding to shorter half-lives was identified.Furthermore,the cluster emitters corresponding to shorter half-lives for different isotopes of a given superheavy element were also identified.Eventually,cluster emissions corresponding to the shortest half-lives T c were identified;these are referred to as cluster-decay half-lives(T c).Thepredicted cluster decay half-livesin the atomic number region 104≤Z≤126 correspond to all the studied cluster emissions,as shown in Fig.3.

    Fig.3 (Color online)Predicted cluster-decay logarithmic halflives in the atomic number range Z=104–126 using AME16 and FRDM95 mass excess values.The hallow bin with different color in each panel shows the cluster emission corresponding to minimum half-lives

    This figure enables us to identify the cluster emission corresponding to the shorter half-lives of a given superheavy element.The half-lives of superheavy nuclei with Z=115–120 against cluster radioactivity are shorter for86Kr than those of the other studied clusters.The superheavy nuclei with Z=104,106,108,110,112,114,124,and 126 have shorter half-lives against96Mo cluster emissions than those of the other studied clusters.The decay half-lives are shorter for thE91Y emission from superheavy nuclei with Z=109,111,113,121,and 123.Similarly,the half-livesof superheavy nucleiwith Z=105 and 107 against cluster radioactivity are shorter for97Tc and101Rh than those of the other studied clusters.

    Cluster radioactivity in thesuperheavy nuclei region has shorter half-lives for cluster neutron numbers 44–48 from parent nuclei with neutron numbers 130–200,as shown in Fig.4.

    Fig.4 (Color online)Map of nucleireflecting thelogarithmic clusterdecay half-lives for neutron number of parent and cluster isotopes of elements with Z=104–126

    The range of cluster decay half-lives for superheavy elements with 104≤Z≤126 is shown in Fig.5.

    Fig.5 (Color online)Heat map showing the variations of lowest logarithmic half lives of clusters with 104<Z<126

    Shorter half-lives are observed for N/Z>1.37068,and larger half-lives are observed for N/Z<1.37068.From the figure,it is clear that up to superheavy nuclei 104≤Z≤115,larger cluster-decay half-livesareobserved,whereasshorter cluster-decay half-livesareobserved in the superheavy region 116≤Z≤126.The inset of Fig.5 on the top-left side shows a magnified portion of the logarithmic half-lives (T c) in the superheavy region 104≤Z≤115,whereas the inset at the right bottom shows a magnified portion of the shorter logarithmic half-lives(T c)in the superheavy region 116≤Z≤126.This figure also shows that some of the superheavy nuclei have lifetimes of the order of ns toμs and exhibit cluster decay.

    The other prominent decay mode that was studied is spontaneous fission,which is also energetically feasible in heavy and superheavy nuclei.It may occur in such nuclei owing to an increase in the Coulomb interactions.References[10,11,38,74–77]report consistentα-decay chains from superheavy nuclei followed by spontaneous fission.The spontaneous fission half-lives are studied using the theory explained in Sect.2.3.The variations of spontaneous fission half-lives in the superheavy region Z=104–126 are shown in Fig.6.

    Fig.6 (Color online)Heat map of the variations of logarithmic halflives for spontaneous fission for 104<Z<126

    The log T SF valuesvary between-50(dark blue region)and 50(dark-red region).For instance,at atomic number Z=104,for isotopes 245–275,the log T1/2(SF)values ranging from-50 to 5 are shown,whereas the half-lives with smaller values are indicated by the color range from navy blue to blue.The half-livesranging from nanoseconds to 105s are indicated by the color range from yellow to light orange.Similarly,in the atomic number range Z=119 and above,larger values of spontaneous-fission logarithmic half-lives are indicated by the red color range.Thus,on either sideof Fig.6,for isotopescorresponding to the atomic number range Z=104–126,smaller half-lives are observed,whereas in the middle region of the figure,larger values of log T1/2are observed up to Z=116.In contrast,smaller half-lives are observed for higher isotopes(Z>116),and larger log T1/2for lower isotopes(Z<116).A similar trend was also observed in a previous study[78],in which the half-lives of nuclei Z=92–104 were compared with experimentally available values.

    A detailed investigation of the Q-values corresponding to β-decay in the superheavy region demonstrates thatβ+-decay is energetically possible with Z=105,107,113,114,115,117,119,121,123,125,and 126,whereasβ--decay is not energetically possible.Furthermore,we also studiedβ-decay half-lives using the formalism explained in Sect.2.2.1.

    The competition between different possible decay modes,namely α-decay,cluster-decay,β-decay and spontaneous fission,enables us to identify the dominant decay modefor superheavy elementsin theatomic number region 104≤Z≤126 of all possible isotopes Fig.7.

    Fig.7 (Color online)Chart of spontaneous fission(purple),alpha decay(brown),β+-decay(cyan),and cluster emitters(yellow)with atomic numbers Z=104–126.The Q-values were calculated using the FRDM95 mass tables

    shows the decay modes of the superheavy nuclei.In the studied superheavy region,we identified around 20 β+emitters,which are presented in Table 6.We also identified 35 cluster emitters,which are presented in Table 4.

    It was demonstrated that the majority of superheavy nuclei undergo α-decay and spontaneous fission.The αemitting superheavy nuclei are listed in Table 5.

    The identified alpha emitters have half-lives of approximatelyμs to 100 s in the superheavy region 104≤Z≤126.Table 4 liststheidentified cluster emissions with the corresponding half-lives.The amount of energy released during cluster emission,cluster emitted,and log T1/2values are presented in the table.The minimum cluster decay half-lives correspond to86Kr,94Zr,91Y,and96Mo for the nuclei292-293Og,298,300122,299123,300124,and306126,respectively.From the available literature,it is also evident that the heavy particle radioactivity of86Kr is observed in the superheavy nucleus Z=118[36,79].In addition,Rb,Sr,Y,Zr,Nb,and Mo cluster emissions[80]were observed for Z=119–124,respectively.As in previous studies,in the present study,shorter half-lives in the superheavy region Z=118,122-124,and 126 were observed,with86Kr,94Zr,91Y,and96Mo cluster emissions,respectively.Similarly,approximately 20 possibleβ+emitters were identified in the superheavy region 105≤Z≤125,and they are presented in Table 6.

    The information provided Table 7 regarding the halflives and branching ratios presents ambiguities in terms of determining a single decay mode.The branching ratios relativeto theminimum half-livesamong thestudied decay modes are obtained,and the second column of the table shows the log T1/2values corresponding to spontaneous-fission,α-decay,β+-decay,and cluster-decay halflives.For instance,the superheavy nucleus263Rf exhibits shorter log T1/2values for spontaneous fission andβ+-decay than for other decay modes.The branching ratio of spontaneous fission andβ+-decay was obtained,and it was found that the branching ratio corresponding to spontaneous fission andβ+was 55%and 45%,respectively.Similarly,we identified the branching ratios for the superheavy region 104≤Z≤126,which are presented in Table 7.

    Finally,Fig.8 shows the lifetimes of the superheavy elements after the competition between different decay modes was studied.

    Fig.8 (Color online)Heat map showing the variations of atomic number,massnumber of parent and logarithmic half-livesof different decay modes(life times)for 104<Z<126

    It can be seen that the lifetime varies from ns to min and decreases asthe atomic number increases.For instance,the average lifetime of a superheavy element with Z=104 is approximately 10 min,whereas that of a hypothetical superheavy element with Z=126 is of the order of ms.

    3 Conclusion

    We systematically investigated all possible decay modes,namely α-decay, β-decay,cluster decay,and spontaneous fission, in the superheavy region 104≤Z≤126.Thefindingsof thisstudy werevalidated by comparison with experiments.Approximately 20β+and 7 heavy particle emitters were found in the superheavy region.Furthermore,the nuclei with almost the same halflives for the two decay modes were also reported,with the corresponding branching ratios.However,an experimental study is necessary to draw definite conclusions.

    Author contributionsAll authors contributed to the study conception and design.Conceived of the original idea,developed the theory formulation of work and writing by HCM and NS.Performed the computations and performed,analytic calculations and graphical representation by PSDG,KNS and AMN,data analysis by LS and SACR.All authors read and approved the final manuscript.

    99视频精品全部免费 在线| 久久午夜综合久久蜜桃| 午夜日本视频在线| 午夜免费男女啪啪视频观看| 国产精品久久久久久精品电影小说| 日本猛色少妇xxxxx猛交久久| 日产精品乱码卡一卡2卡三| 伦精品一区二区三区| 女的被弄到高潮叫床怎么办| 久久国产亚洲av麻豆专区| 亚洲av男天堂| 天堂中文最新版在线下载| 免费少妇av软件| 久久人人爽人人片av| 夫妻午夜视频| 亚洲天堂av无毛| 大香蕉久久网| 国产又爽黄色视频| 成人二区视频| 丰满少妇做爰视频| 2018国产大陆天天弄谢| 最新的欧美精品一区二区| 看十八女毛片水多多多| 夜夜爽夜夜爽视频| 欧美成人午夜精品| 嫩草影院入口| 欧美少妇被猛烈插入视频| 国产 一区精品| 成人免费观看视频高清| 精品亚洲成a人片在线观看| 五月天丁香电影| 久热久热在线精品观看| 色吧在线观看| 国产精品国产三级专区第一集| 成年人免费黄色播放视频| 精品人妻熟女毛片av久久网站| 精品一区二区三卡| 如何舔出高潮| 精品卡一卡二卡四卡免费| 多毛熟女@视频| 美女大奶头黄色视频| 少妇被粗大的猛进出69影院 | 亚洲国产成人一精品久久久| 国产在视频线精品| 亚洲精品久久午夜乱码| 在线观看三级黄色| 大香蕉久久成人网| 看十八女毛片水多多多| 久久婷婷青草| 内地一区二区视频在线| 大香蕉97超碰在线| 日本猛色少妇xxxxx猛交久久| 精品一区二区三区四区五区乱码 | 国产一区有黄有色的免费视频| 蜜桃国产av成人99| 女性被躁到高潮视频| 久久影院123| 美女国产视频在线观看| 天天躁夜夜躁狠狠久久av| 香蕉丝袜av| 精品人妻熟女毛片av久久网站| 成年动漫av网址| 国产亚洲精品第一综合不卡 | 90打野战视频偷拍视频| 午夜免费观看性视频| 丁香六月天网| 欧美3d第一页| 又黄又粗又硬又大视频| 爱豆传媒免费全集在线观看| 超色免费av| 精品熟女少妇av免费看| 九草在线视频观看| 久久婷婷青草| 国产淫语在线视频| 飞空精品影院首页| 欧美精品高潮呻吟av久久| 午夜免费观看性视频| 欧美 日韩 精品 国产| 免费日韩欧美在线观看| 综合色丁香网| 日韩三级伦理在线观看| 久久久久久久久久久久大奶| 永久网站在线| 国产午夜精品一二区理论片| 午夜精品国产一区二区电影| av卡一久久| 在线观看免费高清a一片| 久久99热这里只频精品6学生| 热re99久久国产66热| 一级黄片播放器| 在线观看人妻少妇| 国产欧美亚洲国产| 女性生殖器流出的白浆| 亚洲成色77777| 成年美女黄网站色视频大全免费| 亚洲欧美一区二区三区黑人 | 色5月婷婷丁香| 免费人成在线观看视频色| 久久久久久久国产电影| 成年人午夜在线观看视频| 精品卡一卡二卡四卡免费| 夫妻午夜视频| 美女国产视频在线观看| 亚洲精品久久久久久婷婷小说| 色网站视频免费| 亚洲精品一二三| 新久久久久国产一级毛片| 国产精品久久久久久精品古装| 男女午夜视频在线观看 | av片东京热男人的天堂| 人成视频在线观看免费观看| 国产 一区精品| 波多野结衣一区麻豆| 亚洲精品,欧美精品| 亚洲美女黄色视频免费看| 国产欧美日韩综合在线一区二区| 菩萨蛮人人尽说江南好唐韦庄| 久久久国产一区二区| 大片电影免费在线观看免费| 国产精品久久久久久精品古装| 久久久久久伊人网av| 亚洲天堂av无毛| av免费观看日本| a级毛片在线看网站| 这个男人来自地球电影免费观看 | 亚洲欧洲精品一区二区精品久久久 | 在线 av 中文字幕| 啦啦啦视频在线资源免费观看| 欧美少妇被猛烈插入视频| 97在线视频观看| 搡老乐熟女国产| 亚洲精品色激情综合| 国产精品偷伦视频观看了| 精品国产一区二区三区久久久樱花| 乱人伦中国视频| 女的被弄到高潮叫床怎么办| 日韩免费高清中文字幕av| 久久午夜综合久久蜜桃| 少妇的丰满在线观看| 搡老乐熟女国产| 中文字幕av电影在线播放| 久久精品国产综合久久久 | 在线精品无人区一区二区三| 亚洲中文av在线| 亚洲欧美精品自产自拍| 欧美精品一区二区大全| 亚洲精品视频女| av免费在线看不卡| 亚洲综合色网址| 亚洲精品av麻豆狂野| 国产精品国产三级国产专区5o| 男女无遮挡免费网站观看| 国产白丝娇喘喷水9色精品| 国产精品成人在线| 国产亚洲欧美精品永久| 国产成人av激情在线播放| 中文乱码字字幕精品一区二区三区| 久久99热6这里只有精品| 日本欧美视频一区| 大香蕉久久网| 国产精品成人在线| 国国产精品蜜臀av免费| videossex国产| 久久久久久人人人人人| 成年动漫av网址| 日本黄大片高清| 精品一区二区免费观看| 国产成人a∨麻豆精品| 成年av动漫网址| 夫妻性生交免费视频一级片| 免费观看a级毛片全部| 青春草亚洲视频在线观看| 日日摸夜夜添夜夜爱| 香蕉精品网在线| 久久久国产欧美日韩av| 日韩一区二区视频免费看| 女人久久www免费人成看片| 欧美日韩一区二区视频在线观看视频在线| 午夜av观看不卡| 免费av中文字幕在线| 亚洲精品日本国产第一区| 久久精品国产亚洲av天美| 国产毛片在线视频| 在线看a的网站| 国产成人91sexporn| xxxhd国产人妻xxx| 亚洲欧美中文字幕日韩二区| 97在线人人人人妻| 成人亚洲欧美一区二区av| 91国产中文字幕| 制服人妻中文乱码| 日本欧美视频一区| 日本wwww免费看| 九色亚洲精品在线播放| 熟女电影av网| 国产精品女同一区二区软件| 日本91视频免费播放| 这个男人来自地球电影免费观看 | 国产亚洲午夜精品一区二区久久| 久久精品国产亚洲av天美| 日韩一本色道免费dvd| 在线观看免费高清a一片| 秋霞伦理黄片| 亚洲精品乱久久久久久| 中文字幕亚洲精品专区| 熟女电影av网| 桃花免费在线播放| 999精品在线视频| 看免费成人av毛片| 国产一区二区在线观看av| 欧美日韩av久久| 五月开心婷婷网| 一级a做视频免费观看| 女的被弄到高潮叫床怎么办| 亚洲精品久久久久久婷婷小说| 亚洲国产最新在线播放| 黄色一级大片看看| 一级a做视频免费观看| 久久人人97超碰香蕉20202| 男女啪啪激烈高潮av片| 日韩欧美精品免费久久| 99re6热这里在线精品视频| 精品午夜福利在线看| 香蕉精品网在线| 人成视频在线观看免费观看| 王馨瑶露胸无遮挡在线观看| 色5月婷婷丁香| 五月玫瑰六月丁香| 在线亚洲精品国产二区图片欧美| 丝袜人妻中文字幕| 亚洲国产色片| 亚洲第一区二区三区不卡| 综合色丁香网| 18禁国产床啪视频网站| 免费不卡的大黄色大毛片视频在线观看| 妹子高潮喷水视频| 视频在线观看一区二区三区| 一区二区三区乱码不卡18| 1024视频免费在线观看| 中文字幕av电影在线播放| 久久久精品免费免费高清| 少妇人妻久久综合中文| 久久久久久久大尺度免费视频| 国产永久视频网站| 国产精品嫩草影院av在线观看| 看免费成人av毛片| 天堂8中文在线网| 丝袜喷水一区| 亚洲av成人精品一二三区| 99热6这里只有精品| videosex国产| 亚洲情色 制服丝袜| 久久精品国产鲁丝片午夜精品| 在线观看免费日韩欧美大片| 亚洲av免费高清在线观看| 日本黄色日本黄色录像| 黑人欧美特级aaaaaa片| 国产成人午夜福利电影在线观看| 草草在线视频免费看| 水蜜桃什么品种好| 亚洲欧美精品自产自拍| 亚洲欧美清纯卡通| 国产精品.久久久| 久久热在线av| 宅男免费午夜| 国产片内射在线| 日本91视频免费播放| 中文字幕人妻熟女乱码| 久久人人爽人人片av| 寂寞人妻少妇视频99o| 亚洲欧美日韩卡通动漫| 尾随美女入室| 青春草亚洲视频在线观看| 精品一区二区三卡| 国产精品国产三级国产av玫瑰| 美女视频免费永久观看网站| 国产黄频视频在线观看| 三级国产精品片| 大话2 男鬼变身卡| 亚洲四区av| 岛国毛片在线播放| 9色porny在线观看| 尾随美女入室| 伦理电影免费视频| 青春草视频在线免费观看| 美女国产视频在线观看| 久久精品国产a三级三级三级| 制服诱惑二区| 麻豆精品久久久久久蜜桃| 亚洲欧美一区二区三区黑人 | 免费黄网站久久成人精品| 2022亚洲国产成人精品| 亚洲四区av| 国产欧美亚洲国产| 午夜免费鲁丝| 久久精品国产亚洲av涩爱| 男女边吃奶边做爰视频| 成年女人在线观看亚洲视频| 久久久久视频综合| 国产日韩欧美亚洲二区| 99热网站在线观看| 香蕉精品网在线| 亚洲精品视频女| 国产精品熟女久久久久浪| 90打野战视频偷拍视频| 91在线精品国自产拍蜜月| 久久这里有精品视频免费| www.色视频.com| 韩国高清视频一区二区三区| 精品国产一区二区三区久久久樱花| 高清不卡的av网站| 日韩,欧美,国产一区二区三区| 亚洲av.av天堂| 欧美丝袜亚洲另类| 日韩av在线免费看完整版不卡| 精品国产国语对白av| av有码第一页| 国产成人a∨麻豆精品| 免费看光身美女| av天堂久久9| 久久久精品免费免费高清| 中文精品一卡2卡3卡4更新| 丁香六月天网| 亚洲国产av新网站| 韩国精品一区二区三区 | 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 99久久人妻综合| 久久久久精品久久久久真实原创| 熟妇人妻不卡中文字幕| 草草在线视频免费看| 亚洲精品成人av观看孕妇| 国产白丝娇喘喷水9色精品| 各种免费的搞黄视频| 最后的刺客免费高清国语| videos熟女内射| 成年美女黄网站色视频大全免费| 国产高清国产精品国产三级| 午夜精品国产一区二区电影| 国产免费现黄频在线看| 极品人妻少妇av视频| 性色avwww在线观看| 久热这里只有精品99| 永久免费av网站大全| 精品久久国产蜜桃| 18禁动态无遮挡网站| 欧美日韩一区二区视频在线观看视频在线| 免费av中文字幕在线| 亚洲国产精品成人久久小说| 黄色一级大片看看| 在线观看三级黄色| 麻豆乱淫一区二区| 精品国产一区二区三区久久久樱花| 久久久久人妻精品一区果冻| 国产av一区二区精品久久| 精品久久国产蜜桃| 国产又色又爽无遮挡免| 成人漫画全彩无遮挡| 赤兔流量卡办理| 国国产精品蜜臀av免费| 80岁老熟妇乱子伦牲交| 久久久久人妻精品一区果冻| 人人妻人人爽人人添夜夜欢视频| 亚洲精品久久午夜乱码| 伦理电影大哥的女人| 国产av国产精品国产| 久久鲁丝午夜福利片| 看免费av毛片| 国产精品.久久久| 久久韩国三级中文字幕| 亚洲经典国产精华液单| 国语对白做爰xxxⅹ性视频网站| 国产男女超爽视频在线观看| 一本大道久久a久久精品| 你懂的网址亚洲精品在线观看| 国产在视频线精品| 日本免费在线观看一区| 久久久久网色| 亚洲欧美精品自产自拍| 国产 精品1| 交换朋友夫妻互换小说| 亚洲精品成人av观看孕妇| 丝袜人妻中文字幕| 国产精品久久久久久精品古装| 国产一区二区三区av在线| 蜜桃国产av成人99| 国产国语露脸激情在线看| 狠狠精品人妻久久久久久综合| 美女主播在线视频| 亚洲伊人色综图| 菩萨蛮人人尽说江南好唐韦庄| 亚洲伊人色综图| 亚洲av男天堂| 国产精品久久久av美女十八| 有码 亚洲区| 亚洲丝袜综合中文字幕| 欧美精品av麻豆av| 欧美激情国产日韩精品一区| 国产免费视频播放在线视频| 欧美激情极品国产一区二区三区 | 国产免费福利视频在线观看| 亚洲欧美清纯卡通| 中文天堂在线官网| 建设人人有责人人尽责人人享有的| freevideosex欧美| 亚洲国产欧美日韩在线播放| 精品久久久久久电影网| 在线观看人妻少妇| 国产激情久久老熟女| 在线观看免费视频网站a站| 欧美少妇被猛烈插入视频| 宅男免费午夜| 日日爽夜夜爽网站| 欧美日韩综合久久久久久| 又黄又粗又硬又大视频| 高清黄色对白视频在线免费看| 人妻少妇偷人精品九色| 国产极品天堂在线| 亚洲精品第二区| a级毛片在线看网站| 夫妻性生交免费视频一级片| 国产精品久久久久成人av| 免费观看在线日韩| 欧美日韩av久久| 汤姆久久久久久久影院中文字幕| 热re99久久精品国产66热6| 国产 一区精品| 国产爽快片一区二区三区| 国产麻豆69| 黑人巨大精品欧美一区二区蜜桃 | 少妇的逼水好多| 久久ye,这里只有精品| 色5月婷婷丁香| 久久人人爽av亚洲精品天堂| 国产一区二区三区综合在线观看 | 亚洲第一区二区三区不卡| 五月天丁香电影| 精品人妻偷拍中文字幕| 侵犯人妻中文字幕一二三四区| 插逼视频在线观看| 亚洲精品日本国产第一区| 高清黄色对白视频在线免费看| 丝袜脚勾引网站| 亚洲精品aⅴ在线观看| 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 侵犯人妻中文字幕一二三四区| 一区二区三区精品91| 999精品在线视频| 日韩精品有码人妻一区| 伊人久久国产一区二区| 国产男女内射视频| 久久av网站| av在线播放精品| 欧美最新免费一区二区三区| av不卡在线播放| 国产免费一级a男人的天堂| 久久精品夜色国产| 精品国产乱码久久久久久小说| 丰满饥渴人妻一区二区三| 久久99精品国语久久久| 热99国产精品久久久久久7| 国产成人精品久久久久久| 美女国产高潮福利片在线看| 丝袜人妻中文字幕| 亚洲成人一二三区av| 永久网站在线| 99热6这里只有精品| 亚洲天堂av无毛| 亚洲图色成人| 自拍欧美九色日韩亚洲蝌蚪91| 成年动漫av网址| 欧美另类一区| 欧美+日韩+精品| 久久韩国三级中文字幕| 天堂中文最新版在线下载| 一本久久精品| 日韩av不卡免费在线播放| 妹子高潮喷水视频| 精品一区二区免费观看| 免费少妇av软件| 多毛熟女@视频| 国产在线视频一区二区| 久久国内精品自在自线图片| 亚洲精品视频女| 纯流量卡能插随身wifi吗| 一级毛片 在线播放| 考比视频在线观看| 国产免费视频播放在线视频| 国产成人精品无人区| 国产69精品久久久久777片| 国产麻豆69| 成年动漫av网址| 亚洲精品久久午夜乱码| 免费观看在线日韩| 韩国高清视频一区二区三区| 99热网站在线观看| 国产不卡av网站在线观看| 亚洲欧美日韩卡通动漫| 制服丝袜香蕉在线| 天堂俺去俺来也www色官网| 国产一区二区在线观看日韩| 又黄又爽又刺激的免费视频.| 久久久精品94久久精品| 黄色一级大片看看| 国内精品宾馆在线| 大香蕉久久成人网| 精品亚洲乱码少妇综合久久| 国产亚洲午夜精品一区二区久久| 国产一区二区三区综合在线观看 | 伊人亚洲综合成人网| 国产精品一二三区在线看| 亚洲欧洲国产日韩| 日韩电影二区| 青青草视频在线视频观看| 成人亚洲精品一区在线观看| av在线播放精品| 捣出白浆h1v1| 狂野欧美激情性xxxx在线观看| 国产成人精品一,二区| 日本色播在线视频| 妹子高潮喷水视频| 国产黄频视频在线观看| 国产精品国产三级国产专区5o| 男女啪啪激烈高潮av片| 亚洲人与动物交配视频| 精品人妻偷拍中文字幕| 97人妻天天添夜夜摸| 永久免费av网站大全| av女优亚洲男人天堂| 在线看a的网站| 中文乱码字字幕精品一区二区三区| 男女免费视频国产| 热99国产精品久久久久久7| 永久网站在线| 高清欧美精品videossex| 亚洲第一区二区三区不卡| 久久国内精品自在自线图片| 天堂中文最新版在线下载| 日本欧美国产在线视频| 国产成人精品无人区| 51国产日韩欧美| 街头女战士在线观看网站| 国产无遮挡羞羞视频在线观看| 欧美人与性动交α欧美精品济南到 | 免费大片18禁| 你懂的网址亚洲精品在线观看| 久久精品久久久久久久性| 亚洲精品乱久久久久久| 免费看光身美女| a级毛片黄视频| 久久精品久久久久久噜噜老黄| 日韩欧美一区视频在线观看| 卡戴珊不雅视频在线播放| 日韩av不卡免费在线播放| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 夫妻性生交免费视频一级片| 美女内射精品一级片tv| 亚洲欧美精品自产自拍| 国产成人免费观看mmmm| 欧美国产精品va在线观看不卡| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| 九草在线视频观看| 最近中文字幕高清免费大全6| 九草在线视频观看| 国产1区2区3区精品| 国产毛片在线视频| 日本av手机在线免费观看| 国产午夜精品一二区理论片| 亚洲综合色网址| 成年动漫av网址| 日韩三级伦理在线观看| 国产成人精品一,二区| 18禁裸乳无遮挡动漫免费视频| 久久婷婷青草| 成人漫画全彩无遮挡| 亚洲国产日韩一区二区| 精品少妇内射三级| 免费女性裸体啪啪无遮挡网站| 黄色视频在线播放观看不卡| 两性夫妻黄色片 | 少妇被粗大的猛进出69影院 | 狠狠精品人妻久久久久久综合| 亚洲av男天堂| 男女无遮挡免费网站观看| 精品一区在线观看国产| 色5月婷婷丁香| 热re99久久国产66热| 午夜福利在线观看免费完整高清在| 成年人午夜在线观看视频| 欧美日韩成人在线一区二区| 免费少妇av软件| 男男h啪啪无遮挡| av播播在线观看一区| 成年女人在线观看亚洲视频| 欧美日韩国产mv在线观看视频| 黄片无遮挡物在线观看| 国产又爽黄色视频| 国产成人免费无遮挡视频| 国产精品国产三级专区第一集| 免费少妇av软件| 亚洲高清免费不卡视频| 超色免费av| 免费人成在线观看视频色| 丰满少妇做爰视频| 五月伊人婷婷丁香| 蜜桃在线观看..| 亚洲精品日本国产第一区| 亚洲国产精品专区欧美| 亚洲在久久综合| 侵犯人妻中文字幕一二三四区| av又黄又爽大尺度在线免费看| 女性生殖器流出的白浆|