Emerging concepts in nuclear structure based on the shell model

Some emerging concepts of nuclear structure are overviewed. (1) Background: the many-body quantum structure of atomic nucleus, a complex system comprising protons and neutrons (called nucleons collectively), has been studied largely based on the idea of the quantum liquid (a la Landau), where nucleons are quasiparticles moving in a (mean) potential well, with weak"residual"interactions between nucleons. The potential is rigid in general, although it can be anisotropic. While this view was a good starting point, it is time to look into kaleidoscopic aspects of the nuclear structure brought in by underlying dynamics and nuclear forces. (2) Methods: exotic features as well as classical issues are investigated from fresh viewpoints based on the shell model and nucleon-nucleon interactions. The 70-year progress of the shell-model approach, including effective nucleon-nucleon interactions, enables us to do this. (3) Results: we go beyond the picture of the solid potential well by activating the monopole interactions of the nuclear forces. This produces notable consequences in key features such as the shell/magic structure, the shape deformation, the dripline, etc. These consequences are understood with emerging concepts such as shell evolution (incl. type-II), T-plot, self-organization (for collective bands), triaxial-shape dominance, new dripline mechanism, etc. The resulting predictions and analyses agree with experiment. (4) Conclusion: atomic nuclei are surprisingly richer objects than initially thought.


I. INTRODUCTION
The atomic nucleus is in a unique position in physics that it is an isolated object but comprises many quantum ingredients.Some emerging concepts for the structure of atomic nuclei will be overviewed in this article, focusing on the works in which the author was involved.Obviously, those concepts have been found or clarified thanks to the great progress of nuclear structure physics over 70 years, including the shell model.
In fact, the understanding of nuclear structure is based, to a great extent, on the shell model, which was introduced by Mayer [1] and Jensen [2] in 1949.Since then, the shell model has been developed significantly in many ways: an initial phase as many-body physics was presented, for instance, by Talmi in [3], in contrast to Mayer-Jensen's independentparticle model.The subsequent developments are reviewed, for instance by Caurier et al. in [4] up to 2005, and in this volume up to date.I would like to sketch emerging concepts of nuclear structure based on recent shell-model studies involving the author, as many other studies are to be presented in other articles of the same volume.
The atomic nucleus comprises Z protons and N neutrons.Their sum is called the mass number A = Z + N.Among atomic nuclei, stable nuclei are characterized by their infinite or practically infinite life times, and are characterized by rather balanced Z to N ratios, with N/Z ranging from about 1 up to about 1.5.There are about 300 nuclear species of this category.Other nuclei are called exotic (or unstable) nuclei.The total number of them is unknown, but seems to be between 7000 and 10000, providing a huge show window of various features as well as the paths of nucleosynthesis in the cosmos (see, for instance, [5][6][7]).The exotic nuclei decay, by β (i.e.weak) processes, to other nuclei where Z and N are better balanced, as the β decay alters a neutron to a proton or vice versa.This decay occurs successively, until the process terminates at a stable nucleus.Thus, only stable nuclei exist on earth, while exotic nuclei do not, being exotic literally.Some of the emerging concepts were conceived in the study of exotic nuclei, particularly by looking at the shell structure and magic numbers of them.The obtained concepts were found later not to be limited to exotic nuclei.In this way, after the initial trigger by exotic nuclei, the overall picture of the nuclear shell structure has been renewed, and Sec. 2 of this article is devoted to a sketch of it with two major keywords, the monopole interaction and the shell evolution.
We then focus on the deformation of nuclear surface.The surface deformation from the sphere has been very important subjects since 1950's, as initiated by Rainwater [8], and by Bohr and Mottelson in [9][10][11][12].In particular, the shape coexistence phenomenon is discussed as the crossroad between the shell evolution and the deformation, leading to the concept of type II shell evolution.Although I will not discuss extensively on the methodology of the shell model calculation in this article because of the length limitation, the T-plot of the Monte Carlo Shell Model will be mentioned as an essential theoretical tool for many physics cases of this article.These are main subjects of Sec. 3.
The in-depth clarification of the collective band is connected to the fundamental question on the relation between the single-particle degrees of freedom and the collective motion of nucleons.These two must be connected through nuclear forces.This question has not been clarified enough as also addressed by G. E. Brown [13].I shall focus, in Sec.IV, how this question may be understood more deeply, by intro-ducing the self-organization aspect of the collective bands and by raising the importance of the triaxiality of nuclear shapes including the ground states.
The interplay between the monopole interaction and the quadrupole deformation is shown to be a major mechanism of the determination of the neutron driplines.This approach explains neutron driplines observed recently.We are led to two dripline mechanisms: the traditional one with the singleparticle origin and the present one.The monopole-quadrupole interplay responsible for this new dripline mechanism is explained in detail in Sec.V.As an alternative case, spherical isotopes, such as Ca, Ni, Sn and Pb, are predicted to have longer isotopic chains.
The intention of this article is to show the major flow of basic ideas and related results without going into details.I hope that the reader can grasp this flow, and could become interested in watching further developments.The past 70 years are really great for the shell model, but the coming years look equally or even more brilliant.I apologize for not covering many of major developments in the last 70 years, as such coverage is not possible within this article but the other articles of this volume are expected to help.

