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

    麻豆一二三区av精品| av在线观看视频网站免费| av在线观看视频网站免费| 黄色女人牲交| 十八禁网站免费在线| 成人特级av手机在线观看| 一级毛片久久久久久久久女| 一区二区三区免费毛片| 久9热在线精品视频| 桃色一区二区三区在线观看| 老女人水多毛片| 成人欧美大片| .国产精品久久| 日韩中文字幕欧美一区二区| 少妇被粗大猛烈的视频| 国产av不卡久久| 观看美女的网站| 欧美精品国产亚洲| 嫁个100分男人电影在线观看| 美女 人体艺术 gogo| 久久久久性生活片| 精品人妻一区二区三区麻豆 | 99视频精品全部免费 在线| 悠悠久久av| 亚洲精品成人久久久久久| 亚洲一级一片aⅴ在线观看| 丰满乱子伦码专区| 18禁裸乳无遮挡免费网站照片| bbb黄色大片| 成人亚洲精品av一区二区| 国产大屁股一区二区在线视频| 欧美精品啪啪一区二区三区| 12—13女人毛片做爰片一| 美女大奶头视频| 久久精品国产清高在天天线| 超碰av人人做人人爽久久| 亚洲五月天丁香| 一进一出抽搐动态| 成人二区视频| 波多野结衣巨乳人妻| 免费无遮挡裸体视频| 亚洲一区二区三区色噜噜| 日韩一本色道免费dvd| 国产成年人精品一区二区| 超碰av人人做人人爽久久| 好男人在线观看高清免费视频| 国产精品久久久久久久久免| 国产精品人妻久久久久久| 伦理电影大哥的女人| 不卡一级毛片| 他把我摸到了高潮在线观看| 在线观看美女被高潮喷水网站| 春色校园在线视频观看| 又爽又黄无遮挡网站| 搡老熟女国产l中国老女人| 亚洲中文字幕一区二区三区有码在线看| 国产精品无大码| 久久久久久久久大av| 亚洲七黄色美女视频| 在线免费十八禁| 午夜福利成人在线免费观看| 国产亚洲精品av在线| 嫩草影视91久久| 亚洲综合色惰| 一个人免费在线观看电影| 深夜精品福利| 日本a在线网址| 97热精品久久久久久| 无人区码免费观看不卡| 亚洲美女搞黄在线观看 | 久久99热6这里只有精品| 国产在线精品亚洲第一网站| .国产精品久久| 久久这里只有精品中国| 欧美成人性av电影在线观看| 欧美在线一区亚洲| 免费大片18禁| 欧美zozozo另类| 午夜爱爱视频在线播放| 欧美一级a爱片免费观看看| 99热这里只有是精品50| 动漫黄色视频在线观看| 黄色女人牲交| 欧美xxxx黑人xx丫x性爽| 久久久国产成人精品二区| 成人欧美大片| 男人舔奶头视频| 免费电影在线观看免费观看| 免费搜索国产男女视频| 亚洲久久久久久中文字幕| 搡老熟女国产l中国老女人| 日本黄大片高清| 亚洲精品成人久久久久久| a级毛片a级免费在线| 毛片女人毛片| 人妻丰满熟妇av一区二区三区| 白带黄色成豆腐渣| 日韩精品有码人妻一区| 成人午夜高清在线视频| 免费看日本二区| 成人综合一区亚洲| 小说图片视频综合网站| 一a级毛片在线观看| 国产一区二区在线av高清观看| 男插女下体视频免费在线播放| 国产单亲对白刺激| 国产探花极品一区二区| 精品人妻偷拍中文字幕| 国产精品免费一区二区三区在线| 欧美zozozo另类| 啦啦啦啦在线视频资源| 少妇的逼水好多| 一区二区三区激情视频| 国产精品久久久久久精品电影| 欧美性猛交黑人性爽| 给我免费播放毛片高清在线观看| 日韩欧美在线二视频| 搡老熟女国产l中国老女人| 亚洲专区中文字幕在线| 久久精品人妻少妇| 长腿黑丝高跟| 久久午夜亚洲精品久久| 亚洲第一电影网av| 人人妻人人看人人澡| 精品一区二区三区人妻视频| 国产精品国产三级国产av玫瑰| 麻豆国产97在线/欧美| 欧美性猛交╳xxx乱大交人| 成人无遮挡网站| 99riav亚洲国产免费| 国产单亲对白刺激| 少妇的逼水好多| 国产免费一级a男人的天堂| 丰满的人妻完整版| 中国美女看黄片| 欧美3d第一页| 人妻丰满熟妇av一区二区三区| .国产精品久久| 国产一区二区在线观看日韩| 亚洲最大成人av| 精品欧美国产一区二区三| 亚洲图色成人| 又紧又爽又黄一区二区| 色吧在线观看| 亚洲五月天丁香| 日韩欧美在线二视频| 九色国产91popny在线| 搞女人的毛片| eeuss影院久久| а√天堂www在线а√下载| 精品无人区乱码1区二区| 久久精品国产99精品国产亚洲性色| 真人做人爱边吃奶动态| 一进一出抽搐gif免费好疼| 99久久中文字幕三级久久日本| 女人十人毛片免费观看3o分钟| 少妇人妻精品综合一区二区 | 日本-黄色视频高清免费观看| 色吧在线观看| 一本一本综合久久| 中文在线观看免费www的网站| 黄色配什么色好看| 男女做爰动态图高潮gif福利片| 欧美性猛交╳xxx乱大交人| 亚洲欧美激情综合另类| 亚洲精华国产精华液的使用体验 | 国产伦在线观看视频一区| 国产三级中文精品| 变态另类丝袜制服| 午夜老司机福利剧场| 日本在线视频免费播放| 国产老妇女一区| 国产91精品成人一区二区三区| 日本 欧美在线| 欧美极品一区二区三区四区| av福利片在线观看| 国产精品免费一区二区三区在线| 国产精华一区二区三区| 久久久久久久亚洲中文字幕| 欧美日韩黄片免| 美女黄网站色视频| 久久久国产成人免费| av在线老鸭窝| 身体一侧抽搐| 日韩在线高清观看一区二区三区 | 免费观看精品视频网站| 国产在线精品亚洲第一网站| 日本 欧美在线| 午夜福利视频1000在线观看| 久久6这里有精品| 久久久成人免费电影| 草草在线视频免费看| 成人国产一区最新在线观看| 看十八女毛片水多多多| 一区二区三区激情视频| 成人av一区二区三区在线看| 97人妻精品一区二区三区麻豆| 哪里可以看免费的av片| 日韩高清综合在线| 久久人妻av系列| 国产精品永久免费网站| 免费无遮挡裸体视频| 欧美日韩精品成人综合77777| 12—13女人毛片做爰片一| 日日撸夜夜添| 中文字幕精品亚洲无线码一区| 欧美激情久久久久久爽电影| 久久精品国产自在天天线| 哪里可以看免费的av片| 婷婷色综合大香蕉| 色播亚洲综合网| 校园春色视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲18禁久久av| 在现免费观看毛片| 桃色一区二区三区在线观看| 亚洲av成人精品一区久久| 12—13女人毛片做爰片一| 精品久久久久久久久av| 国产高清有码在线观看视频| 成人鲁丝片一二三区免费| 亚洲中文日韩欧美视频| 尾随美女入室| 日韩欧美国产在线观看| 我要看日韩黄色一级片| 国内毛片毛片毛片毛片毛片| 99在线视频只有这里精品首页| 国产精品久久久久久久久免| 免费黄网站久久成人精品| 精品日产1卡2卡| 日韩 亚洲 欧美在线| 97超级碰碰碰精品色视频在线观看| 人人妻人人看人人澡| 精品人妻1区二区| a级毛片a级免费在线| 夜夜夜夜夜久久久久| 欧美人与善性xxx| www日本黄色视频网| 免费人成视频x8x8入口观看| 久久亚洲真实| 日韩亚洲欧美综合| 国产视频内射| 黄色欧美视频在线观看| 一级a爱片免费观看的视频| 欧美性猛交黑人性爽| 高清在线国产一区| 12—13女人毛片做爰片一| a级毛片免费高清观看在线播放| av女优亚洲男人天堂| 国产在线精品亚洲第一网站| av中文乱码字幕在线| 一个人免费在线观看电影| 少妇猛男粗大的猛烈进出视频 | www.www免费av| av在线老鸭窝| 美女被艹到高潮喷水动态| 精品人妻视频免费看| 女人十人毛片免费观看3o分钟| 在线播放国产精品三级| 91在线观看av| 久久久成人免费电影| 99久久久亚洲精品蜜臀av| 久久人人爽人人爽人人片va| 国产黄a三级三级三级人| or卡值多少钱| 制服丝袜大香蕉在线| 一卡2卡三卡四卡精品乱码亚洲| av福利片在线观看| 国产精品福利在线免费观看| av天堂在线播放| 一区福利在线观看| 精品欧美国产一区二区三| 国产精品伦人一区二区| 国产亚洲精品久久久久久毛片| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清专用| 免费看日本二区| 免费观看精品视频网站| 亚洲精品亚洲一区二区| 成人午夜高清在线视频| 亚洲avbb在线观看| 老熟妇仑乱视频hdxx| 一个人免费在线观看电影| 国产亚洲精品久久久com| 亚洲午夜理论影院| 欧美bdsm另类| 给我免费播放毛片高清在线观看| 18+在线观看网站| 一a级毛片在线观看| 国产亚洲av嫩草精品影院| 久久精品国产亚洲av天美| 一个人看的www免费观看视频| 亚洲av.av天堂| 亚洲成人久久爱视频| 内射极品少妇av片p| 国产蜜桃级精品一区二区三区| 精品99又大又爽又粗少妇毛片 | 97超视频在线观看视频| 免费看日本二区| 精品福利观看| 国国产精品蜜臀av免费| 日本在线视频免费播放| 日本五十路高清| 欧美黑人巨大hd| 国产探花在线观看一区二区| www日本黄色视频网| 亚洲欧美日韩高清专用| 国产伦在线观看视频一区| 高清日韩中文字幕在线| 真人一进一出gif抽搐免费| 成年免费大片在线观看| 在线a可以看的网站| 久久精品久久久久久噜噜老黄 | 亚洲狠狠婷婷综合久久图片| 99久久成人亚洲精品观看| 两个人的视频大全免费| 亚洲精品亚洲一区二区| 亚洲aⅴ乱码一区二区在线播放| 老师上课跳d突然被开到最大视频| 国产伦在线观看视频一区| 成人毛片a级毛片在线播放| 日本免费a在线| 国内揄拍国产精品人妻在线| 久9热在线精品视频| 人妻少妇偷人精品九色| 国产精品国产高清国产av| 国产黄a三级三级三级人| 日本精品一区二区三区蜜桃| 精品日产1卡2卡| 91久久精品电影网| 国产精品久久视频播放| 桃色一区二区三区在线观看| 国产极品精品免费视频能看的| 国产av不卡久久| 中出人妻视频一区二区| 国产高潮美女av| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 日本色播在线视频| 最近在线观看免费完整版| 可以在线观看毛片的网站| 精华霜和精华液先用哪个| 桃红色精品国产亚洲av| 久久精品国产亚洲av涩爱 | 午夜福利在线观看吧| 毛片一级片免费看久久久久 | 国产精品久久久久久久久免| 嫩草影视91久久| 精品午夜福利在线看| 亚洲欧美清纯卡通| 午夜日韩欧美国产| 午夜影院日韩av| 亚洲综合色惰| 国产69精品久久久久777片| 久久久久国产精品人妻aⅴ院| 国产真实伦视频高清在线观看 | 日本三级黄在线观看| 12—13女人毛片做爰片一| 久久6这里有精品| 国产av麻豆久久久久久久| 婷婷六月久久综合丁香| 欧美激情在线99| 十八禁国产超污无遮挡网站| 一级a爱片免费观看的视频| 我要搜黄色片| 欧美最黄视频在线播放免费| 在线观看一区二区三区| 国产淫片久久久久久久久| 久久久久久久久久成人| 久久婷婷人人爽人人干人人爱| 欧美色欧美亚洲另类二区| 日本-黄色视频高清免费观看| 日韩欧美免费精品| 91精品国产九色| 永久网站在线| 亚洲av熟女| 人人妻人人看人人澡| 欧美日韩综合久久久久久 | 99久久精品国产国产毛片| 在线国产一区二区在线| 91麻豆精品激情在线观看国产| 国产一区二区在线av高清观看| 最后的刺客免费高清国语| 亚洲男人的天堂狠狠| 九九在线视频观看精品| 美女黄网站色视频| 日本三级黄在线观看| 国产av在哪里看| 亚洲欧美日韩高清专用| 亚洲人成网站在线播放欧美日韩| 国产黄色小视频在线观看| 国内精品久久久久精免费| 日本 av在线| 99riav亚洲国产免费| 亚洲成a人片在线一区二区| 老司机深夜福利视频在线观看| 免费高清视频大片| 蜜桃久久精品国产亚洲av| 国产三级在线视频| 一本久久中文字幕| 国产成人a区在线观看| 国内精品久久久久精免费| 久久精品人妻少妇| 男女下面进入的视频免费午夜| 黄色视频,在线免费观看| 久久精品国产鲁丝片午夜精品 | АⅤ资源中文在线天堂| 能在线免费观看的黄片| 小说图片视频综合网站| 舔av片在线| 国产在视频线在精品| 中文字幕久久专区| 神马国产精品三级电影在线观看| 99riav亚洲国产免费| 国产毛片a区久久久久| 久久久久久久久久黄片| 女人被狂操c到高潮| 欧美日韩黄片免| 亚洲国产精品成人综合色| 别揉我奶头 嗯啊视频| 最好的美女福利视频网| 久9热在线精品视频| 国产精品人妻久久久久久| 日本精品一区二区三区蜜桃| 精品久久久久久久久久久久久| 国产高清视频在线播放一区| 久久九九热精品免费| h日本视频在线播放| 我要搜黄色片| 别揉我奶头 嗯啊视频| 亚洲欧美日韩无卡精品| 日本黄色片子视频| 精品一区二区三区人妻视频| 欧美在线一区亚洲| 亚洲av一区综合| 欧美潮喷喷水| 日日摸夜夜添夜夜添小说| 国产一区二区在线观看日韩| 亚洲欧美激情综合另类| 国产爱豆传媒在线观看| 俺也久久电影网| 美女被艹到高潮喷水动态| 国产精品国产三级国产av玫瑰| 又爽又黄无遮挡网站| 欧美日韩国产亚洲二区| 国产精品av视频在线免费观看| 18禁在线播放成人免费| 99国产精品一区二区蜜桃av| 成年免费大片在线观看| 日韩欧美国产一区二区入口| 国国产精品蜜臀av免费| 成人一区二区视频在线观看| 精品午夜福利在线看| 国产av一区在线观看免费| 国产精品亚洲一级av第二区| 国产精品久久久久久av不卡| 国内精品一区二区在线观看| 我的老师免费观看完整版| 亚洲精品粉嫩美女一区| 国产成人av教育| 全区人妻精品视频| 国产精品精品国产色婷婷| 人人妻人人看人人澡| 欧美一区二区国产精品久久精品| 欧美高清性xxxxhd video| 美女免费视频网站| 超碰av人人做人人爽久久| 一区二区三区免费毛片| 亚洲经典国产精华液单| 最近最新免费中文字幕在线| 欧美日韩国产亚洲二区| 免费大片18禁| 搡老熟女国产l中国老女人| 久久久久久久午夜电影| 可以在线观看毛片的网站| 精品午夜福利在线看| 国产在线精品亚洲第一网站| 国产一区二区亚洲精品在线观看| 在线国产一区二区在线| 亚洲av第一区精品v没综合| 国产探花极品一区二区| 性色avwww在线观看| 国产精品久久久久久久电影| 可以在线观看的亚洲视频| 午夜福利高清视频| 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 精品无人区乱码1区二区| 午夜免费成人在线视频| 一进一出抽搐动态| 色在线成人网| 麻豆av噜噜一区二区三区| 99国产极品粉嫩在线观看| 欧美高清性xxxxhd video| 男女下面进入的视频免费午夜| 精品无人区乱码1区二区| 午夜福利高清视频| 国产人妻一区二区三区在| 亚洲欧美日韩东京热| 亚洲无线在线观看| 99久久无色码亚洲精品果冻| 欧美日韩国产亚洲二区| 99九九线精品视频在线观看视频| 一级av片app| 五月伊人婷婷丁香| 日韩欧美在线二视频| 男女下面进入的视频免费午夜| 欧美成人免费av一区二区三区| 亚洲成人中文字幕在线播放| 夜夜夜夜夜久久久久| 国产精品乱码一区二三区的特点| 校园人妻丝袜中文字幕| 国产伦一二天堂av在线观看| 亚洲av美国av| 99久久成人亚洲精品观看| 色综合亚洲欧美另类图片| 亚洲第一电影网av| 国产精品乱码一区二三区的特点| 男女做爰动态图高潮gif福利片| 99久国产av精品| 久久久午夜欧美精品| 色视频www国产| 嫩草影视91久久| 桃色一区二区三区在线观看| 亚洲精品日韩av片在线观看| 91在线精品国自产拍蜜月| 色5月婷婷丁香| 国产国拍精品亚洲av在线观看| 91在线观看av| 免费观看人在逋| 亚洲性久久影院| 久久久久久久午夜电影| 亚洲成av人片在线播放无| 免费观看人在逋| 少妇熟女aⅴ在线视频| 日韩强制内射视频| av在线亚洲专区| 亚洲aⅴ乱码一区二区在线播放| 成人高潮视频无遮挡免费网站| 亚洲av第一区精品v没综合| av专区在线播放| 亚洲无线在线观看| 久久久久久久久久成人| 日本精品一区二区三区蜜桃| 色5月婷婷丁香| 亚洲国产精品成人综合色| h日本视频在线播放| 欧美xxxx黑人xx丫x性爽| 真人做人爱边吃奶动态| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 精品人妻偷拍中文字幕| 美女高潮喷水抽搐中文字幕| 99热精品在线国产| 亚洲精品粉嫩美女一区| 国产欧美日韩精品一区二区| 91久久精品国产一区二区成人| 九九在线视频观看精品| 九九久久精品国产亚洲av麻豆| 亚洲人与动物交配视频| 色综合色国产| 午夜免费激情av| 国产熟女欧美一区二区| 最近在线观看免费完整版| 中文字幕av成人在线电影| 在现免费观看毛片| av天堂在线播放| 一本一本综合久久| 亚洲精品乱码久久久v下载方式| 麻豆久久精品国产亚洲av| 久99久视频精品免费| 熟妇人妻久久中文字幕3abv| 乱人视频在线观看| 国产亚洲精品综合一区在线观看| 亚洲av免费在线观看| 成人一区二区视频在线观看| 亚洲最大成人av| a级一级毛片免费在线观看| 久久精品影院6| 午夜免费男女啪啪视频观看 | 变态另类成人亚洲欧美熟女| 99热只有精品国产| 99在线视频只有这里精品首页| 中文字幕av成人在线电影| 我要搜黄色片| 久久久久久久久中文| 超碰av人人做人人爽久久| 熟妇人妻久久中文字幕3abv| 国产精品三级大全| 成人永久免费在线观看视频| 一级a爱片免费观看的视频| 亚洲av第一区精品v没综合| 看免费成人av毛片| 少妇猛男粗大的猛烈进出视频 | 三级男女做爰猛烈吃奶摸视频| 色综合站精品国产| 国产色爽女视频免费观看| АⅤ资源中文在线天堂| 成人一区二区视频在线观看| 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av香蕉五月| 精品日产1卡2卡| 99热网站在线观看| 日韩亚洲欧美综合| 亚洲 国产 在线| 亚洲内射少妇av| 欧美高清成人免费视频www| 91久久精品国产一区二区三区| 国产午夜精品论理片| 亚洲精华国产精华精| 国产探花在线观看一区二区| a级毛片免费高清观看在线播放| 久久精品国产亚洲av涩爱 | 成年版毛片免费区| 深爱激情五月婷婷|