• <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.

    婷婷亚洲欧美| 久久久久九九精品影院| 欧美3d第一页| 可以在线观看毛片的网站| 亚洲专区国产一区二区| 日本 av在线| 日本免费a在线| 三级毛片av免费| 久久久久久久亚洲中文字幕 | 免费高清视频大片| 最近最新免费中文字幕在线| 久久久久精品国产欧美久久久| 又黄又爽又免费观看的视频| 国产主播在线观看一区二区| 久久天躁狠狠躁夜夜2o2o| 成熟少妇高潮喷水视频| 天天一区二区日本电影三级| 深夜a级毛片| 深夜精品福利| 久9热在线精品视频| 久久人妻av系列| av天堂中文字幕网| 亚洲综合色惰| 亚洲欧美清纯卡通| 亚洲人成网站高清观看| 1000部很黄的大片| 亚洲国产欧美人成| 国内精品久久久久久久电影| 好看av亚洲va欧美ⅴa在| 91久久精品电影网| 90打野战视频偷拍视频| 在线看三级毛片| 免费一级毛片在线播放高清视频| 91久久精品电影网| 亚洲性夜色夜夜综合| 可以在线观看的亚洲视频| 怎么达到女性高潮| 亚洲专区中文字幕在线| 国产一区二区激情短视频| 久久午夜亚洲精品久久| 国产探花在线观看一区二区| 悠悠久久av| 长腿黑丝高跟| 亚洲成人久久爱视频| 精品一区二区三区av网在线观看| 日本黄大片高清| 亚洲av美国av| 亚洲aⅴ乱码一区二区在线播放| 国产精品亚洲av一区麻豆| 亚洲成av人片在线播放无| 久久久久久九九精品二区国产| 亚洲av成人av| 狂野欧美白嫩少妇大欣赏| 免费人成在线观看视频色| 1000部很黄的大片| 国产精品三级大全| 女人被狂操c到高潮| 好看av亚洲va欧美ⅴa在| 最好的美女福利视频网| av视频在线观看入口| 国产精品伦人一区二区| 亚洲av二区三区四区| av欧美777| av视频在线观看入口| 欧美极品一区二区三区四区| 久久香蕉精品热| 国产欧美日韩一区二区三| 色精品久久人妻99蜜桃| 高潮久久久久久久久久久不卡| 亚洲精品在线美女| a级毛片免费高清观看在线播放| 精品欧美国产一区二区三| 成人美女网站在线观看视频| 人人妻,人人澡人人爽秒播| a级毛片a级免费在线| 国产成人aa在线观看| 免费一级毛片在线播放高清视频| 国产真实伦视频高清在线观看 | 伦理电影大哥的女人| 国产午夜精品久久久久久一区二区三区 | 国产精品久久电影中文字幕| 日本成人三级电影网站| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| 亚洲av中文字字幕乱码综合| 村上凉子中文字幕在线| 免费人成视频x8x8入口观看| 老女人水多毛片| 欧美日韩黄片免| 亚洲国产欧洲综合997久久,| 一级作爱视频免费观看| 日日摸夜夜添夜夜添小说| 久久精品人妻少妇| 99国产综合亚洲精品| 国产精品亚洲av一区麻豆| 欧美+日韩+精品| 亚洲国产精品sss在线观看| 99国产精品一区二区蜜桃av| 亚洲在线自拍视频| 成人av一区二区三区在线看| 免费观看人在逋| 久久精品国产99精品国产亚洲性色| 日韩国内少妇激情av| avwww免费| 亚洲五月天丁香| 黄色丝袜av网址大全| 全区人妻精品视频| 久久精品综合一区二区三区| 夜夜夜夜夜久久久久| 欧美绝顶高潮抽搐喷水| 我的女老师完整版在线观看| 国产美女午夜福利| 欧美精品啪啪一区二区三区| 高清日韩中文字幕在线| 欧美日韩福利视频一区二区| 国产免费男女视频| 色吧在线观看| 久久久色成人| 午夜视频国产福利| 丁香欧美五月| 免费观看的影片在线观看| 中文字幕久久专区| 日韩欧美精品v在线| 亚洲aⅴ乱码一区二区在线播放| 99久久成人亚洲精品观看| 熟女电影av网| 美女xxoo啪啪120秒动态图 | 女生性感内裤真人,穿戴方法视频| 一级黄片播放器| 国产精品99久久久久久久久| 我的女老师完整版在线观看| 99热6这里只有精品| 亚洲va日本ⅴa欧美va伊人久久| 国产主播在线观看一区二区| 日韩有码中文字幕| 黄色配什么色好看| 色吧在线观看| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 91麻豆精品激情在线观看国产| 精品欧美国产一区二区三| 在线国产一区二区在线| 九九久久精品国产亚洲av麻豆| 欧美xxxx黑人xx丫x性爽| 亚洲在线自拍视频| 精品一区二区三区视频在线观看免费| 久久精品国产自在天天线| 亚洲中文字幕日韩| 女同久久另类99精品国产91| 我的老师免费观看完整版| 最新中文字幕久久久久| 51午夜福利影视在线观看| 别揉我奶头 嗯啊视频| 日本黄色视频三级网站网址| 少妇裸体淫交视频免费看高清| 精品日产1卡2卡| 成人特级黄色片久久久久久久| 麻豆国产97在线/欧美| 久久精品国产99精品国产亚洲性色| 看十八女毛片水多多多| 久久精品综合一区二区三区| 级片在线观看| ponron亚洲| 精品久久久久久久久av| 久9热在线精品视频| 国产又黄又爽又无遮挡在线| 亚洲最大成人中文| 少妇高潮的动态图| 亚洲国产精品999在线| 国产一区二区在线av高清观看| 两个人视频免费观看高清| 一级作爱视频免费观看| 欧美黄色淫秽网站| 脱女人内裤的视频| 又紧又爽又黄一区二区| 国产老妇女一区| 亚洲欧美日韩无卡精品| 三级国产精品欧美在线观看| 搡老岳熟女国产| 麻豆国产97在线/欧美| 欧美在线一区亚洲| 午夜激情福利司机影院| 99热这里只有是精品50| 伦理电影大哥的女人| 亚洲av.av天堂| 狂野欧美白嫩少妇大欣赏| 91在线观看av| 欧美黄色片欧美黄色片| 婷婷丁香在线五月| 亚洲av成人不卡在线观看播放网| 午夜老司机福利剧场| 丁香六月欧美| 国产精品免费一区二区三区在线| 亚洲人成伊人成综合网2020| 国产精品1区2区在线观看.| 国产国拍精品亚洲av在线观看| 女人被狂操c到高潮| 在线天堂最新版资源| 久久久久久大精品| 久久国产精品人妻蜜桃| 中文字幕人成人乱码亚洲影| 午夜精品在线福利| 亚洲精品一卡2卡三卡4卡5卡| 九色成人免费人妻av| 久久精品国产99精品国产亚洲性色| 夜夜看夜夜爽夜夜摸| 免费看美女性在线毛片视频| 国产探花在线观看一区二区| 欧美bdsm另类| 在现免费观看毛片| 国产探花在线观看一区二区| 黄色日韩在线| 中文在线观看免费www的网站| 怎么达到女性高潮| av天堂中文字幕网| 日本与韩国留学比较| 长腿黑丝高跟| 男人和女人高潮做爰伦理| 国产高清三级在线| 亚洲真实伦在线观看| 有码 亚洲区| 熟女人妻精品中文字幕| 国内精品久久久久精免费| av天堂中文字幕网| 天堂网av新在线| 在线免费观看不下载黄p国产 | 国产一区二区在线观看日韩| 久久久精品欧美日韩精品| 超碰av人人做人人爽久久| 欧美日韩乱码在线| 美女高潮的动态| 亚洲国产精品合色在线| 国产精品一及| 亚洲人与动物交配视频| 免费av观看视频| 一a级毛片在线观看| 男插女下体视频免费在线播放| 窝窝影院91人妻| 亚洲欧美日韩无卡精品| 国模一区二区三区四区视频| 亚洲在线自拍视频| 又粗又爽又猛毛片免费看| 舔av片在线| 日本撒尿小便嘘嘘汇集6| 99久久99久久久精品蜜桃| 波多野结衣高清无吗| 国产一区二区在线av高清观看| 亚洲成人久久爱视频| 99久久精品一区二区三区| 日本免费一区二区三区高清不卡| 久久婷婷人人爽人人干人人爱| 99国产精品一区二区蜜桃av| 久久精品国产亚洲av天美| 亚洲国产高清在线一区二区三| 丰满人妻一区二区三区视频av| 最近最新免费中文字幕在线| 黄色丝袜av网址大全| 身体一侧抽搐| 观看美女的网站| 久久国产精品人妻蜜桃| 好男人电影高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 日本五十路高清| 女人十人毛片免费观看3o分钟| 国产在视频线在精品| 精品一区二区三区视频在线观看免费| 国产91精品成人一区二区三区| 岛国在线免费视频观看| 中出人妻视频一区二区| 亚洲一区二区三区不卡视频| 午夜福利在线在线| 露出奶头的视频| 在线看三级毛片| 久久久久国产精品人妻aⅴ院| 97人妻精品一区二区三区麻豆| 久久午夜福利片| 国产精品美女特级片免费视频播放器| 亚洲精品久久国产高清桃花| 中文字幕人妻熟人妻熟丝袜美| 夜夜看夜夜爽夜夜摸| 欧美黄色片欧美黄色片| 午夜免费激情av| 亚洲欧美日韩高清在线视频| 一进一出抽搐gif免费好疼| 麻豆国产97在线/欧美| 无人区码免费观看不卡| 中文字幕av在线有码专区| 此物有八面人人有两片| 亚洲最大成人av| 国产成人av教育| 91狼人影院| 日韩成人在线观看一区二区三区| 两个人视频免费观看高清| 国产成人啪精品午夜网站| 国产视频一区二区在线看| 国产精品一区二区三区四区久久| 欧美日韩福利视频一区二区| 国产在线男女| 国产午夜精品久久久久久一区二区三区 | 国产精品三级大全| 高潮久久久久久久久久久不卡| 九九久久精品国产亚洲av麻豆| 亚洲精品456在线播放app | 国产亚洲欧美在线一区二区| 可以在线观看毛片的网站| 黄片小视频在线播放| 亚洲欧美精品综合久久99| 此物有八面人人有两片| 国产精品嫩草影院av在线观看 | 婷婷精品国产亚洲av| 午夜亚洲福利在线播放| 内射极品少妇av片p| 国产一区二区亚洲精品在线观看| 精品国内亚洲2022精品成人| 天堂动漫精品| 99热这里只有是精品50| 亚洲av免费在线观看| 国产av一区在线观看免费| 午夜影院日韩av| 国产精品美女特级片免费视频播放器| 成人毛片a级毛片在线播放| 18禁裸乳无遮挡免费网站照片| 国产精品国产高清国产av| 午夜亚洲福利在线播放| 天天一区二区日本电影三级| 深爱激情五月婷婷| 91九色精品人成在线观看| 色综合站精品国产| 亚洲av中文字字幕乱码综合| 欧美性猛交黑人性爽| 国产色婷婷99| 亚洲自偷自拍三级| 免费在线观看日本一区| 长腿黑丝高跟| 丁香六月欧美| 九九热线精品视视频播放| 国产精品野战在线观看| 国产一区二区激情短视频| 免费观看精品视频网站| 99在线人妻在线中文字幕| 国内久久婷婷六月综合欲色啪| 人人妻人人澡欧美一区二区| 国产成人影院久久av| 国产精品亚洲美女久久久| 国产伦精品一区二区三区视频9| 成年人黄色毛片网站| 麻豆一二三区av精品| 美女大奶头视频| www.www免费av| 99久国产av精品| 不卡一级毛片| 欧美一区二区亚洲| 国产久久久一区二区三区| 我的老师免费观看完整版| 一个人免费在线观看电影| av在线老鸭窝| av中文乱码字幕在线| 久久亚洲真实| 人人妻人人看人人澡| 男女做爰动态图高潮gif福利片| bbb黄色大片| 欧美激情久久久久久爽电影| 国产av不卡久久| 亚洲欧美精品综合久久99| 国内精品美女久久久久久| 天美传媒精品一区二区| 脱女人内裤的视频| eeuss影院久久| 精品熟女少妇八av免费久了| 看十八女毛片水多多多| 国产亚洲欧美在线一区二区| 色综合欧美亚洲国产小说| 国产极品精品免费视频能看的| 亚洲国产欧洲综合997久久,| 日本撒尿小便嘘嘘汇集6| 99视频精品全部免费 在线| 中文字幕熟女人妻在线| 国产精品永久免费网站| av在线观看视频网站免费| 男插女下体视频免费在线播放| 欧美激情国产日韩精品一区| 深夜精品福利| 色av中文字幕| 我的老师免费观看完整版| 欧美在线一区亚洲| 赤兔流量卡办理| 高清日韩中文字幕在线| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 69av精品久久久久久| 啦啦啦韩国在线观看视频| 一区福利在线观看| 午夜久久久久精精品| 婷婷六月久久综合丁香| 欧美中文日本在线观看视频| 久久香蕉精品热| 久久久久久大精品| 日本在线视频免费播放| 自拍偷自拍亚洲精品老妇| 天天躁日日操中文字幕| 国产大屁股一区二区在线视频| 亚洲av中文字字幕乱码综合| 五月玫瑰六月丁香| 欧美中文日本在线观看视频| 99久久久亚洲精品蜜臀av| 国产伦人伦偷精品视频| 欧美最黄视频在线播放免费| 久久久久久久久大av| 欧美色欧美亚洲另类二区| 亚洲 欧美 日韩 在线 免费| 18禁在线播放成人免费| 国产在线精品亚洲第一网站| 日韩欧美精品v在线| 男人和女人高潮做爰伦理| 内射极品少妇av片p| 亚洲av免费在线观看| 内地一区二区视频在线| 人人妻,人人澡人人爽秒播| 亚洲精华国产精华精| 看免费av毛片| 3wmmmm亚洲av在线观看| 欧美成狂野欧美在线观看| 国产伦在线观看视频一区| 88av欧美| 亚洲不卡免费看| 日本与韩国留学比较| 欧美成狂野欧美在线观看| 欧美精品啪啪一区二区三区| 久久久精品欧美日韩精品| 久久久久久久亚洲中文字幕 | 国产精华一区二区三区| 内射极品少妇av片p| 国产精品日韩av在线免费观看| 国产av一区在线观看免费| 亚洲国产精品999在线| 久99久视频精品免费| 精品不卡国产一区二区三区| 我要搜黄色片| 亚洲av中文字字幕乱码综合| 国产高清视频在线观看网站| 国产久久久一区二区三区| 一区福利在线观看| 亚洲精品456在线播放app | 亚洲专区国产一区二区| 国产精品精品国产色婷婷| 最近最新中文字幕大全电影3| 久久久久亚洲av毛片大全| 亚洲欧美日韩卡通动漫| 欧美高清成人免费视频www| 在线免费观看的www视频| 国产成人欧美在线观看| 亚洲 国产 在线| 日韩亚洲欧美综合| 99视频精品全部免费 在线| 亚洲av熟女| 在线播放无遮挡| 有码 亚洲区| 男女那种视频在线观看| 亚洲激情在线av| 亚洲avbb在线观看| 无遮挡黄片免费观看| 免费在线观看成人毛片| 国产极品精品免费视频能看的| 青草久久国产| 中文亚洲av片在线观看爽| 日韩亚洲欧美综合| 免费看a级黄色片| 亚洲av电影不卡..在线观看| 日韩欧美三级三区| av中文乱码字幕在线| 午夜福利成人在线免费观看| 亚洲自偷自拍三级| 搡老熟女国产l中国老女人| 免费搜索国产男女视频| 国内精品美女久久久久久| 免费在线观看成人毛片| 草草在线视频免费看| 怎么达到女性高潮| 欧美日韩乱码在线| 中文亚洲av片在线观看爽| 亚洲人与动物交配视频| 少妇丰满av| 国产 一区 欧美 日韩| 变态另类丝袜制服| 天天一区二区日本电影三级| 国产亚洲av嫩草精品影院| 久久欧美精品欧美久久欧美| 日韩欧美在线二视频| 亚洲av成人av| 久久精品国产亚洲av涩爱 | 国产毛片a区久久久久| 午夜免费激情av| 97碰自拍视频| 国产久久久一区二区三区| 变态另类丝袜制服| 亚洲最大成人中文| 精品欧美国产一区二区三| 久久人人精品亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 变态另类丝袜制服| 欧美日韩瑟瑟在线播放| 欧美乱妇无乱码| 日本五十路高清| 亚洲国产精品久久男人天堂| 长腿黑丝高跟| 老司机深夜福利视频在线观看| 高清毛片免费观看视频网站| 熟女电影av网| 国产亚洲精品久久久com| 国产精品久久久久久久电影| 十八禁人妻一区二区| 在线a可以看的网站| 色播亚洲综合网| 亚洲国产精品合色在线| 91在线精品国自产拍蜜月| 国产精品99久久久久久久久| or卡值多少钱| 国产在线男女| 亚洲国产精品成人综合色| 每晚都被弄得嗷嗷叫到高潮| 国产精品美女特级片免费视频播放器| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 亚洲av五月六月丁香网| 夜夜爽天天搞| 欧美激情久久久久久爽电影| 亚洲欧美日韩卡通动漫| 一个人看的www免费观看视频| 国产真实乱freesex| 精品久久久久久久久亚洲 | 国产精品亚洲一级av第二区| 97超级碰碰碰精品色视频在线观看| 国产精品久久久久久精品电影| 99久久精品一区二区三区| 69av精品久久久久久| 一进一出抽搐gif免费好疼| 午夜福利成人在线免费观看| 国产成人福利小说| 看十八女毛片水多多多| 少妇人妻精品综合一区二区 | 99久久精品国产亚洲精品| 草草在线视频免费看| 小蜜桃在线观看免费完整版高清| 两个人视频免费观看高清| 每晚都被弄得嗷嗷叫到高潮| 精品99又大又爽又粗少妇毛片 | 国产成人av教育| 久久99热这里只有精品18| 亚洲中文日韩欧美视频| 色精品久久人妻99蜜桃| 国产av在哪里看| 国产 一区 欧美 日韩| 国内毛片毛片毛片毛片毛片| 国产午夜福利久久久久久| 麻豆av噜噜一区二区三区| 美女高潮喷水抽搐中文字幕| 亚洲最大成人中文| 精品熟女少妇八av免费久了| 草草在线视频免费看| 欧美高清成人免费视频www| 最后的刺客免费高清国语| 午夜福利18| 亚洲欧美日韩高清在线视频| 亚洲av熟女| 好男人电影高清在线观看| 在线播放国产精品三级| 舔av片在线| 欧美三级亚洲精品| 婷婷丁香在线五月| 色哟哟·www| 欧美最新免费一区二区三区 | 两个人的视频大全免费| 亚洲av电影不卡..在线观看| a级毛片免费高清观看在线播放| 久久国产精品影院| 亚洲一区高清亚洲精品| АⅤ资源中文在线天堂| 一个人看的www免费观看视频| 免费在线观看成人毛片| 深夜精品福利| 国产毛片a区久久久久| 三级毛片av免费| 丰满人妻一区二区三区视频av| 欧美绝顶高潮抽搐喷水| 国产人妻一区二区三区在| 日本 av在线| 在线观看舔阴道视频| 一本一本综合久久| av欧美777| 国产亚洲精品av在线| 十八禁人妻一区二区| 一级a爱片免费观看的视频| 亚洲av中文字字幕乱码综合| a在线观看视频网站| 18禁黄网站禁片免费观看直播| 成人一区二区视频在线观看| 午夜福利在线观看吧| 久久人人精品亚洲av| 欧美最黄视频在线播放免费| 色精品久久人妻99蜜桃| 亚洲成人精品中文字幕电影| 成人欧美大片| 十八禁人妻一区二区| 3wmmmm亚洲av在线观看| 久久久久久久久中文| 青草久久国产| 精品国产三级普通话版| 搡老熟女国产l中国老女人| 99在线视频只有这里精品首页| 久久99热这里只有精品18| 久久久久久久久久成人| 亚洲成人久久爱视频| 丰满人妻一区二区三区视频av| 欧美激情久久久久久爽电影| 久久久久性生活片|