II. SHELL EVOLUTION DUE TO MONOPOLE INTERACTION
A. Mayer-Jensen's shell model and observed magic numbers Mayer [1] and Jensen [2] proposed, in 1949, the model of the shell structure and magic numbers of atomic nuclei.This model provided major guides for deeper and wider understanding of the structure of atomic nuclei.While this is a similar situation to electrons in atoms, there are some differences.Figure 1 depicts the basic idea and consequences of the Mayer-Jensen's scheme.We start with the nuclear matter composed of protons and neutrons.This matter shows an almost constant density of nucleons (collective name of protons and neutrons) inside the surface which is a sphere as a natural assumption (see Fig. 1 (a)).Because of the short-range character of nuclear forces, this constant density results in a mean potential with a constant depth inside the surface, as shown in Fig. 1 (b).Let's assume that the density distribution is isotropic, producing an isotropic mean potential.Figure 1 (b) also suggests that the Harmonic Oscillator (HO) potential is a good approximation to this mean potential as long as the mean potential shows negative values as a function of r, the radius from the center of the nucleus.We then switch from the mean potential to the HO potential, which is analytically more tractable.Thus, the HO potential can be introduced from the constant density (sometimes referred to as "density saturation") and the short-range attraction due to nuclear forces.
The eigenstates of the HO potential are single-particle states shown in the far-left column of Fig. 1 (c) with associated magic numbers and HO quanta, N.These HO magic numbers do not change by adding the minor correction of the ℓ 2 term, the scalar product of the orbital angular momentum l (see the second column from left in Fig. 1 (c); for details see [11]).
The crucial factor introduced by Mayer and Jensen was the spin-orbit (SO) term, ( l • s), the effect of which is shown in the third column from left in Fig. 1 (c).The two orbits with the same orbital angular momentum, ℓ, and the same HO quanta are denoted as, where 1/2 is due to the spin, s = 1/2.The notation of j > and j < will be used frequently in this article.The spin-orbit term, is added to the HO + ℓ 2 potential, where f is the strength parameter.With f < 0 as is the case for nuclear forces, the j > state is lowered in energy, whereas the j < state is raised.The value of f is known empirically to be about -20A −2/3 MeV (see eq. (2-132) of [11]).
The final pattern of the single-particle energies (SPE) is shown schematically in Fig. 1 (c).The single-particle states are labelled in the standard way up to their j values, and both HO and spin-orbit magic gaps are indicated in black and red, respectively.The magic numbers have been considered to be Z, N = 2, 8, 20, 28, 50, 82 and 126, because the effect of the spin-orbit term becomes stronger as j becomes larger.In fact, the magic numbers 28, 50, 82 and 126 are all due to this effect.Instead, the HO magic numbers beyond 20 were considered to be absent or show only minor effects.We shall look back on them, from modern views of the nuclear structure covering stable and exotic nuclei.
We now investigate to what extent magic gaps in Fig. 1 (c) have been observed.Figure 2 displays the observed excitation energies of the first 2 + states of even-even nuclei as a function of N, where even-even stands for even-Z-even-N.These excitation energies tend to be high at the magic numbers, because excitations across the relevant magic gap are needed.The conventional magic numbers of Mayer and Jensen, N= 2, 8, 20, 28, ... 126 are expected to arise, and we indeed see sharp spikes at these magic numbers in Fig. 2a where the excitation energies are shown for stable and long-lived (i.e.meta stable) nuclei.The panel b includes all measured first 2 + excitation energies as of 2016.In addition to the spikes in panel a, one sees some new ones.One of them is at N=40, which corresponds to 68 Ni 40 , representing a HO magic gap at N=40.There are three others corresponding to the nuclei, 24 O 16 , 52 Ca 32 and 54 Ca 34 , as marked in red in the panel.The 2 + excitation energies of these nuclei are about a factor two higher than the overall trend, suggesting that N=16, 32 and 34 can be magic numbers, although none of them is present in Fig. 1 (c).These new possible magic numbers are consequences of what are missing in the argument for deriving magic gaps in Fig. 1 (c).We now turn to follow some passages along which this subject has been studied.The spin-orbit term, ( l • s), is included further, and magic gaps emerging from this term are shown in red.The single-particle orbits are labeled to the right, including j = l + s.The magic gaps are classified as "HO" and "SO" for the HO potential and spin-orbit origins, respectively.Taken from Fig. 2 of [14], which was adopted from [15].

