arXiv:1709.03355v1 [nucl-ex] 11 Sep 2017 fnmes h rtsre ,8 0–i trbtdto attributed is second – the 20 while 8, potential, 2, (HO) – oscillator series harmonic and the first stable series The in different nuclei. numbers. two know of magic-number of we consist the as nuclei, to numbers, naturally-occurring magic leads shells, nuclear nuclear hand, The the other particu- of are the filling – The on gases noble chemically. the stable filled as larly with known Atoms – shells atoms. of electron to analogy structure in shell collec- structure, electronic shell (known the exhibit neutrons – nucleons) and as protons tively of consisting tems tmcnce h nt unu aybd sys- many-body quantum finite the – nuclei Atomic 10 16 .R Jansen, R. G. aut fPyisadEgneig oCiMn iyUnivers City Minh Chi Ho Engineering, and Physics of Faculty ainlCne o opttoa cecs a ig Nati Ridge Oak Sciences, Computational for Center National ..Tran, D.T. .Sakaguchi, H. ..Geng, L.S. .Momota, S. 5 eateto hsc,Cleeo uaiisadSciences, and Humanities of College Physics, of Department 4 .Geissel, H. eateto hsc n srnm,Uiest fTenness of University Astronomy, and Physics of Department asso ih uli Performing e nuclei. radii, light distribution of masses point-proton of analysis systematic hrlncenncenadtrencenfre,wihar by which explained forces, be three-nucleon can and nucleon-nucleon closure chiral subshell and radii distribution Ooiiae ai ubr(uselcoue tteproto the at closure) 194 pres (subshell we in Here, number force questions. magic open coupling SO-originated still are (SO) origins fingerpr spin-orbit its ext numbers, remarkable with phenomenological associated most a neutrons Its and of un protons our of for forces. numbers’ ‘magic guide nuclear important underlying an the provides potential, average an 17 8 vdnefor Evidence colo hsc n ula nryEgneig ehn U Beihang Engineering, Energy Nuclear and Physics of School h ula hl tutr,wihoiiae ntenearly the in originates which structure, shell nuclear The srnm n hsc eatet an aysUniversity Mary’s Saint Department, Physics and Astronomy 2 nttt fPyis ita cdm fSineadTechno and Science of Academy Vietnam Physics, of Institute 3 18 hsc iiin a ig ainlLbrtr,OkRidge Oak Laboratory, National Ridge Oak Division, Physics ,2 1, 8 12 15 nttt fMdr hsc,CieeAaeyo cecs L Sciences, of Academy Chinese Physics, Modern of Institute 6 3 16, 19 .Terashima, S. S emotznrmfu cwroefrcug D-64291 f¨ur Schwerionenforschung, Helmholtzzentrum GSI aeIooeSinePoet nttt o ai Science, Basic for Institute Project, Science Isotope Rare 1 ..Ong, H.J. 9 2 13 12, .Nagae, D. hmNo hc nvriyo eiie oCiMn 000 V 700000, Minh Chi Ho Medicine, of University Thach Ngoc Pham .Scheidenberger, C. .Kanungo, R. 23 25 24 6 ainlAtooia bevtr fJpn oy 181-858 Tokyo Japan, of Observatory Astronomical National 11 7 eateto hsc,Uiest fTko oy 113-0033 Tokyo Tokyo, of University Physics, of Department eateto hsc,NiaaUiest,Niaa950-21 Niigata 305-8 University, Ibaraki Niigata Tsukuba, Physics, of of University Department Physics, of Institute 26 eateto hsc,KooUiest,Koo6680,J 606-8502, Kyoto University, Kyoto Physics, of Department nvriyo rnne,94 AGoign h Netherland The Groningen, AA 9747 Groningen, of University eateto hsc,OaaUiest,Oaa5004,J 560-0043, Osaka University, Osaka Physics, of Department ..Harakeh, M.N. ylto nttt,TxsA Texas Institute, Cyclotron 19 Z 1, 22 oh nvriyo ehooy oh 8-52 Japan 782-8502, Kochi Technology, of University Kochi 13 20 20 1 21 6‘ai ubr nnurnrc abnisotopes carbon neutron-rich in number’ ‘magic =6 ∗ oy nvriyo cec,Cia2881,Japan 278-8510, Chiba Science, of University Tokyo 14 8 CP sk nvriy sk 6-07 Japan. 567-0047, Osaka University, Osaka RCNP, utsLei nvriy 59 ise,Germany. Giessen, 35392 University, Liebig Justus IE ihn etr atm 5-18 Japan. 351-0198, Saitama Center, Nishina RIKEN ogNiUiest,Dn a 10,Vietnam. 81000, Nai Dong University, Nai Dong .Hagen, G. V etrfrAvne aito Technology, Radiation Advanced for Center KVI ..Nguyen, N.D. .Tanihata, I. 17 .Kawabata, T. 2 13 12, 14 binitio ab Dtd etme ,2017) 4, September (Dated: ,4 3, .Hashimoto, T. .Tanaka, J. ,8 1, .D Morris, D. T. 21 ..Nguyen, T.T. acltoson calculations .Nishimura, D. 7 ..Khiem, L.H. & nvriy ea 74,USA. 77840, Texas University, M 1 .Takechi, M. n he-ulo ocs(Ns.Atog h theoret- the Although (3NFs). atomic forces of three-nucleon structure and the describe to from nuclei difficulty the microscopic unresolved to remained the have due However, force coupling SO Prize. the of Nobel origins the Jensen Goeppert-Mayer and won and numbers, magic the completely fancenadisobtlaglrmmnu,adthe and momentum, spin of angular intrinsic orbital the introduction on its so-called the and depends (SO) was nucleon that spin-orbit force a It the a of to – 1). force due SO is Fig. this – (see 126 and force 82, coupling 50, 28, – one 15 ,4 3, nlLbrtr,OkRde ense 73,USA 37831, Tennessee Ridge, Oak Laboratory, onal ..Hoang, T.H. etoantctasto ae n atomic and rates transition lectromagnetic otdi h unu chromodynamics. quantum the in rooted e t fSine oCiMn iy720 Vietnam. 70250, City Minh Chi Ho Science, of ity ,10 9, n xeietleiec o h smallest the for evidence experimental ent 14 .Aoi, N. 22 h tt-fteatncerter with theory nuclear state-of-the-art the io nvriy oy 5-50 Japan. 156-8550, Tokyo University, Nihon 2 , asaiiy lhuhteintroduction the Although stability. ra 15 ubr6in 6 number n .Otsuka, T. .Ayyad, Y. ..Lin, W.P. j esadn fncersrcueand structure nuclear of derstanding n steeitneo h so-called the of existence the is int epdepantencermagic nuclear the explain helped 9 e nxil,Tnese396 USA 37996, Tennessee Knoxville, ee, ,w hwta h bevdproton observed the that show we C, - needn oino ulosin nucleons of motion independent j 25 binitio ab aia,N 3 C,Canada. 3C3, B3H NS Halifax, , iest,Biig109,China 100191, Beijing niversity, opigshm hthle xli 1 2] [1, explain helped that scheme coupling 1 .Wada, R. oy ao 00,Vietnam. 10000, Hanoi logy, ense 73,USA. 37831, Tennessee , .Suzuki, T. ajo 44,Korea. 34047, Daejeon ,2 1, amtd,Germany. Darmstadt, nhu700,China. 730000, anzhou 18 .Ideguchi, E. 1 ula hois[–]wt w-(NN) two- with [3–5] theories nuclear 23 7,Japan. 571, ..Chan, P.Y. 1 Japan. 81, .Matsuta, K. ,Japan. 8, .Ozawa, A. 13 8 26 18, Japan. , − apan. ,6 5, apan 20 s. ietnam bandfrom obtained C n .Yamamoto T. and .Kanada-En’yo, Y. 1 1 24 .Inoue, A. 11 .Fukuda, M. ..Ren, P.P. .Mihara, M. 1 18 11 1 7 11 2 ical study [6] of the SO splitting of the 1p and 1p 50 1/2 3/2 1g 9/2 15 2p single-particle states in N has suggested possible roles 1/2 2p 1f 5/2 of two-body SO and tensor forces, as well as three-body 1f 2p 3/2 28 forces, the discovery of the so-far elusive SO-type magic 20 1f 7/2 20 number 6 is expected to offer unprecedented opportuni- 2s 1d 3/2 1d ties to understand its origins. 2s 1/2 1d 5/2 In her Nobel lecture, Goeppert-Mayer had mentioned 8 8 1p 1/2 the magic numbers 6 and 14 – which she described as 1p 6 ? 1p 3/2 2 “hardly noticeable” – but surmised that the energy gap 2 between the 1p1/2 and 1p3/2 orbitals due to the SO force 1s 1s 1/2 harmonic spin-orbit is “fairly small” [7]. That the j-j coupling scheme ap- + l 2 + pears to fail in the p-shell light nuclei was discussed and oscillator coupling attributed to their tendency to form clusters of nucle- FIG. 1. Nuclear shell structure for a harmonic oscillator po- ons [8]. Experimental and theoretical studies in recent tential plus a small orbital angular momentum (l2) term. The decades, however, have hinted at the possible existence right diagram shows the splitting of the single-particle orbitals of the magic number 6 in some semimagic unstable nu- by an additional spin-orbit coupling force. clei, each of which has a HO-type magic number of the opposite type of nucleons. For instance, possible sub- shell closures have been suggested in 8He [9, 10], 14C and the observations are supported by the ab-initio nuclear 14O [11]. Whether such subshell closures exist, and are model that employs the nuclear forces derived from the predominently driven by the shell closures of the counter- effective field theory of the quantum chromodynamics. part protons/neutrons are of fundamental importance. Although still not well established, the size of a nu- The isotopic chain of carbon – with six protons and cleus, which can be defined as the root-mean-square consisting of thirteen particle-bound nuclei, provides an (rms) radius of its nucleon distribution, is expected to important platform to study the SO splitting of the provide important insights on the evolution of the magic 1p1/2 and 1p3/2 orbitals. Like other lighter isotopes, numbers. Recently, an unexpectedly large proton rms ra- the isotopes of carbon are known to exhibit both clus- dius (denoted simply as ‘proton radius’ hereafter) was re- tering and single-particle behaviours. While the sec- ported [21], and suggested as a possible counter-evidence ond excited J π=0+ state in 12C – the famous Hoyle for the double shell closure in 52Ca [22]. Attempts to state and important doorway state that helps produce identify any emergence of non-traditional magic numbers 12C in stars – is well understood as a triple-alpha state, based on the analysis of the systematics of the exper- it seems that the effect of the alpha-cluster-breaking imental proton radii have been reported [11, 23]. For 1p3/2 subshell closure is important to reproduce the the 4
R culations without 3NFs underestimate the radii by about 10%, thus suggesting the importance of 3NFs. oxygen boron
neutron number
FIG. 2. Proton radii of carbon, beryllium, boron and oxygen isotopes. The red- and black-filled circles are, respectively, the proton radii from this and our recent work [24–27]. The open diamonds are the data from electron scattering and isotope- shift methods [11]. The small symbols connected with dashed and dotted lines are the predictions from the AMD [13] and RMF [28] models, respectively. The small blue symbols with solid and dash-dotted lines are the results from the ab-initio N isotones coupled-cluster calculations with NNLOsat [29] and the NN- 3 6 9 12 15 4 7 10 13 16 only interaction NNLOopt [30]. 5 8 11 14
For simplicity, we show only the results for 17−19C in 12−16 Table I; for results on C, see ref. [27]. Rp’s are the proton radii extracted using the Glauber model formu- lated in ref. [27]. The values for 17,18C are the weighted (d) mean extracted using σCC’s at the two energies, while the one for 19C was extracted using the high-energy data. In determining the proton radii, we have assumed harmonic- oscillator (HO)-type distributions for the protons in the Glauber calculations. The uncertainties shown in the 12 5 brackets include the statistical uncertainties, the exper- 8 10 imental systematic uncertainties, and the uncertainties 4 15
0 8 6 20 attributed to the choice of functional shapes, i.e. HO or = = 25 20 Z Z 25 Woods-Saxon, assumed in the calculations. 15 neutron number proton number10 5 To get an overview of the isotopic dependence, we com- 0 pare the proton radii of the carbon isotopes with those of cal the neighbouring beryllium, boron and oxygen isotopes. FIG. 3. Evolution of (a) Rp/Rp , (b) B(E2), and (c) (3) Figure 2 shows the proton radii for carbon, beryllium, ∆p (N, Z)−∆p with proton number up to Z=22 for all avail- boron and oxygen isotopes. The red- and black-filled able isotonic chains. Vertical dotted lines and thin-dashed circles are the data for 12−19C, beryllium and boron iso- denote positions of the traditional proton magic numbers and topes extracted in this and our previous work [24–27]. Z=6, respectively. (d) Two dimensional lego plot of (c). 4
It is interesting to note that Rp’s are almost constant where Sp(N,Z) is the one-proton separation energy. In throughout the isotopic chain from 12C to 19C, fluctuat- the absence of many-body correlations such as pair- ing by less than 5%. Whereas this trend is similar to the ings, Sp(N,Z) resembles the single-particle energy, and (3) one observed/predicted in the proton-closed shell oxygen ∆p (N,Z) yields the proton single-particle energy-level isotopes, it is in contrast to those in the beryllium and spacing or shell gap between the last occupied and first boron isotopic chains, where the proton radii change by unoccupied proton orbitals in the nucleus with Z protons as much as 10% (for berylliums) or more (for borons). (and N neutrons). To eliminate the effect of proton- It is also worth noting that most theoretical calculations proton (p-p) pairing, we subtract out the p-p pairing en- −1/2 shown predict almost constant proton radii in carbon and ergies using the empirical formula: ∆p = 12A MeV. oxygen isotopes. The large fluctuations observed in Be (3) Figure 3(c) shows the systematics of ∆p (N,Z) − ∆p for and B isotopes can be attributed to the development of even-Z nuclides. All data were evaluated from the bind- cluster structures, whereas the almost constant Rp’s for 12−19 ing energies [20]. Here, we have omitted odd-Z nuclides C observed in the present work may indicate an in- to avoid odd-even staggering effects. The cusps observed ert proton core, i.e. 1p3/2 proton-subshell closure. at Z=N for all isotonic chains are due to the Wigner ef- Examining the Z dependence of the proton rms radii fect [32]. Apart from the Z=N nuclides, sizable ‘gaps’ et al. along the N=8 isotonic chain, Angeli have pointed (> 5 MeV) are also observed at Z=6 for N=7 – 15, and out [11, 23] a characteristic change of slope (existence of at Z=8 for N=8 – 10 and 12 – 16, although the ‘gap’ ap- a ‘kink’), a feature closely associated with shell closure, pears to decrease from 18C to 22C. For clarity, we show at Z=6. Here, by combining our data with the recent the corresponding two-dimensional lego plot in Fig. 3(d). data [24–27], as well as the data from ref. [11], we plot By requiring a ‘magic’ nucleus to fulfil at least two the experimental Rp’s against proton number. To elimi- nate the smooth mass number dependence of the proton signatures in Fig. 3(a)–(c), we conclude that we have observed a prominent proton subshell closure at Z=6 in rms radii, we normalized all Rp’s by the following mass- 13−20 (3) 12 dependent rms radii [31]: C. Although the empirical ∆p for C is large (∼14 MeV), applying the prescription from ref. [33], we obtain cal −2/3 −4/3 1/3 about 10.7 MeV for the total p-p and p-n pairing energy. Rp = 3/5 1.15+1.80A − 1.20A A fm. p This estimated large pairing energy indicates possible sig- nificant many-body correlations such as cluster correla- Figure 3(a) shows the evolution of R /Rcal with pro- p p tions, which compete and coexist with the single-particle ton number up to Z=22 for all available isotonic chains. 1p subshell closure in 12C. Based on this fact, as well Each isotonic chain is connected by a solid line. For 3/2 as its non-minimum proton radius and relatively large clarity, the data for N=6 – 13 isotones are displayed in B(E2) value, we have excluded 12C from our conclusion. colours. For nuclides with more than one experimental value, we have adopted the weighted mean values. The It is surprising that the systematics of the proton radii, discontinuities observed at Z=10 and Z=18 are due to B(E2) values and the empirical proton ‘subshell gaps’ for the lack of experimental data in the proton-rich region. most of the carbon isotopes are comparable to those for Note the increase/change in the slope at the traditional proton-closed shell oxygen isotopes. To understand the magic numbers Z=8 and 20. Marked ‘kinks’, similar to observed ground-state properties, i.e. the proton radii those observed at Z=20, 28, 50 and 82 [23], are observed and ‘subshell gap’ of the carbon isotopes, we performed at Z=6 for isotonic chains from N=7 to N=13, indicat- ab-initio CC calculations on 14,15C using various state-of- ing a possible major structural change, e.g., emergence the-art chiral interactions. We employed the CC method of a subshell closure, at Z=6. in the singles-and-doubles approximation with perturba- The possible emergence of a proton subshell closure tive triples corrections [Λ-CCSD(T)] [34] to compute the at Z=6 in neutron-rich even-even carbon isotopes is ground-state binding energies and proton radii for the also supported by the small B(E2) values observed in closed-(sub)shell 14C. To compute 15C (1/2+), we used 14−20C [17–19]. Figure 3(b) shows the systematics of the particle-attached equation-of-motion CC (EOM-CC) B(E2) values in Weisskopf unit (W.u.) for even-even nu- method [35], and included up to three-particle-two-hole clei up to Z=22. All data are available in ref. [19]. Nuclei (3p-2h) and two-particle-three-hole (2p-3h) corrections. with shell closures manifest themselves as minima. Be- Figure 4 shows the binding energies as functions of the sides the traditional magic number Z=8, clear minima proton radii for (a) 14C and (b) 15C. The coloured with B(E2) values smaller than 3 W.u. are observed at bands are the experimental values; the binding ener- Z=6 for N=8,10,12 and 14 isotones. gies (red horizontal lines) are taken from ref. [20], while To further examine the possible subshell closure at proton radii are from ref. [27] (orange bands) and the Z=6, we consider the second derivative of binding en- electron scattering data [11] (green band). The filled ergies defined as follows: black symbols are CC predictions with the NN+3NF chi- ral interactions from ref. [36] labeled 2.0/2.0 (EM)(), (3) Z ∆p (N,Z) ≡ (−1) [Sp(N,Z) − Sp(N,Z + 1)], (1) 2.0/2.0 (PWA)(H), 1.8/2.0 (EM)( ), 2.2/2.0 (EM)(), 5
exp SM CC
exp SM CC
(3/2 -)
6 7 8 9 10 11 12 13 14
14 15 FIG. 4. Binding energies of (a) C and (b) C as functions FIG. 5. Empirical one-proton addition (ǫ+p) and removal of their proton radii. The coloured bands and red-horizontal (ǫ−p) energies (horizontal bars) for (a) carbon, and (b) oxy- lines are the experimental values. The green band repre- gen isotopes deduced from one-proton separation energies and − sents the proton radius from the electron scattering. The the excitation energies of the lowest 3/2 states in the odd- filled black symbols are the CC predictions with SRG-evolved even nitrogen isotopes. The dotted bars indicate the adopted NN+3NF chiral effective interactions at different NN/3NF values for the observed excited states in 19,21N, which have − cutoffs and NNLOsat, whereas the open symbols are the pre- been tentatively assigned as 3/2 [43]. Other experimental dictions with the NN-only EM and NNLOsat interactions. See data are taken from refs. [20, 44, 45]. The blue symbols are text for details. the shell model calculations using the YSOX interactions [46]. Results of the CC calculations with and without 3NFs are shown by the red-solid and red-dashed lines, respectively. 2.8/2.0 (EM)(N), and NNLOsat [29] (⋆). Here, the NN interactions are the next-to-next-to-next-to leading- order (N3LO) chiral interaction from ref. [37], evolved to lower cutoffs (1.8/2.0/2.2/2.8 fm−1) via the similarity- energies and the driplines of nitrogen and oxygen isotopes renormalization-group (SRG) method [38], while the 3NF have been suggested in recent theoretical studies. Here, is taken at NNLO with a cutoff of 2.0 fm−1 and ad- to shed light on the role of 3NFs on the observed ‘sub- justed to the triton binding energy and 4He charge ra- shell gap’, i.e. the SO splitting in the carbon isotopes, dius. The error bars are the estimated theoretical uncer- we investigate the evolution of one-proton separation en- tainties due to truncations of the employed method and ergies for carbon and oxygen isotopes. In Fig. 5, the hor- model space [39]. Note that the error bars for the bind- izontal bars represent the experimental one-proton ad- ing energies are smaller than the symbols. Depending on dition (ǫ+p) and removal (ǫ−p) energies for (a) carbon the NN cutoff, the calculated binding energy correlates and (b) oxygen isotopes deduced from one-proton sep- strongly with the calculated proton radius. In addition, aration energies and the excitation energies of the low- we performed the CC calculations with chiral effective est 3/2− states in the odd-even nitrogen isotopes. The interactions without 3NFs, i.e. the NN-only EM inter- dotted bars indicate the adopted values for the observed actions with NN cutoffs at 1.8 (#), 2.0 (2), 2.2 (♦) and excited states in 19,21N, which have been tentatively as- 2.8 fm−1 (△), and the NN-only part of the chiral in- signed as 3/2− [43]. Other experimental data are taken teraction NNLOsat (▽). Overall, most calculations that from refs. [20, 44, 45]. For comparison, we show the one- include 3NFs reproduce the experimental proton radii proton addition and removal energies (blue symbols) cal- well. For the binding energies, the calculations with the culated using the shell model with the YSOX interac- EM(1.8/2.0) and NNLOsat interactions reproduce both tion [46], which was constructed from a monopole-based data very well. It is important to note that without 3NFs universal interaction (VMU ). Because the phenomenolog- the calculated proton radii are about 9 – 15% (18%) ical effective two-body interactions were determined by smaller, while the ground states are overbound by as fitting experimental data, they are expected to partially much as about 24% (26%) for 14C (15C). These results include the three-nucleon effect and thus can reproduce highlight the importance of comparing both experimental relatively well the ground-state energies, drip lines, en- observables to examine the employed interactions. ergy levels, as well as the electric and spin properties of The importance of the Fujita-Miyazawa type [40] or carbon and oxygen isotopes. As shown in Fig. 5, the the chiral NNLO 3NFs [41, 42] in reproducing the binding shell model calculations reproduce the binding energies 6 of nitrogen, oxygen, fluorine, as well as boron and carbon from the measurements of the proton pick-up (d,3He) isotopes around 14C very well, but underbind boron and reaction on 12,13,14C targets [51], consistent with shell carbon isotopes with N ≥ 10 by as much as 4.5 MeV. model predictions. An attempt to study the ground- As mentioned earlier, in the absence of many-body state configurations with protons outside the 1p3/2 or- 14,15 correlations, ǫ±p resemble the proton single-particle en- bital in C has also been reported [52] very recently. ergies, and the gap between them can be taken as the To further investigate the proton subshell closure in the (sub)shell gap. In the following, we consider 14,15C and neutron-rich carbon isotopes, more experiments using the closed-shell 14,16,22O isotopes in more detail. We one-proton transfer and/or knockout reactions induced computed their ground-state binding energies and those by radioactive boron, carbon and nitrogen beams at fa- of their neighbouring isotones 13,14B, 13,15,16,21N and cilities such as ATLAS, FAIR, FRIB, RCNP, RIBF and 15,17,23F. We applied the Λ-CCSD(T) and the particle- SPIRAL2 are anticipated. attached/removed EOM-CC methods to compute the binding energies for the closed-(sub)shell and open-shell Acknowledgements We thank T. Shima, H. Toki, nuclei, respectively. The ground-state binding energies of K. Ogata and H. Horiuchi for discussion, and K. Hebeler 14B (2−) and 16N (2−) were computed using the EOM- for providing matrix elements in Jacobi coordinates for CC method with reference to 14C and 16O employing 3NFs at NNLO. H.J.O. and I.T. acknowledge the sup- the charge-exchange EOM-CC technique [47]. Results port of A. Tohsaki and his spouse. D.T.T. and T.T.N. of the CC calculations on 14,15C and 14,16,22O with and appreciate the support of RCNP Visiting Young Scientist without 3NFs are shown by the red-solid and red-dashed Support Program. This work was supported in part by lines, respectively. Here, we have opted for EM(1.8/2.0 Nishimura and Hirose International Scholarship Founda- fm−1), which yield the smallest chi-square value for the tions, the JSPS-VAST Bilateral Joint Research Project, calculated and experimental binding energies considered, Grand-in-Aid for Scientific Research Nos. 20244030, as the NN+3NF interactions. For the NN-only interac- 20740163, and 23224008 from Japan Monbukagakusho, tion, we show the calculations with EM(2.8 fm−1). The the Office of Nuclear Physics, U.S. Department of Energy, − 22 −1 calculated ǫ−p(3/2 ) for O with EM(2.8 fm ) (and under grants DE-FG02-96ER40963, de-sc0008499 (NU- other NN-only interactions) has an unrealistic positive CLEI SciDAC collaboration), the Field Work Proposal value, and is thus omitted. We found the norms of the ERKBP57 at Oak Ridge National Laboratory (ORNL), wave functions for the one-particle (1p) 1/2− and one- and the Vietnam government under the Program of De- hole (1h) 3/2− states of 14C, and the two corresponding velopment in Physics by 2020. Computer time was pro- 1p and 1h states of 15C (2− states in 14B and 16N) to be vided by the Innovative and Novel Computational Impact almost 90%. The calculations suggest that these states on Theory and Experiment (INCITE) program. This re- can be accurately interpreted by having dominant single- search used resources of RCNP Accelerator Facility and particle structure, and that the gaps between these 1p-1h the Oak Ridge Leadership Computing Facility located states resembles the ‘proton-subshell gaps’. It is obvious at ORNL, which is supported by the Office of Science from the figure that the calculations with the NN+3NF of the Department of Energy under Contract No. DE- 14,15 interactions reproduce the experimental ǫ±p for C AC05-00OR22725. and 14,16,22O very well. Overall, the calculations without 3NFs predict overbound ‘proton states’, and in the case 14,15 of C, much reduced ‘subshell gaps’. The present re- ∗ 14 15 Correspondence: [email protected] sults show that C is a double-magic nucleus, and C a [1] Goeppert Mayer, M. On closed shells in nuclei II. Phys. proton-closed shell nucleus. Rev. 75, 1969 – 1970 (1949). Finally, we would like to point out that an inert 14C [2] Haxel, O., Jensen, J.H.D. and Suess, H.E. On the “magic core, built on the N=8 closed shell, has been postulated numbers” in nuclear structure. Phys. Rev. 75, 1766 to explain several experimental data for 15,16C. For in- (1949). stance, a 14C+n model was successfully applied [48] to ex- [3] Barrett, B.R., et al. Ab initio no core shell model Prog. Part. Nucl. Phys. 69, 131 – 181 (2013). plain the consistency between the measured g-factor and [4] Carlson, J. et al. Quantum Monte Carlo methods for nu- the single-particle-model prediction (the Schmidt value) clear physics. Rev. Mod. Phys. 87, 1067 (2015). + 15 of the excited 5/2 state in C. Wiedeking et al., on the [5] Hagen, G., et al. Coupled-cluster computations of atomic other hand, have explained [17] the small B(E2) value in nuclei Rep. Prog. Phys. 77, 096302 (2014). 16C assuming a 14C+n+n model in the shell-model cal- [6] Pieper, S.C., Pandharipande, V.R. Origins of spin-orbit 15 culation. In terms of spectroscopy studies using trans- splitting in N. Phys. Rev. Lett. 70, 2541 – 2544 (1993). fer reactions, the results from the 14C(d,p)15C [49] and [7] Goeppert Mayer, M. The shell model. (1963) Nobel Lec- tures, Physics, 20 – 37 (1963). 15C(d,p)16C [50] measurements are also consistent with 14 [8] Inglis, D.R. The energy levels and the structure of light the picture of a stable C core. On the proton side, a nuclei. Rev. Mod. Phys. 25, 390 – 450 (1953). possible consolidation of the 1p3/2 proton subshell closure [9] Otsuka, T. et al. Magic numbers in exotic nuclei and when moving from 12C to 14C was reported decades ago spin-isospin properties of the NN interaction. Phys. Rev. 7
Lett. 87, 082502 (2001). Lett. 110, 192502 (2013). [10] Skaza, F. et al. Experimental evidence for subshell clo- [31] Collard, H.R., Elton, L.R.B. and Hofstadter, R. Nuclear sure in 8He and indication of a resonant state in 7He radii. Springer, Berlin Vol.2 (1967). below 1 MeV. Phys. Rev. C 73, 044301 (2006). [32] Wigner, E.P. and Feenberg, E. Symmetry properties of [11] Angeli, I. and Marinova, K.P. Table of experimental nu- nuclear levels. Rep. Prog. Phys. 8, 274 – 317 (1941). clear ground state charge radii: An update. At. Data [33] Chasman, R.R. n-p pairing, Wigner energy, and shell Nucl. Data Tables 99, 69 – 95 (2013). gaps. Phys. Rev. Lett. 99, 082501 (2007). [12] Suhara, T. and Kanada-En’yo, Y. Effects of α-cluster [34] Taube, A.G. and Barlett, R.J. Improving upon 12 breaking on 3α-cluster structures in C. Phys. Rev. C CCSD(T): ΛCCSD(T). I. Potential energy surfaces. J. 91, 024315 (2015). Chem. Phys. 128, 044110 (2008). [13] Kanada-En’yo, Y. Proton radii of Be, B, and C isotopes. [35] Gour, J.R. et al. Coupled-cluster calculations for valence Phys. Rev. C 91, 014315 (2015). systems around 16O. Phys. Rev. C 74, 024310 (2006). [14] Fujimoto, R. Shell model description of light unstable [36] Hebeler, K. et al. Improved nuclear matter calculations nuclei. Ph.D. Thesis, University of Tokyo (2003). from chiral low-momentum interactions. Phys. Rev. C 83, [15] Fujii, S. et al. Microscopic shell-model description of the 031301(R) (2011). 16 exotic nucleus C. Phys. Lett. B 650 9 – 14 (2007). [37] Entem, D.R. and Machleidt, R. Accurate charge- + [16] Forssen, C. et al. Systematics of 2 states in C isotopes dependent nucleon-nucleon potential at fourth order of from the no-core shell model. J. Phys. G: Nucl. Part. chiral perturbation theory. Phys. Rev. C 68, 041001(R) Phys. 40 055105 (2013). (2003). [17] Wiedeking, M. et al. Lifetime measurement of the first [38] Bogner, S.K., Furnstahl, R.J. and Perry, R.J. Similarity + 16 excited 2 state in C. Phys. Rev. Lett. 100, 152501 renormalization group for nucleon-nucleon interactions. (2008). Phys. Rev. C 75, 061001(R) (2007). [18] Ong, H.J. et al. Lifetime measurements of first excited [39] Hagen, G. et al. Neutron and weak-charge distributions 16,18 states in C. Phys. Rev. C 78, 014308 (2008). of the 48Ca nucleus. Nat. Phys. 12, 186 – 190 (2016). [19] Pritychenko, B., Birch, M., Singh, B. and Horoi, M. Ta- [40] Otsuka, T. et al. Three-body forces and the limit of oxy- + bles of E2 transition probabilities from the first 2 states gen isotopes. Phys. Rev. Lett. 105, 032501 (2010). in even-even nuclei. At. Data Nucl. Data Tables 107, 1 – [41] Hagen, G. et al. Continuum effects and three-nucleon 139 (2016). forces in neutron-rich oxygen isotopes. Phys. Rev. Lett. [20] Wang, M. et al. The AME2012 atomic mass evaluation. 108, 242501 (2012). Chinese Phys. C 36, 1603 – 2014 (2012). [42] Cipollone, A., Barbieri, C. and Navratil, P. Isotopic [21] Garcia Ruiz, R.F. et al. Unexpectedly large charge radii chains around oxygen from evolved chiral two- and of neutron-rich calcium isotopes. Nat. Phys. 12, 594 – three-nucleon interactions. Phys. Rev. Lett. 111, 062501 598 (2016). (2013). [22] Wienholtz, F. et al. Masses of exotic calcium isotopes pin [43] Sohler, D. et al. In-beam γ-ray spectroscopy of the down nuclear forces. Nature 498, 346 – 349 (2013). neutron-rich nitrogen isotopes 19−22N. Phys. Rev. C 77, [23] Angeli, I. and Marinova, K.P. Correlations of nuclear 044303 (2008). charge radii with other nuclear observables. J. Phys. G: [44] Ajzenberg-Selove, F. Energy levels of light nuclei A=13– Nucl. Part. Phys. 42, 055108 (2015). 15. Nucl. Phys. A 523, 1 – 196 (1991). 14 [24] Terashima, S. et al. Proton radius of Be from measure- [45] Tilley, D.R., Weller, H.R., Cheves, C.M. Energy levels of ment of charge-changing cross sections. Prog. Theor. Exp. light nuclei A=16–17. Nucl. Phys. A 565, 1 – 184 (1993). Phys. 101D02 (2014). [46] Yuan, C., Suzuki, T., Otsuka, T., Xu, F.-R. and Tsun- 12−17 [25] Estrade, A. et al. Proton radii of B define a thick oda, N. Shell-model study of boron, carbon, nitrogen and 17 neutron surface in B. Phys. Rev. Lett. 113, 132501 oxygen isotopes with a monopole-based universal inter- (2014). action. Phys. Rev. C 85, 064324 (2012). 12−19 [26] Kanungo, R. et al. Proton distribution radii of C [47] Ekstrom, A. et al. Effects of three-nucleon forces and two- illuminate features of neutron halos. Phys. Rev. Lett. 117 body currents on Gamow-Teller strengths. Phys. Rev. 102501 (2016). Lett. 113, 262504 (2014). [27] Tran, D.T. et al. Charge-changing cross-section measure- [48] Hass, M., King, H.T., Ventura, E., Murnick, D.E. Mea- 12−16 ments of C at around 45A MeV and development of surement of the magnetic moment of the first excited a Glauber model for incident energies 10A – 2100A MeV. state of 15C. Phys. Lett. B59, 32 – 34 (1975). Phys. Rev. C 94 064604 (2016). [49] Goss, J.D. et al. Angular distribution measurements for [28] Geng, L.S., Toki, H. and Meng, J. Masses, deformations 14C(d,p)15C and the level structure of 15C. Phys. Rev. C and charge radii – Nuclear ground-state properties in the 12, 1730 – 1738 (1975). relativistic mean field model. Prog. Theor. Phys. 113, 785 [50] Wuosmaa, A.H. et al. 15C(d,p)16C reaction and exotic – 800 (2005). behavior in 16C. Phys. Rev. Lett 105, 132501 (2010). [29] Ekstrom, A. et al. Acccurate nuclear radii and bind- [51] Mairle, G. and Wagner, G.J. The decrease of ground- ing energies from chiral interaction. Phys. Rev. C 91, state correlations from 12C to 14C. Nucl. Phys. A 253, 051301(R) (2015). 253 – 262 (1975). [30] Ekstrom, A. et al. Optimized chiral nucleon-nucleon in- [52] Bedoor, S. et al. Structure of 14C and 14B from the teraction at next-to-next-to-leading order. Phys. Rev. 14,15C(d,3He)13,14B reactions. Phys. Rev. C 93, 044323 (2016).