Importance of the single-particle continuum in BCS pairing with a pseudostate basis

Importance of the single-particle continuum in BCS pairing with a pseudostate basis J. A. Lay1,C.E.Alonso 2,L.Fortunato 1,and A. Vitturi1 1Dipartimento di Fisica e Astronomia “Galileo Galilei”, Universit`a di Padova, and INFN, Sezione di Padova, via Marzolo, 8, I-35131 Padova, Italy 2Departamento de F´ısica At´omica, Molecular y Nuclear, Facultad de F´ısica, Universidad de Sevilla, Apartado 1065, 41080 Sevilla, Spain Abstract In a recent work [arXiv:1510.03185] the use of the Transformed Harmonic Oscillator (THO) basis for the discretization of the singleparticle continuum into a Generalized Bardeen-Cooper-Schrieffer (BCS) formalism was proposed for the description of weakly bound nuclei. We make use of the flexibility of this formalism to study the evolution of the pairing when the nucleus becomes more and more weakly bound. Specifically we focus on the evolution of the occupation of the different partial waves in 22O when the Fermi level approaches zero. Present developments in nuclear facilities allow to study new exotic nuclei, approaching the proton and neutron driplines for both light and medium mass regions of the Segr`e Chart. Nuclei close to the dripline are weakly bound and, therefore, affected by the continuum, situations where the last nucleon or nucleons are no longer bound and can explore distances far from the rest of the nucleus. The effects of the continuum have been traditionally explored for light nuclei where the proximity of the driplines provides a rich variety of weakly 201 ,0601 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0601 epj 5 1177 NN 8 8  © The Authors, published by EDP Sciences SI . This is an open access article distributed under the terms of the Creative Commons Attribution License 4.0 (http://creativecommons.org/licenses/by/4.0/). F - bound nuclei and experimental data. In the medium mass region, several attempts for generalizing mean field approaches have been proposed with different results [1,2]. Regarding the Bardeen-Cooper-Schrieffer (BCS) approximation, the single-particle continuum is first discretized and then treated normally splitting resonant and background contributions or including a continuum single-particle level density [2]. In a recent work [3] we studied the possibility of treating background and resonant parts of the continuum on the same footing. In order to do so, the single-particle continuum is discretized using a square-integrable basis. The pseudostates or eigenstates obtained through the diagonalization of the matrix elements of the single-particle Hamiltonian are expected to give a good representation of the continuum. Then, the pairing Hamiltonian is written in terms of the pair creator and annihilator operators built on top of these pseudostates. The strength of the pairing is determined by a densitydependent delta interaction and matrix elements are calculated separately for each couple of pairs. Once we set a maximum energy, it is possible to solve the BCS equations. The basis chosen for this purpose is the Transformed Harmonic Oscillator (THO) basis. It is constructed by applying a Locale Scale Tranformation to the Harmonic Oscillator basis of the form [4]: s(r)=⎡ ⎣ 1 1 rm+1 γ√rm⎤ ⎦ 1 m ,(1) that depends on the parameters mand γ. Even though we increase the number of parameters with respect to the simpler HO, two main advantages are obtained: the asymptotic behaviour of the functions is now exponential, rather than gaussian as in the HO case, and γcan be used to tailor the basis to the different cases. This is of great help when assessing the importance of the low-energy continuum since, by reducing the value of γ, the pseudostates tend to concentrate at low energies, thus optimizing the calculation. Already in [3], an interesting increase of the occupation of the nonresonant low-energy continuum was found for extremely weakly bound situations. In this work, we will concentrate on the evolution of the total occupation when we vary the strength of the single-particle continuum. Remember that, in the BCS theory, the quantities v2and u2stand for the occupation and deoccupation of the different level. We will choose the same case and parameters chosen in [3], so that we will look to the occupations in 22O when varying the strength of the central single-particle potential V. 201 ,0601 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0601 epj 5 1177 NN 8 8 2 48 50 52 54 56 V (MeV) 10-4 10-3 10-2 10-1 100 6(2j+1)v2 Total d5/2 bound s1/2 bound d5/2 continuum s1/2 continuum d3/2 p1/2 p3/2 Figure 1: Occupation corresponding to the six valence neutrons considered in 22O for the different partial waves. We separate the d5/2(solid line) and s1/2(dashed line) bound contribution from the corresponding d5/2(dotted line) and s1/2(dotted dashed line) continuum contribution. The thick dashed line corresponds to the total occupation as the sum of the different contributions considered which, by construction, is always the number of valence neutrons. In [3], the pairing interaction between the different pairs is calculated by using a density dependent delta interaction (DDDI) fitted to reproduced the neutron separation energy in 18O, and then renormalized along the isotopic chain following a 1/A reduction. In fig. 1 we show the occupations in 22O for all the partial waves considered, calculated for different depths of the central single-particle potential V. We have also considered separately the bound and continuum parts of the occupation in the d5/2−and s1/2−waves. It is commonly stated that the contribution of the continuum is mainly governed by the resonances present in each case. In our case, this fact should be reflected in the dominance of the d3/2contribution where a low-lying resonance can be found. This is true for a well-bound 22O, which in fig. 1 corresponds to the right part of the figure around V= 57 MeV. It is interesting to see how the occupation decreases when the strength of the potential is reduced. However, the occupation of the resonant part of the continuum is still larger by more than one order of magnitude than the non-resonant part of the continuum. This situation changes for really weakly-bound situations. A sudden increase of the s1/2continuum contribution is found, even though only for 201 ,0601 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0601 epj 5 1177 NN 8 8 3 an extreme situation. One has to keep in mind than the Fermi level for V=47.5 MeV is practically zero. On the other hand, the bound part of the s1/2contribution is more or less constant with a slight increase. If we look at the remaining contributions, only the d5/2continuum part is also constantly increasing when reducing the strength V. However, this represents a modest variation if we compare with the behaviour of the s1/2 waves. Nevertheless, one has to remind that we are using a DDDI, which can overestimate the relative importance of the s−wave in detritment of contributions with larger orbital angular momentum. A more realistic treatment of the pairing interaction could lead to a larger increase of the occupation of the d5/2continuum, at least in comparison with the s1/2contribution. The fact that only those continuum components with a bound state of the same orbital and total angular momentum show such an increase is telling us that this effect is not related with the extension of the wave functions as in the case of the formation of a halo nucleus. In such a case, p3/2occupation should show an increase since the extension of p−waves also goes to infinity when the binding energy goes to zero. Instead, they decrease as we can see in fig.1. In conclusion, non-resonant continuum contributions to the ground state of a nucleus cannot be neglected when dealing with weakly bound nuclei. The importance of the resonant configuration decreases consistently when the system becomes more and more weakly bound until a situation is reached where is no longer the dominant part of the single-particle continuum. However, we also show that for not so extreme situations, it is the resonant part that carries almost all the importance of the continuum. Acknowledgments This work has been supported by MIUR research fund PRIN 2009TWL3MX, by the Spanish Ministerio de Econom´ıa y Competitividad under Project FIS2014-53448-C2-1-P by Junta de Andaluc´ıa under group number FQM- 160 and Project P11-FQM-7632, and by the Consolider-Ingenio 2010 Program CPAN (CSD2007-00042). The research leading to these results has received funding from the European Commission, Seventh Framework Programme (FP7/2007-2013) under Grant Agreement no600376. J.A.L. is a Marie Curie Piscopia fellow at the University of Padova. 201 ,0601 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0601 epj 5 1177 NN 8 8 4 References [1] J. Dobaczewski et al.,Phys. Rev. C,53 (1996) 2809. [2] R. Id Betan, Nucl. Phys. A,879 (2012) 14. [3] J.A.Layet al., arXiv:1510.03185 Preprint, 2015. [4] J.A.Layet al.,Phys. Rev. C,82 (2010) 024605. 201 ,0601 (2016) EPJ Web of Conferences DOI: 10.1051/ conf/201611 0601 epj 5 1177 NN 8 8 5