arXiv:1005.3183v1 [nucl-th] 18 May 2010 where eairo h eto rpln.I hsfrua the formula, this In line. energy drip binding neutron nuclear the of behavior nteform the in number sbsdo iwo ulia nopesbespherical incompressible density as nuclei uniform of of drops view liquid a on based is of prediction the line. affects drip turn dependence neutron in density the energy the this symmetry how In the investigate of [3]. will energy the we of symmetry through Letter, equation the matter of the nuclear by dependence asymmetric density shape controlled of and be size (EOS) to the state shown that are remark nuclei We of [2]. stars evolution neutron rotational of and are thermal which observed vortices, the neutron to relevant well superfluid as of and crust dynamics thermal the the in as the matter of control typical properties a to transport are considered electric stars are neutron of and neutral crust example neutron dense the the in in present Nuclei beyond be Even matter. can nuclei [1]. however, nuclei line, respectively.drip light for state, acces- only ground marginally sible is (pro- the line drip neutrons in neutron two the nucleus Experimentally, and a one from remove tons) to separa- required (proton) energy two-neutron and energies, one- (proton) tion the neutron nu- the are superheavy study line to to drip quantities light key from The drip ranging clides. proton chart or nuclear neutron a the against in probe nuclei experimentally atomic to of able stability be might one facilities, a a a sym surf Coul hoeial,aWisakrBtems oml,which Weizs¨acker-Bethe formula, a mass Theoretically, beam ion radioactive of developments recent to Thanks A [( Z N 2 2 / E /A 3 A vol − − 1 n hrenumber charge and stesraeeeg,and energy, surface the is 4 / Z E dacdSineRsac etr aa tmcEeg Agen Energy Atomic Japan Center, Research Science Advanced 3 ASnmes 16.f 21.10.Dr 21.65.Ef, neutron- numbers: of s PACS neutron- property Such surface a the increases. to from energy shift mainly symmetry coming appreciably feature the to mat of nuclear shown derivative asymmetric t is density of us masses state allows these of that from equation model the nuclear on macroscopic dependent a using by energy, = ) B S /A steCuobeeg.Te,teone- the Then, energy. Coulomb the is eivsiaehwtenurndi iei eae oteden the to related is line drip neutron the how investigate We n = a 3 ( eto rpln n h qaino tt fncermatter nuclear of state of equation the and line drip Neutron ] p eateto aua cec,KciUiest,Akebono- University, Kochi Science, Natural of Department 2 vol E ) A vol and A 2 stesmer energy, symmetry the is IE ihn etr IE,Hrsw,Wk,Siaa351-0 Saitama Wako, Hirosawa, RIKEN, Center, Nishina RIKEN + E stevlm energy, volume the is S B E 2 n sym switna ucino mass of function as written is 1 (2 auioOyamatsu Kazuhiro eateto ua nomtc,AciSuuouUniversi Shukutoku Aichi Informatics, Human of Department p aaue aaueco ih-u,Aci4019,Japan 480-1197, Aichi Aichi-gun, Nagakute-cho, Nagakute, ) Z + hc orsodt an to correspond which , n E o eto number neutron (or 0 surf rvdsastandard a provides , + E Coul E Dtd coe 1 2018) 31, October (Dated: , Coul E E 1 , sym surf 2 e Iida Kei , N (1) = = = ) er energy metry hc sdfie as defined is which h ymtyeeg sadniysmer coefficient symmetry density a is of energy dependence symmetry impor- density the The the characterizing seriously. the neutron taken parameter of tant be the dependence to need density predicting energy the In symmetry in uncertainties line, [11–13]). neutron-rich drip Refs. involving collisions (e.g., heavy-ion nuclei Refs. (e.g., and 5]) resonances 10]), [4, dipole [9, Refs. [6–8]), Refs. (e.g., (e.g., masses radii the nuclear and to of relevant is dependence it isospin because attention much attracts ergy eciigtemse n ai fusal ulii a in for nuclei used unstable be of can radii and and masses nuclei the stable known describing the of reproduce to properties as way global a such in constructed was nuclear the on line drip this how chart. neutron address depen- the we Here isospin affects [5]. and dependence tension density surface the the of from dence understood be can S where ie h eie rpln scoeto close is line drip derived The line. hw ohv napeibedpnec on dependence and appreciable model an nuclear have macroscopic to a shown nuclei from neutron-rich calculated extremely were of Masses matter. nuclear n on,btoenest obyn q 2 ytaking den- by nuclear (2) start- the Eq. from of good sity beyond deviation a go nonnegligible be to account may needs into one This but light chart. point, the ing nuclear in the the except between of terms, region competition energy symmetry the and by volume controlled basically is eto eaaineeg a eeautdas evaluated be can energy separation neutron 2 n , 3 eety h est eedneo h ymtyen- symmetry the of dependence density the Recently, ebgnwt arsoi oe fnce 6,which [6], nuclei of model macroscopic a with begin We S n iouiKoura Hiroyuki and , ie mohdbhvo ftenurndrip neutron the of behavior smoothed a gives 0 = n x = ≈ ihnuclei. rich iti lal enfrlgtnce,a nuclei, light for seen clearly is hift = n − aclt ula assi way a in masses nuclear calculate o e.Tenurndi ieobtained line drip neutron The ter. ∂E 0 ihsd nancercata the as chart nuclear a in side rich y oa,Iaai3919,Japan 319-1195, Ibaraki Tokai, cy, ∂N a . Z/A iydpnec ftesymmetry the of dependence sity vol h,Kci7082,Japan 780-8520, Kochi cho, B S
( − Z n stepoo rcin h condition The fraction. proton the is eedn ntedensity the on dependent ) a sym L 3 = (1 9,Japan 198, − 2 n ty, , 4 0 4 ( x dS/dn 2 ) − 3 2 ) A a n = surf 1 n / 3 0 x + ihtesym- the with 0 = a 3 Coul A n . L ,which 3, 4 which , fbulk of / Z 3 2 , (2) L , 2 manner that is dependent on the EOS of nuclear matter. 100 Exp. (CN2004) This model can be summarized as follows: mass-measured (AWT03) (i) We set the bulk energy per nucleon as 80 stable for beta-decay
drip line 2 2 2/3 EOS C 3¯h (3π ) 5/3 5/3 60 EOS G w = (nn + np ) KTUY05 10mnn 2 2 +(1 − α )vs(n)/n + α vn(n)/n, (3) 40 where Proton number 20 3 2 a2n vs = a1n + (4) 1+ a3n 0 0 20 40 60 80 100 120 140 and Neutron number 3 FIG. 1: (Color online) The neutron and proton drip lines 2 b2n obtained from the EOS models C and G by using the macro- vn = b1n + (5) 1+ b3n scopic nuclear model and from a contemporary mass formula [17]. The regions filled with squares correspond to empirically are the potential energy densities for symmetric nuclear known nuclides [18, 19]. matter and pure neutron matter, nn and np are the neutron and proton number densities, n = nn + np, α = (nn − np)/n is the neutron excess, and mn is the is the Coulomb energy, and mp is the proton mass. neutron mass. A set of expressions (3)–(5) is one of the (iii) For simplicity we use the following parametrization simplest that reduces to a usual expansion [14] for the nucleon distributions ni(r) (i = n,p):
3 ti K0 2 L 2 r w = w0 + 2 (n − n0) + S0 + (n − n0) α (6) nin 1 − , r
Z S0 ≈ B + CL, (12) is the gradient energy with adjustable constant F0, where B ≈ 28 MeV and C ≈ 0.075. 2 ′ e ′ n (r)n (r ) We proceed to obtain the neutron and proton drip E = d3r d3r p p (10) C 2 |r − r′| lines from the macroscopic nuclear model. For various Z Z 3
15 120 EOS C 10 neutron drip EOS G 80 neutron drip 5 40 proton drip 0 proton drip
0 neutron number diff. neutron - proton number 0 20 40 60 80 100 20 40 60 80 100 proton number proton number FIG. 2: (Color online) N −Z obtained for nuclides at neutron FIG. 3: (Color online) Differences in the neutron number N and proton drip from the EOS models C and G. obtained for nuclides at neutron and proton drip between the calculations from the EOS models C and G. sets of L and K0, we first evaluate the binding energy 0.35 EB of nuclei in the ground state by minimizing the en- ergy (7) for fixed N and Z. We then draw the neu- 0.30 tron (proton) drip line by identifying nuclides at neu- neutron drip tron (proton) drip with those neighboring to nuclides for 0.25 which Sn = EB (Z,N)−EB(Z,N −1) (Sp = EB(Z,N)− EOS C
proton fraction 0.20 EOS G EB(Z − 1,N)) and S2n = EB(Z,N) − EB (Z,N − 2) (S2p = EB(Z,N) − EB (Z − 2,N)) are positive and be- 0.15 yond which at least one of them is negative. The results 0 20 40 60 80 100 obtained from the two extreme EOS models denoted as proton number EOS C (L = 146 MeV and K0 = 360 MeV) and EOS FIG. 4: (Color online) Proton fraction for nuclides at neutron G (L =5.7 MeV and K0 = 180 MeV) are shown in Fig. drip obtained from the EOS models C and G. 1, together with the empirically known nuclides [18, 19] and the prediction from a contemporary mass formula [17]. We remark that inclusion of the condition for S2n shows a stronger dependence on L for lighter nuclei, as and S2p in addition to Sn and Sp in drawing the drip lines shown in Fig. 4. This is advantageous because heavier ra- makes only a little difference in the case of the present dioisotopes are more difficult to produce in experiments. model calculations, while being essential in the case of the The L dependence of the neutron drip line as obtained prediction from the mass formula because of the Wigner, above can be understood within the framework of a com- shell, and even-odd terms included therein. We remark pressible liquid-drop model in which nuclei in equilibrium that the rms deviations of the calculated masses from the are allowed to have a density different from the saturation measured values [19] are about 3 MeV, which is of order density n0 of symmetric nuclear matter. By following a the deviations obtained from a Weizs¨acker-Bethe mass line of argument of Ref. [5], we first add the surface sym- 2/3 2 formula. metry term, assymA [(N −Z)/A] , to the mass formula We find from Fig. 1 that the obtained neutron drip (1) based on an incompressible liquid-drop model. This lines show an appreciable L dependence, while the pro- surface symmetry term gives rise to additional contribu- ton ones do not. This is reasonable because nuclides at tion, neutron (proton) drip are far away from (near) N = Z 2a δS = − ssym (1 − 2x)(1 + 4x), (13) (see Fig. 2). The neutron drip line does shift to a neutron- n 3A1/3 rich side as L increases, a feature that will be discussed later in terms of a compressible liquid-drop model. We to the neutron separation energy (2). Next, we consider remark that the EOS dependence of the obtained drip the density-dependent surface tension [20], lines comes predominantly from L because of negligible 2 nin − n0 K dependence of the calculated masses [5]. σ(nin, αin)= σ0 1 − Csymα + χ , (14) 0 in n In order to see the L dependence more clearly, we plot 0 in Fig. 3 the difference in the neutron number of nuclides where nin and αin are the density and neutron ex- at neutron and proton drip between the calculations from cess inside a liquid drop, σ0 = σ(n0, 0), Csym is the EOS models C and G. The difference shows only a the surface symmetry energy coefficient, and χ = weak dependence on Z both in the case of neutron and (n0/σ0)∂σ/∂nin|nin=n0,αin=0. By taking a limit of van- 2 2/3 proton drip. This indicates that the L dependence can ishing compressibility, one obtains 4πσ0R = asurf A 2 2/3 be seen more clearly for lighter nuclei. In fact, the cor- and 4πσ0CsymR = −assymA with the liquid-drop ra- responding proton fraction of nuclides at neutron drip dius R. Typically, fitting to the empirical mass data 4
−2 yields σ0 ≈ 1 MeV fm and Csym = 1.5–2.5. As we play an important role in Sn. In fact, the magnitude of shall see below, nonvanishing compressibility effectively even-odd staggering in Sn could be comparable to that of introduces the L dependence into the surface symmetry change in Sn due to uncertainties in L. This implies that term through the parameter χ characterizing the den- some kind of smoothing would be required in deriving sity dependence of the surface tension. The value of χ is information about the EOS from future empirical data poorly known, but likely to be positive [5]. For example, on neutron drip. Shell corrections would further compli- χ =4/3 for the Fermi gas model. cate such derivation, but are intriguing in the context of If one ignores Coulomb and surface corrections, the magicity of nuclei close to the neutron drip line [21]. equilibrium density and energy per nucleon of a liquid- drop, ns and ws, can be evaluated from Eq. (6) as
2 ws = w0 + S0α (15) We are grateful to Dr. A. Kohama for useful discus- sion and late Prof. M. Uno for continuous encourage- and ments. Authors K.O. and K.I. acknowledge the hospital- 3n0L 2 ity of the Yukawa Institute for Theoretical Physics during ns = n0 − α . (16) K0 the workshop “New Frontiers in QCD 2010,” where this work was initiated. This work was supported in part by Strictly speaking, expressions (15) and (16) are applica- Grant-in-Aid for Scientific Research through Grant No. ble only for nearly symmetric nuclear matter. We nev- 19740151, which was provided by the Ministry of Educa- ertheless use these expressions for the purpose of char- tion, Culture, Sports, Science, and Technology of Japan. acterizing the liquid-drop properties because the typical value of α along the neutron drip line is of order 0.35–0.4, considerably smaller than unity. By substituting ns, Eq. (16), into nin in Eq. (14), we thus obtain [1] M. Notani et al., Phys. Rev. C 76, 044605 (2007); T. 3Lχ 2 Baumann et al., Nature (London) 449, 1022 (2007). σ(ns, αin)= σ0 1 − Csym + αin . (17) K0 [2] C.J. Pethick and D.G. Ravenhall, Annu. Rev. Nucl. Part. Sci. 45, 429 (1995). 2/3 75 This result allows one to identify assymA with [3] K. Oyamatsu and K. Iida, Phys. Rev. C , 015801 2 (2007). −4πσ0R (Csym +3Lχ/K0) and hence to conclude that [4] P. Danielewicz and J. Lee, Nucl. Phys. A818, 36 (2009). with increasing L, S increases through the surface sym- n [5] K. Oyamatsu and K. Iida, Phys. Rev. C 81, 054302 metry contribution (13) for x< 1/2. (2010). This conclusion is consistent with the L dependence of [6] K. Oyamatsu and K. Iida, Prog. Theor. Phys. 109, 631 the neutron drip line shown in Fig. 1 because any positive (2003). 67 corrections to Sn act to enhance the stability of nuclei [7] A.R. Bodmer and Q.N. Usmani, Phys. Rev. C , 034305 against neutron emission. It is important to note that (2003). [8] M. Warda, X. Vi˜nas, X. Roca-Maza, and M. Centelles, the bulk symmetry term gives rise to a negative contri- 80 bution to S for x< 1/2, as shown in Eq. (2), and that Phys. Rev. C , 024316 (2009). n [9] A. Klimkiewicz et al., Phys. Rev. C 76, 051603(R) the parameter asym corresponds to the symmetry energy (2007). coefficient S0, which in turn is related to L by the relation [10] L. Trippa, G. Col`o, and E. Vigezzi, Phys. Rev. C 77, (12) obtained from fitting to empirical masses and charge 061304(R) (2008). radii of stable nuclei. Since one obtains a larger asym for [11] D.V. Shetty, S.J. Yennello, and G.A. Souliotis, Phys. larger L, the effect of a tends to decrease S with L Rev. C 76, 024606 (2007). sym n 464 and hence to facilitate neutron drip. However, this effect [12] B.A. Li, L.W. Chen, and C.M. Ko, Phys. Rep. , 113 (2008). is relatively small compared with the above-mentioned [13] M.B. Tsang et al., Phys. Rev. Lett. 102, 122701 (2009). effect of assym. This is consistent with the fact that the [14] J.M. Lattimer, Ann. Rev. Nucl. Part. Sci. 31, 337 (1981). L dependence is clearer for lighter nuclei. [15] K. Oyamatsu, Nucl. Phy. A561, 431 (1993). In summary, we have investigated the influence of the [16] B. Friedman and V.R. Pandharipande, Nucl. Phys. density dependence of the symmetry energy on the drip A361, 502 (1981). lines by using a macroscopic nuclear model that depends [17] H. Koura, T. Tachibana, M. Uno, and M. Yamada, Prog. 113 explicitly on the EOS of nuclear matter. We find that an Theor. Phys. , 305 (2005). [18] Chart of the Nuclides 2004, compiled by T. Horiguchi et L dependence appears appreciably in the neutron drip al., JNDC and Nuclear Data Center, JAERI (2005). line and it is clearer for lighter nuclei, a feature com- [19] G. Audi, A.H. Wapstra, and C. Thibault, Nucl. Phys. ing mainly from the surface property through the den- A729, 337 (2003). sity and neutron excess dependence of the surface ten- [20] K. Iida and K. Oyamatsu, Phys. Rev. C 69, 037301 sion. We note that our calculations do not include even- (2004). odd or shell corrections. The even-odd corrections would [21] A. Ozawa et al., Phys. Rev. Lett. 84, 5493 (2000).