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

    The Gamow shell model with realistic interactions: a theoretical framework for ab initio nuclear structure at drip-lines

    2022-10-22 08:15:56NicolasMichelJianguoLiandFurongXu
    Communications in Theoretical Physics 2022年9期

    Nicolas Michel,Jianguo Li and Furong Xu

    1 CAS Key Laboratory of High Precision Nuclear Spectroscopy,Institute of Modern Physics,Chinese Academy of Sciences,Lanzhou 730000,China

    2 School of Nuclear Science and Technology,University of Chinese Academy of Sciences,Beijing,China

    3 School of Physics,and State Key Laboratory of Nuclear Physics and Technology,Peking University,Beijing 100871,China

    Abstract Ab initio approaches are among the most advanced models to solve the nuclear many-body problem.In particular,the no-core—shell model and many-body perturbation theory have been recently extended to the Gamow shell model framework,where the harmonic oscillator basis is replaced by a basis bearing bound,resonance and scattering states,i.e.the Berggren basis.As continuum coupling is included at basis level and as configuration mixing takes care of internucleon correlations,halo and resonance nuclei can be properly described with the Gamow shell model.The development of the no-core Gamow shell model and the introduction of the -box method in the Gamow shell model,as well as their first ab initio applications,will be reviewed in this paper.Peculiarities compared to models using harmonic oscillator bases will be shortly described.The current power and limitations of ab initio Gamow shell model will also be discussed,as well as its potential for future applications.

    Keywords: binding energy,lifetimes,nucleon—nucleon interactions,shell model,continuum coupling

    1.Introduction

    The most advanced models to describe atomic nuclei are those where all nucleons are active and where all inter-particle correlations are included.These methods are called ab initio,because only the form of the nuclear Hamiltonian is postulated at first so that nuclear many-body states are exact eigensolutions of the associated Schr?dinger equation in principle.In practice,this means that many-body eigenstates can be calculated up to a margin that can be specified during the construction of the Hamiltonian.It has been possible to develop ab initio nuclear models because present theoretical interactions between nucleons are nowadays sufficiently accurate to reproduce the physical observables of few-body systems,where the three-body force can be included [1,2].Indeed,three-body forces are inherent of the nuclear interaction.They are responsible for the saturation of nuclear matter and they have been shown to exist at several orders of any effective field theory for nucleon.

    However,despite the relatively direct and simple theoretical nature of ab initio nuclear models,their applications in numerical calculations are far from obvious.The most important problem is the hard core repulsion between nucleons,which prevents from using an independent-particle wave function as a starting point.A practical solution is to use equivalent lowenergy inter-nucleon interactions,which remove the highenergy components of the nuclear Hamiltonian,while still reproducing the two-nucleon observables up to about 300 MeV[1,2].Nevertheless,model space dimensions grow very rapidly with the number of valence nucleons,as their increase is combinatorial.In fact,only nuclei bearing fewer than about 15 nucleons can be numerically treated with nocore—shell model (NCSM) [3].

    Due to these theoretical and practical conundrums,it has then become clear that another route should be followed if one aims at considering heavier nuclei with ab initio shell model methods.A powerful idea is not to consider the initial Hamiltonian,but,on the contrary,to transform it so that a relatively small number of configurations is necessary to diagonalize it.One can cite for that matter the in-medium similarity renormalization group(IMSRG),which decouples particle-particle and hole-hole configurations so that a simple Hartree—Fock (HF) or Hartree—Fock—Bogoliubov ansatz can give a quasi-exact energy at a given flow parameter with the renormalized Hamiltonian [4—8].Another method for transforming the initial Hamiltonian in order to be able to treat heavier nuclei is the many-body perturbation theory(MBPT),where an effective Hamiltonian is calculated perturbatively from the initial realistic Hamiltonian in small model spaces [9,10].

    Even though the latter numerical approaches allowed to greatly improve the number of nuclei that can be treated with ab initio methods,they still rely on the use of bases of square-integrable states,so that a precise calculation of halo and resonance nuclei is precluded.An alternative ab initio theoretical approach must then be implemented,which includes continuum coupling.Hence,one has generalized the Gamow shell model (GSM) [11]to ab initio frameworks,namely the no-core Gamow shell model (NCGSM)[12]and GSM with-box method [13].Indeed,GSM allows for the description of weakly bound and unbound states,because continuum coupling is included at basis level,while inter-nucleon correlations are present via configuration mixing [11].

    Several calculations have been made in ab initio GSM and it is the object of this review to give a short description of its most important applications.For this,we will shortly present the basics of the one-body Berggren basis expansion in section 2,as it is fundamental to formulate GSM.We will then focus on the numerical methods used to diagonalize the GSM matrix in section 3,as the numerical problems met therein are different from those of the shell model using the harmonic oscillator basis.Then,we will quickly present the theory of NCGSM in section 4,and its first applications afterward,i.e.the unbound nucleus5He in section 5,A=3,4 unbound nuclei in section 6,and the oxygen and fluorine chains in section 7,where the NCGSM Hamiltonian is transformed via MBPT to obtain GSM with-box method.A conclusion will be made in section 8.

    Figure 1.Representation of the complex k-plane,showing the positions of bound states,well bound or of halo character (empty disks) and resonances,narrow or broad (grey filled disks).L+ is the contour representing the non-resonant continuum.The Berggren completeness relation involves the bound states,resonance states lying between L+and the real k-axis,and the scattering states on L+(from [35]).

    2.Berggren basis and completeness relations

    Many-body resonance wave functions are expanded on a basis of complex-energy independent-particle states in GSM.Thus,the completeness relation of complex-energy one-body eigenstates is fundamental and will be firstly depicted.For this,one has to consider the standard real-energy completeness relation generated by a one-body Hermitian Hamiltonian,i.e.the Newton completeness relation [14]:

    where un(r) are the bound states and real-energy scattering states uk(r) run over the positive real-k axis.

    Let us now deform the integration contour of equation(1)into the complex k-plane in order for complex-energy states to explicitly appear,as shown in figure 1.

    Following the complex integration Cauchy theorem,one obtains:

    where knare the poles ofuk(r)uk(r′)lying between the real kaxis and the L+complex contour.The Berggren completeness relation [15]follows immediately:

    where un(r) are now the bound states and resonances present between the real k-axis and the L+contour now consists of complex-energy scattering states uk(r)defined in the complex k-plane.All bound,resonance and scattering states un(r) and uk(r) are generated by the same basis potential.

    Equations (1),(3) are two completeness relations differing by way of their domain of definition.Indeed,equation(1)is similar to the Fourier transform as only bound states can expanded from the Newton completeness relation [14].Conversely,equation (3) allows for the expansion of bound states and resonance states whose complex momentum is situated between the L+contour and the real k-axis [15].

    In order to be used in numerical calculations,equation(3)must be discretized,which is most efficiently performed with the Gaussian quadrature [11]:

    where N is the total number of used resonant and discretized scattering states,uki(r)=un(r)if it is a resonant state anduki(r)=if it is a discretized scattering state of momentum ki=k and Gaussian weight wk.The renormalization of discretized scattering states,depending on Gaussian weight,is necessary for one-body bound,resonance and scattering states to be treated on the same footing [11].

    Equation (4) is formally identical to the standard completeness relation in a discrete basis so that Slater determinants of the many-body Berggren basis are straightforward to generate.However,as the formalism of Gamow states is non-Hermitian,the obtained Hamiltonian matrix H is complex symmetric.New numerical methods then had to be devised to solve the eigenproblem specific to the GSM Hamiltonian matrix.

    3.Numerical implementation of the Gamow shell model

    3.1.Calculation of two-body matrix elements with the Berggren basis

    Two-body matrix elements are straightforward to calculate with a basis of harmonic oscillator states.For this,one firstly evaluates nuclear matrix elements in relative coordinates with harmonic oscillator states,as the nuclear Hamiltonian is defined with the relative coordinates of the interacting nucleons.Then,one calculates the nuclear matrix elements in laboratory coordinates,which are necessary as Slater determinants are defined in the laboratory frame.This is conveniently done using the Talmi-Brody-Moshinsky transformation [16].

    However,one cannot use the Talmi-Brody-Moshinsky transformation with the Berggren basis.In fact,a generalization of the Talmi-Brody-Moshinsky transformation to arbitrary bases exists and consists of the vector brackets method [17—19].Unfortunately,it is of little practical use,as it leads to very time-consuming calculations.The reason for it is that the finite sums entering the Talmi-Brody-Moshinsky transformation become two-dimensional integrals with vector brackets [17—19].Added to this difficulty,Dirac delta and Heaviside functions enter vector brackets,whose non-analytical character in the complex plane even complicates the extension of vector brackets to complex-energy states.

    A new theoretical method clearly had to be devised for the numerical evaluation of nuclear matrix elements with the Berggren basis.For this,one has to rely on the fundamental property of the nuclear interaction in halo and resonance many-body states,which is its localized character in the vicinity of the considered nucleus.Indeed,it implies that an expansion of the nuclear interaction in a basis of harmonic oscillator states will converge rapidly,even though the manybody nuclear states of interest are either very extended in space (halos) or even diverge therein (resonances).

    For simplicity,we will present the harmonic oscillator expansion method for a two-body interaction.Generalization to three-body,four-body,...interactions is straightforward.Let us represent the matrix elements of the two-body nuclear interactionin the Berggren basis with a finite number of harmonic oscillator states [20]:

    whereNmaxis the number of used harmonic oscillator states and where greek (Latin) letters refer to harmonic oscillator(Berggren basis) states.As the action of the nuclear interaction is localized,one can indeed assume thatcan be reliably reproduced with a fairly small value forNmax.The two-body matrix elements appearing on the right-hand side of equation (5) only involve harmonic oscillator states,so their evaluation is straightforward using the Talmi-Brody-Moshinsky transformation.Moreover,Berggren basis states only appear in overlaps of the form 〈∣αa〉,i.e.in one-dimensional integrals of the products of Berggren basis and harmonic oscillator states.As harmonic oscillator states decrease like Gaussians when r→+∞and as Berggren basis states can at most increase exponentially in modulus therein,these overlaps are always well defined and converge rapidly when r→+∞.

    As a demonstration of the effectiveness of the harmonic oscillator expansion method,we considered the18O nucleus with an16O core[20].The18O wave functions then consist of two valence nucleons,so that GSM calculations can be quickly effected.The used nuclear interaction consists of the N3LO realistic interaction,renormalized with the Vlow-kmethod using a cut-off momentum Λ=1.9 fm-1[20].

    The convergence of the energies of the ground and excited states of18O nucleus is depicted in table 1.Convergence occurs for the rather small valuenmax~10,and evenfor resonance states,where one can also see the fast convergence of the imaginary part of energies,proportional to widths.

    Apparently,the calculation of Coulomb two-body matrix elements with the Berggren basis still remains problematic.Indeed,the Coulomb interaction is of infinite range,so that equation (5) cannot be directly applied for that matter.The solution to this problem relies on the fact that the Coulomb interaction must behave as (Z-1)/r at large distances.Thus,one can rewrite the Coulomb interactionas is a given Coulomb potential of charge Z-1.The interest of this decomposition is that(Z-1)is of shortrange character so that it can be expanded using equation(5).The one-body part(Z-1)poses no problem,as it is spherical and hence can be directly included in the generating potential of the Berggren basis.

    Table 1.Convergence of the real and imaginary parts of energies(E)of the 01 +,0 2+,21 +,22 +,41+ and 42 +eigenstates of 18O(in MeV).The first column consists of the number of nodes of the harmonic oscillator basis states used in equation (5) (adapted from [20]).

    As a consequence,the calculation of two-body matrix elements with the Berggren basis is affected by separating its finite-range and infinite-range parts and by introducing a harmonic oscillator expansion afterward (see equation (5)).It only relies on the Talmi-Brody-Moshinsky transformation and quickly converging one-dimensional overlap integrals,so that it is convenient to use in practical applications.

    3.2.Determination of eigenvalues in Gamow shell model: the overlap method.

    In order to diagonalize the standard shell model matrices,the Lanczos method is widely employed,as it allows for the determination of low-energy nuclear states without having to fully diagonalize the Hamiltonian matrix[21].However,even though the Lanczos method can be straightforwardly extended to complex symmetric matrices,it cannot be applied directly in GSM,because resonance A-body states are surrounded by many scattering A-body states.Indeed,the lowenergy spectrum obtained in GSM with the Lanczos method consists of bound,resonance and scattering eigenstates,which cannot be separated out.

    In fact,it is impossible to know whether a GSM eigenstate is a resonance or scattering state from the knowledge of its eigenenergy only.Only bound states can be identified unambiguously from their eigenenergies,which are always negative.Hence,one has to determine the resonant or scattering nature of eigenstates,even though one does not have access to their asymptotic properties.A solution to this problem could be devised with the overlap method[22],which is performed in two steps:

    · The Hamiltonian matrix is diagonalized at pole approximation level,i.e.with a basis of Slater determinants where only resonant one-body states are occupied.One then obtains the∣〉 zeroth-order eigenvector,issued from an exact diagonalization of the Hamiltonian matrix at the pole approximation level.

    · The full Hamiltonian matrix is diagonalized,i.e.with a basis of Slater determinants generated by the complete Berggren basis,with both resonant and scattering onebody states.The sought GSM eigenstate∣〉 is then the eigenstate which maximizes the∣〈〉∣overlap.

    As the occupation of non-resonant basis Slater determinants is small compared to pole configuration components,this method works very well in practice.Note that the Jacobi-Davidson method [23]is used at this level instead of the Lanczos method for Hamiltonian diagonalization.Indeed,as resonance energies are interior eigenvalues,the Jacobi-Davidson method can efficiently provide associated eigenstates,in contrast to the Lanczos method,well suited for extremal eigenvalues [21,24].

    3.3.Natural orbitals

    In order to diminish the size of Hamiltonian matrices,it is possible to use optimized one-body states.While the HF method provides by definition the best basis-generating potential [25],the obtained one-body basis states are not optimal because internucleon correlations are neglected in the HF potential.In fact,the optimal one-body states are those with which nucleon occupation is maximized in the considered many-body eigenstate,which are called natural orbitals [26].

    Natural orbitals are the eigenvectors of the scalar onebody density matrix [26],whose matrix elements read:

    where m,n are the indices of the one-body states of the used Berggren basis and∣〉is the sought many-body state.As∣〉is not known,one uses instead an approximate eigenstate,calculated in a truncated model space generated by the Berggren basis,in which numerical calculations are tractable.Note that,contrary to well bound systems,natural orbitals in unbound many-body states are unbound as well.Added to that,they contain sizable scattering components,issued from the continuum coupling present in∣〉.Consequently,the use of natural orbitals is not restricted to bound NCGSM eigenstates,but is also of practical interest for resonance NCGSM eigenstates.

    As natural orbitals are defined with respect to the sought many-body state,only a few natural orbitals are necessary to accurately represent a given partial wave.One has seen in practice that typically 5—7 natural orbitals per partial wave have to be used,which is to be compared to the 30—50 discretized scattering states necessary when using the Berggren basis [27].Due to the quickly increasing matrix dimensions with the number of basis one-body states,the size of matrices can be decreased by several orders of magnitude when replacing Berggren basis states by natural orbitals.Consequently,calculations impossible to do with the Berggren basis can be performed without truncations with natural orbitals.

    However,if the number of valence nucleons and/or valence states becomes too large,even the use of natural orbitals is not sufficient.In fact,it is necessary to resort to other diagonalization methods as the Jacobi-Davidson method[23]in this case.Hence,one has developed for that purpose the density matrix renormalization group method (DMRG),where the matrix-vector applications involving the full Hamiltonian matrix are avoided.

    3.4.Diagonalization of very large Gamow shell model matrices with the density matrix renormalization group approach

    DMRG is an iterative method that allows diagonalizing Hamiltonians without having to explicitly store their configuration interaction matrix.For this,Hamiltonians are instead represented by bases where particles are increasingly correlated through a renormalization process,so that matrix dimensions to consider remain small,while their eigenvectors converge to those of the initial Hamiltonians.In fact,DMRG has become a standard method to solve the many-body problem in matter physics[28].DMRG has been tested a few years ago in nuclear physics with the harmonic oscillator shell model,where it only found little success,however (see [29—32]).Convergence was noticed to be rather slow,as a large part of the total shell model space was needed in order for eigenenergies to converge.The reason for this situation is the large nucleon—nucleon correlations present in atomic nuclei.Indeed,DMRG converges the fastest when correlations are weak,even if numerous,as is the case in spin chains,where DMRG was introduced [33,34].

    Similarly,the nuclear many-body states studied in GSM are very good candidates for DMRG.Indeed,contrary to wellbound nuclei,correlations involving the numerous non-resonant configurations are weak,while the strong nuclear interaction is mainly acting among a few pole configurations.Consequently,nuclear states in GSM can be expected to converge more rapidly than well-bound nuclei with DMRG.However,DMRG was initially implemented for Hermitian Hamiltonians,so it was necessary to extend it to the complex symmetric case.

    For that matter,the GSM valence model space is separated into two spaces: the first one,the A reference space,is generated by the one-body resonant states of the Berggren basis,while the second one,the B space,is built in the DMRG process,by adding one after the other the remaining scattering states of the Berggren basis (see figure 2).When dimensions become too large in the non-resonant space in DMRG,the Hamiltonian matrix is diagonalized,so that it provides an approximate nuclear eigenstate by using the overlap method (see section 3.2).At this level,the reduced density matrix is calculated and diagonalized,as its eigenvectors have been shown to be the closest to exact Hamiltonian eigenvectors[33,34].A new optimal basis space is then calculated from the Nopteigenvectors of the reduced density matrix bearing the largest eigenvalues.All suboperators associated with the Hamiltonian must then be recalculated,which then forms a DMRG step.Several‘sweeps’,i.e.inclusions of the same one-body basis states by enumerating several times all scattering basis states,are necessary in order to recapture the inter-nucleon correlations neglected in the previous renormalization processes.

    Figure 2.Simplified illustration of the DMRG method.The {kA}states come from the A reference space and the αB states belong to the B space.The shell (lj)s is added to the B space,which generates the new basis states {kA ?{αB ?(lj)s}}J (from [35]).

    The first application of DMRG in GSM was done in [35],where only valence neutrons were considered.The more general case of model spaces containing both valence protons and neutrons was introduced a few years later[36].DMRG has been applied in light nuclei in the core+valence nucleon picture.Theground state and thefrist excited state of7He were considered in [35],and the ground states of the7Li and8Li nuclei were studied afterwards in [36].The fast convergence of DMRG in GSM is depicted in figure 3 for the7He nucleus.

    4.No-core Gamow shell model

    4.1.Hamiltonian and interaction

    In ab initio approaches,and hence in NCGSM,the aim is to solve the A-body Schr?dinger equation,where H is the intrinsic nuclear Hamiltonian:

    Figure 3.Convergence of the real and imaginary parts of the energy of ground state of 7He(upper and lower parts,respectively),as a function of the number of steps Nstep for several values for the number of retained eigenvectors of the reduced density matrix Nopt(from [35]).

    In NCGSM,local or non-local interactions can be included,contrary to the Green’s Function Monte Carlo approach,for example,where non-local potentials cannot be handled as Green’s Function Monte Carlo are solved in coordinate space.Thus,a local interaction such as the Argonne υ18potential[1],or a non-local interaction such as the CD-Bonn 2000[2],as well as various chiral interactions,can be considered.

    Softening techniques such as Vlow-k[37],similarity renormalization group [38,39],or G-matrix approaches[9,40],must usually be applied in order to be able to use bare nuclear interactions in small model spaces.This will be the case in the examples discussed in the next sections,where one uses the nuclear chiral interactions [41],renormalized with the Vlow-ktechnique,using a sharp momentum cutoff Λ~2 fm-1in order to decouple the degrees of freedom of low and high momenta [42].Note that,in principle,the microscopic potential and chiral interactions derived from effective field theory are different mathematical objects.However,as we do not deal with the nuclear interaction itself,but only with its application in GSM,we will not insist on this point.

    4.2.Center of mass treatment

    As in the harmonic oscillator shell model,a spurious center of mass excitations can appear in NCGSM.This arises because the number of physical degrees of freedom of the nucleus is 3A-3,where as the many-body wave function depends on A nucleons,and hence on 3A space coordinates.Thus,spurious center-of-mass motion arises due to the three redundant degrees of freedom.

    The solution to this problem in the harmonic oscillator shell model has been known for a long time and consists of the application of the Lawson method[43].For this,by using shell model spaces defined in full N?ω space,the total wave function factorizes into∣ψrel〉 ?∣ψCM〉,where ∣ψCM〉 is an eigenstate of the harmonic oscillator center of mass Hamiltonian [43].However,the Lawson method cannot be extended in general to the NCGSM.Indeed,as∣ψCM〉 is localized in space,the former decomposition of the total wave function in relative and center of mass parts only applies to well-bound many-body states.In fact,an exact separation in relative and center of mass parts of NCGSM eigenstates cannot occur except in infinite model spaces.

    Fortunately,it is not necessary for the NCGSM ground state wave function to be separated in relative and center of mass parts to reliably represent the exact many-body wave function in infinite space,which is due to the generalized variational principle [44,45].Indeed,similarly to the variational principle governing bound states,one can show that the eigenenergy of NCGSM ground states is stationary with respect to small changes of NCGSM space[44].Thus,due to the intrinsic character of the Hamiltonian,eigenenergies converge to their exact value in infinite space when NCGSM model spaces are becoming larger,even in the absence of an exact separation in relative and center of mass parts.The difference between NCGSM eigenstates in truncated and infinite space can similarly be assumed to be small if the model space used is sufficiently large.The fundamental advantage of this method is that it can be applied to many-body unbound states.This will be the method used for the next applications,dedicated to the study of the resonance ground states of5He and3,4n.

    Note,however,that the generalized variational principle does not apply to excited states,as we have no guarantee that an obtained excited state is not a repetition of the ground state with a different center of mass motion.In fact,how to correct in the general case for center-of-mass motion in NCGSM is still an open question.

    5.Ab initio description of the 5He resonance nucleus

    5.1.5He resonance properties

    Due to the well-unbound character of both its ground and first excited states,the ab initio description of5He is a difficult task.Its resonance properties had been evaluated along with those of other resonance many-body states with the variational Monte-Carlo method in [46].The associated nucleonalpha scattering reaction had also been studied with the quantum Monte-Carlo method in [47].As both halo and resonance many-body states can be calculated in NCGSM,the study of5He is well-suited in the frame of NCGSM.

    For that matter,as the width of5He mainly comes from the neutron p3/2partial wave,one will use the Berggren basis to expand this partial wave,where 0p3/2is a one-body resonance state,of energy 1.193 MeV and width 1.267 MeV.The neutron s1/2and p1/2partial waves are chosen to be real,where the p1/2partial wave consists of a real-energy contour only,while the s1/2contour must be complemented by the well-bound 0s1/2pole state,of energy equal to —24 MeV.As other partial waves add binding energy to the5He system without changing its asymptotic behavior,they are represented with harmonic oscillator basis states.

    Due to the very large dimensions occurring in the NCGSM description of5He,one has used DMRG to diagonalize the NCGSM Hamiltonian (see section 3.4).Note that the dimension of the DMRG Hamiltonian matrix to diagonalize is 105,which is much smaller than the NCGSM Hamiltonian matrix in a basis of Slater determinant,which is close to 3×109.

    Results are shown in the upper part of figure 4,for which one has used Λ=1.9 fm-1.The obtained energy for the ground state of5He is R(E)=-26.31 MeV,hence one MeV too high with respect to experimental data [48].From the NCGSM energies of the ground states of4He and5He,one can determine the5He resonance to be at neutron separation energy 1.17 MeV,which is less than 500 keV away from the experimental energy of 0.735 MeV.Presumably,the binding energy lacking in4He and5He could be attributed to the missing three-nucleon forces.The chiral interaction used might also play a role.This would explain the better reproduction of the neutron separation energy of5He compared to its total binding energy.

    Resonance properties are extracted differently according to the model used [49—51].One can either use the Lane and Thomas prescription to determine the energy and width of resonance from the conventional R-matrix approach on the real axis[52],or use the extended R-matrix approach[50,51],where the complex poles of the S-matrix provide the energy and width of the resonance from their real and imaginary parts.As NCGSM is a numerical method where many-body S-matrix poles are calculated,the definition of resonances of the extended R-matrix approach is closest to that of NCGSM.The latter extraction method approach determines the position of the 3/2-ground state of5He at an energy 0.798 MeV above the neutron threshold,while its width Γ is equal to 648 keV.The width of the ground state of5He is provided by the imaginary part of the5He ground state energy (see the lower part of figure 4).The obtained imaginary part of the energy of the ground state of5He is I(E)=-0.2 MeV,i.e.Γ=400 keV.As for total binding energy,but to a lesser extent,the difference between the widths arising from NCGSM and experimental data might be due to missing many-body forces,or to the chiral interaction itself.

    Figure 4.Convergence of the real and imaginary parts (upper and lower parts,respectively)of 5He ground state energy,calculated with the density matrix renormalization group,is shown as a function of the number of iterations,also called sweeps.Experimental data are also provided [48](from [12]).

    In order to verify these assumptions,the 3/2-ground state of5He has been recalculated in NCGSM using slightly different Hamiltonian parameters.For this,one has used two different Vlow-krenormalizing momenta,which are Λ=1.9 fm-1and 2.1 fm-1.For this comparative study,a smaller model space has been utilized,where only the neutron p3/2partial wave is represented by Berggren basis states,other one-body basis states being harmonic oscillator states.Added to that,one truncated the NCGSM space so that a maximum of four nucleons in the non-resonant continuum is accepted in the basis of Slater determinants.The NCGSM is then sufficiently small to be able to use the Jacobi-Davidson diagonalization method [11,23].

    With Λ=1.9 fm-1,the obtained energy for4He is equal to-27.386 MeV,while that of5He is E=-25.825 MeV,with a width Γ=370 keV.These results are close to the DMRG results and hence justify a posteriori the use of truncated model spaces.When the Λ momentum has the value 2.1 fm-1,one obtains instead-26.06 MeV for the binding energy of4He,where as the total binding energy and width of5He are E=-23.90 MeV and Γ=591 keV,respectively.Thus,the neutron separation of the5He ground state is 2.15 MeV above the α+n decay threshold.Hence,clearly,the resonance properties of the ground state of5He largely depend on the interaction used,which explains the discrepancies encountered with respect to experimental data.

    Table 2.The DMRG results for the energy and width of the 5He ground state obtained in NCGSM are compared to experimental data,determined from several R-matrix procedures.Shown energies are provided with respect to the α+n decay threshold (adapted from [12]).

    The position of the energy of the5He resonance ground state and its associated width,obtained from R-matrix theory,are listed in table 2.One can see that the values obtained are rather different from each other,which reflects the broad character of the5He resonance ground state.

    5.2.Radial overlap integral and asymptotic normalization coefficient

    The reproduction of the energy and width of a many-body resonance state is not sufficient to assess the quality of the many-body wave function.The asymptotic behavior of a many-body resonance state can be quantitatively studied by considering the overlap function and the associated asymptotic normalization coefficient (ANC) (see [53]for definition and details).The overlap function is a radial function providing information about nuclear structure in a many-body state[46,54].As we are considering the ab initio description of5He,we will show the overlap function in the context of the reaction4He0++n→:

    The radial overlap integral for+n→is shown in figure 5,where the NCGSM calculation is performed using the chiral N3LO interaction [41]with Λ=1.9 fm-1.One can determine the value of the ANC(denoted by C in the following) by fitting the radial overlap integral in the asymptotic region by a Hankel function (see figure 5).Indeed,the ANC is the multiplicative coefficient of the fitting Hankel function and has been found equal to C=0.197 fm-1/2.

    Figure 5.The radial overlap integral for 4He 0++n→ .The asymptotic region is fitted by a Hankel function multiplied by the ANC (from [12]).

    Table 3.NCGSM results calculations of the one-neutron separation

    energy S1nand spectroscopic amplitudeas a function of the Vlow-kcut-off parameter Λ (adapted from [12]).

    Note that results heavily depend on the cutoff parameter Λ,because S1n=-1.56 MeV and S1n=-2.15 MeV for Λ=1.9 fm-1and 2.1 fm-1,respectively(see table 3).Indeed,associated ANCs are respectively equal to 0.197 fm-1/2and 0.255 fm-1/2in the latter cases.

    A convenient test of the accuracy of the ANC is the relation existing between the ANC and the resonance width Γ(see [46,53,55]and references cited therein):

    where μ is the effective mass and k is the linear momentum(k=)of the one-neutron separation energy S1n.Note that,even though k is a complex number in general,its imaginary part can be neglected if Γ ?|S1n|.

    Using equation (9),one obtains ΓANC=311 keV and ΓANC=570 keV for Λ=1.9 fm-1and 2.1 fm-1,respectively.These values for the widths are consistent with those obtained from the imaginary part of energies in NCGSM (see section 5),as the 10%—15% difference between NCGSM widths and ΓANCcan be traced back to the rather broad character of the ground state width of5He,implying that equation (9) is only approximate therein.

    The spectroscopic factor is the norm of the square of the radial overlap integral I?j(r)2:

    While S?jis not an observable [56],its value nevertheless provides important information about nuclear structure (see[54,57—61]for spectroscopic factor calculations in various situations).Indeed,a large value for the spectroscopic factor indicates that the nuclear many-body function has a singleparticle character,where as a small value,in this case,would depict a highly correlated structure.

    An interesting tendency of spectroscopic factors as a function of the momentum cut-off Λ has been found inHe pairs calculated in NCGSM(see table 3).Indeed,by decreasing the momentum cut-off Λ,or,equivalently,by softening the interaction,one can see thatdeparts more and more from unity.This shows that one cannot infer a priori any single-particle character of the ground state of5He,i.e.of one neutron occupying a 0p3/2resonance state in the mean field potential generated by the4He core.On the contrary,the5He wave function content heavily depends on the Λ value used in the Hamiltonian.

    6.A=3,4 unbound nuclei

    6.1.Trineutron and tetraneutron in the auxiliary potential method

    The trineutron and tetraneutron systems are not well known despite their small number of nucleons (see [62]and references therein).The study of multi-neutron systems,which are broad resonances,is difficult with standard ab initio methods,as the latter are tailored for well-bound atomic nuclei[63,64].

    Hence,NCGSM is the tool of choice to study multineutron systems,as continuum coupling and inter-nucleon correlations can be both included in this framework.This has been done firstly in [65]for tetraneutron,where ,unfortunately,calculations were still incomplete,as they could be performed only in truncated model spaces or with overbinding interactions so that one could not make any definite conclusion.Consequently,multi-neutron systems have been reconsidered in [62]in the context of NCGSM,where converged results could be obtained in physical situations,albeit in truncated spaces.

    The main issue to solve in the first place is the identification of the resonance states of multi-neutron systems.Indeed,multi-neutron systems are so dilute that the overlap method presented in section 3.2 is no longer sufficient to extract resonances from the many-body spectrum made almost exclusively of scattering states of complex energy.As a consequence,we introduced auxiliary potentials to confine neutrons in an external trap.This method was initially introduced in the quantum Monte-Carlo framework [64,66].For this,one adds to the ab initio Hamiltonian of equation(7)an auxiliary Woods—Saxon potential of potential depth Vaux,which will vary from a finite negative value to zero.As multineutron systems become well bound for sufficiently large negative values of Vaux,the calculation of Hamiltonian eigenstates is stable therein when applying the overlap method (see section 3.2).From this starting point,one can follow Hamiltonian eigenstates adiabatically when Vaux→0.It has been checked that this method provides nearly exact results for the dineutron system [67].

    The ground states of the trineutron and tetraneutron will be calculated using the auxiliary potential method.The interaction used is the chiral two-body interaction N3LO[68]renormalized with Λ=2.1 fm-1.The single-particle Berggren basis,possessing bound,resonance and scattering states,is generated by a standard Woods—Saxon potential [68].One imposes either 3p-3h and 4p-4h truncations in NCGSM model spaces,i.e.only three or four neutrons can occupy nonresonant basis states in Slater determinants,respectively.

    One also uses the so-called ‘large’ model space,where basis Slater determinants are generated from s1/2,p3/2Berggren basis states,as well as p1/2,f5/2,f7/2,g9/2harmonic oscillator states.Alternatively,in the so-called ‘small’ model space,Slater determinants are built using the same partial waves as in the ‘large’ model space,but with the p1/2,d3/2,d5/2and f5/2,f7/2,g9/2partial waves represented by only two and one harmonic shells,respectively.Eigenenergies have been checked to vary smoothly vary when Vaux→0,i.e.that no multi-neutron channel opens when Vauxdecreases.Results for both the trineutron and tetraneutron in various model spaces are shown in figure 6.

    However,despite the progress made compared to [65],tetraneutron energies could not be reached in the ‘large’model space beyond the value of Vaux=-1.5 MeV.In fact,the physical region could be attained either in NCGSM truncated spaces,in ‘small’ model spaces,denoted as 'small’in figure 6,or using 3p-3h truncations.Therefore,it is still necessary to rely on extrapolations to have an estimate of the energy and width of tetraneutron with realistic interactions.

    The auxiliary potential method,used with the Jacobi-Davidson diagonalization technique (see figure 6) [23],and DMRG (see section 3.4) [62],both predict similar tetraneutron energy,which is close to 2.6 MeV.This energy lies within the range provided by experimental errors [62].The width is estimated to be equal to 2.38 MeV from the extrapolation of the results obtained in the ‘large’ model space with Vaux<-1.5 MeV.Conversely,it was possible to obtain the trineutron ground state without model space truncation in the physical region.One obtained therein E3n=1.29 and Γ3n=0.91 MeV.One should recall,nevertheless,that these results are tentative,and hence not definitive and still subject to debate.For example,the results obtained with Faddeev—Yakubovsky techniques in [69]indicate that no n4resonance can be produced from realistic nuclear interactions unless they are unphysically modified.

    Interestingly,the energy and width of the trineutron are smaller than that of the tetraneutron (see figure 6),its energy still being larger than that of the dineutron.Consequently,it might be that the trineutron resonance might not be as broad as tetraneutron,which should be confirmed experimentally if it is the case.

    Figure 6.Energies and widths of the trineutron and tetraneutron systems as a function of the depth of the auxiliary Woods—Saxon potential.For tetraneutron calculations,NCGSM spaces with three particles at most in the non-resonant continuum,and without model space truncation,have been considered.Berggren or natural orbital(n.o.)bases have been used,while‘small’or‘large’denotes a model space where the number of basis shells has been reduced or not(see text for details).Tetraneutron width was decreased by a factor 3 for readability.Trineutron results are always obtained without model space truncation.The experimental data for tetraneutron (see [62])are also shown (adapted from [68]).

    6.2.A=4 unbound isospin multiplet

    Nuclei bearing four nucleons possess resonance states in their spectrum.Hence,they form good laboratories to test realistic nuclear interactions in nuclei far from stability with NCGSM.Of particular interest is the A=4 T=1 isospin multiplet formed by the first 2-states of4H,4He and4Li.This isospin multiplet has unique features,because continuum coupling is maximal in unbound states,so that one can assess the influence of the Coulomb force in broad resonances and then quantitatively measure isospin symmetry breaking in this case.

    The isospin symmetry in the considered A=4 T=1 multiplet states is indeed expected to be well broken.This is especially due to the Coulomb force[70—72],which acts very differently in4H,4He and4Li.Indeed,both proton—proton and Coulomb forces are absent in4H,where as4Li has no neutron-neutron interaction.Consequently,we studied the unbound spectrum formed by the first 2-states of the spectrum of4H,4He and4Li with NCGSM in [45].

    Figure 7.Unbound 2-states of T=1 isospin multiplet of the A=4 nuclei,calculated using NCGSM and the N3LO-SRG2.0 interaction.Experimental data come from [71,90].Dotted lines connect the T=1 isobaric triplet states of A=4 nuclei to guide the eye,while the obtained energies for A=3 are also given for comparison.The experimental energy of the T=1 Jπ=2- state in 4He is deduced from the R-matrix analysis in [71](from [45]).

    Results are shown in figure 7,where one presents manybody states in A=3,4 nuclei,along with the isospin triplet states calculated with NCGSM using the realistic interaction N3LO-SRG2.0,where SRG stands for similarity renormalization group.Both the ground states of4H and4Li are unbound and bear large particle-emission widths.The eigenstates of A=3,4 nuclei can be seen to be well described in figure 7.The broad character of the T=1 isospin triplet states of A=4 nuclei is reproduced[70].In particular,the calculated energy of the T=1 Jπ=2-excited state in4He is close to that arising from an R-matrix analysis [71].Indeed,the4He 2-excited state is found to be 22.03 MeV above the4He ground state.This result agrees with the experimental value of[70],which is equal to 22.5±0.3 MeV.Nevertheless,the width of the4He 2-excited state obtained in the R-matrix analysis is about three times larger than in our calculation.

    The calculations illustrated in figure 7 allow us to exhibit the isospin symmetry breaking occurring in the T=1 isospin triplet states of A=4 nuclei.One can see that the 2-state of4H,which is also its ground state,is the lowest state in energy.The4H isotope is more bound by about 0.5 and 1.7 MeV when it is compared to the 2-states of the same multiplet in4He and4Li,respectively.As could be expected,the width of the 2-excited state in4He,equal to about 1.5 MeV,has a value in between those of4H and4Li.

    In order to have a well-defined estimate of isospin symmetry breaking in the T=1 triplet states,the isospin expectation value of their 2-analog states has been calculated.While the value of isospin in the ground states of4Li and4H ground states is equal to 1,that for the 2-excited state of4He is equal to about 0.71.This suggests that isospin symmetry breaking is large in this state.Note that the isospin mixing occurring in the 2-excited state of4He is consistent with experimental data [73,74].

    The presented NCGSM calculations of the A=4 T=1 isospin multiplet have clearly exhibited isospin symmetry breaking.While it originates from the Coulomb interaction and continuum coupling,it is difficult to disentangle both these effects,however.

    7.Renormalization of realistic interactions in a finite model space using -box method

    In order to consider heavier systems than in NCGSM,while still dealing with realistic Hamiltonians,it would be useful to extend IMSRG[4—8]or MBPT[9,10]to GSM,as one would reduce model space dimensions.Therefore,the extension of the-box method has been implemented within the GSM framework,similarly to the harmonic oscillator shell model [75].

    However,the generalization of the IMSRG and MBPT methods to GSM is not straightforward.Indeed,a caveat appears with the optimization of the single-particle basis,which does not exist in the harmonic oscillator shell model.It had been noticed that the many-body perturbation expansion entering the-box method poorly converges with a Berggren basis arising from a Woods—Saxon potential [75].The solution to this problem relies on the use of the HF potential of the Hamiltonian instead of a Woods—Saxon potential to generate the one-body basis.As many diagrams present in MBPT cancel out when a HF basis potential is utilized,the convergence of the-box method is much faster,so that IMSRG and MBPT can be efficiently applied in GSM.

    MBPT is based on the extended Kuo-Krenciglowa(EKK)method,which provides a quickly converging series of Hamiltonian powers [76].However,it is derived with degenerate basis states,so it can hardly be used with the Berggren basis,as its one-body basis states all have different energies.Consequently,it is more convenient to use the EKK method generalized to nondegenerate basis states in order to obtain the effective Hamiltonian[10,77,78].

    Figure 8.Ground-state and one-neutron separation energies of oxygen isotopes (upper panel) and fluorine isotopes (lower panel)compared with experimental data [90](The extrapolated values issued from AME2016 are used for 27,28O and 30,31F).For other calculations,one also depicts results obtained with the coupledcluster approach with continuum (CC) [85],GSM (GSM) [91],valence-space IMSRG (VM-IMSRG) [84],shell model with twonucleon forces and three-nucleon forces (SM NN+3NF) [80],shell model with the USDB interaction(SM USDB)[83]and shell model with the SDPF-M interaction (SM SDPF-M) [92](from [13]).

    The GSM ground-state and single-neutron separation energies,denoted as Egsand Sn,respectively,are illustrated in figure 8 and compared to experimental data and theoretical calculations issued from other models.One can see that the neutron dripline is reproduced for oxygen isotopes in shell model calculations when three-nucleon forces are included[80](labeled ‘SM NN+3NF’ in figure 8).However,the quality of results strongly depends on the model used beyond the neutron dripline.This happens to be the case because the coupling to the decay channels must be included for the neutron dripline to be correctly predicted [81,82].Open decay channels and particle continuum are indeed not included in the shell model with two-nucleon and three-nucleon forces [80],shell model USDB [83]and valence-space IMSRG calculations [84](see respectively curves labeled‘SM NN+3NF’,‘SM USDB’and‘VM-IMSRG’in figure 8).Conversely,GSM utilized with MBPT and two-body interactions(labeled‘This work(continuum)’in figure 8)correctly provides the location of the neutron dripline for oxygen isotopes.As a consequence,one can obtain a good description of the25,26O isotopes with GSM,even though they are situated beyond the neutron dripline(see the upper panel in figure 8).However,the effects arising from three-nucleon forces and couplings to the scattering continuum cannot be disentangled.

    In figure 8,we also presented results obtained with the coupled-cluster model using the Berggren basis and a schematic three-nucleon force (labeled ‘CC (continuum)’ in figure 8)[85].Even though the neutron dripline is reproduced in this approach,results depart from experimental data as a whole.The last theoretical calculation shown in figure 8 is that of GSM (see curves labeled ‘GSM (continuum)’ in figure 8),where a two-body interaction fitted to the bound and resonance states of23-26O is utilized with a22O core.Energies are well described,with the28O ground state predicted as very weakly bound or unbound.The28O ground state is also more bound than27O and slightly more unbound than26O.

    All25-28O isotopes are unbound,as can be seen from their one-neutron separation energies (see figure 8).In fact,a possible tetraneutron decay in28O might occur if S4nis negative and Sn,S2nand S3nare all positive [86].However,that decay channel is excluded in GSM using-box interactions as S2n<0 therein.Nevertheless,a sequential emission of two dineutrons cannot be ruled out.

    The ground state wave function of28O is unbound and dominated by theconfiguration.Hence,the four neutrons occupying the single-particle resonance ν0d3/2play a major role in the width value of the28O ground state,which is about several tens of keV.In order to assess coupling to the continuum at the end of the oxygen chain,we did calculations in the presence and absence of d3/2scattering contour (see figure 8).The influence of non-resonant basis states is not negligible as the gain in binding energy is close to 0.5 MeV for26O and 1 MeV for28O.

    Let us now consider the results obtained for fluorine isotopes using GSM with the NNLOoptinteraction (see the lower part of figure 8).As for the oxygen chain,we use a22O core and energies are given with respect to the22O ground state.A good description of ground state energies is obtained for the23-28F isotopes.Nevertheless,GSM provides the energy of31F to be lower than that of30F,even though31F is located at the neutron dripline experimentally [87,88].The imposed model space truncations,which are necessary as Hamiltonian matrix dimensions quickly increase with the number of valence nucleons,might be at the origin of this discrepancy.To give an idea of dimensions in untruncated spaces,the29F Hamiltonian matrix has a dimension close to 2×1010using an sdpf7/2valence space.As one cannot diagonalize matrices that large in GSM at present,2p-2h truncations had to be considered in the non-resonant continuum,i.e.only two neutrons can occupy scattering basis states in Slater determinants.

    The analysis of effective single-particle energies can partially explain the different positions of driplines obtained for oxygen and fluorine isotopes.While the neutron 0d3/2effective single-particle state is found to be unbound in the oxygen isotopes,it is,on the contrary,bound in fluorine isotopes,as close to 1 MeV of binding energy is obtained from the proton-neutron interaction.

    The previous treatment of ab initio Gamow shell model framework has then been seen to be very powerful,as one could consider the oxygen and fluorine chains with realistic interactions,without having to solve the full NCGSM problem.In fact,it has become clear that the use of IMSRG and MBPT in GSM must be generalized to other systems as very light nuclei,as NCGSM can only be applied to the latter.The main problem to tackle will be to generate effective interactions for spectra involving several unbound many-body states.An interesting example for that matter is the mirror nuclei16N and16F,as the spectrum of16F is unbound and has sizeable widths,where as16N is well bound.The first calculations done for these nuclei are encouraging [89].

    8.Conclusion

    The description of atomic nuclei from ab initio frameworks is one of the most advanced nuclear methods,as it takes,in principle,all nuclear degrees of freedom into account.However,its numerical application is difficult in no-core approaches,firstly from a theoretical point of view,because of the presence of a hard core in the bare nuclear interaction,and from a practical point of view,because of the quickly increasing model space dimensions with the number of nucleons.While NCSM cannot be applied for nuclei bearing more than about 15 nucleons,NCGSM can hardly be utilized for nuclei with A≥6,7 due to the inclusion of continuum coupling.

    NCGSM has been firstly presented in the context of5He.This nucleus is well suited for an NCGSM study,as it is unbound at ground state level and bears only a few nucleons.It could be successfully described within NCGSM,where it was also shown that its supposed single-particle character depends,in fact,on the Hamiltonian used and is not a general feature of its wave functions with various realistic interactions.

    Then,we presented calculations related to A=3,4 unbound nuclear states.These nuclei are interesting because they are not well known theoretically and experimentally and are broad resonances.Due to their small number of nucleons,they can also be considered in an NCGSM approach.Obtained results are consistent with available experimental data.However,even though tetraneutron can be assumed to have a width slightly larger than 2 MeV,results are still nonconclusive,because this estimate still relies on extrapolations.We could also see that trineutron has a surprisingly small width in NCGSM calculations,of about 1 MeV.It would then be interesting to confirm this value experimentally.The large isospin symmetry breaking in the isospin multiplet of T=1 A=4 unbound nuclear states could also be unambiguously exhibited from NCGSM calculations,which might be of interest for future experiments in this region.

    The oxygen and fluorine isotopic chains at neutron dripline have also been calculated in the frame of GSM using realistic interactions within thebox method based on MBPT.These chains are interesting because they contain part of the heaviest nuclei with which the neutron dripline can be attained experimentally.A full NCGSM calculation of these isotopes is out of reach so a renormalization scheme such as the-box method is necessary to study such nuclei with ab initio formalisms.Results reproduce experimental data satisfactorily,where one could show that both a realistic Hamiltonian and the presence of continuum coupling are necessary for that matter.

    Thus,GSM has been shown to be a reliable tool to study dripline nuclei with realistic Hamiltonians.While nuclei bearing only 3—5 nucleons can be treated almost exactly with NCGSM,nuclei with more than a few nucleons have to be considered with transformed Hamiltonians via IMSRG or MBPT,for example,with which very precise eigenenergies can be obtained.The halo and resonance nuclei at driplines can then be studied with realistic interactions within the same level of accuracy as well bound nuclei.

    It is then considered to study a larger number of nuclei with GSM extended to ab initio frameworks.For example,the two-nucleon halo nucleus6He is currently being studied in this context but is at the limit of the computational capabilities of NCGSM.Conversely,possibilities are more numerous when using ab initio GSM with a core.An interesting study concerns the depiction of the charge-symmetry breaking features of the nuclear interaction.Indeed,as their effects on nuclear energy spectra depend on the different asymptotes present in mirror pairs of nuclei,GSM is very well suited for this type of quantitative study.Other calculations in the same context involve medium drip-line nuclei,as both proton and neutron drip-line can be attained up to A~40.As deformation becomes prominent in this zone of the nuclear chart,ab initio GSM can provide useful information about the interplay of continuum and deformation degrees of freedom when coupled by a realistic interaction.

    Acknowledgments

    This work has been supported by the National Natural Science Foundation of China under Grants No.12175281,11835001,11921006,12035001,and 11975282;the State Key Laboratory of Nuclear Physics and Technology,Peking University under Grant No.NPT2020ZZ01 and NPT2020KFY13;by the National Key R&D Program of China under Grant No.2018YFA0404401;the Strategic Priority Research Program of Chinese Academy of Sciences under Grant No.XDB34000000;the Key Research Program of the Chinese Academy of Sciences under Grant No.XDPB15;and the CUSTIPEN (China-U.S.Theory Institute for Physics with Exotic Nuclei) funded by the U.S.Department of Energy,Office of Science under Grant No.desc0009971.The High-Performance Computing Platform of Peking University and Gansu Advanced Computing Center are acknowledged.

    ORCID iDs

    国产伦一二天堂av在线观看| 免费人成在线观看视频色| 小蜜桃在线观看免费完整版高清| 日本黄大片高清| 91av网一区二区| 国产一区二区亚洲精品在线观看| 午夜福利视频1000在线观看| av黄色大香蕉| 一个人看的www免费观看视频| 中文欧美无线码| 国产av不卡久久| 美女内射精品一级片tv| 亚洲不卡免费看| 国产精品乱码一区二三区的特点| eeuss影院久久| 久久久久久国产a免费观看| 午夜福利视频1000在线观看| 国产精品麻豆人妻色哟哟久久 | 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 国产精品,欧美在线| 99久久无色码亚洲精品果冻| 日本在线视频免费播放| 国产老妇女一区| 日韩精品有码人妻一区| 好男人视频免费观看在线| 麻豆精品久久久久久蜜桃| 久久这里有精品视频免费| 午夜福利在线观看免费完整高清在 | 国产片特级美女逼逼视频| 国产综合懂色| 国内精品美女久久久久久| 亚洲av.av天堂| 欧美成人a在线观看| 一卡2卡三卡四卡精品乱码亚洲| 伦精品一区二区三区| 男人舔奶头视频| 变态另类丝袜制服| 99热6这里只有精品| 99久久精品一区二区三区| 1000部很黄的大片| 成人特级av手机在线观看| 亚洲欧美中文字幕日韩二区| 午夜视频国产福利| 99热这里只有是精品在线观看| 天堂av国产一区二区熟女人妻| 婷婷精品国产亚洲av| 免费av观看视频| 亚洲国产欧美在线一区| 五月伊人婷婷丁香| 亚洲美女视频黄频| 久久精品国产自在天天线| 久久99蜜桃精品久久| 三级男女做爰猛烈吃奶摸视频| 国内精品一区二区在线观看| 亚洲国产日韩欧美精品在线观看| 久久久欧美国产精品| 成年免费大片在线观看| 久久久久网色| 亚洲欧美清纯卡通| 国产成人影院久久av| av在线播放精品| 男女视频在线观看网站免费| 最近2019中文字幕mv第一页| 国产午夜精品论理片| 日产精品乱码卡一卡2卡三| 久久久精品欧美日韩精品| 日本撒尿小便嘘嘘汇集6| kizo精华| 成人三级黄色视频| 好男人视频免费观看在线| 国产精品永久免费网站| 夜夜看夜夜爽夜夜摸| 国产精品免费一区二区三区在线| 久久久久网色| 亚洲,欧美,日韩| av在线观看视频网站免费| 精品久久久久久久久av| 精品无人区乱码1区二区| a级一级毛片免费在线观看| av.在线天堂| 丰满人妻一区二区三区视频av| 高清毛片免费观看视频网站| 久久久精品大字幕| 男女啪啪激烈高潮av片| 日本免费a在线| 久久草成人影院| 亚洲最大成人中文| av免费在线看不卡| 简卡轻食公司| 久久这里有精品视频免费| 国产高潮美女av| 国产精品三级大全| 在线免费观看不下载黄p国产| 精品久久久久久久久久久久久| 久久草成人影院| 亚洲精品久久久久久婷婷小说 | 1024手机看黄色片| 国产一区亚洲一区在线观看| 最新中文字幕久久久久| 欧美区成人在线视频| 免费av不卡在线播放| 亚洲成人av在线免费| 女同久久另类99精品国产91| 国产亚洲精品久久久久久毛片| 亚洲欧美精品自产自拍| 2021天堂中文幕一二区在线观| 精品国产三级普通话版| 中国美女看黄片| 国产 一区 欧美 日韩| 美女脱内裤让男人舔精品视频 | 啦啦啦啦在线视频资源| 特级一级黄色大片| 热99在线观看视频| 国内精品宾馆在线| 国产精品野战在线观看| 亚洲成人中文字幕在线播放| 国产成人影院久久av| 大又大粗又爽又黄少妇毛片口| 国内少妇人妻偷人精品xxx网站| 丰满乱子伦码专区| 成年av动漫网址| 成人亚洲欧美一区二区av| 国产成人精品婷婷| 性欧美人与动物交配| 黄色配什么色好看| АⅤ资源中文在线天堂| 最近2019中文字幕mv第一页| 99久久精品一区二区三区| 国产伦精品一区二区三区四那| 菩萨蛮人人尽说江南好唐韦庄 | 中文字幕免费在线视频6| 国产在线精品亚洲第一网站| 成人一区二区视频在线观看| 中文字幕av在线有码专区| 日本在线视频免费播放| 又爽又黄无遮挡网站| 一本一本综合久久| 日韩av在线大香蕉| 免费搜索国产男女视频| 一边摸一边抽搐一进一小说| 一本久久精品| 国产片特级美女逼逼视频| 午夜爱爱视频在线播放| 亚洲国产精品成人久久小说 | 亚洲av二区三区四区| 国产高清不卡午夜福利| 天堂√8在线中文| 国产麻豆成人av免费视频| 波多野结衣巨乳人妻| 99久久九九国产精品国产免费| 在线国产一区二区在线| 不卡视频在线观看欧美| 久久久久国产网址| 亚洲无线在线观看| 丰满的人妻完整版| 美女xxoo啪啪120秒动态图| 一边亲一边摸免费视频| 91在线精品国自产拍蜜月| 中国美女看黄片| 18禁在线无遮挡免费观看视频| 欧美成人a在线观看| 可以在线观看的亚洲视频| 1024手机看黄色片| 1024手机看黄色片| 国产 一区 欧美 日韩| 亚洲在久久综合| 99国产极品粉嫩在线观看| 国产高清视频在线观看网站| 最近中文字幕高清免费大全6| 久99久视频精品免费| 国产成人a区在线观看| 三级经典国产精品| 久久久午夜欧美精品| 日韩强制内射视频| 亚洲欧美清纯卡通| 日韩欧美国产在线观看| kizo精华| 日本成人三级电影网站| 精品一区二区三区人妻视频| 黄片wwwwww| 久久人人爽人人爽人人片va| 美女xxoo啪啪120秒动态图| 校园人妻丝袜中文字幕| 偷拍熟女少妇极品色| 亚洲最大成人手机在线| 精品人妻熟女av久视频| 日本撒尿小便嘘嘘汇集6| 老司机福利观看| 国产成人精品一,二区 | 日日摸夜夜添夜夜爱| 亚洲欧美日韩高清在线视频| 最好的美女福利视频网| 日韩亚洲欧美综合| 不卡视频在线观看欧美| a级毛片免费高清观看在线播放| 欧美成人精品欧美一级黄| 永久网站在线| 在线观看美女被高潮喷水网站| 观看免费一级毛片| avwww免费| 欧美色视频一区免费| 少妇熟女aⅴ在线视频| 亚洲av成人av| 国产中年淑女户外野战色| 日日摸夜夜添夜夜添av毛片| 国产精品免费一区二区三区在线| 欧美精品国产亚洲| av在线天堂中文字幕| 欧美xxxx黑人xx丫x性爽| 在线免费十八禁| 精品久久久久久久久久免费视频| 亚洲精品乱码久久久久久按摩| 卡戴珊不雅视频在线播放| 国产欧美日韩精品一区二区| 自拍偷自拍亚洲精品老妇| 一个人免费在线观看电影| 国产真实伦视频高清在线观看| 日日啪夜夜撸| 国产精品免费一区二区三区在线| 国产精品久久久久久av不卡| 美女高潮的动态| 联通29元200g的流量卡| 亚洲人成网站在线播放欧美日韩| 国产精品,欧美在线| videossex国产| 久久精品国产99精品国产亚洲性色| 国产精品麻豆人妻色哟哟久久 | 欧美日韩精品成人综合77777| 国产av不卡久久| 成人无遮挡网站| 婷婷色av中文字幕| 91狼人影院| 又粗又爽又猛毛片免费看| 爱豆传媒免费全集在线观看| 最近手机中文字幕大全| 免费人成在线观看视频色| 国产免费一级a男人的天堂| 99在线人妻在线中文字幕| 久久综合国产亚洲精品| 成人性生交大片免费视频hd| 久久精品国产亚洲av天美| 国产精品人妻久久久久久| 色综合色国产| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| 人妻系列 视频| 天堂中文最新版在线下载 | 久久久精品大字幕| 一区福利在线观看| 欧美成人精品欧美一级黄| 男的添女的下面高潮视频| 国产精品久久视频播放| 亚洲,欧美,日韩| 变态另类丝袜制服| 国产成人精品一,二区 | 中文资源天堂在线| 麻豆国产97在线/欧美| 内地一区二区视频在线| av天堂在线播放| 亚洲四区av| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 国产综合懂色| 久久久久久伊人网av| 久久99热6这里只有精品| 精品一区二区三区人妻视频| 午夜福利视频1000在线观看| 在线观看免费视频日本深夜| 少妇猛男粗大的猛烈进出视频 | 18禁在线无遮挡免费观看视频| 亚洲性久久影院| 亚洲内射少妇av| 直男gayav资源| 人人妻人人澡欧美一区二区| 一进一出抽搐动态| 我的老师免费观看完整版| av卡一久久| 夫妻性生交免费视频一级片| 丰满乱子伦码专区| 天堂av国产一区二区熟女人妻| 国产私拍福利视频在线观看| 久久草成人影院| 欧美又色又爽又黄视频| 最近的中文字幕免费完整| 亚洲人成网站在线播放欧美日韩| 男女边吃奶边做爰视频| 国产亚洲欧美98| 欧美激情国产日韩精品一区| 日韩欧美 国产精品| 亚洲欧洲日产国产| 亚洲精品456在线播放app| 亚洲人成网站在线观看播放| 欧美+亚洲+日韩+国产| 少妇的逼好多水| 高清毛片免费观看视频网站| 国产精品永久免费网站| 色综合色国产| 亚洲精品色激情综合| 亚洲最大成人手机在线| 亚洲第一电影网av| 日本黄色视频三级网站网址| 99久久精品热视频| 亚洲精品国产av成人精品| 99热这里只有是精品50| www.av在线官网国产| 久久亚洲精品不卡| 欧美成人一区二区免费高清观看| 蜜桃久久精品国产亚洲av| 国产精品.久久久| 国产精品国产三级国产av玫瑰| 久久久久性生活片| 网址你懂的国产日韩在线| 国产精品三级大全| 一级av片app| 亚洲av二区三区四区| 国产一区二区激情短视频| 在现免费观看毛片| 人妻制服诱惑在线中文字幕| 精品一区二区免费观看| 99久国产av精品国产电影| 日韩一区二区视频免费看| avwww免费| 美女高潮的动态| 99国产极品粉嫩在线观看| 可以在线观看的亚洲视频| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| av福利片在线观看| 日本av手机在线免费观看| 中文在线观看免费www的网站| .国产精品久久| 亚洲av成人av| 欧美成人免费av一区二区三区| 亚洲国产高清在线一区二区三| 国产精品一区二区三区四区久久| 免费看av在线观看网站| 亚洲国产高清在线一区二区三| 国产黄片视频在线免费观看| 久久人妻av系列| 日韩av不卡免费在线播放| 久久久久久大精品| 丝袜喷水一区| 国产一区二区三区av在线 | 中文资源天堂在线| 亚洲真实伦在线观看| 国产精品久久久久久av不卡| 欧美高清成人免费视频www| 深夜a级毛片| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久噜噜| 永久网站在线| 天美传媒精品一区二区| 直男gayav资源| 久久久午夜欧美精品| 国产黄片视频在线免费观看| 级片在线观看| 夜夜看夜夜爽夜夜摸| 色综合站精品国产| 成年免费大片在线观看| 在线a可以看的网站| 美女黄网站色视频| 欧美色视频一区免费| 成年版毛片免费区| 麻豆国产97在线/欧美| 成年免费大片在线观看| 欧美最新免费一区二区三区| 国产高清视频在线观看网站| 亚洲成a人片在线一区二区| 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 99热这里只有精品一区| 欧美日韩一区二区视频在线观看视频在线 | 成人一区二区视频在线观看| 午夜福利高清视频| 偷拍熟女少妇极品色| 人妻久久中文字幕网| 你懂的网址亚洲精品在线观看 | 乱人视频在线观看| 亚洲天堂国产精品一区在线| 精品人妻偷拍中文字幕| 校园人妻丝袜中文字幕| 日韩一区二区视频免费看| 蜜臀久久99精品久久宅男| 欧美xxxx黑人xx丫x性爽| 国产成人精品一,二区 | 桃色一区二区三区在线观看| 国产熟女欧美一区二区| 久久久精品94久久精品| 亚洲久久久久久中文字幕| 中文欧美无线码| 久久精品国产99精品国产亚洲性色| 精品人妻视频免费看| 欧美成人免费av一区二区三区| 在线观看一区二区三区| 国产乱人偷精品视频| 国产精品免费一区二区三区在线| 欧美最黄视频在线播放免费| 日韩强制内射视频| 日本熟妇午夜| 男女视频在线观看网站免费| 欧美三级亚洲精品| 国产黄色视频一区二区在线观看 | 亚洲欧洲日产国产| 亚洲无线观看免费| 免费人成视频x8x8入口观看| 国产精品一区二区在线观看99 | 亚洲自偷自拍三级| 波多野结衣巨乳人妻| 中国国产av一级| 在线播放无遮挡| 国产一区二区在线av高清观看| 亚洲国产精品成人综合色| 麻豆国产av国片精品| 日韩,欧美,国产一区二区三区 | 日本色播在线视频| 亚洲av熟女| 国产精品电影一区二区三区| 亚洲精品日韩av片在线观看| 大型黄色视频在线免费观看| 精品免费久久久久久久清纯| 18+在线观看网站| 有码 亚洲区| 国产精品麻豆人妻色哟哟久久 | 亚洲丝袜综合中文字幕| 丰满的人妻完整版| 校园春色视频在线观看| 亚洲真实伦在线观看| 三级毛片av免费| 国产熟女欧美一区二区| videossex国产| 黄色日韩在线| 丝袜美腿在线中文| 伦理电影大哥的女人| 成人漫画全彩无遮挡| 日韩一区二区视频免费看| 亚洲欧洲国产日韩| 亚洲人成网站高清观看| 亚洲最大成人av| 免费看光身美女| 亚洲婷婷狠狠爱综合网| 亚洲成人久久性| 日本撒尿小便嘘嘘汇集6| av在线观看视频网站免费| 最近2019中文字幕mv第一页| 亚洲人成网站高清观看| 看免费成人av毛片| 黑人高潮一二区| 久久精品夜色国产| 精品不卡国产一区二区三区| 69av精品久久久久久| 精品久久久噜噜| 国产一区二区在线观看日韩| 狂野欧美白嫩少妇大欣赏| .国产精品久久| 乱人视频在线观看| 天堂影院成人在线观看| 男人舔奶头视频| 日韩亚洲欧美综合| 99久久精品热视频| 麻豆成人午夜福利视频| 日本黄色视频三级网站网址| av卡一久久| www.av在线官网国产| 亚洲内射少妇av| 亚洲国产精品久久男人天堂| 国产精品综合久久久久久久免费| 给我免费播放毛片高清在线观看| 一个人看的www免费观看视频| av视频在线观看入口| 欧美人与善性xxx| 激情 狠狠 欧美| 国产伦精品一区二区三区四那| 九色成人免费人妻av| 精品久久久久久久久久免费视频| 欧美bdsm另类| 狂野欧美白嫩少妇大欣赏| 三级男女做爰猛烈吃奶摸视频| 亚洲av熟女| 国产探花极品一区二区| 亚洲无线在线观看| a级一级毛片免费在线观看| 中文精品一卡2卡3卡4更新| 中文在线观看免费www的网站| 欧美bdsm另类| 欧美3d第一页| 亚洲成人中文字幕在线播放| 久久久国产成人精品二区| 国产av在哪里看| 精品一区二区三区视频在线| 国产伦在线观看视频一区| 国产精品一区二区三区四区久久| 欧美色视频一区免费| 亚洲成av人片在线播放无| 国产男人的电影天堂91| 日韩一区二区三区影片| 自拍偷自拍亚洲精品老妇| 亚洲国产精品国产精品| 综合色av麻豆| 夫妻性生交免费视频一级片| 美女 人体艺术 gogo| 91av网一区二区| 国产在线精品亚洲第一网站| 一级毛片久久久久久久久女| 欧美日韩一区二区视频在线观看视频在线 | 国产v大片淫在线免费观看| 久久精品国产清高在天天线| 欧美高清成人免费视频www| 天天躁日日操中文字幕| 男女啪啪激烈高潮av片| 亚洲不卡免费看| 亚洲国产欧美人成| h日本视频在线播放| 能在线免费观看的黄片| 久久精品久久久久久噜噜老黄 | 国内少妇人妻偷人精品xxx网站| 特大巨黑吊av在线直播| av又黄又爽大尺度在线免费看 | 女同久久另类99精品国产91| 秋霞在线观看毛片| 日日撸夜夜添| 中国美白少妇内射xxxbb| 亚洲精品影视一区二区三区av| 精品人妻偷拍中文字幕| 久久精品国产亚洲av天美| 国产精品日韩av在线免费观看| 国产精品.久久久| 久99久视频精品免费| 久久精品夜夜夜夜夜久久蜜豆| 美女被艹到高潮喷水动态| 91久久精品国产一区二区成人| 一进一出抽搐gif免费好疼| 精品99又大又爽又粗少妇毛片| 亚洲av成人av| 婷婷亚洲欧美| 最近的中文字幕免费完整| 久久精品国产亚洲av香蕉五月| 亚洲人成网站在线观看播放| 中国美女看黄片| 日韩精品青青久久久久久| 亚洲国产精品国产精品| 国产精品精品国产色婷婷| 国产熟女欧美一区二区| 久久久国产成人精品二区| 亚洲一级一片aⅴ在线观看| 中文资源天堂在线| 亚洲最大成人手机在线| 亚洲成人中文字幕在线播放| 一级毛片aaaaaa免费看小| 99久久无色码亚洲精品果冻| 男人舔奶头视频| 午夜福利成人在线免费观看| 色噜噜av男人的天堂激情| 又爽又黄a免费视频| 99久久精品热视频| 最近最新中文字幕大全电影3| 日韩三级伦理在线观看| 日韩视频在线欧美| 天堂√8在线中文| 啦啦啦观看免费观看视频高清| 久久精品影院6| 免费观看精品视频网站| 久久久久久久久久黄片| 老师上课跳d突然被开到最大视频| 精华霜和精华液先用哪个| 欧美成人免费av一区二区三区| 少妇丰满av| av黄色大香蕉| 亚州av有码| 黄色欧美视频在线观看| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 九色成人免费人妻av| 国产激情偷乱视频一区二区| 淫秽高清视频在线观看| 在线天堂最新版资源| 免费看av在线观看网站| 哪里可以看免费的av片| 国产精品一区二区在线观看99 | 国产精品爽爽va在线观看网站| 你懂的网址亚洲精品在线观看 | 99久久九九国产精品国产免费| 三级毛片av免费| 国产精品国产三级国产av玫瑰| 国产成人freesex在线| 国产成人aa在线观看| 91麻豆精品激情在线观看国产| 亚洲在久久综合| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 免费看光身美女| 狂野欧美白嫩少妇大欣赏| 国产精品国产三级国产av玫瑰| 久久鲁丝午夜福利片| 少妇人妻精品综合一区二区 | 精品久久久久久成人av| 丰满乱子伦码专区| av.在线天堂| 校园春色视频在线观看| 乱人视频在线观看| 1000部很黄的大片| 爱豆传媒免费全集在线观看| 看免费成人av毛片| 韩国av在线不卡| 国产精品一区二区三区四区久久| h日本视频在线播放| 欧美3d第一页| 国产精品日韩av在线免费观看| 熟女电影av网| 欧美极品一区二区三区四区| 乱人视频在线观看| 91久久精品电影网| 国内久久婷婷六月综合欲色啪|