B. Monopole interaction
The change from the Mayer-Jensen scheme is discussed from the viewpoint of the NN interaction.The Hamiltonian is written as, where Ĥ0 denotes the one-body term given by and V stands for the NN interaction.Here, np,n j means the proton-or neutron-number operator for the orbit j, and ǫ p,n 0; j implies proton or neutron single-particle energy (SPE) of the orbit j.This SPE is composed of the kinetic energy of the orbit j and the binding energy on the orbit j generated by all nucleons in the inert core.
The interaction V can be decomposed, in general, into the two components: monopole and multipole interactions [16], irrespectively of its origin, derivation or parameters.The monopole interaction, denoted as Vmono , is expressed in terms of the monopole matrix element, which is defined for singleparticle orbits j and j ′ as, where m and m ′ are magnetic substates of j and j ′ , respectively, and the summation over m, m ′ is taken for all ordered pairs allowed by the Pauli principle.The monopole matrix element represents, as displayed schematically in Fig. 3, an orientation average for two nucleons in the orbits j and j ′ .See [14] for more detailed descriptions.The monopole interaction between two neutrons is given as The monopole interaction between two protons is given similarly.The monopole interaction between a proton and a neutron is given as where Ṽmono pn ( j, j ′ ) stands for the monopole matrix element of the proton-neutron interaction slightly changed from the one given in eq. ( 5): changes are due to the isospin symmetry and are made only for j= j ′ (see Sec. III A of [14] for details).
The functional forms in eqs.(6,7) appear to be in accordance with the intuition from the averaging over all orientations: no dependences on angular properties (e.g.coupled J values) between the two interacting nucleons, and the sole dependence on the number of particles in those orbits.The (total) monopole interaction is written as and the monopole Hamiltonian is defined as, The multipole interaction is introduced as and the (total) Hamiltonian is written as Ĥ = Ĥmono + Vmulti .The multipole interaction becomes crucial in many aspects of nuclear structure, for instance, the shape deformation, as touched upon in later sections of this article.The monopole interaction has been studied over decades with many works, for example, [16][17][18][19] (see [14] for more details).We define the effective SPE (ESPE) of the proton (neutron) orbit j, denoted by ǫ p j (ǫ n j ), as the change of the monopole Hamiltonian, Ĥmono in eq. ( 9), due to the addition of one proton (neutron) into the orbit j.This change is nothing but the difference, when n p,n j is replaced by n p,n j +1.For instance, the first term on the right-hand-side (RHS) of eq. ( 9) contributes to ǫ p j by a constant, ǫ p 0; j .As another example, the RHS of eq. ( 7) contributes by Σ j ′ Ṽmono pn ( j, j ′ ){(n p j + 1) nn Combining all terms, the ESPE of the proton orbit j is given as, The second and third terms on the RHS are obviously contributions from valence protons and neutrons, respectively.The neutron ESPE is expressed similarly as In many practical cases, an appropriate expectation value of the ESPE operator is also called the ESPE with an implicit reference to some state characterizing the structure, e.g. the ground state.
The ESPE as an expectation value is often discussed in terms of the difference between two states, e.g.Ψ and Ψ ′ .We here show the formulas for this difference.First we introduce the symbol ∆O for an operator Ô implying the difference, and If Ψ ′ is a doubly closed shell and Ψ is an eigenstate with some valence protons and neutrons on top of this closed shell, these quantities stand for the evolution of ESPE's as functions of Z and N. One can thus see various physics cases represented by Ψ and Ψ ′ .Such ESPE's can provide picturesque prospects and great helps in intuitive understanding without resorting to complicated numerical calculations.The notion of the ESPE has been well utilized, for instance, in empirical studies in [6,21], in certain ways related to the present article.
The interaction V can be decomposed into several parts according to some classifications.The discussions in this subsection can then be applied to each part separately: the monopole interaction of a particular part of V can be extracted, and its resulting ESPE's can be evaluated.Examples are presented in the subsequent subsections.
We note that the definition of the ESPE can have certain variants with similar consequences, for instance, the combination of n p,n j − 1/2 and n p,n j + 1/2 instead of n p,n j and n p,n j + 1. Appendix A shows a note on the relation to Baranger's ESPE.

C. Monopole interaction of the central force
With these formulations, we can discuss a variety of subjects ranging from the shell structure, to the collective bands, and to the driplines.Let us start with the shell structure.While the discussions in Sec.II A are based on basic nuclear properties, some aspects are missing.One of them is the orbital dependences of the monopole matrix element.This dependence generally appears, but shows up more crucially in certain cases.We first look into the monopole interaction of the central-force component of NN interactions.Because the NN interaction keeps the character of short-range attraction after modifications or renormalizations, the monopole matrix elements gain large magnitudes with negative sign (i.e.attractive), if radial wave functions of the single-particle orbits, j number of matrix elements in the summation ✁ ✁ ✁ ✁ : magnetic substates of orbit j ✁ ✁ ✁ ✁ : magnetic substates of orbit j' Monopole matrix element between orbits j and j' 3. Schematic illustration of the monopole matrix element.Taken from Fig. 7 of [14].
wave function of relative motion spin Monopole interaction of the tensor force.Taken from [22].and j ′ in eq.( 5), are similar to each other.This similarity is visible, if these orbits are spin-orbit partners ( j = j > and j ′ = j < ) with the identical radial wave functions (see eq. ( 1)), for instance 1 f 7/2 and 1 f 5/2 .Another example is the coupling between unique-parity orbits, such as 1g 9/2 and 1h 11/2 , for which the radial wave functions are similar because of no radial node.These cases were pointed out by Federman and Pittel [20], as the total effect of the 3 S 1 channel of the NN interaction without the reference to the monopole interaction.

D. Monopole interaction of the tensor force
Another important source of the monopole interaction with strong orbital dependences is the tensor force.The tensor force produces very unique effects on the ESPE.This is shown in Fig. 4: the intuitive argument in [14,22] proves that the monopole interaction of the tensor force is attractive between a nucleon in an orbit j < and another nucleon in an orbit j ′ > , whereas it becomes repulsive for combinations, ( j > , j ′ > ) or ( j < , j ′ < ).The magnitude of such monopole interaction varies also.For example, it is strong in magnitude between spinorbit partners, or between unique-parity orbits, etc. [14].
The ESPE is shifted in very specific ways as exemplified in Fig. 5 (b): if neutrons occupy a j ′ > orbit, the ESPE of the proton orbit j > is raised, whereas that of the proton orbit j < is lowered.This is nothing but a reduction of a proton spin-orbit splitting due to a specific neutron configuration.The amount of the shift is proportional to the number of neutrons in this configuration, as shown in eq. ( 13) and in Fig. 5(c).Other cases follow the same rule shown in Fig. 4.These general features have been pointed out in [22] with an analytic formula and an intuitive description of its origin.

E. Monopole-interaction effects from the central and tensor forces combined
The combined effects of the central and tensor forces were discussed in [23] in terms of realistic shell-model interactions, USD [24] and GXPF1A [25], referring to their roots in microscopic G-matrix NN interactions proposed initially by Kuo and Brown [26,27].Many other valuable shell-model interactions, for instance, (so-called) KB3 [16], Kuo-Herling [28], sn100pn [29] and LNPS [30] interactions, have been constructed from the G-matrix interactions sometimes with refinements like monopole adjustments.The M3Y interaction [31] is related to the G-matrix, too.We appreciate the original contribution of the so-called G-matrix approach to the effective NN interaction [26,27].
The V MU interaction was then introduced as a general and simple shell-model NN interaction.Its central part consists of Gaussian interactions with spin/isospin dependences, and their strength parameters are determined so as to simulate the overall features of the monopole matrix elements of the central part of USD [24] and GXPF1A [25] interactions.Its tensor part is taken from the standard πand ρ-meson exchange potentials [22,32,33].Thus, the V MU interaction is defined as a function of the relative distance of two nucleons with spin/isospin dependences, which enables us to use it in a variety of regions of the nuclear chart, as we shall see.A wide model space, typically a HO shell or more, is required in order to obtain reasonable results, though.
Figure 6 depicts some examples: panel (a) displays the transition from a standard (a la Mayer-Jensen) N=20 magic gap to an exotic N=16 magic gap by plotting ǫn j within the filling scheme (see eq. ( 12)), as Z decreases from 20 to 8. The tensor monopole interaction between the proton d 5/2 and the neutron d 3/2 orbits plays an important role.The small N=20 magic gap for Z=8-12 is consistent with the island of inversion picture (see reviews, e.g.[4,14]).Panel (b) depicts the inversion between the proton f 5/2 and p 3/2 orbits as N increases in Ni  isotopes, by showing ǫ p j (see eq. ( 11)).The figure exhibits exotically ordered single-particle orbits for N > 44.The tensor monopole interactions between the proton f 7/2,5/2 and the neutron g 9/2 orbits produce crucial effects.Panel (c) shows significant changes of the neutron single-particle levels from 90 Zr to 100 Sn, in terms of ǫn j .Without the tensor force, the degeneracy of g 7/2 and d 5/2 orbits in 100 Sn does not show up.
These changes of the shell structure as a function of Z and/or N were collectively called shell evolution in [22].The splitting between proton g 7/2 and h 11/2 in Sb isotopes shows a substantial widening as N increases from 64 to 82 as pointed out by Schiffer et al. [43], which was one of the first experimental supports to the shell evolution partly because this was not explained otherwise.Note that while the origin of the shell evolution can be any part of the NN interaction, its appearance is exemplified graphically in Fig. 5(a,b,c) for the tensor force.The shell-evolution trend depicted in Fig. 6 ap-pear to be consistent with experiment [14,23,[44][45][46][47][48].The monopole properties discussed in this subsection are consistent with the results shown in [49] obtained by the spin-tensor decomposition [14].
F. N=34 new magic number as a consequence of the shell evolution Among various cases of shell evolution, a notable impact was made by predicting a new magic number N=34. Figure 7 displays the shell evolution of some neutron orbits from Ni back to Ca isotopes, as Z decreases from 28 to 20.The 1 f 5/2 orbit is between the 2p 3/2 and 2p 1/2 in Mayer-Jensen's shell model (see Fig. 1).By loosing eight protons in Ni isotopes, this canonical shell structure is destroyed as the 1 f 5/2 orbit moves up above the 2p Proton number, Z creates the N=32 gap as a byproduct.The energy shift of the 1 f 5/2 orbit is due to the central and tensor forces by almost equal amounts: the N=34 magic gap does not appear without this shift, if the Mayer-Jensen scheme holds in Ni isotopes.

G. Monopole interaction of the 2-body spin-orbit force
It is a natural question what effect can be expected from the 2-body spin-orbit force of the NN interaction.This force can be well described by the M3Y interaction, and its monopole effects were described in detail in [14], particularly in its supplementary document.Although the monopole effects of this force contributes to the spin-orbit splitting [14], the effect is much weaker than the tensor force in most cases, as also discussed in the article by Utsuno in this volume.
An interesting case is found in the coupling between an s orbits and the p 3/2,1/2 orbits.There is no monopole effect from the tensor force, if an s orbit is involved.Instead, the s-p coupling due to 2-body spin-orbit force can be exceptionally strong as intuitively stressed in [14].Figure 8 shows that the possible significant change of the neutron 2p 3/2 -2p 1/2 gap between 35 Si and 37 S is explained to a good extent by the shell evolution due to the 2-body spin-orbit force.

H. Monopole interaction from the three-nucleon force
The three-nucleon force (3NF) is currently of intense interest.Among various aspects, we showed [34] the characteristic feature of the monopole interaction of the effective NN interaction derived from the Fujita-Miyazawa 3NF. Figure 9 (a  39 Ar [50] and 41 Ca [51] and [50].The horizontal bars are the energy differences between relevant highest peaks [52], [50].Shell evolution predictions are shown by blue closed symbols and solid line connecting them.The loose binding effect for 35 Si is included in the open circle.The calculation with Woods-Saxon potential with parameters adjusted are shown by the yellowish shaded area [50].(b) Neutron 2p 1/2 single-particle energy (blue solid line) by a Woods-Saxon potential [11] for varying depth parameter, V.The linear dependence of the deeply bound region is linearly extrapolated (blue dashed line) and is compared to the curved dependence that results from the proximity of the continuum.The dashed line is for the 2p 3/2 orbit, and the loose-binding contribution to the present splitting appears to be 0.06 MeV against 1.5 MeV splitting itself.Taken from Fig. 8 of Supplementary Material of [14].In this argument, the 3NF produces a repulsive monopole NN interaction in the valence space, after the summation over the hole states of the inert core (see the right lower panel, labelled (d)).The right upper panel indicates an example of the repulsive effect on the ground-state energy of oxygen isotopes, locating the oxygen dripline at the right place or solving the oxygen anomaly [34].This is rather strong repulsive monopole interaction, which is a consequence of inert core.This means that the present case is irrelevant to the no-core shell model or other many-body approaches without the inert core (e.g.GFMC) [35].This feature has caused some confusions in the past, but the difference is clear.This repulsive effect is much stronger than the other effects of the 3NF [36].I note that this repulsive effect was stressed by Talmi in 1960's [3].

I. Short summary of Sec. II
The shell evolution phenomena are seen in many isotopic and isotonic chains, and sometimes result in the formation of new magic gaps or the vanishing of old ones.Figure 2(b) displays the emergence of such new magic numbers N=16, 32 and 34, whereas the lowering of some 2 + levels can mean the weakening of some magic numbers.More changes may appear in the future studies.Thus, the characteristic monopole features of the central, tensor, 2body-LS and 3NF-based NN interactions and the resulting shell evolution are among the emerging concepts of the nuclear structure.Interestingly, these findings are neither isolated nor limited to particular aspects, but are related to other aspects of the nuclear structure.We now move on to such a case.

III. TYPE-II SHELL EVOLUTION AND SHAPE COEXISTENCE A. Type-II Shell Evolution
The shell evolution shown in Fig. 5(b) and (c) are, respectively, due to the addition of two or four neutrons into the orbit j ′ > .Instead of adding, one can put neutrons into the orbit j ′ > by taking the neutrons from some orbits below j ′ > , or equivalently by creating holes there, as shown in panel (d).If such lower orbit happens to be the j ′′ < orbit as in panel (d), its monopole matrix elements show just the opposite trends compared the j ′ > orbit.However, because holes are created in j ′′ < , the sign of the monopole-interaction effect is reversed, and the final effect FIG. 10.Level scheme of 68 Ni.Taken from Fig. 2 of [53].
has the same sign as the monopole effect form the orbit j ′ > (see panel (d)).Thus, the particle-hole (ph) excitation of the two neutrons in panel (d) reduces the proton j >j < splitting even more than in panel (b).This reduction becomes stronger with the ph excitations of four neutrons, as depicted in panel (d).Such strong reduction of the spin-orbit splitting produces interesting consequences beyond shell-structure changes.This type of the shell structure change within the same nucleus is called Type II Shell Evolution.

B. A doubly-closed nucleus 68 Ni
The Type-II shell evolution was first discussed in [53] for 68 Ni as an example.Figure 10 shows its theoretical and experimental energy levels.The theoretical results were obtained for the A3DA-m interaction by the Monte Carlo Shell Model, which is a powerful methodology for the shell model calculation but is not discussed here due to the length limitation [54][55][56][57].
Because Z=28 and N=40 of 68 Ni are both magic (see Fig. 1), the ground state is primarily a doubly closed shell.Indeed, in the theoretical ground state, the occupation of the neutron g 9/2 orbit is negligibly small.In contrast, the 0 + 3 state located at the excitation energy, Ex∼ 3 MeV, is the band head of a rotational band of an ellipsoidal shape, and its neutron g 9/2 occupation number is as large as ∼4.The mechanism shown in Fig. 5(e) is then switched on, reducing the proton f 5/2f 7/2 splitting.A reduced splitting facilitates more configuration mixing between these two orbits, which can produce notable effects on the quadrupole deformation as stated below.

C. Coexistence between spherical and deformed shapes
We here quickly overview the quadrupole deformation, or the shape deformation from a sphere to an ellipsoid [12].The quadrupole deformation is driven by the quadrupole interaction, a part of the multipole interaction in eq.(10).The quadrupole interaction is a somewhat vague idea because of a certain mathematical complication, but its main effects can be simulated by the (scalar) coupling of the quadrupole moment operators.If the quadrupole moments are larger, i.e., a stronger quadrupole deformation occurs, the nucleus gains more binding energy from the quadrupole interaction.This is a very general phenomena, and because of this the ground states of many nuclei are deformed, although 68 Ni is not among them.
The energy of 68 Ni (intrinsic state) is graphically illustrated in Fig. 11 (left) for various ellipsoidal shapes, spherical, prolate, oblate and in between (called triaxial).The energy is obtained by the constraint Hartree-Fock (CHF) calculation with the same shell-model Hamiltonian as in Fig. 10.The imposed constraints are given by the quadrupole moments in the intrinsic (body-fixed) frame.This plot is usually called the Potential Energy Surface (PES).The minimum energy occurs at the spherical shape (red sphere).The constraints are changed to a more prolate deformed ellipsoid (blue object) along the upper-right axis ("prolate deformation" in the figure).The energy relative to the minimum energy climbs up by 6 MeV first.This is because protons and neutrons must be excited across the magic gaps from the doubly closed shell in order to create states of deformed shapes (see Fig. 1).The energy then start to come down, as the quadrupole moments increase, thanks to the quadrupole interaction.It is lowered by 3 MeV from the local peak to the local minimum.Beyond the localminimum area, the effect of the quadrupole interaction is saturated, and it cannot compete the energy needed for exciting more protons and neutrons across the gaps required by the constraints.This energy variation appears as the basin in the 3-dimensional PES.This is the usual explanation of the local deformed minimum.The appearance of two (or more) different shapes with rather small energy difference is one of the phenomena frequently seen, and is called the shape coexistence [60].The quadrupole interaction is undoubtedly among the essential factors of the shape coexistence.But, this may not be a full story.
Figure 11 (right) exhibits the same energy along the axis lines of the left panel.The positive (intrinsic) quadrupole moments imply prolate shapes (blue object in the left panel), whereas the negative ones oblate shapes (green object).The red solid line shows the CHF results of the full Hamiltonian, whereas for the dashed line, the tensor monopole interactions between the neutron g 9/2 orbit and the proton f 5/2,7/2 orbits are practically removed.This removal means no effects depicted in Fig. 5 (d,e).The dashed line displays a less pronounced prolate local minimum at weaker deformation with much higher excitation energy.The significant difference between the solid and dashed lines suggests that the monopole effects are crucial to lower this local minimum and stabilize it.We now discuss the mechanism for this difference.With the tensor monopole interaction, once sufficient neutrons are in g 9/2 , the proton f 5/2f 7/2 splitting is reduced, and this reduced splitting facilitates the mixing between these two orbits driven by the quadrupole interaction.The resulting deformation is stronger compared to no tensor-force case.The tensor monopole interaction involving the neutron g 9/2 orbit produces extra binding energy, if more protons are in f 5/2 and less are in f 7/2 .This extra binding energy lowers the deformed states, otherwise they are high in energy because of the energy cost for promoting neutrons from the p f shell to g 9/2 .Thus, a strong interplay emerges between the monopole interaction and the quadrupole interaction, and Type II shell evolution materializes this interplay in the present case.It enhances the deformation and lowers the energy of deformed states.Without this interplay, as indicated by blue dashed line in Fig. 11 (right), the rotational band corresponding the local minimum is pushed up by 2 MeV, and may be dissolved into the sea of many other states.It is obvious that this interplay mechanism works self-consistently.

D. Tplot analysis
The T-plot was introduced in the same reference [53], so that what shapes are more relevant to individual eigenstates of the shell-model calculation.Let us take an example.Figure 12 [59] depicts the PES of 66 Ni with the same Hamiltonian as in Fig. 10.The small circles on the PES are the T-plot.The Tplot is obtained from MCSM eigenstate.We therefore briefly explain the MCSM eigenstate.An MCSM eigenstate, Ψ, is written, with the ortho-normalization, as where f k denotes amplitude, PJ π means the projection operator on to the spin/parity J π (this part is more complicated in practice), and φ k stands for Slater determinant called (kth) MCSM basis vector: Here, |0 is the inert core (closed shell), c (k) † i refers to a superposition of usual single-particle states, for instance, those of the HO potential, with a † n being the creation operator of a usual single-particle state, for instance, that of the HO potential, and D (k) i,n denoting a matrix element.By choosing a favorable matrix D (k) , we can select φ k so that such φ k better contributes to the lowering of the corresponding energy eigenvalue.Thus, the determination of D (k) is the core of the MCSM calculation.The index k runs up to 50-100, but sometimes to 300 at maximum.These are much smaller than the dimension of the many-body Hilbert space.
Each φ k has intrinsic quadrupole moments ( φ k | Q0 |φ k and φ k | Q2 |φ k ), and its T-plot circle is placed according to those values on the PES with its area proportional to the overlap probability with the corresponding eigenstate, i.e.Ψ in eq. ( 15).Such T-plot circles are shown in Fig. 12.The white circles represent the MCSM basis vectors for the ground state, while the red circles indicate the MCSM basis vectors for the 0 + 4 state, which is strongly deformed.Although there is no local minimum for oblate shape, the 0 + 2 state is shown to be moderately oblate deformed.The T-plot can thus give partial labelling to fully correlated eigenstates for mean values as well as fluctuations with respect to their quadrupole shapes.The advantages of mean-field approaches are now nicely incorporated into the shell model.

E. Short summary of Sec. III
Type II shell evolution occurs in various cases, especially in a number of shape coexistence cases, providing deformed states with stronger deformation, lower excitation energies and more stabilities.It is an appearance of the monopolequadrupole interplay, and plays crucial roles in various phenomena including the first order phase transition (Zr isotopes [61][62][63]), the second order phase transition (Sn isotopes [64]), the multiple even-odd phase transitions (Hg isotopes [65]) as well as the raising of the intruder band due to the suppression of the type II shell evolution (lighter Ni isotopes [59,66]).As the involvement of the monopole interaction in this manner had not been recognized, Type II shell evolution appears to be among the emerging concepts of nuclear structure.The Type II shell evolution has been clarified by the T-plot in many cases.Including other contributions, the T-plot is undoubtedly one of the emerging concepts of nuclear structure, apart from its impact on the computational methodology.

IV. SELF-ORGANIZATION AND COLLECTIVE BANDS IN HEAVY NUCLEI
We now proceed to more general cases of the monopolequadrupole interplay.This interplay leads to unexpected consequences in the under-lying mechanism of collective bands of heavy nuclei [67], beyond the standard textbooks.
The MCSM has become powerful enough [57] to reproduce collective bands of heavy nuclei such as 154 Sm and 166 Er, with one and half HO major shells [67].We sketch the new findings by using the results of such most advanced MCSM calculations.

A. Shape coexistence in 154 Sm
Figure 13 shows low-lying energy levels of 154 Sm.The present MCSM calculation can describe the four low-lying bands including negative-parity one.The agreement between the theoretical levels in panel a and the experimental levels in panel b is rather good.Although the importance of the quadrupole interaction is evident for the formation of deformed rotational bands, one can investigate to what extent the monopole interaction is involved.The monopole interaction here was obtained from the shell-model interactions, comprising the central, tensor and other components.
The monopole interaction is an operator, but we "freeze" it now: its ESPE expectation values ǫ p,n j are calculated for the state to be specified, and the obtained values are adopted as the SPEs, ǫ p,n 0; j , with the monopole interaction removed.We then perform the shell model calculation and draw the PES.This toy game is called the "monopole frozen" analysis [67], as the monopole properties are included only through the specified state.Panel c of Fig. 13 exhibits the energy levels obtained by the monopole-frozen analysis referring to the ground state.The band built on the 0 + 2 state (often called the β band) is lifted up by 0.5 MeV (∼50% of the original excitation energy), suggesting that the active monopole interaction produces a substantial lowering of this state.Panel d shows the monopolefrozen analysis referring to the spherical HF state: the ground state is no longer prolate, but triaxial, with the wave function close to the 0 + 2 state of the original Hamiltonian.Thus, the crucial effect of the monopole interaction is verified.
Right panel of Fig. 13 shows the actual values of ǫ p,n j for the 0 + 1 and 0 + 2 states.This figure demonstrates the significant differences between two sets of the ESPE values.The occupation numbers are also different: there are more half-filled orbits for the 0 + 2 state, which is indicative of its triaxial nature.The smaller occupation numbers of unique-parity orbits are also consistent with the tendency away from the prolate shape.calc.original  Panels c and d depict a valley of the PES with a local minimum around γ = 15 • .Similar valleys are seen in the PES obtained by the mean-field calculations [68,69], implying that this valley likely has a common origin.On the other hand, one can state that the present monopole effect results in the local minimum, not simply the valley.It is of interest to refine the monopole interaction in mean-field models.
Regarding the so-called β vibration picture of the the 0 + 2 state, the present view is opposed to such a conventional view.The triaxial deformation is shared by the members not only of the 0 + 2 band but also of the 2 + 3 band (usually called γ band), as can be verified by their T-plots.Namely, the 0 + 2 state is the "ground" state of the triaxial states to which both the 0 + 2 and 2 + 3 bands belong.In short, this is a shape coexistence between the prolate and triaxial shapes assisted by the interplay between the monopole interaction and the quadrupole deformation.

B. Collective bands and γ vibration in 166 Er
The features of the collective motion in 166 Er have been studied by the MCSM similarly well (see Fig. 15 a).Among so-called rotational nuclei, 166 Er is characterized by particularly low-lying 2 + 2 state and the so-called γ band built on it.Aage Bohr stressed that this 2 + 2 state was a γ vibration from the prolate ground state [9,10,12].The relatively strong 2 + 2 → 0 + 1 E2 transition (B(E2)∼ 5 e fm 2 , see panel a) was ascribed to the annihilation of one γ phonon in the 2 + 2 state.This was one of the major points of the Nobel lecture by Aage Bohr, and has been a common sense as stated in many textbooks of nuclear physics.We now challenge this traditional belief, by utilizing the recent MCSM calculation.It is reminded that no firm experimental evidence to uniquely pin down the γvibration nature of 166 Er has been reported, and also that in a systematic calculation of many heavy nuclei [70], the excitation energies of the 2 + 2 states in the γ band appeared to be about twice higher than the observed values, despite much better description of those of the 2 + 1 state in the ground band.
Figure 15 b shows the calculated PES, which shows the minimum not at γ = 0 • but around γ = 9 • (see also [71]).The T-plot is shown for the 0 + 1 and 2 + 2 states in panels c and d, respectively.The patterns of the T-plot circles are nearly identical between these two panels.This is consistent with a (rigid) triaxial interpretation, and indeed E2 transition strengths follow the predictions of Davidov triaxial model [72,73] with γ = 9 • .Certainly, a pure rigid triaxiality is not the correct picture, and there are quantum fluctuations, as evident from panels c and d [71].After all, the displacement from the γ = 0 • is ovbious.Note that the so-called K π =4 + double γphonon states are also described well by the present work [71], whereas its reproduction within the phonon model is known to be very difficult.The triaxiality of 166 Er is also suggested by the triaxial projected shell model, although the rigid-triaxiality is not an outcome but an assumption [74,75].
The monopole-frozen analysis referring to the spherical CHF state shows that the ground state moves to γ = 0 • , confirming the important role of the monopole interaction activated.The triaxial ground states are now shown to appear in a large number of nuclei in the nuclear chart, besides the known triaxial domain [76].The collective bands in heavy nuclei have traditionally been understood in terms of the ground band with axially symmetric prolate shape and the side bands with the β or γ vibrational excitations from the ground state.This picture is consistent with the Nilsson model [78], and was confirmed by the Pairing + Quadrupole-Quadrupole (P+QQ) model [79,80], where the monopole interaction is not included, however.It has been shown in this section that the monopole interaction is crucial also for the collective bands in heavy nuclei.We just note that in lighter nuclei, the situation can be different mainly because of small model spaces comprising single or a few active orbits, where the rotational motion has been nicely described by symmetry-based approaches, e.g.SU(3) model of Elliott for the sd shell [81,82], and by realistic calculations, e.g. on 48 Cr [83].
Regarding heavy nuclei, for individual rotational bands, the monopole interaction contributes differently, and the intrinsic structure is determined not only by the quadrupole interaction but also by the monopole interaction, as verified by the monopole-frozen analyses.Thus, the monopole-quadrupole interplay arises.The monopole interaction does not directly drive the deformation, but optimizes the ESPEs so that more binding energy is gained.This gain is state-dependent, and even can alter the ordering of bands as mentioned above.The present monopole-quadrupole interplay can be described also from the viewpoint of the self-organization [67]: the nucleus is changed from a disorder (original SPEs) to an order (ESPEs tailored to the shape of interest) by activating the monopole in-teraction.As this occurs "purposely" towards certain shapes with positive feedback particularly between the monopole and quadrupole effects, the whole picture fits well the (quantal) self-organization [67].The self-organization for collective bands is among the emerging concepts of nuclear structure, showing novel consequences.For example, the dominant fraction of the ground states of heavy nuclei are expected to show triaxial shapes, as another emerging concept of nuclear structure, in contrast to the traditional view of the prolate shape dominance in those states.
Appendix B presents a possible extension or generalization of the current idea to "many-ingredient" systems outside nuclear physics.An isotopic chain is a horizontal belt, and its neutron-rich end is called neutron dripline.The location of the dripline in the nuclear chart implies the extent of the isotopes, and is of fundamental importance to nuclear science.The experimental determination of the dripline is a very difficult task.Very recently, as shown by red squares in Fig. 16 a, the driplines of F and Ne isotopes and its candidate of Na isotope were reported [84].The traditional view of the dripline is shown in Fig. 16 b: all bound single-particle orbits are occupied, and the next neutron goes away.It is an open question whether this view is valid for all nuclei or not.We look into this question now [36].

V. DRIPLINE MECHANISM A. Traditional view
The structure of neutron-rich exotic isotopes of F, Ne, Na and Mg can be well described by the shell model calculation with the full sd+p f shells and the EEdf1 interaction [85].This interaction was derived from the chiral EFT interaction of Machleidt and Entem [86], first processed by the V low−k method [87,88], and then processed by the EKK (Extended  17.Ground-state energies of even-N isotopes of F, Ne, Na and Mg, relative to the 16 O value.Colored segments exhibit decompositions into various effects from the monopole (green), pairing (blue) and rest (such as quadrupole) (red) components of the effective nucleon-nucleon interaction as well as those from Coulomb interaction (black) and single-particle energies (bare SPE; grey).The monopole effect grows steadily as a function of N in all cases, as highlighted by straight lines: dashed (Ne), dotted (Na) and solid (Mg).The experimental values are indicated by black circles [48].The theoretical driplines indicated by red arrows.Modified from Fig. 4 of [36].Krenciglowa-Kuo) method [89][90][91].The V low−k method is used to transform the nuclear forces in the free space into a tractable form for further treatments.The V low−k method has been adopted for the derivation of other modern shell-model interactions, for instance, the one by Corraggio et al., for Sn and Cr-Fe regions [92,93].
The present work is unique in the usage of the EKK method, which is an extension of the many-body perturbation theory (MBPT) [27].The MBPT produced the G-matrix interactions in its early formulations [26], from which many useful shell-model interactions have been constructed (see subsect.II E).However, the G-matrix formulation has a limitation that if two major shells are merged, the results may diverge [91].As the gap between two shells often vanishes or becomes smaller in exotic nuclei, this difficulty can be fatal there.The EKK method nicely avoids this difficulty besides other merits.
The EEdf1 interaction has thus been derived in an ab initio way by the V low−k and EKK methods, and describes the prop-erties of the ground and low-lying states of F, Ne, Na and Mg isotopes quite well [36,85].

B. Monopole-quadrupole interplay for the driplines
Figure 17 shows the ground-state energies of F, Ne, Na and Mg isotopes as functions of the neutron number N. These energies are decomposed into several pieces according to their origins: SPE (on top of the 16 O inert core), monopole, pairing and rest terms.The Coulomb contribution is ignored in the following discussion, because of virtually no relevance.Here, the multipole interaction is divided into the pairing and rest terms.The pairing is the BCS-type pairing interaction acting on two neutrons coupled to J π = 0 + and on two protons coupled to J π = 0 + .The rest term means the multipole interaction subtracted by the pairing term.Although the rest term contains many different pieces, its major effects in the present discussion is simulated by the quadrupole interaction.This is the reason why the rest term is associated with "(quadrupole etc)" in the figure.
The lower edges of the red areas exhibit the ground-state energies as functions of the neutron number N, while only even N values are taken.These values show a good agreement with measured values shown by black dots.As long as the ground-state energy becomes lower as N increases, the isotope gains more binding energy by having more neutrons, and the isotope chain is stretched.However, if the ground-state energy is not lowered, there is no gain of the binding energy by having these extra neutrons; these extra neutrons are emitted, and the neutron dripline implies the nucleus with the lowest groundstate energy.The driplines obtained by the present calculation are shown by red arrows for each isotopic chain, reproducing experimental driplines for F, Ne and Na isotopes [84].
We focus on the lower edge of the green areas in Fig. 17.This represents the monopole contributions comprising the SPE and the monopole interaction.For Ne, Na and Mg isotopes, this edge is lowered almost linearly as N increases from N=16 to each dripline.We then fit the edge with pink dashed, purple dotted and black solid lines for F, Ne and Mg isotopes, respectively.The lines of Ne and Na isotopes are copied to the panel for Mg, with their positions adjusted at a certain N. It is evident that the lines become steeper almost linearly as Z increases.This edge is almost flat for F isotopes for N ≥ 16, and this feature will be discussed later.
Figure 17 indicates that the effect of the pairing term shows small variations.In contrast, the rest terms changes more, and are largely due to the quadrupole interaction.Figure 18(a) schematically indicates the variation of the effect of the quadrupole interaction: The effect is small at the far left position with a spherical shape.As some neutrons are added, the shape is deformed, and the ground-state energy is lowered due to the quadrupole interaction.This trend continues, but becomes its maximum at a certain value of N (red object in the figure).However, the dripline is not determined just by this maximum point.18(a), with the maximum effect at N=22 in all four chains.However, the driplines are different among these four.This is due to the monopole interaction.Let me explain it, by taking the Mg isotopes as an example.The black straight line of the monopole effect in Fig. 17(d) depicts about 3 MeV lowering per additional neutron, implying about 6 MeV for additional two neutrons.After N=22, the rest effect loses its magnitude.If the loss is less than the monopole gain (∼6 MeV), this loss is compensated by the monopole effect.However, the loss becomes larger for N larger, and at a certain point, the loss exceeds the monopole compensation.The dripline thus arises with the "monopole displacement" from N=22 to N=30 as shown in Fig. 18(b) (and also in Fig. 18(a) schematically).
The monopole effect depends directly on the number of protons, as visualized by three straight lines Fig. 17.Consequently, the monopole displacement is ∆N=2 (6) for Ne (Na) isotopes.For F isotopes, the monopole effect is negligibly small for N ≥16, and the dripline is located at the maximum rest (quadrupole etc) effect.

C. Stability of spherical isotopes and the monopole-quadrupole interplay
An immediate lemma of the present dripline mechanism is that the driplines of spherical nuclei, such as Ca, Sn and Pb isotopes, may be further away from the stability line than other elements.One can assume a basically constant pairing contribution and a minor rest-term contribution.These two are thus irrelevant to the driplines of these isotopes.The remaining monopole effect gradually changes, pushing the driplines away.The present new dripline mechanism [36] involves the monopole-quadrupole interplay, and is one of the emerging concepts.It definitely differs from the traditional mechanism of the single-particle origin, where neutron halo arises at extremes [94,95].In the new mechanism, the coupling to continuum may be visible if the monopole effect vanishes like heavy F isotopes [36].As Z changes, two dripline mechanisms may appear alternatively, but the present one may be more relevant to heavier nuclei where the deformation develops more.

VI. PROSPECT
As this article is a kind of summary, I am afraid that a summary section may be redundant.I state some prospects.First of all, ab initio no-core Monte Carlo Shell Model calculations became feasible recently up to 12 C and beyond [96], and as an example, we can look into α clustering in light nuclei, e.g. the Hoyle state, with correlations produced by nuclear forces [97].This direction will produce a major outcome from the shell model.This includes clarifications of α decay, α knockout, etc.Another major frontier is the quest for fission dynamics and superheavy elements, with (almost) full inclusion of the correlations due to nuclear forces.
Although more computer powers and further advancements in computational methodology are needed also, the perspectives of the shell model look unlimited, to me.May the (nuclear) force be with you.

FIG. 1 .
FIG. 1. Schematic illustration of (a) density distribution of nucleons in atomic nuclei, (b) a mean potential (solid line) produced by nucleons in atomic nuclei and an approximation by a Harmonic Oscillator (HO) potential (dashed line).In (a) and (b), the horizontal line denotes the radius from the center of the nucleus.(c) The shell structure produced with resulting magic numbers in circles.(left column) Only the HO potential is taken with HO quanta shown as N=0, N=1, ... (N here does not mean the neutron number, N.) (middle column) The so-called ℓ 2 term is aded to the HO potential, where the magic gaps are shown in circles.The single-particle orbits are labeled in the standard way to the left.(right column)The spin-orbit term, ( l • s), is included further, and magic gaps emerging from this term are shown in red.The single-particle orbits are labeled to the right, including j = l + s.The magic gaps are classified as "HO" and "SO" for the HO potential and spin-orbit origins, respectively.Taken from Fig.2of[14], which was adopted from[15].

FIG. 5 .
FIG. 5. Schematic picture of the ESPE change (i.e.shell evolution) due to the monopole interaction of the tensor force.(a) SPE's with no neutrons in in the orbit j ′ > .(b) The shifts of the proton ESPE's due to two (valence) neutrons in the orbit j ′ > .(c) Same as (b) except for four neutrons.(d,e)Type-II shell evolution due to neutron particle-hole excitations.Taken from[37].

FIG. 6 .
FIG.6.Effective single-particle energies calculated by the VMU interaction.The dashed lines are obtained by the central force only, while the solid lines include both the central-force and the tensor-force contributions.Taken from[23].

FIG. 7 .
FIG. 7. (left) Schematic illustration of the shell evolution from Ni back to Ca for neutron orbits.Blue circles denote protons.The wavy line is the interaction between the proton 1 f 7/2 orbit and the neutron 1 f 5/2 orbit.The numbers in circles indicate magic numbers.(right) Observed excitation energies of the 2 + 1 states.Taken from Fig. 3 of [37] and from Fig. 2 c of [40].

FIG. 8 .
FIG. 8. (a) Neutron 2p 3/2 -2p 1/2 splitting for N = 21 isotones.The symbols are the centroids for 37 S [48],39 Ar [50] and41 Ca[51] and[50].The horizontal bars are the energy differences between relevant highest peaks[52],[50].Shell evolution predictions are shown by blue closed symbols and solid line connecting them.The loose binding effect for35 Si is included in the open circle.The calculation with Woods-Saxon potential with parameters adjusted are shown by the yellowish shaded area[50].(b) Neutron 2p 1/2 single-particle energy (blue solid line) by a Woods-Saxon potential[11] for varying depth parameter, V.The linear dependence of the deeply bound region is linearly extrapolated (blue dashed line) and is compared to the curved dependence that results from the proximity of the continuum.The dashed line is for the 2p 3/2 orbit, and the loose-binding contribution to the present splitting appears to be 0.06 MeV against 1.5 MeV splitting itself.Taken from Fig.8of Supplementary Material of[14].
FIG. 11. (left) Potential energy surface (PES) of68 Ni.Taken from Fig.5of[58].(right) PES of68 Ni with axially symmetric shapes.The solid line shows the PES of the full Hamiltonian, whereas the dashed line is the PES with practically no tensor-force contribution.Taken from Fig.6of[58].

Figure 14 c
Figure 14 c and d show the T-plot for the original interaction, while Fig. 14 e and f depict the T-plot for the the monopole-frozen interaction obtained with the spherical HF state.The T-plot patterns are consistent with the above features suggested by the shell-model diagonalization.The cut of the PES shown in panel b suggests that the local minimum is raised by the monopole-frozen process with the ground state.Panels c and d depict a valley of the PES with a local minimum around γ = 15 • .Similar valleys are seen in the PES obtained by the mean-field calculations[68,69], implying that this valley likely has a common origin.On the other hand, one can state that the present monopole effect results in the local minimum, not simply the valley.It is of interest to refine the monopole interaction in mean-field models.Regarding the so-called β vibration picture of the the 0 +

6 FIG. 14 .
FIG. 14. Properties of the 0 + 1,2 states of 154 Sm.(a) Deformation parameters and shapes.(b) Lowest values of PES for a given γ value for the original, prolate and spherical monopole-frozen cases.(c-f) Three-dimensional T-plot in the original or spherical monopole-frozen case.Based on Fig. 3 of [67].

Figure 16 a
Figure 16 a shows the lower-left part of the nuclear chart (Segrè chart) for Z ≤ 16.The black squares represent stable nuclei, while the orange ones exotic nuclei (see Sec. I).An isotopic chain is a horizontal belt, and its neutron-rich end is called neutron dripline.The location of the dripline in the nuclear chart implies the extent of the isotopes, and is of fundamental importance to nuclear science.The experimental determination of the dripline is a very difficult task.Very recently, as shown by red squares in Fig.16a, the driplines of F and Ne isotopes and its candidate of Na isotope were reported[84].

FIG. 18 .
FIG. 18. a Presently proposed mechanism based on shape evolution and the resulting change of the binding energy.b The rest-term contribution to the ground-state energies for F, Ne, Na and Mg isotopes.Dashed arrows indicate the monopole displacement.Modified from Figs. 2 and 6 of [36].

Figure 18 (
Figure18(b) depicts the actual effect of the rest term.It follows the trend illustrated in Fig.18(a), with the maximum effect at N=22 in all four chains.However, the driplines are different among these four.This is due to the monopole interaction.Let me explain it, by taking the Mg isotopes as an example.The black straight line of the monopole effect in Fig.17(d) depicts about 3 MeV lowering per additional neutron, implying about 6 MeV for additional two neutrons.After N=22, the rest effect loses its magnitude.If the loss is less than the monopole gain (∼6 MeV), this loss is compensated by the monopole effect.However, the loss becomes larger for N larger, and at a certain point, the loss exceeds the monopole compensation.The dripline thus arises with the "monopole displacement" from N=22 to N=30 as shown in Fig.18(b)(and also in Fig.18(a) schematically).