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

    Decay Analysis of Ge Isotopes Formed in Reactions Induced by Tightly and Loosely Bound Projectiles?

    2018-11-19 02:22:54AmandeepKaurKirandeepSandhuandManojSharma
    Communications in Theoretical Physics 2018年11期

    Amandeep Kaur,Kirandeep Sandhu,and Manoj K.Sharma

    1School of Physics and Materials Science,Thapar Institute of Engineering and Technology,Patiala 147004,Punjab,India

    2Department of Physics,GSSDGS Khalsa College,Patiala 147001,Punjab,India

    AbstractThe dynamical cluster-decay model(DCM)is employed to investigate the decay of68,70Ge?compound nuclei formed respectively via tightly(4He)and loosely(6He)bound projectiles,using64Zn target.The study is carried out over a wide energy range(Ec.m.~5 MeV to 16 MeV)by including the quadrupole deformations(β2i)and optimum orientationsof the decaying fragments.The fusion cross-sections,obtained by adding various evaporation channels show nice agreement with the experimental data for4He+64Zn reaction.The contribution from competing compound inelastic scattering channel is also analyzed particularly for68Ge?nucleus at above barrier energies.On the other hand,the decrement in the fusion cross-sections of70Ge?nuclear system is addressed by presuming that65Zn ER is formed via two different modes:(i)the αn evaporation of70Ge? nucleus,and(ii)1n-evaporation of66Zn? nuclear system,formed via breakup and 2n-transfer channels due to halo structure of the6He projectile.Besides this,the suppression in 2np evaporation cross-sections suggests the contribution of another breakup and transfer process of6He i.e.4He+64Zn.The contribution of breakup+transfer channels for6He+64Zn reaction is duly addressed by applying relevant energy corrections due to the breakup of“6He” projectile into 2n and4He.In addition to this,the barrier lowering,angular momentum and energy dependence effects are also explored in view of the dynamics of chosen reactions.

    Key words:tightly and loosely bound projectiles,evaporation residue,breakup and transfer

    1 Introduction

    Over the last few decades,the reactions induced by tightly and loosely bound projectiles have catched the extensive attention of the nuclear physics community.It is so because the nature of projectile plays very important role in the formation and decay process of the compound nucleus(CN).Generally,the CN formed via tightly bound projectile disintegrates through light particle emission or heavy mass fragments such as fusion- fission( ff)channels.The probability of each decay channel depends upon the mass of CN,excitation energy(E?CN)and related angular momentum(?)etc.However,when a loosely bound projectile strikes the target,the breakup and transfer processes also start contributing towards the fusion component.The loosely bound halo projectile gets disengaged in the field of target into two parts,from which one part of the projectile is transferred to the target to form an excited compound system,while the other escapes with relatively unchanged velocity.This breakup+transfer process leads to loss of flux from the fusion channel and thus suppress the cross-section,especially at energies in the vicinity of the Coulomb barrier.[1?2]In fact,the larger radius of loosely bound projectile decreases the Coulomb barrier,which leads to enhancement in the fusion cross-sections at near and sub barrier energies.On the contrary,at above barrier energies the breakup of loosely bound halo nucleons due to low binding energy,give rise to the reduction in complete fusion(CF)cross-sections,known as fusion suppression.[3?5]In the recent years,the reactions induced by unstable nuclei around the neutron drip line with mass A<20 are vigorously used to investigate various reaction mechanisms such as fusion,breakup and transfer etc.[6?9]Beside this,the total reaction cross-sections(σreac)are also studied extensively within the application of tightly bound projectiles,bombarding with wide range of target nuclei such as9Be to208Pb.[10?11]The decay of CN formed in such reactions,leads to various reaction channels of light particles(LPs)emission and heavier fragments such as fusion- fission.The possibility of the compound inelastic scattering may also be perceived for these reactions.[12]Interestingly,it is observed that the total sum of all open channel cross-sections(σreac)is identical to the cross-section calculated from analysis of elastic scattering angular distributions.[13]Therefore,it would be of interest to understand the structure effects and dynamics involved in reactions induced by both tightly and loosely bound nuclei.

    In view of the above discussion,the present work is carried out to analyze different decay mechanisms governed via tightly(4He)and loosely(6He)bound projectiles using same target nucleus(64Zn)within the framework of dynamical cluster-decay model(DCM).[14?19]In the recent experiments,[13,20?21]the reaction cross-sections have been measured for the decay of68Ge?and70Ge?compound nuclei formed respectively in4He+64Zn and6He+64Zn reactions,over the wide range of center-of-mass energies Ec.m.~(5–16)MeV.The total cross-sections for these reactions are considered to have the contribution from the complete fusion and other processes such as compound inelastic scattering,breakup,transfer etc.Note that,in both reactions the fusion process leads to the production of evaporation residues(ERs),hence σfus≈ σER,as the contribution of the ff-channel is not significant for4He+64Zn and6He+64Zn reactions.

    The68Ge?nucleus formed via4He beam,mainly deexcites by one-and two-particle evaporation channels like γ,n,p,2p,nα,pα,and 2α. Interestingly,oneparticle emission is the dominant for68Ge?nuclear system,whereas the contribution of nα,pα,and 2α is negligible.[13,20]The experimental excitation functions of these evaporation channels are addressed individually by optimizing the only parameter of the model,known as neck-length parameter?R.Moreover,another competing channel i.e.compound inelastic scattering for4He+65Zn reaction is also investigated,in particular at the above barrier energies.In addition to this,the fusion cross-sections are addressed for70Ge?nucleus by summing all the excitation functions of(2n)68Ge,(np)68Ga,(2np)67Ga and(nα)65Zn evaporation residues.The comparison between the calculated and experimental fusion cross-sections for6He+64Zn→70Ge?reaction suggests the contribution of breakup and transfer events arise due to the halo structure of the6He projectile.It may be noted that,65Zn ER is also formed via breakup and two neutron transfer process,in addition to the αn evaporation channel.[21]In addition to this,the suppression in 2np cross-sections is associated with the incomplete fusion of other part of6He projectile with the target i.e.6He+64Zn.Thus,an effort is made to estimate the breakup+transfer contribution by including the relevant energy corrections in reference to the breakup of6He projectile into4He and 2n,which is responsible for the enhancement of the total cross-sections in comparison to the ones calculated for4He+64Zn→68Ge?reaction.

    The outline of the paper is as follows:Section 2 gives a brief description of the dynamical cluster-decay model(DCM).The DCM based calculations for the4He+64Zn and6He+64Zn reactions at various centre-of-mass energies(Ec.m.)are given in Sec.3. The calculations are made for the choice of deformed fragments having(β2i)quadrupole deformations with“optimum”orientations taken from.[22]Finally,the results are summarized in Sec.4.

    2 Theoretical Framework

    Based on quantum mechanical fragmentation theory(QMFT),[23?25]the dynamical cluster-decay model(DCM)is used here to study the decay mechanism of hot and rotating nuclear systems.This model is worked out in terms of collective coordinates of mass and charge asymmetriesand ηZ=(Z1?,the relative separation R,the multipole deformations βλi(λ=2,3,4)and orientations θi,where i=1 and 2 stand for heavy and light fragments.For the decay of a excited nucleus,the η and R coordinates refer,respectively,to nucleon division or exchange between outgoing fragments and the transfer of kinetic energy of the incident channel(Ec.m.)to the internal excitation of the outgoing channel.

    The stationary Schr?dinger equation in η coordinate,at a fixed R=Ra(the first turning point of barrier penetration shown in Fig.1),is written as

    with ν =0,1,2,3,...referring to ground state(ν=0)and the excited statesolutions.The mass parameter Bηηare the smooth hydrodynamical masses of Kr?ger and Scheid.[26]By solving Eq.(1),we get the preformation probability P0,which imparts the important structure information of the decaying nucleus and reads as

    The fragmentation potential VR(η,T)used in stationary Schr?dinger Eq.(1)in η,at fixed R is defined as

    Here the VR(η,T)calculated in macro-microscopic method of Strutinsky,where VLDMis the temperature dependent liquid drop energy of Davidson et al.[27]with its constants refitted at T=0,[28]give the experimental binding energies of Audi et al.[29]or theoretical estimates of M?ller et al.[30]wherever not available in Ref.[29].The “empirical”shell correction(δU)of Myers and Swiatecki,[31]also made T-dependent to vanish exponentially at T0=1.5 MeV.[32]The other terms VC,VP,and V?are,respectively,the Coulomb,nuclear proximity,and centrifugal potentials for deformed and oriented nuclei.The Coulomb potential for multipole-multipole interaction between two separated nuclei is given by

    and the nuclear proximity potential for deformed,oriented nuclei employing temperature effects,reads as

    where b(T)is the nuclear surface thickness,γ is the surface energy constant,ˉR(T)is the mean curvature radius and Φ(s0(T))is the universal function(for more details see Refs.[33–34]).The angular momentum effects are included via the following expression

    with moment of inertia I(T)in completes ticking limit,i.e.,or,alternatively,the one calculated in non-sticking limit INS=μR2.Note that the alternative use of non-sticking limit corresponds to the supposition of prompt emission of fragments,generally used in the experimental determination of the ?-value.We find that the sticking limit(IS)used in the present work is more appropriate for the proximity potential,which inturn result in higher ?-values.Hence,the angular momentum values are larger for sticking moment of inertia and correspond to smaller separation distance?R,while the process is vice versa for INS.

    In terms of collective coordinates,the compound nucleus decay cross sections,for ?-partial waves is given by

    where the preformation probability P0(given in Eq.(2))refers to η-motion and the penetrability P defined later to R-motion.?maxis the maximum angular momentum,fixed here for the light particle cross-section approaching zero.The tunneling probability P is given by the Wentzel-Kramers-Brillouin(WKB)integral

    with V(Ra,T)=V(Rb,T)=TKE(T)=Qeff(T)for the two turning points.V(Ra,T)acts like an effective Q value of decay,Qeff(T),and TKE(T)as the total kinetic energy of decaying fragments.The turning point Raof barrier penetration is given as

    with radius vectors(i=1,2)

    and T-dependent nuclear radii R0iof the equivalent spherical nuclei read as,[35]

    The temperature T is related to the incoming center-of mass energy Ec.m.or the compound nucleus(CN)excitation energy E?CNvia the entrance channel Qinvalue asIn Eq.(9),the“?R”is the separation distance between the surfaces of decaying fragments,known to assimilate the neck formation effects,and hence referred as the neck-length parameter.This choice of neck criteria is very much similar to the one used in scission point[36?37]and saddle point[38?40]statistical models for addressing evaporation residue and fission decay modes.It is worth noting that the neck-length increases as a function of center-of-mass energy(Ec.m.)(see Fig.5(a)).The neck values lie below 2.5 fm due to the use of proximity potential.The necklength ?R or,equivalently,the first turning point Rais the only parameter of the model,and inturn to define the effective “barrier-lowering” parameter ?VB(?)as the difference between the actually used barrier V(Ra,?),and the calculate barrier VB(?),as

    ?VB(?)represents a negative quantity,owing to the fact that penetration can occur below the barrier only.

    3 Results and Discussions

    In the present work,we intend to study various decay paths of68,70Ge?nuclear systems using the dynamical cluster decay model(DCM)[14?19]at center-of-mass energies(Ec.m.~5 MeV to 16 MeV)lying in the vicinity of the Coulomb barrier.First,in Subsec.3.1,the calculations are made for the evaporation residues of the68Ge?compound nucleus formed via interaction of tightly bound4He beam on64Zn target.In the subsequent analysis,the dynamics of6He+64Zn→70Ge?reaction is explored in Subsec.3.2,by investigating the fusion,breakup and transfer channels,where breakup+transfer events arise due to the halo structure of6He projectile.Finally,a comparative study of tightly and loosely bound projectile induced nuclear reactions is carried out to analyze the influence of halo structure of the projectile on the considered reaction dynamics.

    3.1 Decay of68Ge?Nucleus Formed Using Tightly Bound4He Beam

    The interaction of alpha particle with64Zn target leads to the formation of68Ge?compound system,which deexcites via various evaporation channels such as γ,p,n,2p,nα,pα,and 2α.[20]It is suggested in Refs.[13,20]that within Ec.m.~ (5–16)MeV energy range,the one-particle(p and n)evaporation channels are dominant,whereas negligible contribution of the two-or multi-particle decay modes is observed.Moreover,the γ-channel is also not included in the present calculations due to negligible cross-sections,contributing towards the total reaction cross-sections.[13]Apart from these decay events,the contribution from another competing channel,i.e.,compound inelastic scattering is also investigated particularly at higher centre-of-mass energies.Hence in the present section,an attempt is made to address the dynamics of4He+64Zn reaction by understanding various decay modes of68Ge?nucleus.

    In order to understand the dynamics involved through the one particle emission,the barrier characteristics of68Ge?→ p+67Ga and68Ge?→ n+67Ge decay channels are plotted in Fig.1(a)at highest incident energy Ec.m.=15.17 MeV and extreme angular momentum states ?=10~and 41~.Interestingly,both the decay channels show different barrier characteristics at different?-states.Clearly,at ?min=10~,the barrier height(VB)is higher for proton-emission in comparison to the neutrondecay,and is maintained upto ?=22~as depicted from Fig.1(b).However,these results show reverse trend at higher angular momentum states up to ?max=41~.This may be due to the interplay between the various potentials contributing towards the barrier height.In order to explore this,Fig.2(a)is plotted,which shows the variation of the centrifugal potential for proton and neutron emissions at different angular momentum states. The centrifugal potential remains higher for neutron-decay for majority of ?-states.The same is true for the proximity potential,but not shown here.However,with the addition of the Coulomb term in Fig.2(b),the magnitude of the potential for the neutron emission decreases for lower ?-values,leading to higher barrier height for the proton decay.Elaborating this,the larger barrier height for the proton emission at lowest ?-state is mainly due to the dominance of Coulomb term at lower angular momentum states,where the suppression in the magnitude of the centrifugal potential is simultaneously observed.At lower ?-states, ?(?+1) ~ ?2term in the centrifugal potential,given by Eq.(6)does not contribute much,hence the dominance of Coulomb factor can be felt for the proton decay.However,the Coulomb potential remains zero for the neutron emission resulting higher barrier for the proton at ?min=10~state.In contrast to this,the magnitude of the centrifugal potential competes the Coulomb term and gradually becomes larger in magnitude for the higher ?-values.Due to less influence of the Coulomb term at higher ?-states,the barrier height gets suppressed for the proton emission(see Fig.1(b)).In other words the Coulomb and centrifugal potentials play an important role to define the barrier characteristics of the proton and neutron decays.Note that the emergence of the flip due to different barrier characteristics for proton and neutron emission channels also effect the penetration path(see Fig.2(c)).In addition to the barrier height,barrier width also plays an important role to describe the penetration of the fragment through the barrier,as described by Eq.(8)of Sec.2.Depending on these factors,the penetrability P shows different behavior for p-and n-decay channels,plotted in Fig.2(c).Specifically,the variation of decay barrier widthcan be depicted from the inset of Fig.2(c)to illustrate the fact that penetrability increases exponentially with decrease in barrier width and vice versa.However,we notice that P increases with increase in ? for both decay channels and in turn sets the minimum value of angular momentum,slightly different for both particles.

    Fig.1 (a)Scattering potential V(R)for68Ge? → p+67Ga and68Ge? → n+67Ge decay channels at Ec.m.=15.17 MeV,calculated for two extreme ?-values, and ?max=41~.(b)Variation of DCM-calculated decay barrier height(VB)as a function of angular momentum ?(~)for both emissions.

    In the next step,to investigate the fragmentation path of68Ge?nucleus via p-and n-particle evaporation channels,the mass-fragmentation potential V(A2)is shown in Figs.3(a)and 3(b)at Ec.m.=15.17 MeV.Note that the fragmentation potential goes as an input to the Schr?dinger equation,which is used to estimate the preformation probability P0,and hence the excitation functions of the nuclei.The results are obtained for both the decay channels using different separation distances(?R)chosen to fit the respective experimental data.In order to see the comparative contribution of energetically feasible fragments,the fragmentation paths are plotted for?minand ?maxvalues.One may notice from Figs.3(a)and 3(b)that the structure of potential energy surfaces(PES)are almost independent at two extreme ?-values for the fission fragments(FF),but the characteristic behavior of fragmentation potential for light particles(LPs)is significantly different.At lower ?-values,the emission of light particles is favorable that means contribution of pand n-particles is prominent,whereas fission fragments(marked in figure as FF)start competing the light particles at higher ?-states.Following this,an attempt is made to predict the fission cross-sections for the decay of68Ge?nucleus within the framework of DCM,which may be verified by future experiment(s)(refer Table 2).

    Fig.2 (a)Variation of the centrifugal potential V?for the proton and neutron emissions with angular momentum ?(~)at highest center-of-mass energy Ec.m.=15.17 MeV.(b)Same as that of(a),but plotted for Coulomb+centrifugal potentials.(c)Variation of DCM-calculated barrier penetrability(P)as a function of angular momentum ?(~)for both emissions.The inset of(c)shows the behavior of barrier width varying with angular momentum.

    After analyzing the ?-dependence of fragmentation decay path for both decay channels,the effect of center-of mass energy on the fragmentation potential(V),preformation probability(P0)and barrier penetrability(P)is illustrated in Figs.4(a),4(b),and 4(c)respectively,at two common ?=15~and 41~angular momentum states.These calculations are done at common neck-length parameter?R,chosen to fit the p-evaporation channel crosssections.Clearly,the magnitude of fragmentation potential increases slightly with center-of-mass energies Ec.m.for p-and n-particle emissions.Interestingly,the fragmentation potentials for both particles overlaps at the ?=15~angular momentum state,indicating the fact that neutron and proton emissions are equally probable at lower angular momentum states.On the other hand,for higher angular momentum the fragmentation behavior of the pand n-decay channels differs significantly.For example at ?=41~,the fragmentation potential is much higher for the neutron-emission as compared to proton-decay.Lower fragmentation potential for the proton-decay seems to suggest that the proton emission is more favorable towards decay at highest angular momentum state.This result is in accordance with the observation of Fig.1(a),where lower barrier is reported for the proton emission at ?max=41~.This difference in fragmentation potential of n-and p-decay channels at higher ?-values enters primarily via rotational energy term due to difference in moment of inertia,as different daughter nuclei are produced in these two channels.It may be noted that67Ga(the daughter nucleus in p-channel)is prolate deformed with β2=0.228,whereas67Ge(the daughter nucleus in n-channel)is oblately deformed with β2= ?0.299.As their radii are changed significantly,so moment of inertia changes accordingly.It is important to mention here that the minima in the fragmentation potential corresponds to maxima in preformation probability and vice versa.Thus,the preformation factor P0shows opposite trend than that of the fragmentation potential in reference to the ?-value.One can see from Fig.4(b),at lower ?-value,both n-and p-particles have almost same and maximum preformation probability,which in turn,represents that LPs are favorable at lower ?-states for all energies.Whereas,at higher?,the P0has much lower magnitude,this is in agreement with Fig.3(a)for the fragmentation case.Next,the variation of penetrability P is defined via Fig.4(c),it shows similar results as shown by fragmentation potential,owing maximum magnitude at higher ?-state,and indicating the different barrier characteristics for both the particle emissions.Also,barrier penetrability(P)shows increment with increase in center-of-mass energy particularly at lower angular momentum.

    Fig.3 Mass-fragmentation potential V(A2)at centre-of-mass energy Ec.m.=15.17 MeV for(a)64Zn(4He,p)67Ga and(b)64Zn(4He,n)67Ge decay channels,using best- fit ?R’s,plotted for ?minand ?maxvalues.

    Fig.4 The energy(Ec.m.)dependence of(a)fragmentation potential(V),(b)preformation probability(P0)and(c)penetrability(P)for proton and neutron evaporation channels at two ?-values,15~and 41~.

    Using the combined effect of P0and P values,the ER cross-sections for p-and n-decay channels are calculated within the framework of DCM,presented in Tables 1–2,and compared with experimental data.[13,20]The DCM calculated results show nice agreement with the experimental measurements.The only parameter of the model is the neck-length(?R),which decides the first turning point of the penetration path,obtained for best- fit of the measured data,listed in the tables and plotted in Fig.5(a).The optimized neck-length parameter“?R”is observed to have larger magnitude for n-channel as compared to that of p-channel.It may be due to the reason that68Ge appears below the β stability line and hence seems more prone for p-emission,which is in agreement with the result depicted in Fig.4(a).As discussed earlier,the fission fragments(FF)may also contribute in the decay of68Ge?nucleus,particularly at higher ?-states.Thus,in this work an effort is made to predict the fission contribution pre-sented in Table 2,using?RFF= ?Rn?0.5 fm.It is relevant to mention that lighter fragments are generally emitted first from the hot and rotating nuclei,establishing the fact that light particle emission occurs prior to the fission process.Keeping this in mind,higher?R is considered for ER formation as compared to the fission decay.Similar analysis is carried out in the previous work[41]of DCM,based on the decay of Xe-nucleus.At different neck-length parameters“?RFF”,the predicted excitation functions of fission fragments(A2=28–34 and their complementary fragments)are found to be negligibly small,indicating the fact that Ge-nucleus decays primarily via light particles,and fission fragments contribution is negligibly small.Moreover,if we calculate FF cross-sections at?RFF=?Rp?0.5 fm i.e.at the neck-length parameter of proton emission,the fission contribution reduces further.Hence the collective clusterization does not reinforce the fragmentation to follow the fission path,instead ER formation is preferred.In addition to this,the“barrier modification” (?VB)another observable of DCM,which can be related to the?R(equivalently,Ra),is examined in Fig.5(b)as a function of Ec.m..Apparently,?VB(in magnitude)for each particle emission is maximum for lowest Ec.m.and minimum at the highest one.It is also evident from figure that larger barrier modification is needed for the case of proton-emission in comparison with neutron one.

    Table 1 The calculated decay cross-sections for p-evaporation channels using DCM at various Ec.m.for4He+64Zn reaction and compared with experimental data.[13,20]The other related observables like neck-length parameter(?R)and temperature(T)are also listed.

    One-particle emission is dominant decay mode of68Ge?nucleus,as explained in Sec.1 of the present work.However,experiments[13,20]also suggest that two-particle emission channels play significant role at energies above 20 MeV due to the large threshold of multi-nucleon clusters.Although,in this work we are dealing with energies below 20 MeV,but an effort is made here to check the relative contribution of two-particle emission for the decay of68Ge?nucleus formed in4He+64Zn reaction.The cross-sections of68Ge?→2p+66Zn and68Ge?→ αn+63Zn decays are calculated at two extreme energies Ec.m.=5.79 MeV and 15.17 MeV.The DCM calculated crosssections are comes out to be σ2p=1.77 × 10?4mb and σαn=1.35 × 10?2mb at lowest Ec.m.=5.79 MeV,which increases with increase in centre-of-mass energy and read as σ2p=1.27 mb and σαn=15.17 mb at highest Ec.m.=15.17 MeV.These cross-sections are calculated at same?R as observed for p-emission,because the previous study[18]of DCM indicates almost similar magnitude of?R for one-particle and multi-nucleon emissions.Hence it is evident from above discussion that the estimated crosssections for 2p and αn do not contribute much towards the total reaction cross-sections,as reported in the experimental work.[20]

    Table 2 Same as Table 1 but for n-evaporation channel.In addition to this,the predicted fission contributions(σFF)are also listed.

    Fig.5 Variation of(a)the neck-length parameter(?R)and(b)the barrier modification(?VB)as a function of Ec.m.for the decay of68Ge?compound system via neutron(solid squares)and proton(open circles)decay channels.

    Apart from one-and two-particle evaporation channels,the compound inelastic scattering is another competing channel emerging in the decay of68Ge?compound system especially at above barrier energies.In this process,the projectile is absorbed by the target nucleus and leads to formation of compound nucleus(CN),which further re-emits the incident particle with low kinetic energy,thus original nucleus remains in the excited state.Here,we have studied the compound inelastic scattering for4He+64Zn →68Ge?→4He′+64Zn reaction.For inelastic scattering process,the entrance channel keeps its identity,and hence the preformation probability is taken P0=1 for(4He,4He′)decay channel,and penetrability is calculated using WKB approximation and presented in Fig.6.It is evident from Fig.6 that larger penetration probability is obtained at highest center-of-mass energy,indicating easier penetration of the projectile like fragments.Besides this,the penetrability corresponding to each ?-state increases with increase in angular momentum.Using P0and P,the measured[13]cross-sections of(81±13.8)mb and(201±37.1)mb are best fitted at?R=0.699 fm and 0.827 fm values respectively,for two highest energies Ec.m.=11.29 MeV and 15.17 MeV,at which signature of the inelastic scattering is observed.It is relevant to mention that,if energy loss up to 30%is considered in the inelastic scattering process then the neck-length parameter“?R”and angular momentum effect remain almost unchanged.Hence,negligible change in ?R and ?maxvalues is observed even with the reduction of 30%in the energy of incident flux.

    Fig.6 Barrier penetrability P as a function of ? for 64Zn(4He,4He′)64Zn reaction plotted at two center-of mass energies.

    3.2 Decay of70Ge?Formed via Loosely Bound Projectile6He

    After analyzing the dynamics of tightly bound4Heinduced reaction,we intend to explore the decay of excited70Ge?nuclear system formed via interaction of two neutron halo “6He” projectile and64Zn target.Interestingly,in recent experiment,[21]6He+64Zn→70Ge?reaction has been studied at various energies across the barrier ranging Ec.m.=(7.54 to 16.1)MeV.The total fusion cross-sections(σfus)have been obtained by measuring the excitation-functions of various evaporation residues(ERs)(2n)68Ge,(np)68Ga,(2np)67Ga,and(nα)65Zn,as light charged particle emission is the primary decay mode of70Ge?nuclear system. Apart from ERs,the breakup and transfer events also contribute to the total crosssections due to loosely bound nature of the6He projectile.Hence in the present section,the evaporation residue,breakup+transfer cross-sections(or equivalently incomplete fusion cross-sections)are addressed using DCM,in order to have overall understanding of the dynamics of6He+64Zn→70Ge?reaction so that loosely bound nature of6He nucleus be explored.The detailed analysis of the above work is carried out in the subsequent discussion.

    In order to investigate the decay of70Ge?nucleus,the fragmentation potential V(?)and preformation probability P0(?)are presented in Fig.7 for the evaporation residues of70Ge?nucleus,corresponding to energy Ec.m.=12.7 MeV,above the Coulomb barrier.The different ER decay modes,which are shown in figure represent multi-nucleon clusters and heavy residual nuclei,i.e.,X+Y→2n+68Ge(solid squares),np+68Ga(open squares),2np+67Ga(solid triangles)and nα+65Zn(open triangles).Clearly,the fragmentation potentials for all ER channels show increment with increase in angular momentum.The result reflects that light charged particle emission is probable at small values of angular momentum states,due to lower fragmentation potential and larger preformation probability,as depicted from Figs.7(a)and 7(b)respectively.It is relevant to mention that similar results are extracted for the decay of68Ge?compound nucleus,where the evaporation residue formation dominates at lower angular momentum states.In addition to this,Fig.7(a)shows that at lower ?-values,the fragmentation potential increases with increase in mass of the multi-particle cluster,i.e.,maximum for nα and minimum for 2n cluster,whereas reverse trend is observed at higher angular momentum values.As explained earlier in Subsec.3.1,this happens because of the relative contribution of Coulomb,proximity and centrifugal potentials to the total fragmentation potential.Due to this,an flip in the magnitude of the fragmentation potential arises around ?=30~.Similar results are obtained from Fig.7(b),plotted for the preformation probabilities of the decay fragments.The preformation probability shows increment up to certain ?-states,then sharply decreases near the ?max-values of the cluster emission.The maximum values of the angular momentum can be estimated from Fig.7(b),designated as a point where the contribution of the preformation probability is negligible.Above the?max-values the total cross-sections becomes significantly small.

    In reference to above discussion,the angular momentum window ?min≤ ?≤ ?maxis fixed and the excitation functions for the evaporation residues of70Ge?nucleus are estimated using Eq.(7).The ?minfor all the decay fragments lies around 0~,and ?max-values are pointed in Fig.7(b),which are different for different clusters.The fitted cross-sections are tabulated in Table 3 for all the mentioned clusters.The best- fit values of neck-length parameters(?R)for individual ER channels along with temperature(T)and center-of-mass energies(Ec.m.)are also listed in Table 3.Different?R values indicate unique reaction time scales for individual evaporation residue channels.It is of interest to note that the total calculated fusion crosssectionsof70Ge?compound system,i.e.,sum of individual ER excitation functions,deviate from the experimental measured fusion cross-sectionsIt is clearly depicted from Fig.8(a)and 8(b)that the 2n and np clusters show nice agreement with the experimental data.However,the strong mismatch of the experimental and our DCM based cross-sections is obtained for the(nα)65Zn-cluster(see Fig.8(c)),which is responsible for the decrement in the fusion cross-sections and signify the large contribution of breakup and transfer channels.Additionally,the cross-sections of 2np-cluster are also underestimated at above barrier energies as shown in Fig.8(d).The results remain intact even at the maximum allowed neck-length parameters.Hence the wide range of?R could not address the experimental crosssections upto~ 2.5 fm.The suppression in the nα and 2np cross-sections is governed mainly through the loosely bound nature of the6He projectile,responsible for the breakup and transfer events.The loosely bound6He nucleus breaks up into two parts4He and 2n,which act as new projectiles and interact with the64Zn target nucleus individually.It may be noted that65Zn nucleus is not formed only via nα evaporation channel,but also find its genesis in the breakup and two-neutron transfer channels of6He projectile,which is addressed in the subsequent discussion,by employing the relevant energy correction.Also,the suppression in 2np cross-sections is associated with the incomplete fusion of other part of6He projectile with the target i.e.4He+64Zn.In reference to the above discussion,the contribution of breakup and transfer channels is included in the present analysis as the difference between estimated and experimental cross-sections,which are read as

    Note that the efficacy of DCM in reference to breakup of loosely bound projectile and interaction of its one part with the target has been analysed in previous studies.[18?19]The 2n and4He now act as new projectiles and interact with the64Zn target nucleus individually,and lead to the formation of66Zn?and68Ge?nuclear systems,which further decays via n+65Zn and 2p+66Zn channels respectively.In the above process,the energy correction is required to address the breakup+transfer process,because the interaction energy must be modified according to the charge and mass number of the new projectiles.The corrected energy(Ec)of new projectile is obtained as[19]

    Here Apis the mass of projectile before breakup and Aeis mass of ejectile,which is an uninteracted outgoing fragment. Hence by including the relevant energy corrections,we look for the breakup and transfer processes and subsequently the decay of66Zn?and68Ge?nuclear systems.The modified center-of-mass energiesand respective excitation energiesare listed in Table 4.The corrected energies for breakup channels are smaller in magnitude as compared to the fusion process.

    Fig.7 The fragmentation potential(V)and preformation probability P0vs.? plotted for evaporation residues(2n)68Ge,(np)68Ga,(2np)67Ga and(nα)65Zn for the de-excitation of70Ge? system formed via using loosely bound nucleus as a projectile at energy 12.7 MeV.

    Accordingly,the cross-sections are calculated for 1n decay of66Zn?and 2p decay of68Ge?via breakup and transfer reactions.The temperature(T)and best- fit?R values are reported in Table 4 and both increase with increase in energy.The angular momentum contributing to the breakup+transfer cross-sections has the limiting valuesand.The total sum of DCM calculated breakup+transfer and fusion cross-sections show nice agreement with the experimental fusion data(refer Fig.9).Broadly speaking,the suppression in the total fusion cross-sections of70Ge?CN due to nα and 2np channels is compensated through the breakup+transfer events leading to66Zn?and68Ge?nuclear systems,followed by 1n and 2p emission respectively.Conclusively,the loosely bound nature of the6He projectile influences the dynamics of6He+64Zn reaction as the total contribution of breakup+transfer channel is reasonably higher than that of the fusion cross-sections.

    Table 3 The DCM calculated cross-sections for individual ER 2n,np,2np,and nα channels with best- fit values of?R for the decay of70Ge?at various centre-of-mass energies.obtained by summing all ERs cross-sections and compared with experimental fusion data.[21]

    Table 3 The DCM calculated cross-sections for individual ER 2n,np,2np,and nα channels with best- fit values of?R for the decay of70Ge?at various centre-of-mass energies.obtained by summing all ERs cross-sections and compared with experimental fusion data.[21]

    Ec.m.E?CNT?R2nσ2n(68Ge)?Rnpσnp(68Ga)?R2npσ2np(67Ga)?Rnασnα(65Zn)σDCM fusσExpt.fus/MeV /MeV /MeV /fm /mb /fm /mb /fm /mb /fm /mb /mb /mb 7.54 30.20 2.036 1.592 9.97 - - - - 1.349 14.2 24.17 128.43 9.15 31.81 2.088 1.860 44.6 - - 1.433 3.37 1.447 18.7 66.67 283.07 10.60 33.26 2.133 1.980 70.1 1.890 106 1.700 16.2 1.650 32.8 225.1 501.08 11.83 34.49 2.171 2.140 116 1.980 139 1.910 34.0 2.020 53.8 342.8 553.85 12.70 35.36 2.197 2.300 165 2.070 174 2.120 55.1 2.070 51.7 445.8 809.01 13.80 36.46 2.230 2.306 164 2.030 150 2.350 72.9 2.150 49.0 435.9 817.81 15.00 37.66 2.266 2.300 153 2.170 209 2.400 72.2 2.200 45.9 480.1 863.29 16.10 38.76 2.298 2.260 141 1.965 114 2.430 71.5 2.250 42.0 368.5 781.04

    Fig.8 Comparison of experimental cross-sections with calculated ones for(a)2n,(b)np,(c)nα,and(d)2np evaporation channels at various Ec.m.’s for the6He+64Zn reaction.

    Finally,the comparison of4He+64Zn and6He+64Zn reactions is investigated in Fig.9.As mentioned above,the fusion cross-sections of nucleus formed in4He+64Zn reaction matches the experimental data nicely.On the other hand,suppression in the DCM calculated fusion cross-sections are observed for6He+64Zn reaction channel,due to mismatch in the cross-sections of nα and 2npclusters. If one compares the fusion cross-sections for4He+64Zn and6He+64Zn reactions,sub-barrier enhancement in the fusion cross-sections of6He+64Zn reaction is observed.However,the DCM fitted σfusfor6He+64Zn reaction bear lower magnitude at above barrier energies than the same calculated for4He+64Zn.These two results can be explained on the basis of structure properties of loosely bound projectile.The larger radius of the loosely bound projectile reduces the Coulomb barrier,as a result one can expect the fusion enhancement,as shown in Fig.9 at sub-barrier energies.

    On the other hand,the low binding energy of the outer two neutrons of the6He nucleus implies large contribution of breakup+transfer channels,which could remove flux from the complete fusion and thus suppress the fusion cross-sections.The DCM results support that at the above-barrier energies,the loss of flux due to the breakup+transfer channels dominates and the resulting fusion cross-section is suppressed.Note that similar kind of results have been analyzed for the reactions of loosely bound projectile with heavy target(e.g.11Be+208Pb[7]).Hence from the above discussion the loosely bound nature of the6He can be easily illustrated by concluding following results:(i)The loosely bound character of the projectile is responsible for the breakup+transfer channels,which are addressed here through the two different breakup and transfer channels of6He i.e.2n+64Zn and4He+64Zn.(ii)The comparison of the tightly and loosely bound projectile induced reactions reveals the peculiar dynamics of6He+64Zn interacting pair.It has been emphasized that the loosely bound nature of6He is responsible for the sub-barrier fusion enhancement and above barrier fusion suppression.

    Table 4 The calculated breakup+transfer cross-sections using DCM for 1n decay channel for 66Zn? excited system formed in the 2n+64Zn reaction,and 2p decay channel for68Ge? nucleus formed in4He+64Zn reaction at modified energies

    Table 4 The calculated breakup+transfer cross-sections using DCM for 1n decay channel for 66Zn? excited system formed in the 2n+64Zn reaction,and 2p decay channel for68Ge? nucleus formed in4He+64Zn reaction at modified energies

    Ecc.m./MeVE?CN/MeVT/MeV?R/fmσ/mb 2n+64Zn→66Zn?→1n+65Zn 2.670 22.08 1.805 1.265 89 3.229 22.64 1.827 1.340 162 3.746 23.16 1.847 1.398 241 3.824 23.24 1.850 1.366 177 4.489 23.90 1.875 1.462 349 4.878 24.92 1.889 1.472 345 5.304 24.72 1.905 1.473 320 5.690 25.10 1.920 1.475 305 4He+64Zn→68Ge? → 2p+66Zn 9.470 12.82 1.370 1.570 17.5 10.30 13.65 1.412 2.060 76.0 11.05 14.40 1.448 2.450 112.0

    Fig.9 Variation of cross-sections as a function of centreof-mass energies for4He+64Zn and6He+64Zn reactions and compared with experimental data.[13,20?21]

    4 Summary

    In this work,we have investigated the decay of68Ge?and70Ge?nuclear systems formed via tightly bound“4He” and loosely bound “6He” projectiles over a wide range of center-of-mass energies Ec.m..The DCM calculated cross-sections find nice agreement with the experimental data for all the decay channels of68Ge?nucleus.The p-and n-emissions,which are dominant decay modes of68Ge?nucleus,show different barrier characteristics with respect to angular momentum(?).Also,the Ec.m.-dependence and ?-dependence of fragmentation potential,preformation probability and barrier penetrability are duly addressed.Moreover,the fission contribution for the decay of68Ge?has been predicted,which comes out to be negligibly small possibly due to lighter mass and lower excitation energies of compound system.Apart from this,another competing channel i.e.compound inelastic scattering is also analyzed at above barrier energies.

    The excitation functions of various evaporation residues(ERs)(2n)68Ge,(np)68Ga,(2np)67Gaand(nα)65Zn for the decay of70Ge?system are estimated and found that the σERfor(nα)65Zn and(2np)67Ga channel deviates from the experimental data.This deviation suggests the breakup of6He projectile(i.e.6He→2n+4He)and transfer of 2n and4He leads to the formation of66Zn?and68Ge?nuclei.The breakup+transfer cross-sections have been calculated by applying relevant energy corrections. The total sum of fusion and breakup+transfer cross-sections show nice agreement with experimentally measured fusion data.

    Lastly,in order to understand the influence of the halostructure of6He,a comparative study of tightly bound“4He” and loosely bound “6He”projectile induced reactions has been carried out.A lower Coulomb barrier and smaller binding energy of loosely bound nucleus in turn influences the fusion cross-sections and related reaction mechanisms.Due to the loosely bound nature of6He nucleus,the sub-barrier enhancement and above-barrier suppression are observed in the fusion cross-sections of6He+64Zn reaction.

    亚洲精华国产精华精| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 精品国内亚洲2022精品成人| 草草在线视频免费看| 中出人妻视频一区二区| 村上凉子中文字幕在线| 久久久精品国产亚洲av高清涩受| 欧美激情久久久久久爽电影| 久久香蕉国产精品| 亚洲精品av麻豆狂野| 在线免费观看的www视频| 色av中文字幕| 国产精品久久久久久亚洲av鲁大| 亚洲狠狠婷婷综合久久图片| 老司机午夜福利在线观看视频| 999久久久精品免费观看国产| 俄罗斯特黄特色一大片| 久久婷婷人人爽人人干人人爱| 18禁观看日本| 又爽又黄无遮挡网站| 欧美日韩亚洲综合一区二区三区_| 一级黄色大片毛片| 在线观看午夜福利视频| 91老司机精品| 日本 av在线| 人妻久久中文字幕网| 久久久久亚洲av毛片大全| 国产真实乱freesex| 精品国内亚洲2022精品成人| 日韩欧美国产一区二区入口| 97超级碰碰碰精品色视频在线观看| 少妇被粗大的猛进出69影院| 看黄色毛片网站| 欧美一级毛片孕妇| cao死你这个sao货| 我要搜黄色片| 国产乱人伦免费视频| 99久久无色码亚洲精品果冻| 岛国在线观看网站| 午夜福利在线观看吧| 熟女电影av网| 麻豆国产av国片精品| 在线观看日韩欧美| 欧美成狂野欧美在线观看| 亚洲色图 男人天堂 中文字幕| 欧美色视频一区免费| 国产麻豆成人av免费视频| 一进一出抽搐gif免费好疼| 1024手机看黄色片| 亚洲精品在线美女| 成人av在线播放网站| 最好的美女福利视频网| 两人在一起打扑克的视频| 久久久久久人人人人人| 国产亚洲av嫩草精品影院| 12—13女人毛片做爰片一| 在线观看日韩欧美| 久久精品91蜜桃| www.www免费av| 午夜视频精品福利| 在线观看免费午夜福利视频| 欧美一区二区精品小视频在线| 首页视频小说图片口味搜索| 国产成人系列免费观看| 99久久综合精品五月天人人| 亚洲片人在线观看| www.熟女人妻精品国产| 久久香蕉国产精品| 国产1区2区3区精品| 国产精品爽爽va在线观看网站| 一本大道久久a久久精品| 亚洲精华国产精华精| 亚洲欧美日韩东京热| 人妻夜夜爽99麻豆av| 欧美日韩福利视频一区二区| 91字幕亚洲| 精品久久久久久成人av| 久久婷婷人人爽人人干人人爱| 免费看十八禁软件| 黄色 视频免费看| 欧美3d第一页| 国产三级在线视频| 国产精品久久久久久人妻精品电影| 两个人看的免费小视频| 一区二区三区激情视频| 后天国语完整版免费观看| 天天躁夜夜躁狠狠躁躁| 桃色一区二区三区在线观看| 日本 av在线| 日韩欧美在线乱码| 亚洲国产欧美人成| 成人欧美大片| 一个人免费在线观看电影 | 色综合亚洲欧美另类图片| 中出人妻视频一区二区| 宅男免费午夜| 欧美高清成人免费视频www| 精品不卡国产一区二区三区| 久久天躁狠狠躁夜夜2o2o| 在线观看免费视频日本深夜| 久久午夜亚洲精品久久| 国产真人三级小视频在线观看| 久久精品aⅴ一区二区三区四区| 亚洲精品国产精品久久久不卡| av超薄肉色丝袜交足视频| 国模一区二区三区四区视频 | 久久久久国内视频| 性欧美人与动物交配| 免费在线观看日本一区| 国产精品亚洲美女久久久| 国产一区二区三区视频了| 女人高潮潮喷娇喘18禁视频| 五月伊人婷婷丁香| 久久精品国产清高在天天线| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品成人综合色| 国产97色在线日韩免费| 国产激情久久老熟女| 国产精品电影一区二区三区| 精品福利观看| 一级黄色大片毛片| 美女午夜性视频免费| 国产高清视频在线观看网站| av在线天堂中文字幕| 免费在线观看完整版高清| 一区福利在线观看| 亚洲精品粉嫩美女一区| 一本久久中文字幕| 91大片在线观看| 久久久久久九九精品二区国产 | 国产精品av视频在线免费观看| 国产三级黄色录像| a级毛片在线看网站| 熟女少妇亚洲综合色aaa.| 国产久久久一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲人成网站在线播放欧美日韩| 啦啦啦免费观看视频1| 亚洲男人天堂网一区| 村上凉子中文字幕在线| 欧美日韩一级在线毛片| 久久香蕉激情| 国产精品 欧美亚洲| 午夜日韩欧美国产| 又黄又粗又硬又大视频| 欧美日韩瑟瑟在线播放| 色噜噜av男人的天堂激情| 他把我摸到了高潮在线观看| 香蕉av资源在线| 国产免费av片在线观看野外av| 91字幕亚洲| 成人国产一区最新在线观看| 午夜视频精品福利| 在线免费观看的www视频| 久久久国产成人免费| 成人亚洲精品av一区二区| 91老司机精品| 黑人巨大精品欧美一区二区mp4| 久久久久久人人人人人| 亚洲欧美精品综合一区二区三区| 天堂影院成人在线观看| 丰满人妻熟妇乱又伦精品不卡| 黑人操中国人逼视频| 亚洲成人久久性| 丰满人妻一区二区三区视频av | 国内久久婷婷六月综合欲色啪| 亚洲国产精品合色在线| 午夜精品久久久久久毛片777| 可以免费在线观看a视频的电影网站| 两性夫妻黄色片| 亚洲人成伊人成综合网2020| 天堂av国产一区二区熟女人妻 | 在线视频色国产色| 99精品欧美一区二区三区四区| 日本精品一区二区三区蜜桃| 极品教师在线免费播放| av福利片在线| 成人精品一区二区免费| 国内少妇人妻偷人精品xxx网站 | 97碰自拍视频| 久久精品人妻少妇| 18禁国产床啪视频网站| 欧美成人性av电影在线观看| av福利片在线| 国产激情偷乱视频一区二区| 精品国产超薄肉色丝袜足j| 色在线成人网| 国产蜜桃级精品一区二区三区| 亚洲国产欧洲综合997久久,| 成人午夜高清在线视频| 别揉我奶头~嗯~啊~动态视频| 久9热在线精品视频| 国产av一区在线观看免费| 免费在线观看视频国产中文字幕亚洲| 曰老女人黄片| 国产三级在线视频| 久久亚洲精品不卡| 久久人妻福利社区极品人妻图片| 男女视频在线观看网站免费 | 国产精品亚洲一级av第二区| av中文乱码字幕在线| 久久精品国产清高在天天线| 欧美一区二区国产精品久久精品 | 国产激情久久老熟女| 精品久久久久久久人妻蜜臀av| 亚洲专区国产一区二区| 国产精品98久久久久久宅男小说| 国产精品影院久久| 一区福利在线观看| 国产精品亚洲一级av第二区| 男人舔女人下体高潮全视频| 在线免费观看的www视频| 999精品在线视频| av在线播放免费不卡| 久久久水蜜桃国产精品网| 男插女下体视频免费在线播放| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| 色哟哟哟哟哟哟| 91麻豆精品激情在线观看国产| 麻豆国产av国片精品| 久久久久久久久免费视频了| 老汉色∧v一级毛片| 床上黄色一级片| 成人手机av| 国产久久久一区二区三区| 777久久人妻少妇嫩草av网站| 久久精品aⅴ一区二区三区四区| 欧美黑人巨大hd| 欧美人与性动交α欧美精品济南到| 99热这里只有精品一区 | 久久久久性生活片| 国产人伦9x9x在线观看| 午夜两性在线视频| 香蕉av资源在线| a级毛片在线看网站| 精品午夜福利视频在线观看一区| 免费在线观看视频国产中文字幕亚洲| 久久精品人妻少妇| 一个人免费在线观看的高清视频| 俺也久久电影网| 午夜老司机福利片| 波多野结衣巨乳人妻| 午夜福利在线在线| 成在线人永久免费视频| 无人区码免费观看不卡| 国产1区2区3区精品| 国产精品影院久久| 视频区欧美日本亚洲| 天天添夜夜摸| 久99久视频精品免费| 欧美一级a爱片免费观看看 | 午夜福利免费观看在线| 午夜免费激情av| 日韩欧美在线二视频| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品影院| 每晚都被弄得嗷嗷叫到高潮| 午夜影院日韩av| 国产午夜精品久久久久久| 这个男人来自地球电影免费观看| 美女 人体艺术 gogo| 国产又黄又爽又无遮挡在线| 免费在线观看日本一区| 亚洲午夜精品一区,二区,三区| 午夜免费成人在线视频| 久久中文字幕人妻熟女| 国产欧美日韩一区二区精品| 在线观看免费视频日本深夜| 午夜福利免费观看在线| 91老司机精品| 日韩精品免费视频一区二区三区| 国产97色在线日韩免费| 国产精品免费视频内射| 一夜夜www| 50天的宝宝边吃奶边哭怎么回事| 亚洲美女视频黄频| 欧美成狂野欧美在线观看| 看免费av毛片| 男女床上黄色一级片免费看| 亚洲色图 男人天堂 中文字幕| 99久久99久久久精品蜜桃| 亚洲精品一区av在线观看| 久久久久九九精品影院| www日本在线高清视频| 99热只有精品国产| 观看免费一级毛片| 99re在线观看精品视频| av片东京热男人的天堂| 老熟妇乱子伦视频在线观看| 黄频高清免费视频| 日韩大尺度精品在线看网址| 天天躁夜夜躁狠狠躁躁| 制服诱惑二区| 777久久人妻少妇嫩草av网站| 午夜福利免费观看在线| 国产av麻豆久久久久久久| 亚洲精品久久成人aⅴ小说| 欧美日韩一级在线毛片| 高清毛片免费观看视频网站| 国内精品久久久久精免费| 午夜a级毛片| 精品一区二区三区四区五区乱码| 日韩免费av在线播放| 麻豆av在线久日| 香蕉久久夜色| 日本精品一区二区三区蜜桃| √禁漫天堂资源中文www| 在线观看www视频免费| 五月玫瑰六月丁香| 中文字幕最新亚洲高清| 日韩成人在线观看一区二区三区| 高清毛片免费观看视频网站| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| 色综合欧美亚洲国产小说| 亚洲国产高清在线一区二区三| 亚洲精品中文字幕一二三四区| 亚洲人成网站高清观看| 国产一区二区在线观看日韩 | 亚洲国产精品sss在线观看| 欧美一区二区国产精品久久精品 | 欧美一区二区精品小视频在线| 国产亚洲精品久久久久5区| 亚洲第一电影网av| 久久久久久九九精品二区国产 | 国产精品98久久久久久宅男小说| 露出奶头的视频| 久久久精品大字幕| 亚洲欧美精品综合一区二区三区| 亚洲精品中文字幕一二三四区| 久久 成人 亚洲| 日韩精品免费视频一区二区三区| 欧美一级毛片孕妇| 久久精品91蜜桃| 欧美日韩乱码在线| 久久午夜综合久久蜜桃| 国产免费av片在线观看野外av| 成年免费大片在线观看| 日本黄大片高清| 最近最新免费中文字幕在线| 免费在线观看视频国产中文字幕亚洲| 一个人免费在线观看的高清视频| 天堂av国产一区二区熟女人妻 | 成年免费大片在线观看| 精品人妻1区二区| 18禁裸乳无遮挡免费网站照片| 欧美 亚洲 国产 日韩一| 露出奶头的视频| 久久国产精品人妻蜜桃| 国产高清激情床上av| 国产亚洲欧美98| 久久久久久免费高清国产稀缺| 免费人成视频x8x8入口观看| 欧美日韩国产亚洲二区| 99在线人妻在线中文字幕| 成年免费大片在线观看| 波多野结衣高清无吗| 这个男人来自地球电影免费观看| 亚洲精品av麻豆狂野| 在线永久观看黄色视频| 亚洲av电影不卡..在线观看| 久久久久久九九精品二区国产 | 麻豆久久精品国产亚洲av| 亚洲国产欧美人成| 亚洲国产精品久久男人天堂| 国产野战对白在线观看| 国产麻豆成人av免费视频| 又黄又爽又免费观看的视频| 18禁美女被吸乳视频| 两个人看的免费小视频| 午夜福利18| 男女视频在线观看网站免费 | 免费搜索国产男女视频| 亚洲欧美日韩东京热| 草草在线视频免费看| www.熟女人妻精品国产| 男女做爰动态图高潮gif福利片| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 国产精品 国内视频| 亚洲最大成人中文| 久久中文字幕一级| 国产精品野战在线观看| 久久久久久久午夜电影| 一级毛片女人18水好多| 九色成人免费人妻av| 久久精品国产99精品国产亚洲性色| 午夜福利18| 国产视频内射| 日韩欧美国产一区二区入口| 非洲黑人性xxxx精品又粗又长| 国产成+人综合+亚洲专区| 亚洲美女黄片视频| 久久草成人影院| 又粗又爽又猛毛片免费看| 久久久国产成人免费| 久久久久免费精品人妻一区二区| 91麻豆精品激情在线观看国产| www国产在线视频色| 国产高清视频在线播放一区| 99国产综合亚洲精品| 欧美大码av| 美女免费视频网站| 国产成人啪精品午夜网站| 久久这里只有精品中国| 国产成人av激情在线播放| 十八禁网站免费在线| 日日夜夜操网爽| 老汉色av国产亚洲站长工具| 波多野结衣巨乳人妻| 九九热线精品视视频播放| 91九色精品人成在线观看| 露出奶头的视频| 成人永久免费在线观看视频| av在线天堂中文字幕| 一区二区三区国产精品乱码| 露出奶头的视频| 亚洲av片天天在线观看| 舔av片在线| 男男h啪啪无遮挡| 视频区欧美日本亚洲| 亚洲午夜理论影院| 亚洲av电影在线进入| 成人精品一区二区免费| 99国产综合亚洲精品| 国产免费男女视频| 一夜夜www| 国产不卡一卡二| 这个男人来自地球电影免费观看| 激情在线观看视频在线高清| 99riav亚洲国产免费| av免费在线观看网站| 黄片大片在线免费观看| 国产乱人伦免费视频| 中亚洲国语对白在线视频| 久久中文字幕一级| 精品欧美一区二区三区在线| av国产免费在线观看| 狂野欧美白嫩少妇大欣赏| 成人国产一区最新在线观看| 免费在线观看成人毛片| ponron亚洲| 欧美人与性动交α欧美精品济南到| 香蕉丝袜av| 亚洲欧美日韩无卡精品| 舔av片在线| 亚洲国产精品久久男人天堂| 日韩欧美国产一区二区入口| 两性夫妻黄色片| 日韩精品中文字幕看吧| 19禁男女啪啪无遮挡网站| 成人av在线播放网站| 国产片内射在线| 久久人人精品亚洲av| 高清毛片免费观看视频网站| 亚洲成av人片在线播放无| 久久精品国产清高在天天线| 国产精品亚洲av一区麻豆| 悠悠久久av| 久久精品国产亚洲av高清一级| 一本精品99久久精品77| 久久久久久九九精品二区国产 | 人妻丰满熟妇av一区二区三区| 亚洲精品在线观看二区| 成人国语在线视频| 美女免费视频网站| www.熟女人妻精品国产| 国产精品一区二区三区四区久久| 欧洲精品卡2卡3卡4卡5卡区| 香蕉av资源在线| 久久热在线av| 黄色视频不卡| netflix在线观看网站| 日日干狠狠操夜夜爽| 国产黄片美女视频| 高清毛片免费观看视频网站| 人人妻,人人澡人人爽秒播| 亚洲欧美激情综合另类| 国产av一区二区精品久久| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费一区二区三区在线| 免费在线观看日本一区| 亚洲无线在线观看| 国产亚洲精品综合一区在线观看 | 亚洲色图av天堂| 精品久久久久久久末码| 手机成人av网站| 高清毛片免费观看视频网站| 少妇熟女aⅴ在线视频| 色播亚洲综合网| 夜夜看夜夜爽夜夜摸| 别揉我奶头~嗯~啊~动态视频| 日韩免费av在线播放| 国产成人一区二区三区免费视频网站| 成人三级黄色视频| www国产在线视频色| 黄片小视频在线播放| 正在播放国产对白刺激| 亚洲av电影在线进入| 亚洲精品国产一区二区精华液| 久久天堂一区二区三区四区| 亚洲av中文字字幕乱码综合| 免费看美女性在线毛片视频| 小说图片视频综合网站| 最近最新中文字幕大全免费视频| 日韩中文字幕欧美一区二区| 两性夫妻黄色片| 国产乱人伦免费视频| 男插女下体视频免费在线播放| 精品人妻1区二区| 亚洲一卡2卡3卡4卡5卡精品中文| av中文乱码字幕在线| 久久久精品国产亚洲av高清涩受| 无遮挡黄片免费观看| 欧美极品一区二区三区四区| 宅男免费午夜| 一区二区三区国产精品乱码| av视频在线观看入口| 国产精品1区2区在线观看.| 久久性视频一级片| 悠悠久久av| 国产av又大| 日韩欧美 国产精品| 操出白浆在线播放| 亚洲自拍偷在线| 禁无遮挡网站| 婷婷精品国产亚洲av在线| 亚洲欧美激情综合另类| 国产精品免费一区二区三区在线| 男女那种视频在线观看| 18禁国产床啪视频网站| 国产又黄又爽又无遮挡在线| av有码第一页| 777久久人妻少妇嫩草av网站| 精品欧美国产一区二区三| 日韩成人在线观看一区二区三区| 成人18禁在线播放| 国产一区二区激情短视频| 国产伦在线观看视频一区| 国产精品久久久久久人妻精品电影| 免费在线观看影片大全网站| 国产真实乱freesex| 亚洲精品久久成人aⅴ小说| 啦啦啦观看免费观看视频高清| 国产精品自产拍在线观看55亚洲| 99国产精品一区二区三区| 男插女下体视频免费在线播放| 一进一出抽搐gif免费好疼| 国产成人欧美在线观看| 欧美黄色片欧美黄色片| 国产精品亚洲美女久久久| 欧美zozozo另类| 妹子高潮喷水视频| 国内精品久久久久久久电影| 国产精品一及| 亚洲av中文字字幕乱码综合| av超薄肉色丝袜交足视频| 午夜福利高清视频| 久久伊人香网站| 欧美黄色片欧美黄色片| 可以在线观看毛片的网站| 制服丝袜大香蕉在线| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| 欧美在线黄色| 日韩大尺度精品在线看网址| 久久人妻福利社区极品人妻图片| 变态另类成人亚洲欧美熟女| 亚洲人与动物交配视频| 国产精品久久久久久精品电影| 国产激情偷乱视频一区二区| 人成视频在线观看免费观看| 亚洲中文字幕日韩| 久久精品国产99精品国产亚洲性色| 精品国产美女av久久久久小说| 午夜日韩欧美国产| 一个人观看的视频www高清免费观看 | 国内毛片毛片毛片毛片毛片| 国产精品乱码一区二三区的特点| 麻豆成人午夜福利视频| 国产精品一区二区精品视频观看| 欧美成狂野欧美在线观看| 伊人久久大香线蕉亚洲五| 女警被强在线播放| 变态另类成人亚洲欧美熟女| 欧美黑人欧美精品刺激| 国产午夜精品论理片| 成在线人永久免费视频| 亚洲欧美日韩东京热| 国产精品一及| 丁香六月欧美| 两个人视频免费观看高清| 免费在线观看视频国产中文字幕亚洲| 最近在线观看免费完整版| ponron亚洲| 亚洲avbb在线观看| 久久人人精品亚洲av| 狂野欧美激情性xxxx| 国产精品,欧美在线| 国产成人精品久久二区二区免费| 99国产精品一区二区蜜桃av| 久久久久久人人人人人| 丰满人妻熟妇乱又伦精品不卡| 免费观看人在逋| 在线a可以看的网站| 女人被狂操c到高潮| 琪琪午夜伦伦电影理论片6080| 精品少妇一区二区三区视频日本电影| 久久久久久亚洲精品国产蜜桃av| 亚洲精品粉嫩美女一区| 欧美精品亚洲一区二区| 成人国语在线视频| 欧美一区二区国产精品久久精品 |