The limits of the nuclear landscape explored by the relativistic continuum HatreeBogoliubov theory
Abstract
The groundstate properties of nuclei with 8\leqslant Z\leqslant 120 from the proton drip line to the neutron drip line have been investigated using the spherical relativistic continuum HartreeBogoliubov (RCHB) theory with the relativistic density functional PCPK1. With the effects of the continuum included, there are totally 9035 nuclei predicted to be bound, which largely extends the existing nuclear landscapes predicted with other methods. The calculated binding energies, separation energies, neutron and proton Fermi surfaces, rootmeansquare (rms) radii of neutron, proton, matter, and charge distributions, groundstate spins and parities are tabulated. The extension of the nuclear landscape obtained with RCHB is discussed in detail, in particular for the neutronrich side, in comparison with the relativistic mean field calculations without pairing correlations and also other predicted landscapes. It is found that the coupling between the bound states and the continuum due to the pairing correlations plays an essential role in extending the nuclear landscape. The systematics of the separation energies, radii, densities, potentials and pairing energies of the RCHB calculations are also discussed. In addition, the \alphadecay energies and proton emitters based on the RCHB calculations are investigated.
Contents
 1 Introduction
 2 Theoretical framework
 3 Numerical details
 4 Results and discussions
 5 Summary
1 Introduction
Nuclear landscape or nuclear chart is a diagrammatic representation of atomic nuclei consisting of a bound collection of Z protons and N neutrons using Z and N as axes in the plot. Along or near to the diagonal N=Z line, the stable nuclei inhabit, forming the valley of stability. Away from this valley by adding or removing nucleons, there emerge the vast territory of shortlived radioactive nuclei, which disintegrate by emitting particles or other transmutation. Sequentially adding or removing more nucleons, the boundary of nuclear landscape, i.e., the drip line, is reached. Exploring the limit of nuclear binding has always been a priority in nuclear physics and investigating the properties of all nuclei in the nuclear landscape is crucial for understanding the nuclear force, and the formation of elements in astrophysics, etc. [1, 2].
To date, the existence of about 3000 nuclides has been confirmed experimentally [2, 1]. The protonrich boundary of the nuclear territory has been experimentally determined up to protactinium [1, 2]. However, the neutronrich boundary is known only up to oxygen [2, 1], and most of the neutronrich nuclei far from the valley of stability seem still beyond the experimental capability in the foreseeable future.
The nuclear mass or binding energy is of crucial importance not only to nuclear structure, but also to the weak interactions [3, 4] and the origin of elements [5]. Despite of the considerable achievements in the precise mass measurements [6], the vast territory of neutronrich nuclei are still beyond the experimental capability. Therefore, the development of a reliable theoretical nuclear mass model is very important to grasp a complete understanding of the nature.
Many theoretical efforts have been made to predict nuclear binding energies and to explore the great unknowns of the nuclear landscape [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. With macroscopicmicroscopic models, precise descriptions of nuclear masses have been achieved [7, 8, 9, 10, 11, 12, 13, 14]. Based on the nonrelativistic density functional theory (DFT), a series of HartreeFockBogoliubov mass models has been developed [15, 16, 17, 18, 19]. Recently, to explore the limit of the nuclear landscape and to estimate theoretical uncertainties in drip lines, the SkyrmeDFT calculations with six different parametrizations have been performed, and approximately 7000 bound nuclei with Z\leq120 were predicted [20]. However, the continuum effects are absent in these investigations.
The covariant density functional theory (CDFT) has attracted wide attentions for its successful description of many nuclear phenomena [29, 30, 31, 32, 33, 34]. It includes naturally the nucleonic spin degree of freedom and automatically results in the nuclear spinorbit potential with the empirical strength in a covariant way. It can reproduce well the isotopic shifts in the Pb region [35], and give naturally the origin of the pseudospin symmetry [36] and the spin symmetry in the antinucleon spectrum [37, 38]. Furthermore, it can include the nuclear magnetism [39], that is, a consistent description of currents and timeodd fields, which plays an important role in nuclear magnetic moments [40, 41, 42, 43, 44] and nuclear rotations [33, 45, 46, 47, 48, 49, 50, 51].
Therefore, it is natural to investigate the nuclear landscape based on the CDFT. The masses of about 2000 eveneven nuclei with 8\leq Z\leq 120 up to the proton and neutron drip lines were investigated using relativistic mean field (RMF) theory without pairing correlations in 1997 [21]. Later on, by including pairing correlations with BardeenCooperSchrieffer (BCS) method, the groundstate properties of 1315 eveneven nuclei with 10\leq Z\leq 98 were calculated [22]. By employing the statedependent BCS method with a zerorange pairing force, the first systematic study of the groundstate properties for over 7000 nuclei ranging from the proton drip line to the neutron drip line was performed [23]. To explore the location of the proton and neutron drip lines, a systematic investigation has been performed for eveneven nuclei by employing relativistic HartreeBogoliubov (RHB) theory with four sets of parametrizations [24, 25, 26]. Very recently, the binding energies for eveneven nuclei with proton number ranging from Z=8 to Z=108 were calculated and the dynamical correlation energies which include the rotational correction energies and quadrupole vibrational correction energies were investigated [27, 28].
It is well known that continuum and pairing correlations play a critical role in exotic nuclei. For the exotic nuclei close to the nucleon drip lines, where the Fermi levels are very close to the continuum threshold, pairing correlations can scatter the valence nucleons between the bound states and continuum, and therefore provide a significant coupling between them. As a result, some unbound nuclei predicted without pairing correlations can exist as a bound state. For example, it is found that after taking into account pairing correlations and the contribution from continuum, the neutronrich nuclei {}^{6272}Ca predicted unbound without pairing correlations are found to be bound [52]. Therefore, the couplings between the bound states and the continuum due to the pairing correlations can strongly influence the position of the drip line.
A consistent description of neutronrich nuclei requires a unified and selfconsistent treatment of continuum, meanfield potentials and pairing correlations. By extending the relativistic mean field theory with the Bogoliubov transformation in the coordinate representation, the relativistic continuum HartreeBogoliubov (RCHB) theory was developed [53, 54] and it provides a proper treatment of pairing correlations and meanfield potentials in the presence of the continuum. With the RCHB theory, the first microscopic selfconsistent description of halo in {}^{11}Li has been provided [53] and the giant halos in light and mediumheavy nuclei have been predicted [55, 52, 56]. The RCHB theory has been generalized to treat the odd nucleon system [57], and together with the Glauber model, the chargechanging cross sections from carbon (C) to fluorine (F) isotopes on a carbon target have been well reproduced [58]. For deformed nuclei, much effort has been made to develop a deformed RHB theory in continuum and an interesting shape decoupling between the core and the halo was predicted [59, 60]. Later, the deformed RHB theory in continuum has been extended to incorporate the densitydependent mesonnucleon couplings [61], and the blocking effect which is required for the description of oddnucleon systems [62].
As one of the most successful relativistic energy density functionals, PCPK1 [63], which is fitted to the binding energies, charge radii, and empirical pairing gaps of 60 selected spherical nuclei, has been used successfully in describing not only nuclear groundstate properties [64, 27, 28] but also various excited state properties [65, 66, 67, 68, 69, 70, 71, 72]. In particular, PCPK1 provides a good description for the isospin dependence of binding energy along either the isotopic or the isotonic chain, which makes it more reliable for describing exotic nuclei [63, 64].
Focusing on the continuum effects on the limits of the nuclear landscape, systematic spherical calculations over the nuclear landscape in the framework of RCHB theory with PCPK1 will be performed. As a first attempt, the results from O to Ti isotopes were reported in Ref. [73], and the nuclear landscape from O to Ti is remarkably extended. Here, the results for all isotopes from Z=8 and Z=120 are reported. The paper is organized as follows: In Section 2, the theoretical framework is introduced briefly. The numerical details are given in Section 3. Extensive results are compiled in Section 4, including the binding energies, onenucleon and twonucleon separation energies, neutron and proton Fermi surfaces, rootmeansquare (rms) radii of neutron, proton, matter, and charge distributions, groundstate spins and parities. Finally, a brief summary is presented.
2 Theoretical framework
The covariant density functional theory is constructed with either the finiterange mesonexchange interation or the contact interaction in the pointcoupling representation between nucleons [29, 31, 34, 74, 75, 63]. For the former, the nucleus is described as a system of Dirac nucleons that interact with each other via the exchange of mesons (e.g., scalar meson \sigma, vector meson \omega, and isovector meson \rho). For the latter, the meson exchange in each channel (scalarisoscalar, vectorisoscalar, scalarisovector, and vectorisovector) is replaced by the corresponding local fourpoint (contact) interaction between nucleons.
Following the pointcoupling representation in Ref. [63], one starts with the Lagrangian density as
\displaystyle{\cal L}  \displaystyle=  \displaystyle\bar{\psi}(i\gamma_{\mu}\partial^{\mu}M)\psi\frac{1}{2}\alpha_{% S}(\bar{\psi}\psi)(\bar{\psi}\psi)\frac{1}{2}\alpha_{V}(\bar{\psi}\gamma_{\mu% }\psi)(\bar{\psi}\gamma^{\mu}\psi)\frac{1}{2}\alpha_{TV}(\bar{\psi}\vec{\tau}% \gamma_{\mu}\psi)(\bar{\psi}\vec{\tau}\gamma^{\mu}\psi)  (1)  
\displaystyle\frac{1}{2}\alpha_{TS}(\bar{\psi}\vec{\tau}\psi)(\bar{\psi}\vec{% \tau}\psi)\frac{1}{3}\beta_{S}(\bar{\psi}\psi)^{3}\frac{1}{4}\gamma_{S}(\bar% {\psi}\psi)^{4}\frac{1}{4}\gamma_{V}[(\bar{\psi}\gamma_{\mu}\psi)(\bar{\psi}% \gamma^{\mu}\psi)]^{2}  
\displaystyle\frac{1}{2}\delta_{S}\partial_{\nu}(\bar{\psi}\psi)\partial^{\nu% }(\bar{\psi}\psi)\frac{1}{2}\delta_{V}\partial_{\nu}(\bar{\psi}\gamma_{\mu}% \psi)\partial^{\nu}(\bar{\psi}\gamma^{\mu}\psi)\frac{1}{2}\delta_{TV}\partial% _{\nu}(\bar{\psi}\vec{\tau}\gamma_{\mu}\psi)\partial^{\nu}(\bar{\psi}\vec{\tau% }\gamma_{\mu}\psi)  
\displaystyle\frac{1}{2}\delta_{TS}\partial_{\nu}(\bar{\psi}\vec{\tau}\psi)% \partial^{\nu}(\bar{\psi}\vec{\tau}\psi)\frac{1}{4}F^{\mu\nu}F_{\mu\nu}e% \frac{1\tau_{3}}{2}\bar{\psi}\gamma^{\mu}\psi A_{\mu}, 
where M is the nucleon mass; A_{\mu} and F_{\mu\nu} are respectively the fourvector potential and field strength tensor of the electromagnetic field. Here \alpha_{S},~{}\alpha_{V},~{}\alpha_{TS} and \alpha_{TV} represent the coupling constants for fourfermion pointcoulping terms, \beta_{S},~{}\gamma_{S} and \gamma_{V} are those for the higherorder terms which are responsible for the effects of medium dependence, and \delta_{S},~{}\delta_{V},~{}\delta_{TS} and \delta_{TV} refer to those for the gradient terms which are included to simulate the finite range effects. The subscripts S, V and T stand for scalar, vector and isovector, respectively.
The Hamiltonian density can be obtained by the Legendre transformation,
\displaystyle\mathcal{H}=\frac{\partial\mathcal{L}}{\partial(\partial_{0}\phi)% }\partial_{0}\phi\mathcal{L},  (2) 
where \phi represents the nucleon or photon field. Then, the total Hamiltonian reads,
\displaystyle H=  \displaystyle~{}\int d^{3}r\mathcal{H}  
\displaystyle=  \displaystyle~{}\int d^{3}r\{\bar{\psi}(i\gamma^{i}\partial_{i}+M)\psi  
\displaystyle+\frac{1}{2}\alpha_{S}(\bar{\psi}\psi)(\bar{\psi}\psi)+\frac{1}{2% }\alpha_{V}(\bar{\psi}\gamma_{\mu}\psi)(\bar{\psi}\gamma^{\mu}\psi)+\frac{1}{2% }\alpha_{TV}(\bar{\psi}\vec{\tau}\gamma_{\mu}\psi)(\bar{\psi}\vec{\tau}\gamma^% {\mu}\psi)  
\displaystyle+\frac{1}{3}\beta_{S}(\bar{\psi}\psi)^{3}+\frac{1}{4}\gamma_{S}(% \bar{\psi}\psi)^{4}+\frac{1}{4}\gamma_{V}[(\bar{\psi}\gamma_{\mu}\psi)(\bar{% \psi}\gamma^{\mu}\psi)]^{2}  
\displaystyle\frac{1}{2}\delta_{S}[\partial_{0}(\bar{\psi}\psi)\partial^{0}(% \bar{\psi}\psi)+\nabla(\bar{\psi}\psi)\cdot\nabla(\bar{\psi}\psi)]  
\displaystyle\frac{1}{2}\delta_{V}[\partial_{0}(\bar{\psi}\gamma_{\mu}\psi)% \partial^{0}(\bar{\psi}\gamma^{\mu}\psi)+\nabla(\bar{\psi}\gamma_{\mu}\psi)% \cdot\nabla(\bar{\psi}\gamma^{\mu}\psi)]  
\displaystyle\frac{1}{2}\delta_{TV}[\partial_{0}(\bar{\psi}\vec{\tau}\gamma_{% \mu}\psi)\partial^{0}(\bar{\psi}\vec{\tau}\gamma^{\mu}\psi)+\nabla(\bar{\psi}% \vec{\tau}\gamma_{\mu}\psi)\cdot\nabla(\bar{\psi}\vec{\tau}\gamma^{\mu}\psi)]  
\displaystyle+\frac{1}{4}F^{\mu\nu}F_{\mu\nu}+e\bar{\psi}\gamma^{\mu}\frac{1% \tau_{3}}{2}A_{\mu}\psiF^{0\mu}\partial^{0}A_{\mu}\}.  (3) 
With the meanfield and the nosea approximations, the nucleon field operators can be written as
\displaystyle\psi(x)=\sum_{k}\psi_{k}(x)c_{k},  
\displaystyle\psi^{{\dagger}}(x)=\sum_{k}\psi^{{\dagger}}_{k}(x)c^{{\dagger}}_% {k},  (4) 
where c_{k} is the annihilation operator for a nucleon in the state k and \psi_{k} is the corresponding singleparticle wave function. The operator c_{k} and its conjugate c_{k}^{{\dagger}} satisfy the anticommutation relation for fermions,
\displaystyle\{c_{k},c^{\dagger}_{k^{\prime}}\}=\delta_{kk^{\prime}},  
\displaystyle\{c_{k},c_{k^{\prime}}\}=\{c^{\dagger}_{k},c^{\dagger}_{k^{\prime% }}\}=0.  (5) 
The ground state of a nucleus can be written as
\displaystyle\Phi\rangle=\prod_{i}^{A}c_{i}^{\dagger}0\rangle~{}~{}{\rm with% }~{}~{}\langle\Phi\Phi\rangle=1.  (6) 
The energy density functional for the nucleus system can be represented as
\displaystyle E_{DF}  \displaystyle=\langle\PhiH\Phi\rangle  
\displaystyle=\int d^{3}r\{\sum_{k}\psi^{{\dagger}}_{k}(\bm{\alpha}\cdot\bm{p}% +\beta M)\psi_{k}+\frac{1}{2}\alpha_{S}\rho_{S}^{2}+\frac{1}{2}\alpha_{V}j_{% \mu}j^{\mu}+\frac{1}{2}\alpha_{TV}(\vec{j}_{TV})_{\mu}\vec{j}^{\mu}_{TV}  
\displaystyle~{}~{}+\frac{1}{3}\beta_{S}\rho_{S}^{3}+\frac{1}{4}\gamma_{S}\rho% _{S}^{4}+\frac{1}{4}\gamma_{V}(j_{\mu}j^{\mu})^{2}\frac{1}{2}\delta_{S}[% \partial_{0}\rho_{S}\partial^{0}\rho_{S}+\nabla\rho_{S}\cdot\nabla\rho_{S}]  
\displaystyle~{}~{}\frac{1}{2}\delta_{V}[\partial_{0}j_{\mu}\partial^{0}j^{% \mu}+\nabla j_{\mu}\cdot\nabla j^{\mu}]\frac{1}{2}\delta_{TV}[\partial_{0}(% \vec{j}_{TV})_{\mu}\partial^{0}\vec{j}^{\mu}_{TV}+\nabla(\vec{j}_{TV})_{\mu}% \cdot\nabla\vec{j}^{\mu}_{TV}]  
\displaystyle~{}~{}+eA_{\mu}j_{p}^{\mu}F^{0\mu}\partial^{0}A_{\mu}+\frac{1}{4% }F_{\mu\nu}F^{\mu\nu}\},  (7) 
where the local densities and currents are given by
\displaystyle\rho_{S}  \displaystyle=\langle\Phi:\bar{\psi}\psi:\Phi\rangle=\sum_{k}^{A}\bar{\psi}_% {k}(x)\psi_{k}(x),  (8a)  
\displaystyle j^{\mu}  \displaystyle=\langle\Phi:\bar{\psi}\gamma^{\mu}\psi:\Phi\rangle=\sum_{k}^{A% }\bar{\psi}_{k}(x)\gamma^{\mu}\psi_{k}(x),  (8b)  
\displaystyle\vec{j}^{\mu}_{TV}  \displaystyle=\langle\Phi:\bar{\psi}\gamma^{\mu}\vec{\tau}\psi:\Phi\rangle=% \sum_{k}^{A}\bar{\psi}_{k}(x)\gamma^{\mu}\vec{\tau}\psi_{k}(x),  (8c)  
\displaystyle j^{\mu}_{p}  \displaystyle=\langle\Phi:\bar{\psi}\gamma^{\mu}\frac{1\tau_{3}}{2}\psi:% \Phi\rangle=\sum_{k}^{A}\bar{\psi}_{k}(x)\gamma^{\mu}\frac{1\tau_{3}}{2}\psi_% {k}(x).  (8d) 
By minimizing the energy density functional with respect to the densities, one obtains the Dirac equation for the nucleons:
\displaystyle[\bm{\alpha}\cdot(\bm{pV})+V^{0}+\beta(M+S)]\psi_{k}=\varepsilon% \psi_{k},  (9) 
with the local scalar S(\bm{r}) and vector V^{\mu}(\bm{r}) potentials
\displaystyle S(\mathbf{r})  \displaystyle=  \displaystyle\alpha_{S}\rho_{S}+\beta_{S}\rho_{S}^{2}+\gamma_{S}\rho_{S}^{3}+% \delta_{S}\triangle\rho_{S},  (10)  
\displaystyle V^{\mu}(\mathbf{r})  \displaystyle=  \displaystyle\alpha_{V}j^{\mu}+\gamma_{V}(j_{\mu}j^{\mu})j^{\mu}+\delta_{V}% \triangle j^{\mu}+eA^{\mu}+\alpha_{TV}\tau_{3}\vec{j}^{\mu}_{TV}+\delta_{TV}% \tau_{3}\triangle\vec{j}^{\mu}_{TV}.  (11) 
For a system with time reversal invariance, the spacelike components of the current and the vector potential vanish. Furthermore, one can assume that the nucleon singleparticle states do not mix isospin, i.e. the singleparticle states are eigenstates of \tau_{3}. Therefore, the energy density functional can be written as
\displaystyle E_{\rm DF}  \displaystyle=\langle\PhiH\Phi\rangle  
\displaystyle=\int d^{3}r\{\sum_{k}\psi^{{\dagger}}_{k}(\bm{\alpha}\cdot\bm{p}% +\beta M)\psi_{k}+\frac{1}{2}\alpha_{S}\rho_{S}^{2}+\frac{1}{2}\alpha_{V}\rho_% {V}^{2}+\frac{1}{2}\alpha_{TV}(\rho_{TV})^{2}  
\displaystyle~{}~{}+\frac{1}{3}\beta_{S}\rho_{S}^{3}+\frac{1}{4}\gamma_{S}\rho% _{S}^{4}+\frac{1}{4}\gamma_{V}(\rho_{V})^{4}+\frac{1}{2}\delta_{S}\rho_{S}% \Delta\rho_{S}  
\displaystyle~{}~{}+\frac{1}{2}\delta_{V}\rho_{V}\Delta\rho_{V}+\frac{1}{2}% \delta_{TV}\rho_{TV}\Delta\rho_{TV}+\frac{1}{2}eA_{0}\rho_{p}\}.  (12) 
Pairing correlations are crucial in the description of openshell nuclei. The relativistic HartreeBogoliubov model provides a unified description of particlehole (ph) and particleparticle (pp) correlations on a mean field level by using two average potentials: the selfconsistent mean field \hat{h} and a pairing field \hat{\Delta}. In analogy to Eq. (6) the ground state of a nucleus is described by a generalized Slater determinant \Phi\rangle that represents the vacuum with respect to quasiparticles. The quasiparticle operators are defined by the unitary Bogoliubov transformation of the singlenucleon creation and annihilation operators
\displaystyle\alpha_{k}^{{\dagger}}=\sum_{n}U_{nk}c^{{\dagger}}_{n}+V_{nk}c_{n},  (13) 
and ground state wave function can be written as
\displaystyle\Phi\rangle=\prod_{k}\alpha_{k}\rangle,  (14) 
where the index n refers to the original basis, e.g. an oscillator basis or the coordinates in space, spin and isospin (\bm{r},s,t). The matrix U and V are the HartreeBogoliubov wave functions determined by the variational principle. In the presence of pairing, the singleparticle density matrix is generalized to two different densities: the normal density \hat{\rho} and the pairing tensor \hat{\kappa},
\displaystyle\hat{\rho}_{nn^{\prime}}=\langle\Phic^{{\dagger}}_{n^{\prime}}c_% {n}\Phi\rangle,  
\displaystyle\hat{\kappa}_{nn^{\prime}}=\langle\Phic_{n^{\prime}}c_{n}\Phi\rangle.  (15) 
The RHB energy density functional thus depends on both \hat{\rho} and \hat{\kappa}
\displaystyle E_{\rm RHB}[\hat{\rho},\hat{\kappa}]=E_{\rm DF}[\hat{\rho}]+E_{% \rm pair}[\hat{\kappa}],  (16) 
where E_{DF}[\hat{\rho}] is the usual relativistic meanfield functional defined in Eq. (7) or (12), and the pairing part of the RHB functional reads
\displaystyle E_{\rm pair}[\hat{\kappa}]=\frac{1}{4}\sum_{n_{1}n_{1}^{\prime}}% \sum_{n_{2}n_{2}^{\prime}}\hat{\kappa}^{*}_{n_{1}n_{1}^{\prime}}\langle n_{1}n% _{1}^{\prime}V^{pp}n_{2}n_{2}^{\prime}\rangle\hat{\kappa}_{n_{2}n_{2}^{% \prime}},  (17) 
where \langle n_{1}n_{1}^{\prime}V^{pp}n_{2}n_{2}^{\prime}\rangle represent the matrix elements of the pairing interaction. The RHB equation is obtained by the variational principle,
\displaystyle\left(\begin{array}[]{cc}\hat{h}_{D}\lambda&\hat{\Delta}\\ \hat{\Delta}^{*}&\hat{h}_{D}+\lambda\\ \end{array}\right)\left({U_{k}\atop V_{k}}\right)  \displaystyle=  \displaystyle E_{k}\left({U_{k}\atop V_{k}}\right),  (18) 
where E_{k} is the quasiparticle energy, \lambda the chemical potential or the Fermi surface, and h_{D} the Dirac Hamiltonian in Eq. (9), in which the densities can be constructed by quasiparticle wave functions,
\displaystyle\rho_{S}(\mathbf{r})  \displaystyle=  \displaystyle\sum_{k>0}\bar{V}_{k}(\mathbf{r})V_{k}(\mathbf{r}),  
\displaystyle\rho_{V}(\mathbf{r})  \displaystyle=  \displaystyle\sum_{k>0}V_{k}^{\dagger}(\mathbf{r})V_{k}(\mathbf{r}),  (19)  
\displaystyle\rho_{TV}(\mathbf{r})  \displaystyle=  \displaystyle\sum_{k>0}V_{k}^{\dagger}(\mathbf{r})\tau_{3}V_{k}(\mathbf{r}). 
The pairing potential is determined by
\displaystyle\hat{\Delta}_{n_{1}n_{1}^{\prime}}=\frac{1}{2}\sum_{n_{2}n_{2}^{% \prime}}\langle n_{1}n_{1}^{\prime}V^{pp}n_{2}n_{2}^{\prime}\rangle\hat{% \kappa}_{n_{2}n_{2}^{\prime}},  (20) 
which depends on the pairing tensor \kappa=U^{*}V^{T} and pairing interaction V^{pp} in the particleparticle channel. There have been several types of effective pairing forces V^{pp} in the literatures. The finiterange Gogny force can provide a good treatment of the coupling to the highly excited states but requires more numerical costs [76]. The separable pairing force is also a finiterange force but has a separable form, which reduces significantly the computational costs [77]. The densitydependent zerorange force has a simple form and a realistic densitydependent behavior
\displaystyle V^{pp}(\mathbf{r_{1}},\mathbf{r_{2}})=V_{0}\delta(\mathbf{r_{1}}% \mathbf{r_{2}})\frac{1}{2}(1P^{\sigma})(1\eta\frac{\rho(\mathbf{r_{1}})}{% \rho_{0}}).  (21) 
Here, V_{0} is the interaction strength and \rho_{0} is the saturation density of nuclear matter. Note that \eta here can be either 0 or 1. For \eta=1, it is called as a surface pairing force, because the force is dramatically suppressed in the nuclear interior and has its largest contribution in the surface. Conversely, \eta=0 corresponds to a volume pairing force. For a zerorange pairing force, one has to introduce an energy cutoff to avoid the divergence of the pairing energy and, thus, the interaction strength V_{0} has to be properly justified in a given cutoff.
With spherical symmetry, the quasiparticle wave function in the coordinate space can be written as
\displaystyle U_{k}=\frac{1}{r}\left(\begin{array}[]{c}iG_{U}^{k}(r)Y_{jm}^{l}% (\theta,\phi)\\ F_{U}^{k}(r)(\boldsymbol{\sigma\cdot\hat{r}})Y_{jm}^{l}(\theta,\phi)\end{array% }\right)\chi_{t}(t),~{}~{}V_{k}=\frac{1}{r}\left(\begin{array}[]{c}iG_{V}^{k}(% r)Y_{jm}^{l}(\theta,\phi)\\ F_{V}^{k}(r)(\boldsymbol{\sigma\cdot\hat{r}})Y_{jm}^{l}(\theta,\phi)\end{array% }\right)\chi_{t}(t).  (22) 
The corresponding RHB equation can be expressed as the following radial integraldifferential equations in coordinate space [54]:
\displaystyle\frac{dG_{U}}{dr}+\frac{\kappa}{r}G_{U}(r)(E+\lambdaV(r)+S(r))F% _{U}(r)+r\int r^{\prime}dr^{\prime}\Delta_{F}(r,r^{\prime})F_{V}(r^{\prime})=0,  
\displaystyle\frac{dF_{U}}{dr}\frac{\kappa}{r}F_{U}(r)+(E+\lambdaV(r)S(r))G% _{U}(r)+r\int r^{\prime}dr^{\prime}\Delta_{G}(r,r^{\prime})G_{V}(r^{\prime})=0,  
\displaystyle\frac{dG_{V}}{r}+\frac{\kappa}{r}G_{V}(r)+(E\lambda+V(r)S(r))F_% {V}(r)+r\int r^{\prime}dr^{\prime}\Delta_{F}(r,r^{\prime})F_{U}(r^{\prime})=0,  
\displaystyle\frac{dF_{V}}{r}\frac{\kappa}{r}F_{V}(r)(E\lambda+V(r)+S(r))G_% {V}(r)+r\int r^{\prime}dr^{\prime}\Delta_{G}(r,r^{\prime})G_{U}(r^{\prime})=0.  (23) 
If the zerorange pairing force is applied, the above coupled integraldifferential equations can be reduced to differential ones, which can be solved in coordinate space using the shooting method with RungeKutta algorithms [54]. After the solution, new densities and fields are obtained, which are iterated in the differential equations until convergence is achieved.
Finally, one can calculate the total energy of a nucleus by [54]
\displaystyle E_{\rm RHB}  \displaystyle=\sum_{k}(\lambdaE_{k})V^{2}_{k}E_{\rm pair}\int d^{3}r\{\frac% {1}{2}\alpha_{S}\rho_{S}^{2}+\frac{1}{2}\alpha_{V}\rho_{V}^{2}+\frac{1}{2}% \alpha_{TV}(\rho_{TV})^{2}  
\displaystyle~{}~{}+\frac{2}{3}\beta_{S}\rho_{S}^{3}+\frac{3}{4}\gamma_{S}\rho% _{S}^{4}+\frac{3}{4}\gamma_{V}(\rho_{V})^{4}+\frac{1}{2}\delta_{S}\rho_{S}% \Delta\rho_{S}  
\displaystyle~{}~{}+\frac{1}{2}\delta_{V}\rho_{V}\Delta\rho_{V}+\frac{1}{2}% \delta_{TV}\rho_{TV}\Delta\rho_{TV}+\frac{1}{2}eA_{0}\rho_{p}\}\frac{1}{2MA}% \langle\hat{P}^{2}_{cm}\rangle,  (24) 
and the rootmeansquare (rms) radius by
\displaystyle R_{\rm{rms}}=\langle r^{2}\rangle^{1/2}=\biggl{\{}\int d^{3}\bm{% r}[r^{2}\rho_{V}(\bm{r})]\biggr{\}}^{1/2}.  (25) 
Note that for oddA and oddodd nuclei, the blocking effects of the unpaired nucleon(s) have to be taken into account. Generally, the ground state of an odd system can be written as an onequasiparticle state,
\displaystyle\Phi_{1}\rangle=\alpha^{{\dagger}}_{1}\prod_{k}\alpha_{k}\rangle,  (26) 
where \alpha^{{\dagger}}_{1} corresponds to the quasiparticle state which is blocked. As explained in great detail in Ref. [78], this state \Phi_{1}\rangle can be regarded as the vacuum with respect to a new set of quasiparticle operators (\alpha^{\prime}_{1},...\alpha^{\prime}_{N}) with
\displaystyle\alpha^{\prime}_{1}=\alpha^{{\dagger}}_{1},~{}\alpha^{\prime}_{2}% =\alpha_{1},~{}...,~{}\alpha^{\prime}_{N}=\alpha^{\prime}_{N}.  (27) 
In such a way, the oddnucleon systems can be simply treated just like an eveneven system. The only difference is the exchange of the quasiparticle operator \alpha_{1}^{{\dagger}}\leftrightarrow\alpha_{1}, which corresponds to the exchange of the column (U_{n1},V_{n1})\leftrightarrow(V_{n1}^{*},U_{n1}^{*}) in the Bogoliubov transform matrix.
3 Numerical details
To explore the nuclear landscape extended by the continuum effects, the systematic calculations for nuclei with 8\leqslant Z\leqslant 120 from proton drip line to neutron drip line are performed by using the relativistic continuum HartreeBogoliubov theory with PCPK1. In order to describe the continuum and its coupling to the bound states properly, the RCHB equation (23) is solved in coordinate space using the shooting method. Each atomic nucleus is described within a box by imposing the boundary condition U_{k}(r=R_{\rm box})=V_{k}(r=R_{\rm box})=0.
As shown in Fig. 1, the convergence of the RCHB solutions with respect to the box size, the mesh size and the angular momentum cutoff has been examined for nuclei {}^{304}120. In Fig. 1, the box size R_{\rm box}=20 fm leads to a relative accuracy of 0.002\% in total energy for {}^{304}120 in comparison with R_{\rm box}=27.5 fm. The energy difference between the calculations with \Delta r=0.1 fm and \Delta r=0.05 fm is smaller than 0.0035 MeV for {}^{304}120, which is about 0.0002\% of the total energy. Similar convergence has also been confirmed for the angular momentum cutoff J_{\rm max}=19/2~{}\hbar, as shown in Fig. 1. Note that the pairing correlation is neglected in the convergence check for angular momentum cutoff in order to avoid renormalizing the strength of the zerorange pairing force to the corresponding model space. In the following calculations, we fix the box size R_{\rm box}=20 fm, the mesh size \Delta r=0.1 fm, and the angular momentum cutoff J_{\rm max}=19/2~{}\hbar.
3.1 Relativistic density functional PCPK1
In the present systematic calculations, the relativistic density functional PCPK1 [63] listed in Table 1 is adopted for the particlehole channel. PCPK1 has turned out to be very successful in providing good descriptions for the isospin dependence of the binding energy along either the isotopic or the isotonic chain. In Ref. [64], a crucial test of the predictive power of PCPK1 was performed by the comparison with the new and accurate experimental masses in the element range from 50 to 91 [79], and it was found that the experimental data can be reproduced quite well with a rms deviation 0.859 MeV. PCPK1 has also been applied successfully to investigate nuclear lowlying excited states [65, 66, 67, 68, 69, 70, 71, 72], magnetic rotation [50, 80, 81, 82, 83], antimagnetic rotation [49, 51, 84, 85, 86], chiral rotation modes [87], fission barrier [88], and superheavy nuclei [89], etc.
Coupling constant  Value  Dimension 

\alpha_{S}  3.96291\times 10^{4}  MeV{}^{2} 
\beta_{S}  8.6653\times 10^{11}  MeV{}^{5} 
\gamma_{S}  3.80724\times 10^{17}  MeV{}^{8} 
\delta_{S}  1.09108\times 10^{10}  MeV{}^{4} 
\alpha_{V}  2.6904\times 10^{4}  MeV{}^{2} 
\gamma_{V}  3.64219\times 10^{18}  MeV{}^{8} 
\delta_{V}  4.32619\times 10^{10}  MeV{}^{4} 
\alpha_{TV}  2.95018\times 10^{5}  MeV{}^{2} 
\delta_{TV}  4.11112\times 10^{10}  MeV{}^{4} 
3.2 Pairing strength
For the particleparticle channel, the surface pairing force with \eta=1 in Eq. (21) is used. The parameters in the pairing force include the saturation density \rho_{0} and pairing strength V_{0}. The empirical saturation density of nuclear matter \rho_{0} is 0.152 fm{}^{3}. The pairing strength V_{0} is fitted to the oddeven mass differences for given energy cutoff E_{\rm cut}. In the present RCHB calculations, the pairing strength V_{0}=342.5 MeV fm{}^{3} for both neutron and proton is fixed by reproducing the experimental oddeven mass differences in the Ca, Sn, Pb and U isotopic chains, as well as the N = 20, 50 isotonic chains with a cutoff energy of 100 MeV for the pairing window. For the oddeven mass difference, a threepoint formula is utilized,
\Delta=\frac{(1)^{N}}{2}\left[\ E_{b}(Z,N+1)2E_{b}(Z,N)+E_{b}(Z,N1)\right]\,.  (28) 
The corresponding comparison between the theoretical and experimental oddeven mass differences is shown in Fig. 2.
3.3 Blocking effects for oddA and oddodd nuclei
For the oddA or oddodd nuclei, one has to take into account the blocking effect of the odd nucleon(s) [78]. The procedure to determine the ground state of an oddA nucleus with even Z and odd N in RCHB calculations is as follows:
(i) Perform the calculation for its neighboring eveneven nucleus (Z, N1).
(ii) From the single particle levels in the eveneven nucleus (Z, N1), choose two neutron levels just below the Fermi surface (labeled as E_{1} and E_{2}) and two neutron levels just above the Fermi surface (labeled as E_{3} and E_{4}).
(iii) Perform the calculations for the oddA nucleus (Z, N) with blocking levels E_{1}, E_{2}, E_{3} and E_{4}, respectively, and choose the one with the lowest total energy as the ground state of the oddA nucleus (Z, N).
The blocking procedure for oddodd nucleus (Z, N) is similar to the oddA nuclei, but requires blocking for both the proton and neutron levels at the same time. From the single particle levels in the neighboring eveneven nucleus (Z1, N1), one chooses the blocking proton and neutron levels (two just below and two just above the Fermi surfaces). By performing the calculations for the oddodd nucleus (Z, N) with the blocking proton and neutron levels, the one with the lowest total energy is chosen as the ground state of the oddodd nucleus (Z, N).
4 Results and discussions
4.1 Nuclear mass
We perform systematic calculations for all nuclei from Z=8 to Z=120 by varying the neutron number from the proton drip line to neutron one. In Table 4, the groundstate properties of these nuclei are summarized. The mass number A, neutron number N, binding energy E_{\rm b}, neutron Fermi surface \lambda_{n}, proton Fermi surface \lambda_{p}, matter rootmeansquare (rms) radius R_{m}, neutron rms radius R_{n}, proton rms radius R_{p}, charge radius R_{\mathrm{c}}, groundstate spin and the parity of proton j^{\pi}(P) and neutron j^{\pi}(N) are listed. The binding energy per nucleon E_{\rm b}/A, twoneutron S_{2n}, twoproton S_{2p}, oneneutron S_{n} and oneproton separation energy S_{p} for each nucleus are also provided. In addition, the available experimental binding energies [6] and charge radii [90] are shown for comparison. There are totally 9035 nuclei from O(Z=8) to Z=120 which are predicted to be bound by the RCHB theory with the relativistic density functional PCPK1. For guidance, one or two unbound nuclei just outside the drip lines are listed and underlined in the table as well.
In Fig. 3, the nuclear landscape from O (Z=8) to Z=120 explored by the RCHB theory with PCPK1 is shown, where the squares represent the bound nuclei. The squares with cross denote the 288 stable or practically stable (that is, have halflives longer than the expected life of the Solar System) nuclei existing in nature. Among these 9035 bound nuclei, the masses of 2284 nuclei have been measured experimentally [6]. The relative binding energy differences (E^{\mathrm{Exp.}}_{\mathrm{b}}E^{\mathrm{Cal.}}_{\mathrm{b}})/E^{\mathrm{Exp% .}}_{\mathrm{b}} for these measured nuclei scaled by colors are plotted. It is found that most of the deviations are between 0.75% and 0.75%, and the root of the relative square (rrs) deviation \sigma_{r}=\sqrt{\sum_{i}^{n}\frac{(E^{\rm Exp.}_{\rm b}E^{\rm Cal.}_{\rm b})% ^{2}/(E^{\rm Exp.}_{\rm b})^{2}}{n}} for these 2284 nuclei is around 0.710\%. It is noted that the relative differences increase when the neutron (proton) number goes away from the magic number, due to the spherical symmetry adopted here.
Therefore, it is interesting to focus the comparisons on the nuclei with either neutron or proton magic number. For the nuclei considered here (Z=8, 20, 28, 50, 82 or N=8, 20, 28, 50, 82, 126), the corresponding rrs deviation is reduced to 0.583\%. The rootmeansquare (rms) deviation of the binding energy \sigma=\sqrt{\sum_{i}^{n}(E_{i}^{\mathrm{Exp.}}E_{i}^{\mathrm{Cal.}})^{2}/n} for these either neutron or proton magic nuclei is 2.157 MeV, which is remarkably smaller than the rms deviation of 7.960 MeV for the 2284 nuclei with masses measured. Therefore, we can infer that the deviation between the present RCHB calculations and experimental data comes mainly from the absence of the deformation effects.
To investigate the deformation effects, the rms deviation of binding energy \sigma in RCHB calculations is compared with the average quadrupole deformation parameter \langle\beta_{2}\rangle=(\sum^{n^{\prime}}_{i}\beta_{2}^{i})/n^{\prime} for each isotopic chain in Fig. 4, where the deformation parameters for eveneven nuclei are from 3DRHB [91] calculations with PCPK1. It can be seen that, for isotopic chains, the rms deviations are small at the closed shell Z=20, 28, 50, 82, and large around the middle of two closed shells. Furthermore, a similar tendency for the average deformations is also found in the lower panel of Fig. 4(a). Analogously, for isotonic chains, the variation of \sigma is similar to that of \langle\beta_{2}\rangle as well. These correspondences between the \sigma and \langle\beta_{2}\rangle indicate that taking the deformation effects into account in the future will improve the agreement with the experimental data.
It should be emphasized that it is not the aim of the present RCHB calculations to reproduce the experimental binding energy to a high level of accuracy, but rather to investigate how the continuum states influence the positions of the drip lines, i.e., to which extent the nuclear landscape will be extended by considering the coupling to the continuum. In Fig. 3, the drip lines determined from the mass table TMA [23] , HFB21 [18], FRDM [7], WS3 [11], UNEDF [20], NL3* [24, 25], and also PCPK1 without pairing correlations are shown. The comparisons of the nuclear landscape predicted by different models will be discussed in the following sections in detail.
4.2 The limits of the nuclear landscape
4.2.1 Onenucleon and twonucleon separation energies
The oneneutron separation energies S_{n} and oneproton separation energies S_{p} are respectively defined as
\displaystyle S_{n}(Z,N)=E_{\rm b}(Z,N)E_{\rm b}(Z,N1),  
\displaystyle S_{p}(Z,N)=E_{\rm b}(Z,N)E_{\rm b}(Z1,N).  (29) 
The twoneutron separation energies S_{2n} and twoproton separation energies S_{2p} are respectively defined as
\displaystyle S_{2n}(Z,N)=E_{\rm b}(Z,N)E_{\rm b}(Z,N2),  
\displaystyle S_{2p}(Z,N)=E_{\rm b}(Z,N)E_{\rm b}(Z2,N).  (30) 
These quantities provide information on whether a nucleus is stable against one or two nucleon emissions, and thus define the nucleon drip lines. In the present work, we consider a nucleus is bound only if both the one and twonucleon separation energies of this nucleus are positive. For each isotopic chain, the last bound nucleus in the neutronrich side is the neutron dripline nucleus and that in the neutrondeficient side is the proton dripline nucleus.
In Fig. 5, the twoneutron separation energies S_{2n} of the bound nuclei predicted by the RCHB theory with PCPK1 are shown. From a global view, S_{2n} is large at proton dripline side and close to zero at neutron dripline side. For a given isotopic (isotonic) chain, S_{2n} decreases with the increasing neutron number, while increases with the proton number. It is found that there are 87 nuclei with the predicted twoneutron separation energies larger than 30 MeV, and most of them are located at the protonrich side of nuclear landscape with Z\leqslant 50. There are 388 nuclei with S_{2n} in the range 21–30 MeV, mainly located in the protonrich region of the nuclear chart; 1967 nuclei with S_{2n} in the range of 12–21 MeV, most of which are near the valley of stability; 4715 nuclei with S_{2n} in the range 312 MeV which lie on the neutronrich region mostly. In addition, there are 1828 nuclei with their S_{2n} less than 3 MeV, and most of them are located in the region far from the stability line, and are even approaching to the neutron drip line.
It should be noted that, there are 586 weakly bound nuclei with S_{2n}\leqslant 1 MeV in the RCHB calculations. They are the extremely neutronrich nuclei predicted by RCHB, and many of them lie even beyond the neutron drip lines predicted by the other nuclear mass models. For these weakly bound nuclei, as the neutron Fermi surface is close to the continuum threshold, pairing correlations could scatter the nucleons from bound states to the continuum, thus provide a significant coupling between the continuum and bound states [73]. The RCHB theory allows a proper treatment of the continuum and the coupling to the bound states, so it predicts a more extended neutron drip line than the other models. In addition, the nearly vanishing S_{2n} around the neutron drip line might be regarded as a sign of the neutron giant halo [55, 52], which invokes the further analysis of neutron radii and single particle levels.
S_{2n} is a widely used probe of the neutron shell structure. Generally, for a given isotopic chain S_{2n} decreases smoothly with the neutron number N, except at a magic number where S_{2n} drops significantly. An abrupt decline of S_{2n} indicates the occurrence of neutron shell closure. It can be seen in Fig. 5 that, the significant drops exist at the traditional magic numbers N=20,~{}28,~{}50,~{}82 and 126, which demonstrates that these shell closures are well reproduced in the RCHB theory. Apart from this, a dramatic decline of S_{2n} can be found at N=184, which indicates that neutron number N=184 may be a new magic number [92].
In Fig. 6, the oneneutron separation energies S_{n} of bound nuclei predicted by the RCHB theory are shown with different colors. Analogous to S_{2n}, S_{n} decreases (increases) with the increase of the neutron (proton) number for a given isotopic (isotonic) chain. It is found that 95 nuclei have their S_{n} larger than 15 MeV, 214 nuclei within 12–15 MeV, 610 nuclei within 9–12 MeV, 1516 nuclei within 6–9 MeV, 2583 nuclei within 3–6 MeV, and 4006 nuclei less than 3 MeV. Note that there are 1096 nuclei with S_{n}\leqslant 1 MeV and 466 nuclei with S_{n}\leqslant 0.5 MeV. Due to the pairing correlations, which make the nuclei with even nucleon number more bound than their neighbors with odd nucleon number, S_{n} has a ragged evolution pattern with the variation of neutron number that zigzags between odd and evenparticle nuclei.
In addition to the neutron separation energies S_{2n} and S_{n}, the twoproton S_{2p} and oneproton separation energies S_{p} of the bound nuclei are scaled by colors in Figs. 7 and 8. The proton separation energies increase with the neutron number for a given isotopic chain, while decrease with the proton number for a given isotonic chain. In the present calculations, it is found that there are 520 nuclei with S_{2p}\geqslant 40 MeV, 594 nuclei with 36\leqslant S_{2p}<40 MeV, 802 nuclei with 32\leqslant S_{2p}<36 MeV, 977 nuclei with 28\leqslant S_{2p}<32 MeV, 1150 nuclei with 24\leqslant S_{2p}<28 MeV, 1088 nuclei with 20\leqslant S_{2p}<24 MeV, 957 nuclei with 16\leqslant S_{2p}<20 MeV, 857 nuclei with 12\leqslant S_{2p}<16 MeV, 781 nuclei with 8\leqslant S_{2p}<12 MeV, 731 nuclei with 4\leqslant S_{2p}<8 MeV and 381 nuclei with 0\leqslant S_{2p}<4 MeV. It is noted that the proton drip line predicted by the RCHB calculations is close to the proton drip lines predicted by the other mass tables. At the proton magic number 20, 28, 50 and 82, S_{2p} changes abruptly. This indicates that the RCHB theory reproduces the traditional proton closed shells quite well. Note that a sudden change of the S_{2p} can be also found at Z=92, while this has been considered as a pseudo shell in the previous relativistic mean field calculations [93].
In Fig. 8, there are 589 nuclei with S_{p}\geqslant 20 MeV, 1348 nuclei with 16\leqslant S_{p}<20 MeV, 2083 nuclei with 12\leqslant S_{p}<16 MeV, 2075 nuclei with 8\leqslant S_{p}<12 MeV, 1648 nuclei with 4\leqslant S_{p}<8 MeV and 1227 nuclei with 0\leqslant S_{p}<4 MeV. Due to the pairing correlations, S_{p} is staggering between odd and evenparticle nuclei for a given isotonic chain. It should be noted that the number of the nuclei with S_{2p}\leqslant 1 MeV is 176, which is much less than the number of the nuclei with S_{2n}\leqslant 1 MeV. This is consistent with the number of halo nuclei observed in the protonrich side is much less than that in the neutronrich side, which should be attributed to the existence of the Coulomb barrier. Moreover, although the proton Fermi surface is approaching zero for the neutrondeficient nuclei, the contributions from the continuum are suppressed due to the Coulomb barrier in comparison with that for the neutronrich nuclei. This also explains why the proton drip line obtained by the RCHB calculations is very close to the ones predicted by the other mass models.
4.2.2 Neutron drip line
In Fig. 3, the nucleon drip lines predicted by the mass tables TMA [23], HFB21 [18], FRDM [7], WS3 [11], UNEDF [20], and NL3* [24, 25] have been plotted for comparison, together with the drip line given by the PCPK1 calculations without taking pairing correlations into account. At present, the experimental protonrich border of the nuclear territory has been reached up to protactinium [1] (Z=91). Due to the Coulomb repulsive interaction among protons, the proton drip line does not lie so far away from the valley of stability. Moreover, the proton continuum is effectively shifted up in energy as a result of the confining effect of the Coulomb barrier. Therefore, the proton drip line obtained by the present calculations is very close to those from other mass tables, which are roughly the same with the experimental observations.
On the neutronrich side, however, the neutronrich boundary is known only up to oxygen (Z=8) experimentally. In contrast to the proton drip line, the locations of neutron drip lines from the various mass tables obviously differ with each other, especially for the heavy mass region. The differences of neutron drip lines increase with the proton number. The neutron drip line predicted by the RCHB theory is most extended compared with other models. For example, the last bound neutronrich nucleus of the Kr (Z=36) isotope is {}^{118}Kr in the calculations of FRDM, WS3, HFB21, TMA and the PCPK1 calculations without pairing correlations, {}^{122}Kr in UNDEF and {}^{128}Kr in NL3*, while the present RCHB calculations predict {}^{136}Kr as the last bound nucleus of Kr isotopes in the neutron rich side.
Element  Neutron number  Element  Neutron number  Element  Neutron number  

(Z)  No pairing  RCHB  (Z)  No pairing  RCHB  (Z)  No pairing  RCHB 
O ( 8)  20  20  Pd (46)  112  118  Po ( 84)  184  184 
Ne (10)  20  32  Cd (48)  112  126  Rn ( 86)  184  184 
Mg (12)  34  34  Sn (50)  126  126  Ra ( 88)  184  258 
Si (14)  34  38  Te (52)  126  126  Th ( 90)  184  258 
S (16)  40  40  Xe (54)  126  126  U ( 92)  234  258 
Ar (18)  40  44  Ba (56)  126  126  Pu ( 94)  258  258 
Ca (20)  40  60  Ce (58)  126  126  Cm ( 96)  258  258 
Ti (22)  54  62  Nd (60)  126  168  Cf ( 98)  258  258 
Cr (24)  62  68  Sm (62)  126  168  Fm (100)  258  258 
Fe (26)  62  68  Gd (64)  142  176  No (102)  258  258 
Ni (28)  70  70  Dy (66)  168  184  Rf (104)  258  258 
Zn (30)  70  72  Er (68)  184  184  Sg (106)  258  258 
Ge (32)  70  82  Yb (70)  184  184  Hs (108)  258  258 
Se (34)  82  94  Hf (72)  184  184  Ds (110)  258  258 
Kr (36)  82  100  W (74)  184  184  Cn (112)  258  258 
Sr (38)  102  110  Os (76)  184  184  Fl (114)  258  258 
Zr (40)  112  112  Pt (78)  184  184  Lv (116)  258  288 
Mo (42)  112  112  Hg (80)  184  184  Og (118)  258  288 
Ru (44)  112  112  Pb (82)  184  184  Z=120  288  288 

Element  \Delta N  Element  \Delta N  Element  \Delta N  Element  \Delta N  Element  \Delta N 
(Z)  (Z)  (Z)  (Z)  (Z)  
O (8)  2  Ni (28)  6  Rh (45)  11  Dy (66)  25  Fr (87)  29 
Ne (10)  10  Cu (29)  6  Pd (46)  13  Ho (67)  23  Ra (88)  28 
Na (11)  10  Zn (30)  1  Ag (47)  9  Er (68)  24  Ac (89)  28 
Mg (12)  6  Ga (31)  1  Cd (48)  6  Tm (69)  26  Th (90)  37 
Al (13)  8  Ge (32)  2  In (49)  4  Yb (70)  20  Pa (91)  43 
Si (14)  6  As (33)  1  Sn (50)  6  Lu (71)  21  U (92)  39 
P (15)  8  Se (34)  6  Sb (51)  8  Hf (72)  21  Np (93)  38 
S (16)  6  Br (35)  7  Te (52)  2  Ta (73)  19  Pu (94)  44 
K (19)  14  Kr (36)  9  I (53)  1  W (74)  13  Am (95)  45 
Ca (20)  10  Rb (37)  10  Xe (54)  1  Re (75)  16  Cm (96)  38 
Sc (21)  10  Sr (38)  14  Pr (59)  20  Os (76)  6  Bk (97)  38 
Ti (22)  12  Y (39)  21  Nd (60)  16  Ir (77)  6  Cf (98)  37 
V (23)  10  Zr (40)  19  Pm (61)  15  Pt (78)  6  Es (99)  35 
Cr (24)  10  Nb (41)  12  Sm (62)  19  Au (79)  2  Fm (100)  36 
Mn (25)  8  Mo (42)  12  Eu (63)  22  Hg (80)  2  Md (101)  28 
Fe (26)  5  Tc (43)  13  Gd (64)  25  Tl (81)  2  No (102)  28 
Co (27)  4  Ru (44)  10  Tb (65)  27  Rn (86)  31  Lr (103)  28 

In order to explore where and how the continuum extends the neutron drip line, we compare the neutron drip lines obtained from the RCHB calculations with those from the calculations without pairing correlation. The most neutronrich eveneven nuclei predicted to be bound in the RCHB theory and the calculations without pairing is listed in Table 2 for different isotopic chains. Marked at subshell or shell closures N=40, 70, 82, 112, 126, 184 and 258, the extensions are dramastic when the next subshell or shell is occupied.
A similar comparison is performed for the RCHB theory and the FRDM. The number of neutronrich nuclei predicted to be bound by the RCHB, but unbound in the FRDM are listed in Table 3. It can be found that the neutron drip line of the RCHB results is much farther away from the \betastability line than that of the FRDM for almost all isotopic chains. For example, from the isotopic chains Se (Z=34) to Zr (Z=40), the FRDM gives the drip line at N=82. The RCHB theory predicts 6 to 19 more bound neutronrich nuclei from the isotopic chains Se (Z=34) to Zr (Z=40).
Comparing to the drip line of RHB calculations in harmonic oscillator (HO) basis with NL3* [24, 25], the present RCHB calculations show considerable difference from it. The considerable difference of the drip line between RHB in HO calculations and present RCHB calculations are due to different functionals (PCPK1 or NL3*), different pairing forces used, and the treatment of continuum in coordinate or HO space. Future investigations with same density functional and pairing force in both coordinate space and in harmonic oscillator basis are needed.
4.2.3 Continuum effects
To illustrate the extension of the neutron drip line in RCHB theory, by taking {}^{60,66,80}Ca, {}^{118,124,136}Kr and {}^{188,194,230}Sm as examples, we compare the neutron single particle levels in canonical basis from the RCHB calculations with those of {}^{60}Ca, {}^{118}Kr and {}^{188}Sm from the calculations without pairing, respectively in Figs. 911.
For Ca isotopes in Fig. 9, {}^{60}Ca is the neutron dripline nucleus in the calculations without pairing, while with 20 more neutrons, {}^{80}Ca is the neutron dripline nucleus in the RCHB calculations. In the calculations without pairing, the neutron Fermi surface of {}^{60}Ca is 1f_{5/2}, the last bound level. The occupation of the next level 1g_{9/2} which has a positive single particle energy will lead to the unbound nucleus {}^{62}Ca. In RCHB theory, due to the pairing correlations, the neutron Fermi surface locates between 1f_{5/2} and 1g_{9/2}, and the continuum in the canonical basis also have small occupation probabilities in {}^{60}Ca. Adding more neutrons, see {}^{66}Ca in the figure, although the occupation probabilities of the continuum in the canonical basis become larger, the neutron Fermi surface remains negative. This means that the nucleus is still bound. Moreover, with more neutrons added, the level 1g_{9/2} becomes negative and the neutron Fermi surface remains negative up to {}^{80}Ca, which is turned out to be the last bound nucleus in the RCHB calculations.
4.2.4 Discussion on deformation effects
It should be noted that the present survey is performed with the assumption of spherical symmetry. A natural question is whether and how the deformation will influence the position of the neutron drip line.
To explore the deformation effects on the position of neutron drip line, taking argon isotopic chain as an example, the calculations with the deformed relativistic HartreeBogoliubov theory in continuum (DRHBc) [59, 60] are performed, where the axial deformation, continuum and pairing correlation are treated selfconsistently. In the DRHBc calculations, the relativistic density functional PCPK1 is employed. We use a box size R_{\rm box}= 20 fm, a mesh size \Delta r= 0.1 fm, a cutoff energy of WoodsSaxon basis E^{+}_{\rm cut}= 160 MeV and a cutoff of angular moment J_{\rm max}=25/2\hbar. For the particleparticle channel, the pairing force and pairing parameters are the same as Ref. [60].
In Fig. 12, the twoneutron separation energies from RCHB and DRHBc calculations for argon (Ar) isotopes and the corresponding groundstate deformations from the DRHBc calculations are shown. It is can be seen that most argon isotopes are oblate in ground states, and the deformation parameter \beta_{2} increases with the increase of neutron number for N>40. In Fig. 12(a), S_{2n} from the DRHBc calculations are close to those from the RCHB calculations globally. However, the deviation exists at N>44 where the large deformation occurs. In the RCHB calculations, the S_{2n} becomes negative after N=44 and still decreases with the increase of neutron number, while S_{2n} from DRHBc calculations keeps positive for 44<N<52 and becomes negative for N=54, i.e. the neutron dripline nuclei of argon isotopes are {}^{62}Ar in the RCHB and {}^{70}Ar in the DRHBc calculations, respectively.
In Fig. 13, the pairing energies of argon isotopes from the DRHBc and RCHB calculations are shown. It can be seen that proton pairing energies of DRHBc calculations are close to zero and deviate from the ones of RCHB calculations for {}^{36}Ar,{}^{4456}Ar and {}^{6470}Ar which may be attributed to the occurrence of large deformation there. On the other hand, neutron pairing energies from these two calculations present similar evolution trend. For N>44, the pairing energy in RCHB calculations increases slowly with the neutron number, while the pairing energy in DRHBc calculation increases rapidly.
Fig. 14 shows the neutron single particle energy from DRHBc calculations for {}^{5870}Ar. It can be seen that, with the evolution of neutron number, the single particle levels around the Fermi surface are getting denser and denser, and the pairing correlation is enhanced. The deformation, continuum and pairing correlation extend the neutron dripline nucleus from {}^{62}Ar in the RCHB to {}^{70}Ar in the DRHBc calculations.
From above discussions, one can find that the deformation would affect the position of neutron drip line. For argon isotopes, selfconsistent treatment of the deformation and continuum extends the neutron drip line to more neutronrich region.
In Fig. 3, the deformed RHB calculations with NL3* in harmonic oscillator basis [24] predicted less neutronrich dripline nuclei than the RCHB calculations for most of the isotopic chains. Although quite timeconsuming and numerically challenging, considering the deformation and continuum effects simultaneously is essential to determine the drip line and should be investigated in the future.
4.3 Radii of nucleon distributions
4.3.1 Charge radii
The nuclear charge radius is one of the most important nuclear observables that estimates the size of nuclei. In Fig. 15, the charge radii calculated by the RCHB theory with PCPK1 are shown as a function of neutron number in comparison with the data available. The charge radii in the RCHB calculations reproduce the experiment quite well for the nuclei with data available [90], and the charge radii of the nuclei without experimental data are predicted.
In Fig. 16, the available experimental charge radii, the RCHB calculated charge radii with PCPK1, and the difference between the RCHB calculated charge radii and the experimental data [90] are shown. One can see that most deviations are in the range of 0.05\sim 0.05 fm, and the rms deviation \sigma of the RCHB calculations from the data is 0.0358 fm, which is in comparison with the results from the RHB calculations for NL3* \sigma=0.0407 fm, for DDME2 \sigma=0.0376 fm, for DDME\delta \sigma=0.0412 fm, for DDPC1 \sigma=0.0402 fm [25], the HFB calculations \sigma=0.026 fm [94], and the formula based on the WeizsackerSkyrme mass model \sigma=0.022 fm [95].
In Fig. 16, larger deviations are located between N=28 and 50, N=50 and 82, as well as N=82 and 126. These deviations are expected to be reduced when the deformation effects are included.
4.3.2 Neutron radii
In Fig. 17, the calculated neutron rms radii of the neutron density distribution for eveneven nuclei with 8\leqslant Z\leqslant 120 are shown as a function of neutron number. In addition, the empirical formula R_{n}=r_{0}N^{1/3} is shown for guidance with r_{0}=1.140 fm determined by R_{n} of {}^{208}Pb. Except for extremely neutronrich nuclei, the systematic trend of the neutron radii follows the simple empirical formula quite well. Pronounced deviations of RCHB calculations from the empirical formula can be found in some extremely neutronrich nuclei near the drip line. Such deviations are regarded as signals of the socalled giant halo predicted in Ca [52] and Zr [55].
To give a deeper insight of the deviations from the empirical formula, the differences of neutron radii between the RCHB calculations and the empirical formula are shown in Fig. 18 for the eveneven nuclei with 8\leqslant Z\leqslant 120, in which the smallest ratio R_{n}/N^{1/3} in each isotopic chain is chosen as r_{0} to ensure nonnegative differences.
From Fig. 18, one can see the pronounced deviations in some isotopic chains, such as Ne, Ca, Kr, Nd and Ra, which indicates the possible existence of the giant halos. The deviations show two different behaviors when approaching the neutron drip line. One increases with neutron number monotonically. The other increases with neutron number at the beginning then bends down at next stage. The different behavior is connected to the shell closure and the occupation of levels with high angular momentum.
4.4 Neutron density distributions
To investigate the evolution of the density distribution in each isotopic chain, in Fig. 19, the neutron density distributions of eveneven Ne, Ca, Ni, Kr, Sr and Nd isotopes are shown as examples, where the density distributions of proton and neutron dripline nuclei have been labeled, and the increase of the neutron number is marked.
There are several features to be noted in this figure. Firstly, it is a global trend that the neutron density distributions are extended further with neutron number. The surface expands outward rapidly, and the internal density distribution changes slightly. Secondly, for the Ne, Ca and Kr isotopes, the tails of the density distributions extend with the neutron number monotonically until the neutron drip line. For the Ni, Sr and Nd isotopes, however, the tails of the density distributions reach a maximum at certain neutronrich nucleus, then saturate or even decrease with the neutron number till the neutron drip line. Thirdly, the shell structures influence the density distribution significantly, which can be seen from the dramatic change at the tail of the neutron density distribution when the neutron number N in a nucleus crosses the magic numbers 20, 28, 50, 82 and 126.
4.5 Neutron potential and diffuseness
To examine the isospin dependence of the meanfield potentials, we have investigated the neutron vector plus the scalar potentials V(r)+S(r). The neutron potentials V(r)+S(r) for eveneven Ne, Ca, Ni, Kr, Sr and Nd isotopes from the proton dripline nuclei to the neutron dripline nuclei in the RCHB calculations are shown in Fig. 20. The increase of the neutron number is marked. Generally, the depths of the potentials rise with the neutron number, except for some fluctuation due to the shell structure. At the surface, the potentials extend outward and the diffuseness increases significantly with the neutron number. As a result, the potentials for nuclei near the neutron drip line become highly diffused. Such highly diffused potential will dramatically influence the lowl orbits near the threshold, which may lead to the level crossing and the halos in dripline nuclei.
4.6 Pairing energies
To examine the pairing correlation globally, we investigate the pairing energies E_{\rm pair} for the eveneven nuclei over the nuclear landscape.
In the RCHB calculations, the pairing energy is defined as
\displaystyle E_{\rm pair}=\frac{1}{2}{\rm Tr}(\Delta\kappa).  (31) 
In Fig. 21, the neutron pairing energies E_{\rm pair}^{N} of eveneven nuclei with 8\leqslant Z\leqslant 120 from the RCHB calculations are shown. One can see that the neutron pairing energies are approaching to zero or even vanish for the nuclei near the closed shells N= 8, 20, 28, 50, 82 and 126, and they have maximum values for nuclei in the middle of the shells. Furthermore, the neutron pairing energies vanish at N= 184 and 258, which agrees with the predicted possible magic number N= 184 in Ref. [92]. For the isotopic chains with Z around 10, 20 and 38, the pairing energies for nuclei near the drip line do not vanish, which suggest the disappearance of the neutron magic numbers 28, 50 and 82 in these nuclei.
In Fig. 22, the proton pairing energies E_{\rm pair}^{P} of eveneven nuclei with 8\leqslant Z\leqslant 120 from the RCHB calculations are shown. Similar to the neutron, the proton pairing energies are also approaching to zero or even vanish at the closed shells Z= 8, 20, 28, 50 and 82, and have the maximum value in the middle of the shells. In comparison with the neutron case, the proton pairing energies in each isotopic chain remain almost the same or change modestly with the neutron number. Finally, the total pairing energies of eveneven nuclei with 8\leqslant Z\leqslant 120 from the RCHB calculations are shown in Fig. 23, where the shell structure effects appear again.
4.7 Alpha decay energies
By using the binding energies provided in the RCHB theory with PCPK1, we can extract the Q_{\alpha} values of the bound nuclei with 10\leq Z\leq 120 [96],
\displaystyle Q_{\alpha}=E_{B}(Z2,N2)+E_{B}(2,2)E_{B}(Z,N),  (32) 
where E_{B}(Z,N) is the binding energy for the nucleus with proton number Z and neutron number N. In Fig. 24, the Q_{\alpha} obtained both from the RCHB calculations and the data in AME2012 [6] are shown, where the nuclei observed with \alphadecay radioactivity experimentally [1] are marked with green crosses.
From the RCHB calculations, the Q_{\alpha} values of 3703 nuclei are positive, and they are scaled by colors in Fig. 24(a). Experimentally, as shown in Fig. 24(b), the existing data give 1067 nuclei with their Q_{\alpha} values positive. Comparing the two panels of Fig. 24, we see the following features. Firstly, the positive \alphadecay energies deduced from the RCHB results present a similar pattern with those from available experimental values, although spherical symmetry is assumed. Secondly, in Fig. 24(a), 719 nuclei observed experimentally with the radioactivity of \alphadecay are marked with green crosses. Most of these nuclei are in the heavy and superheavy regions predicted to have Q_{\alpha} values larger than 4 MeV. Thirdly, in the nuclear region with N>184 and Z>92, most nuclei in the upleft corner have Q_{\alpha}> 4 MeV. This indicates a possibility of \alpharadioactivity for these nuclei.
4.8 Proton emitters
The proton radioactivity from neutrondeficient nuclei beyond the proton drip line provides not only spectroscopic information on nuclear structure but also valuable clues on the nuclear force beyond the proton drip line. These proton emitters are also interesting open systems for quantum mechanical tunneling as the emitted proton goes through the barrier of Coulomb and centrifugal potentials.
Based on the RCHB mass table, the properties of (near) spherical proton emitters in the region from I (Z=53) to Bi (Z=83) isotopes are systematically studied [97]. In Fig. 25, the nuclear chart from the I (Z=53) to Bi (Z=83) isotopes is shown, where the experimentally known proton emitters [98] are highlighted.
The proton emitter candidates are located outside the proton drip line with positive Q values. Accordingly, the proton emitter candidates in the calculations satisfy the Fermi surface \lambda_{p}>0, and oneproton separation energy S_{p}<0 or twoproton separation energy S_{2p}<0. The proton emitters observed in this region are well reproduced by the selfconsistent RCHB calculations, except {}^{170,171}Au, {}^{176,177}Tl and {}^{185}Bi. The halflives of the proton emitters are calculated in Ref. [97] by combining the RCHB calculations and WentzelKramersBrillouin (WKB) method, and good agreements with the data are obtained.
5 Summary
In summary, we have performed systematic spherical calculations for all nuclei from Z=8 to Z=120 by using the RCHB theory with the relativistic density functional PCPK1. The calculated binding energies, separation energies, neutron and proton Fermi surfaces, rootmeansquare (rms) radii of neutron, proton, matter, and charge distributions, groundstate spins and parities are tabulated.
With the effects of the continuum included, there are totally 9035 nuclei predicted to be bound from Z=8 to Z=120. Comparing with the calculation without pairing, the RCHB theory predicts more bound nuclei in the neutronrich side and extends the neutron drip line to more neutronrich region. For the nuclei close to the neutron drip line, as the neutron Fermi surface is close to the continuum threshold, pairing correlations could scatter the nucleons from bound states to the continuum, thus provide a significant coupling between the continuum and bound states. As a result, some unbound nuclei predicted without pairing correlations can exist as bound ones. Therefore, the RCHB theory which allows a proper treatment of the continuum and the coupling to the bound states predicts a more extended neutron drip line than the other models. The considerable difference of the drip line between RHB in HO calculations with NL3* and present RCHB calculations are due to different functionals (PCPK1 or NL3*), different pairing forces used, and the treatment of continuum in coordinate or HO space. Future investigations with same density functional and pairing force in both coordinate space and in harmonic oscillator basis are needed. The separation energy evolutions for twoneutron S_{2n}, oneneutron S_{n}, twoproton S_{2p} and oneproton S_{p} are presented and discussed.
The deformation effects on the neutron drip line are investigated with the DRHBc theory by taking Ar isotopes as examples. It is found that the deformation would affect the position of neutron drip line. For Ar isotopic chain, selfconsistent treatment of the deformation and continuum extends the neutron dripline nucleus from {}^{62}Ar in the RCHB to {}^{70}Ar in the DRHBc calculations.
The RCHB calculated charge radii R_{\rm c} with PCPK1 are compared with the experimental values, the rms deviation \sigma of the RCHB calculations from the data is 0.0358 fm. The comparison of the neutron rms radii of eveneven nuclei with 8\leqslant Z\leqslant 120 between the RCHB calculations and the empirical formula are performed. Except for extremely neutronrich nuclei, the systematic trend of the neutron radii in the RCHB calculations follows the simple empirical formula quite well.
The neutron density distributions of eveneven Ne, Ca, Ni, Kr, Sr and Nd isotopes are discussed. For the Ne, Ca and Kr isotopes, the tails of the density distributions extend with the neutron number monotonically until the neutron drip line. For the Ni, Sr and Nd isotopes, the tails of the density distributions reach a maximum at certain neutronrich nucleus, then saturate or even decrease with the neutron number till the neutron drip line.
The neutron potentials of eveneven Ne, Ca, Ni, Kr, Sr and Nd isotopes are investigated. The depths of potentials generally rise with the neutron number, except for some fluctuation due to the shell structure. At the surface, the potentials extend outward and the diffuseness increases significantly with the neutron number.
The pairing energies of the eveneven nuclei over the nuclear landscape are discussed. The pairing energies are approaching to zero or even vanish for the nuclei near the closed shells, and they have maximum values for nuclei in the middle of the shells. The possible magic numbers N= 184 and 258 in superheavy are predicted.
In addition, the \alphadecay energies and proton emitters based on the RCHB calculations are investigated. The \alphadecay energy pattern deduced from the RCHB calculations is similar to the one from available experimental values, and most of the proton emitters observed in the region from I (Z=53) to Bi (Z=83) isotopes are well reproduced by the selfconsistent RCHB calculations.
The successful exploration of the nuclear chart by using the RCHB theory with the relativistic density functional PCPK1 demonstrated the important effects of the continuum on the nuclear landscape. It is quite encouraging to consider the deformation and continuum effects simultaneously in determining the drip line in the future although it is quite timeconsuming and numerically challenging.
References
References
 htt [lgov] Evaluated Nuclear Structure Data File (ENSDF), http://www.nndc.bnl.gov/ensdf.
 Thoennessen [2004] M. Thoennessen, Rep. Prog. Phys. 67 (2004) 1187.
 Lunney et al. [2003] D. Lunney, J. M. Pearson, C. Thibault, Rev. Mod. Phys. 75 (2003) 1021.
 Blaum [2006] K. Blaum, Phys. Rep. 425 (2006) 1.
 Schatz Blaum [2006] H. Schatz K. Blaum, Europhysics News 37 (2006) 16.
 Wang et al. [2012] M. Wang, G. Audi, A. H. Wapstra, F. G. Kondev, M. MacCormick, X. Xu, B. Pfeiffer, Chin. Phys. C 36 (2012) 1603.
 Möller et al. [1995] P. Möller, J. R. Nix, W. Myers, W. Swiatecki, At. Data Nucl. Data Tables 59 (1995) 185.
 Möller et al. [2012] P. Möller, W. D. Myers, H. Sagawa, and S. Yoshida, Phys. Rev. Lett. 108 (2012) 052501.
 Möller et al. [2016] P. Möller, A. Sierk, T. Ichikawa, H. Sagawa, At. Data Nucl. Data Tables 109 (2016) 1.
 Aboussir et al. [1995] Y. Aboussir, J. Pearson, A. Dutta, F. Tondeur, At. Data Nucl. Data Tables 61 (1995) 127.
 Wang et al. [2010] N. Wang, Z. Liang, M. Liu, X. Wu, Phys. Rev. C 82 (2010) 044304.
 Liu et al. [2011] M. Liu, N. Wang, Y. Deng, X. Wu, Phys. Rev. C 84 (2011) 014333.
 Wang et al. [2014] N. Wang, M. Liu, X. Wu, J. Meng, Phys. Lett. B 734 (2014) 215.
 Zhang et al. [2014a] H. Zhang, J. Dong, N. Ma, G. Royer, J. Li, H. Zhang, Nucl. Phys. A 929 (2014a) 38.
 Goriely et al. [2002] S. Goriely, M. Samyn, P.H. Heenen, J. M. Pearson, F. Tondeur, Phys. Rev. C 66 (2002) 024326.
 Goriely et al. [2007] S. Goriely, M. Samyn, and J. Pearson, Phys. Rev. C 75 (2007) 064312.
 Goriely et al. [2009] S. Goriely, S. Hilaire, M. Girod, S. Péru, Phys. Rev. Lett. 102 (2009) 242501.
 Goriely et al. [2010] S. Goriely, N. Chamel, and J. Pearson, Phys. Rev. C 82 (2010) 035804.
 Goriely et al. [2013a] S. Goriely, N. Chamel, and J. Pearson, Phys. Rev. C 88 (2013a) 024308.
 Erler et al. [2012] J. Erler, N. Birge, M. Kortelainen, W. Nazarewicz, E. Olsen, A. M. Perhac, M. Stoitsov, Nature 486 (2012) 509.
 Hirata et al. [1997] D. Hirata, K. Sumiyoshi, I. Tanihata, Y. Sugahara, T. Tachibana, H. Toki, Nucl. Phys. A 616 (1997) 438.
 Lalazissis et al. [1999] G. Lalazissis, S. Raman, and P. Ring, At. Data Nucl. Data Tables 71 (1999) 1.
 Geng et al. [2005] L. S. Geng, H. Toki, J. Meng, Prog. Theor. Phys. 113 (2005) 785.
 Afanasjev et al. [2013] A. V. Afanasjev, S. E. Agbemava, D. Ray, and P. Ring, Phys. Lett. B 726 (2013) 680.
 Agbemava et al. [2014] S. E. Agbemava, A. V. Afanasjev, D. Ray, and P. Ring, Phys. Rev. C 89 (2014) 054320.
 Afanasjev et al. [2015] A. V. Afanasjev, S. E. Agbemava, D. Ray, and P. Ring, Phys. Rev. C 91 (2015) 014324.
 Zhang et al. [2014b] Q.S. Zhang, Z. M. Niu, Z. P. Li, J. M. Yao, J. Meng, Front. Phys. 9 (2014b) 529.
 Lu et al. [2015] K. Q. Lu, Z. X. Li, Z. P. Li, J. M. Yao, J. Meng, Phys. Rev. C 91 (2015) 027304.
 Ring [1996] P. Ring, Prog. Part. Nucl. Phys. 37 (1996) 193.
 Vretenar et al. [2005] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, P. Ring, Phys. Rep. 409 (2005) 101.
 Meng et al. [2006] J. Meng, H. Toki, S.G. Zhou, S. Q. Zhang, W. H. Long, L. S. Geng, Prog. Part. Nucl. Phys. 57 (2006) 470 .
 Nikšić et al. [2011] T. Nikšić, D. Vretenar, P. Ring, Prog. Part. Nucl. Phys. 66 (2011) 519.
 Meng et al. [2013] J. Meng, J. Peng, S. Q. Zhang, P. W. Zhao, Front. Phys. 8 (2013) 55.
 Meng Zhou [2015] J. Meng S.G. Zhou, J. Phys. G: Nucl. Part. Phys. 42 (2015) 093101.
 Sharma et al. [1993] M. Sharma, G. Lalazissis, P. Ring, Phys. Lett. B 317 (1993) 9.
 Liang et al. [2015] H. Z. Liang, J. Meng, S.G. Zhou, Phys. Rep. 570 (2015) 1 .
 Zhou et al. [2003] S.G. Zhou, J. Meng, P. Ring, Phys. Rev. Lett. 91 (2003) 262501.
 Liang et al. [2010] H. Liang, W. H. Long, J. Meng, N. Van Giai, Eur. Phys. J. A 44 (2010) 119.
 Koepf Ring [1989] W. Koepf P. Ring, Nucl. Phys. A 493 (1989) 61.
 Yao et al. [2006] J. M. Yao, H. Chen, J. Meng, Phys. Rev. C 74 (2006) 024307.
 Arima [2011] A. Arima, Sci. China Phys., Mech. Astron. 54 (2011) 188.
 Li et al. [2011a] J. Li, J. Meng, P. Ring, J. M. Yao, A. Arima, Sci. China Phys., Mech. Astron. 54 (2011a) 204.
 Li et al. [2011b] J. Li, J. M. Yao, J. Meng, A. Arima, Prog. Theor. Phys. 125 (2011b) 1185.
 Wei et al. [2012] J. X. Wei, J. Li, J. Meng, Prog. Theor. Phys. Suppl. 196 (2012) 400.
 König Ring [1993] J. König P. Ring, Phys. Rev. Lett. 71 (1993) 3079.
 Afanasjev et al. [2000] A. V. Afanasjev, P. Ring, and J. König, Nucl. Phys. A 676 (2000) 196.
 Afanasjev Ring [2000] A. V. Afanasjev P. Ring, Phys. Rev. C 62 (2000) 031302 .
 Afanasjev Abusara [2010] A. V. Afanasjev H. Abusara, Phys. Rev. C 82 (2010) 034329 .
 Zhao et al. [2011a] P. W. Zhao, J. Peng, H. Z. Liang, P. Ring, J. Meng, Phys. Rev. Lett. 107 (2011a) 122501.
 Zhao et al. [2011b] P. W. Zhao, S. Q. Zhang, J. Peng, H. Z. Liang, P. Ring, J. Meng, Phys. Lett. B 699 (2011b) 181 .
 Zhao et al. [2012a] P. W. Zhao, J. Peng, H. Z. Liang, P. Ring, J. Meng, Phys. Rev. C 85 (2012a) 054310.
 Meng et al. [2002a] J. Meng, H. Toki, J. Y. Zeng, S. Q. Zhang, S.G. Zhou, Phys. Rev. C 65 (2002a) 041302.
 Meng Ring [1996] J. Meng P. Ring, Phys. Rev. Lett. 77 (1996) 3963.
 Meng [1998] J. Meng, Nucl. Phys. A 635 (1998) 3.
 Meng Ring [1998] J. Meng P. Ring, Phys. Rev. Lett. 80 (1998) 460.
 Zhang et al. [2003] S. Q. Zhang, J. Meng, S.G. Zhou, Science in China (Series G) 46 (2003) 632.
 Meng et al. [1998] J. Meng, I. Tanihata, and S. Yamaji, Phys. Lett. B 419 (1998) 1 .
 Meng et al. [2002b] J. Meng, S.G. Zhou, and I. Tanihata, Phys. Lett. B 532 (2002b) 209.
 Zhou et al. [2010] S.G. Zhou, J. Meng, P. Ring, E.G. Zhao, Phys. Rev. C 82 (2010) 011301.
 Li et al. [2012a] L. L. Li, J. Meng, P. Ring, E.G. Zhao, S.G. Zhou, Phys. Rev. C 85 (2012a) 024312.
 Chen et al. [2012a] Y. Chen, L. L. Li, H. Z. Liang, J. Meng, Phys. Rev. C 85 (2012a) 067301.
 Li et al. [2012b] L. L. Li, J. Meng, P. Ring, E.G. Zhao, S.G. Zhou, Chin. Phys. Lett. 29 (2012b) 042101.
 Zhao et al. [2010] P. W. Zhao, Z. P. Li, J. M. Yao, J. Meng, Phys. Rev. C 82 (2010) 054319.
 Zhao et al. [2012b] P. W. Zhao, L. S. Song, B. Sun, H. Geissel, J. Meng, Phys. Rev. C 86 (2012b) 064324.
 Yao et al. [2013] J. M. Yao, H. Mei, Z. P. Li, Phys. Lett. B 723 (2013) 459.
 Yao et al. [2014] J. M. Yao, K. Hagino, Z. P. Li, J. Meng, P. Ring, Phys. Rev. C 89 (2014) 054306.
 Wu et al. [2014] X. Y. Wu, J. M. Yao, Z. P. Li, Phys. Rev. C 89 (2014) 017304.
 Li et al. [2012c] Z. P. Li, C. Y. Li, J. Xiang, J. M. Yao, J. Meng, Phys. Lett. B 717 (2012c) 470.
 Fu et al. [2013] Y. Fu, H. Mei, J. Xiang, Z. P. Li, J. M. Yao, J. Meng, Phys. Rev. C 87 (2013) 054305.
 Li et al. [2013a] Z. P. Li, B. Y. Song, J. M. Yao, D. Vretenar, J. Meng, Phys. Lett. B 726 (2013a) 866.
 Xiang et al. [2013] J. Xiang, Z. P. Li, J. M. Yao, W. H. Long, P. Ring, J. Meng, Phys. Rev. C 88 (2013) 057301.
 Wang et al. [2015] Z. H. Wang, J. Xiang, W. H. Long, Z. P. Li, J. Phys. G: Nucl. Part. Phys. 42 (2015) 045108.
 Qu et al. [2013] X. Y. Qu, Y. Chen, S. Q. Zhang, P. W. Zhao, I. J. Shin, Y. Lim, Y. Kim, J. Meng, Sci. China Phys. Mech. Astron. 56 (2013) 2031.
 Nikolaus et al. [1992] B. Nikolaus, T. Hoch, and D. Madland, Phys. Rev. C 46 (1992) 1757.
 Bürvenich et al. [2002] T. Bürvenich, D. G. Madland, J. A. Maruhn, P.G. Reinhard, Phys. Rev. C 65 (2002) 044308.
 Berger et al. [1984] J. Berger, M. Girod, and D. Gogny, Nucl. Phys. A 428 (1984) 23.
 Tian et al. [2009] Y. Tian, Z. Y. Ma, P. Ring, Phys. Lett. B 676 (2009) 44.
 Ring Schuck [2004] P. Ring P. Schuck, The nuclear manybody problem (Springer Science & Business Media, 2004).
 Chen et al. [2012b] L. Chen, W. Plass, H. Geissel, R. Knöbel, C. Kozhuharov, Y. A. Litvinov, Z. Patyk, C. Scheidenberger, K. SiegieńIwaniuk, B. Sun, et al., Nucl. Phys. A 882 (2012b) 71.
 Steppenbeck et al. [2012] D. Steppenbeck, R. V. F. Janssens, S. J. Freeman, M. P. Carpenter, P. Chowdhury, A. N. Deacon, M. Honma, H. Jin, T. Lauritsen, C. J. Lister, J. Meng, J. Peng, D. Seweryniak, J. F. Smith, Y. Sun, S. L. Tabor, B. J. Varley, Y.C. Yang, S. Q. Zhang, P. W. Zhao, S. Zhu, Phys. Rev. C 85 (2012) 044316.
 Yu et al. [2012] L. F. Yu, P. W. Zhao, S. Q. Zhang, P. Ring, J. Meng, Phys. Rev. C 85 (2012) 024318.
 Li et al. [2013b] J. Li, C. Y. He, Y. Zheng, C. B. Li, K. Y. Ma, J. B. Lu, Phys. Rev. C 88 (2013b) 014317.
 Zhao et al. [2015] P. W. Zhao, S. Q. Zhang, and J. Meng, Phys. Rev. C 92 (2015) 034319.
 Li et al. [2012d] X. W. Li, J. Li, J. B. Lu, K. Y. Ma, Y. H. Wu, L. H. Zhu, C. Y. He, X. Q. Li, Y. Zheng, G. S. Li, X. G. Wu, Y. J. Ma, Y. Z. Liu, Phys. Rev. C 86 (2012d) 057305.
 Zhang et al. [2014c] P. Zhang, B. Qi, S. Y. Wang, Phys. Rev. C 89 (2014c) 047302.
 Peng Zhao [2015] J. Peng P. W. Zhao, Phys. Rev. C 91 (2015) 044329 .
 Kuti et al. [2014] I. Kuti, Q. B. Chen, J. Timár, D. Sohler, S. Q. Zhang, Z. H. Zhang, P. W. Zhao, J. Meng, K. Starosta, T. Koike, E. S. Paul, D. B. Fossan, C. Vaman, Phys. Rev. Lett. 113 (2014) 032501 .
 Lu et al. [2012] B.N. Lu, E.G. Zhao, S.G. Zhou, Phys. Rev. C 85 (2012) 011301.
 Agbemava et al. [2015] S. E. Agbemava, A. V. Afanasjev, T. Nakatsukasa, P. Ring, Phys. Rev. C 92 (2015) 054310.
 Angeli Marinova [2013] I. Angeli K. Marinova, At. Data Nucl. Data Tables 99 (2013) 69.
 Nikšić et al. [2010] T. Nikšić, P. Ring, D. Vretenar, Y. Tian, Z. Y. Ma, Phys. Rev. C 81 (2010) 054318.
 Zhang et al. [2005] W. Zhang, J. Meng, S. Q. Zhang, L. S. Geng, H. Toki, Nucl. Phys. A 753 (2005) 106.
 Geng et al. [2006] L. S. Geng, J. Meng, H. Toki, W. H. Long, G. Shen, Chin. Phys. Lett. 23 (2006) 1139.
 Goriely et al. [2013b] S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. C 88 (2013b) 024308.
 Wang Li [2013] N. Wang T. Li, Phys. Rev. C 88 (2013) 011301.
 Zhang Xia [2016] L. F. Zhang X. W. Xia, Chin. Phys. C 40 (2016) 54102 .
 Lim et al. [2016] Y. Lim, X. W. Xia, Y. Kim, Phys. Rev. C 93 (2016) 014314.
 Blank Borge [2008] B. Blank M. Borge, Prog. Part. Nucl. Phys. 60 (2008) 403..
Table 1. Groundstate properties of nuclei calculated by RCHB theory with relativistic density functional PCPK1
Z  Proton number 
N  Neutron number 
A  Mass number 
E^{\mathrm{Cal.}}_{\mathrm{b}}  Binding energy from RCHB calculations 
E^{\mathrm{Exp.}}_{\mathrm{b}}  Binding energy from experimental data 
E^{\mathrm{Cal.}}_{\mathrm{b}}/A  Per nucleon binding energy from RCHB calculations 
E^{\mathrm{Exp.}}_{\mathrm{b}}/A  Per nucleon binding energy from experimental data 
S_{2n}  Twoneutron separation energy 
S_{2p}  Twoproton separation energy 
S_{n}  Oneneutron separation energy 
S_{p}  Oneproton separation energy 
\lambda_{n}  Neutron Fermi surface 
\lambda_{p}  Proton Fermi surface 
R_{m}  Matter rootmeansquare radius 
R_{n}  Neutron rootmeansquare radius 
R_{p}  Proton rootmeansquare radius 
R^{\mathrm{Cal.}}_{\mathrm{c}}  Charge radius from RCHB calculations 
R^{\mathrm{Exp.}}_{\mathrm{c}}  Charge radius from experimental data 
j^{\pi}(P)  Ground state spin j and parity \pi of proton 
j^{\pi}(N)  Ground state spin j and parity \pi of neutron 
\sigma  rms deviations of binding energies and charge radii for each isotopic chain. 
A  N  E^{\mathrm{Cal.}}_{\mathrm{b}}  E^{\mathrm{Exp.}}_{\mathrm{b}}  E^{\mathrm{Cal.}}_{\mathrm{b}}/A  E^{\mathrm{Exp.}}_{\mathrm{b}}/A  S_{2n}  S_{2p}  S_{n}  S_{p}  \lambda_{n}  \lambda_{p}  R_{m}  R_{n}  R_{p}  R^{\mathrm{Cal.}}_{\mathrm{c}}  R^{\mathrm{Exp.}}_{\mathrm{c}}  j^{\pi}(P)  j^{\pi}(N) 

(\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm MeV)  (\rm fm)  (\rm fm)  (\rm fm)  (\rm fm)  (\rm fm)  
Z=8~{}~{}(\rm{O})  
12  4  60.44  58.68  5.04  4.88  19.45  0.44  2.738  2.335  2.919  3.027  0^{+}  0^{+}  
13  5  77.81  75.56  5.99  5.81  17.38  17.01  0.99  2.619  2.373  2.762  2.876  0^{+}  3/2^{}  
14  6  101.86  98.73  7.28  7.05  41.42  24.05  19.96  1.34  2.465  2.238  2.621  2.741  0^{+}  0^{+}  
15  7  112.81  111.96  7.52  7.46  34.99  10.95  15.48  5.43  2.599  2.537  2.652  2.770  0^{+}  1/2^{}  
16  8  127.29  127.62  7.96  7.98  25.43  14.48  11.64  7.83  2.638  2.626  2.649  2.768  2.701  0^{+}  0^{+}  
17  9  132.48  131.76  7.79  7.75  19.67  5.19  13.49  9.79  2.690  2.733  2.642  2.760  2.695  0^{+}  5/2^{+}  
18  10  141.63  139.81  7.87  7.77  14.34  9.15  6.94  11.79  2.736  2.807  2.644  2.763  2.775  0^{+}  0^{+}  
19  11  145.27  143.76  7.65  7.57  12.79  3.64  6.18  13.50  2.792  2.896  2.642  2.760  0^{+}  5/2^{+}  
20  12  153.16  151.37  7.66  7.57  11.54  7.89  5.76  15.10  2.842  2.964  2.649  2.767  0^{+}  0^{+}  
21  13  156.04  155.18  7.43  7.39  10.77  2.88  5.74  16.16  2.935  3.096  2.653  2.771  0^{+}  1/2^{+}  
22  14  162.91  162.03  7.41  7.36  9.75  6.87  4.79  18.54  2.955  3.111  2.659  2.777  0^{+}  0^{+}  
23  15  166.49  164.77  7.24  7.16  10.45  3.58  3.63  19.61  3.001  3.173  2.650  2.768  0^{+}  1/2^{+}  
24  16  170.90  168.95  7.12  7.02  7.99  4.41  3.39  20.52  3.082  3.268  2.672  2.789  0^{+}  0^{+}  
25  17  172.02  168.18  6.88  6.73  5.53  1.12  3.41  23.79  3.169  3.370  2.692  2.809  0^{+}  3/2^{+}  
26  18  175.10  168.86  6.73  6.49  4.20  3.08  2.14  22.36  3.230  3.428  2.732  2.847  0^{+}  0^{+}  
27  19  175.57  6.50  3.56  0.47  0.89  23.36  3.310  3.516  2.760  2.874  0^{+}  3/2^{+}  
28  20  178.02  6.36  2.92  2.45  0.89  23.36  3.370  3.576  2.790  2.903  0^{+}  0^{+}  
\sigma  2.55  0.054  
Z=9~{}~{}(\rm{F})  
17  8  128.76  128.22  7.57  7.54  28.40  16.64  1.47  13.36  8.41  2.704  2.623  2.774  2.887  5/2^{+}  0^{+}  
18  9  135.84  137.37  7.55  7.63  23.72  7.08  3.36  12.04  8.74  2.734  2.720  2.749  2.863  5/2^{+}  5/2^{+}  
19  10  146.83  147.80  7.73  7.78  18.07  10.99  5.20  8.65  11.89  2.766  2.789  2.741  2.855  2.898  5/2^{+}  0^{+}  
20  11  152.29  154.40  7.61  7.72  16.45  5.46  7.02  7.85  12.81  2.805  2.866  2.729  2.844  5/2^{+}  5/2^{+}  
21  12  161.83  162.50  7.71  7.74  15.00  9.54  8.67  7.33  14.38  2.845  2.928  2.730  2.844  5/2^{+}  0^{+}  
22  13  166.10  167.73  7.55  7.62  13.81  4.27  10.05  6.28  15.07  2.888  2.997  2.723  2.838  5/2^{+}  5/2^{+}  
23  14  174.70  175.29  7.60  7.62  12.87  8.60  11.79  6.12  17.77  2.933  3.057  2.729  2.844  5/2^{+}  0^{+}  
24  15  179.70  179.11  7.49  7.46  13.60  5.00  13.21  4.87  22.62  2.978  3.122  2.720  2.836  5/2^{+}  1/2^{+}  
25  16  184.97  183.38  7.40  7.34  10.28  5.28  14.07  4.39  21.13  3.051  3.213  2.742  2.856  5/2^{+}  0^{+}  
26  17  186.99  184.15  7.19  7.10  7.30  2.02  14.97  4.37  23.63  3.130  3.307  2.766  2.879  5/2^{+}  3/2^{+}  
27  18  191.04  186.25  7.08  6.89  6.06  4.05  15.94  3.09  22.83  3.192  3.369  2.805  2.917  5/2^{+}  0^{+}  
28  19  192.47  186.03  6.87  6.64  5.48  1.43  16.90  1.61  23.41  3.265  3.450  2.834  2.945  5/2^{+}  3/2^{+}  
29  20  195.92  6.76  4.88  3.45  17.90  1.65  24.73  3.321  3.507  2.865  2.975  5/2^{+}  0^{+}  
30  21  195.25  6.51  2.78  0.67  18.02  1.67  24.39  3.654  3.943  2.867  2.977  5/2^{+}  3/2^{}  
31  22  196.42  6.34  0.50  1.17  18.79  0.49  25.40  3.487  3.704  2.890  2.998  5/2^{+}  0^{+}  
32  23  195.57  6.11  0.32  0.85  0.39  25.31  3.810  4.114  2.893  3.002  5/2^{+}  3/2^{}  
33  24  196.34  5.95  0.08  0.77  0.27  25.61  3.662  3.907  2.911  3.019  5/2^{+}  0^{+}  
\sigma  2.70  0.043  
Z=10~{}~{}(\rm{Ne})  
16  6  100.87  97.33  6.30  6.08  0.99  0.51  22.43  0.37  2.735  2.280  2.975  3.081  0^{+}  0^{+}  
17  7  115.25  112.89  6.78  6.64  2.45  14.38  3.13  21.08  1.02  2.757  2.553  2.891  2.999  3.043  0^{+}  1/2^{}  
18  8  133.74  132.14  7.43  7.34  32.87  6.45  18.48  4.98  15.60  2.73  2.754  2.630  2.850  2.960  2.971  0^{+}  0^{+} 
19  9  142.63  143.78  7.51  7.57  27.38  10.15  8.90  6.79  16.73  4.40  2.774  2.718  2.823  2.934  3.008  0^{+}  5/2^{+} 
20  10  155.45  160.64  7.77  8.03  21.71  13.82  12.82  8.62  10.32  6.05  2.798  2.784  2.811  2.923  3.005  0^{+}  0^{+} 
21  11  162.69  167.41  7.75  7.97  20.06  17.42  7.24  10.40  9.49  7.75  2.824  2.851  2.795  2.907  2.970  0^{+}  5/2^{+} 
22  12  173.93  177.77  7.91  8.08  18.48  20.77  11.24  12.10  8.91  9.35  2.856  2.908  2.791  2.904  2.953  0^{+}  0^{+} 
23  13  180.01  182.97  7.83  7.96  17.32  23.97  6.08  13.91  7.65  11.03  2.883  2.961  2.777  2.890  2.910  0^{+}  5/2^{+} 
24  14  190.09  191.84  7.92  7.99  16.16  27.17  10.08  15.39  7.46  12.48  2.922  3.019  2.781  2.894  2.901  0^{+}  0^{+} 
25  15  196.36  196.02  7.85  7.84  16.35  29.87  6.27  16.66  6.74  13.83  2.966  3.086  2.777  2.890  2.932  0^{+}  1/2^{+} 
26  16  202.63  201.55  7.79  7.75  12.54  31.73  6.27  17.65  5.41  14.75  3.035  3.174  2.799  2.911  2.925  0^{+}  0^{+} 
27  17  205.67  203.07  7.62  7.52  9.31  33.65  3.04  18.68  5.39  15.78  3.107  3.262  2.825  2.936  0^{+}  3/2^{+}  
28  18  210.78  206.88  7.53  7.39  8.15  35.68  5.11  19.74  4.09  16.81  3.169  3.327  2.863  2.973  2.964  0^{+}  0^{+} 
29  19  213.25  207.84  7.35  7.18  7.58  37.67  2.47  20.78  2.69  17.82  3.235  3.402  2.892  3.001  0^{+}  3/2^{+}  
30  20  217.80  211.28  7.26  7.04  7.02  39.77  4.55  21.87  2.49  18.83  3.290  3.459  2.923  3.031  0^{+}  0^{+}  
31  21  217.26  211.57  7.01  6.82  4.01  40.03  0.54  22.01  2.56  19.00  3.568  3.836  2.926  3.033  0^{+}  3/2^{}  
32  22  219.51  6.86  1.72  41.88  2.26  23.09  1.05  20.25  3.423  3.617  2.951  3.058  0^{+}  0^{+}  
33  23  218.86  6.63  1.60  0.65  23.29  0.98  20.48  3.680  3.954  2.955  3.061  0^{+}  3/2^{}  
34  24  220.49  6.49  0.98  1.63  24.16  0.78  21.46  3.565  3.784  2.975  3.081  0^{+}  0^{+}  
35  25  219.73  6.28  0.86  0.77  0.69  21.74  3.810  4.095  2.981  3.086  0^{+}  3/2^{}  
36  26  221.17  6.14  0.68  1.45  0.62  22.56  3.703  3.941  2.998  3.103  0^{+}  0^{+}  
37  27  220.38  5.96  0.65  0.79  0.60  22.75  3.962  4.264  3.001  3.106  0^{+}  1/2^{}  
38  28  221.66  5.83  0.49  1.28  0.49  23.59  3.832  4.083  3.019  3.123  0^{+}  0^{+}  
39  29  220.86  5.66  0.48  0.80  0.41  23.86  4.060  4.360  3.024  3.128  0^{+}  1/2^{}  
40  30  221.97  5.55  0.31  1.11  0.34  24.59  3.950  4.209  3.041  3.144  0^{+}  0^{+}  
41  31  221.02  5.39  0.16  0.95  0.20  24.91  4.184  4.489  3.047  3.151  0^{+}  1/2^{}  
42  32  222.04  5.29  0.06  1.02  0.12  25.58  4.057  4.320  3.063  3.166  0^{+}  0^{+}  
43  33  220.79  5.13  0.23  1.25  0.03  25.80  4.365  4.688  3.069  3.172  0^{+}  1/2^{}  
44  34  221.59  5.04  0.45  0.79  0.26  26.50  4.164  4.431  3.089  3.191  0^{+}  0^{+}  
\sigma  3.56  0.046  
Z=11~{}~{}(\rm{Na})  
17  6  98.28  5.78  2.08  2.59  23.41  0.74  3.851  2.284  4.480  4.551  1/2^{+}  0^{+}  
18  7  113.49  111.64  6.30  6.25  1.36  15.21  1.77  17.86  0.61  3.139  2.557  3.458  3.550  1/2^{+}  1/2^{}  
19  8  133.04  131.82  7.00  6.94  34.76  4.28  19.55  0.70  16.92  1.71  2.830  2.633  2.966  3.072  5/2^{+}  0^{+}  
20  9  143.69  145.97  7.18  7.30  30.20  7.85  10.65  1.05  16.07  3.28  2.826  2.713  2.915  3.023  2.972  5/2^{+}  5/2^{+} 
21  10  158.26  163.08  7.54  7.77  25.22  11.43  14.57  2.81  12.00  4.88  2.836  2.775  2.890  2.999  3.014  5/2^{+}  0^{+} 
22  11  167.31  174.15  7.61  7.92  23.63  15.02  9.05  4.62  11.19  6.56  2.848  2.833  2.863  2.972  2.985  5/2^{+}  5/2^{+} 
23  12  180.26  186.56  7.84  8.11  22.00  18.43  12.95  6.33  10.57  8.18  2.868  2.884  2.850  2.960  2.994  5/2^{+}  0^{+} 
24  13  188.50  193.52  7.85  8.06  21.19  22.41  8.24  8.49  9.08  10.03  2.873  2.917  2.820  2.932  2.974  5/2^{+}  5/2^{+} 
25  14  200.03  202.53  8.00  8.10  19.77  25.33  11.53  9.94  8.87  11.53  2.904  2.970  2.819  2.931  2.977  5/2^{+}  0^{+} 
26  15  207.26  208.11  7.97  8.00  18.76  27.56  7.23  10.90  7.26  12.69  2.957  3.050  2.825  2.936  2.993  5/2^{+}  1/2^{+} 
27  16  214.37  214.84  7.94  7.96  14.34  29.39  7.11  11.74  6.38  13.61  3.022  3.136  2.848  2.958  3.013  5/2^{+}  0^{+} 
28  17  218.37  218.38  7.80  7.80  11.11  31.37  4.00  12.70  6.32  14.63  3.090  3.221  2.875  2.985  3.040  5/2^{+}  3/2^{+} 
29  18  224.49  222.78  7.74  7.68  10.13  33.45  6.13  13.71  5.10  15.68  3.150  3.287  2.913  3.021  3.092  5/2^{+}  0^{+} 
30  19  227.96  225.06  7.60  7.51  9.59  35.49  3.47  14.71  3.37  16.70  3.212  3.358  2.941  3.048  3.118  5/2^{+}  3/2^{+} 
31  20  233.54  229.34  7.53  7.39  9.05  37.62  5.58  15.75  3.43  17.73  3.264  3.414  2.972  3.078  3.170  5/2^{+}  0^{+} 
32  21  233.16  230.87  7.29  7.21  5.20  37.91  0.38  15.91  3.53  17.95  3.488  3.728  2.976  3.081  5/2^{+}  3/2^{}  
33  22  236.73  7.17  3.19  40.31  3.57  17.22  1.74  19.34  3.377  3.551  3.001  3.106  5/2^{+}  0^{+}  
34  23  236.32  6.95  3.16  40.76  0.40  17.46  1.69  19.62  3.574  3.816  3.006  3.110  5/2^{+}  3/2^{}  
35  24  238.96  6.83  2.23  42.62  2.63  18.46  1.39  20.71  3.498  3.694  3.027  3.131  5/2^{+}  0^{+}  
36  25  238.50  6.62  2.17  0.46  18.77  1.32  21.06  3.673  3.921  3.033  3.137  5/2^{+}  3/2^{}  
37  26  240.76  6.51  1.81  2.27  19.59  1.19  21.93  3.619  3.834  3.050  3.154  5/2^{+}  0^{+}  
38  27  240.19  6.32  1.70  0.57  19.81  1.08  22.33  3.782  4.040  3.058  3.161  5/2^{+}  3/2^{}  
39  28  242.32  6.21  1.56  2.13  20.66  1.03  23.04  3.736  3.967  3.072  3.175  5/2^{+}  0^{+}  
40  29  241.80  6.04  1.60  0.53  20.94  0.98  23.36  3.900  4.170  3.077  3.180  5/2^{+}  1/2^{}  
41  30  243.69  5.94  1.37  1.89  21.72  0.86  24.09  3.847  4.089  3.093  3.194  5/2^{+}  0^{+}  
42  31  243.17  5.79  1.38  0.52  22.15  0.73  24.50  3.990  4.262  3.099  3.201  5/2^{+}  1/2^{}  
43  32  244.81  5.69  1.12  1.64  22.78  0.63  25.10  3.953  4.203  3.113  3.214  5/2^{+}  0^{+}  
44  33  243.91  5.54  0.74  0.90  23.12  0.36  25.46  4.160  4.453  3.120  3.221  5/2^{+}  1/2^{}  
45  34  245.37  5.45  0.56  1.46  23.79  0.15  26.04  4.052  4.307  3.136  3.236  5/2^{+}  0^{+}  
46  35  243.97  5.30  0.06  1.40  0.03  26.13  4.412  4.742  3.141  3.241  5/2^{+}  1/2^{}  
47  36  244.66  5.21  0.71  0.69  0.34  26.72  4.174  4.436  3.171  3.270  5/2^{+}  0^{+}  
\sigma  3.86  0.055  
Z=12~{}~{}(\rm{Mg})  
19  7  114.54  112.14  6.03  5.84  0.71  1.05  23.65  0.36  2.935  2.573  3.127  3.228  0^{+}  1/2^{}  
20  8  136.23  134.47  6.81  6.72  2.49  21.69  3.19  18.59  1.01  2.886  2.643  3.037  3.141  0^{+}  0^{+}  
21  9  148.42  149.20  7.07  7.10  33.88  5.79  12.19  4.74  19.93  2.47  2.874  2.719  2.985  3.090  0^{+}  5/2^{+}  
22  10  164.63  168.58  7.48  7.66  28.40  9.18  16.21  6.37  13.56  3.97  2.875  2.777  2.954  3.061  0^{+}  0^{+}  
23  11  175.35  181.73  7.62  7.90  26.93  12.66  10.72  8.04  12.79  5.59  2.876  2.828  2.920  3.028  0^{+}  5/2^{+}  
24  12  189.98  198.26  7.92  8.26  25.35  16.05  14.62  9.72  12.18  7.18  2.887  2.873  2.901  3.010  3.057  0^{+}  0^{+} 
25  13  200.11  205.59  8.00  8.22  24.76  20.10  10.13  11.61  10.49  9.03  2.881  2.897  2.864  2.974  3.029  0^{+}  5/2^{+} 
26  14  213.56  216.68  8.21  8.33  23.59  23.48  13.45  13.53  10.30  10.77  2.890  2.927  2.845  2.955  3.034  0^{+}  0^{+} 
27  15  221.50  223.12  8.20  8.26  21.39  25.14  7.93  14.24  9.68  11.69  2.955  3.024  2.867  2.976  0^{+}  1/2^{+}  
28  16  229.57  231.63  8.20  8.27  16.01  26.94  8.07  15.20  7.41  12.61  3.019  3.110  2.893  3.002  0^{+}  0^{+}  
29  17  234.63  235.28  8.09  8.11  13.14  28.97  5.07  16.27  7.37  13.64  3.083  3.193  2.921  3.028  0^{+}  3/2^{+}  
30  18  241.87  241.64  8.06  8.06  12.30  31.09  7.23  17.37  6.15  14.71  3.141  3.258  2.957  3.063  0^{+}  0^{+}  
31  19  246.40  243.95  7.95  7.87  11.77  33.16  4.54  18.44  4.81  15.75  3.199  3.327  2.985  3.091  0^{+}  3/2^{+}  
32  20  253.06  249.72  7.91  7.81  11.19  35.26  6.66  19.52  4.42  16.79  3.249  3.382  3.015  3.119  0^{+}  0^{+}  
33  21  253.45  252.00  7.68  7.64  7.05  36.19  0.39  20.29  4.63  17.77  3.303  3.450  3.029  3.133  0^{+}  7/2^{}  
34  22  257.96  256.71  7.59  7.54  4.90  38.44  4.51  21.23  2.51  18.55  3.349  3.504  3.043  3.147  0^{+}  0^{+}  
35  23  257.86  257.47  7.37  7.36  4.41  39.00  0.10  21.53  2.48  18.88  3.502  3.717  3.049  3.152  0^{+}  3/2^{}  
36  24  261.63  260.80  7.27  7.24  3.67  41.14  3.77  22.67  2.06  20.07  3.454  3.631  3.070  3.173  0^{+}  0^{+}  
37  25  261.55  7.07  3.70  41.83  0.08  23.05  2.00  20.48  3.585  3.805  3.077  3.179  0^{+}  3/2^{}  
38  26  264.70  6.97  3.07  43.53  3.15  23.94  1.78  21.41  3.562  3.758  3.094  3.196  0^{+}  0^{+}  
39  27  264.58  6.78  3.03  44.20  0.12  24.39  1.69  21.88  3.676  3.904  3.102  3.203  0^{+}  3/2^{}  
40  28  267.42  6.69  2.72  45.76  2.84  25.09  1.57  22.60  3.669  3.882  3.116  3.217  0^{+}  0^{+}  
41  29  267.22  6.52  2.64  46.36  0.20  25.42  1.56  22.98  3.793  4.038  3.121  3.222  0^{+}  1/2^{}  
42  30  269.87  6.43  2.45  47.90  2.65  26.18  1.37  23.69  3.774  4.001  3.135  3.236  0^{+}  0^{+}  
43  31  269.84  6.28  2.62  48.82  0.03  26.67  1.29  24.15  3.873  4.122  3.141  3.241  0^{+}  1/2^{}  
44  32  272.04  6.18  2.17  50.00  2.20  27.22  1.12  24.71  3.877  4.116  3.154  3.253  0^{+}  0^{+}  
45  33  271.72  6.04  1.88  50.93  0.32  27.81  0.68  25.22  3.984  4.244  3.161  3.261  0^{+}  1/2^{}  
46  34  273.60  5.95  1.56  1.87  0.56  25.64  3.978  4.225  3.174  3.273  0^{+}  0^{+}  
47  35  272.80  5.80  1.08  0.80  0.57  25.65  4.300  4.623  3.173  3.273  0^{+}  1/2^{+}  
48  36  273.55  5.70  0.05  0.75  0.04  26.36  4.076  4.325  3.214  3.312  0^{+}  0^{+}  
\sigma  3.51  0.062  
Z=13~{}~{}(\rm{Al})  
20  7  112.14  5.61  1.35  2.40  18.81  0.73  3.645  2.577  4.107  4.184  1/2^{+}  1/2^{}  
21  8  134.70  6.00  1.66  22.56  1.53  19.28  0.62  3.079  2.648  3.317  3.412  1/2^{+}  0^{+}  
22  9  147.87  6.72  35.73  4.18  13.71  0.55  18.08  2.18  2.991  2.724  3.163  3.263  1/2^{+}  5/2^{+}  
23  10  165.51  168.72  7.20  7.33  30.81  7.25  17.64  0.88  15.11  2.53  2.919  2.771  3.027  3.131  5/2^{+}  0^{+}  
24  11  178.22  183.59  7.43  7.65  30.35  10.91  12.71  2.86  14.55  3.85  2.897  2.809  2.970  3.076  5/2^{+}  5/2^{+}  
25  12  194.86  200.53  7.79  8.02  29.35  14.60  16.64  4.88  14.04  5.20  2.892  2.844  2.935  3.042  5/2^{+}  0^{+}  
26  13  207.00  211.89  7.96  8.15  28.78  18.50  12.14  6.89  11.99  6.69  2.880  2.866  2.894  3.003  5/2^{+}  5/2^{+}  
27  14  222.44  224.95  8.24  8.33  27.57  22.41  15.44  8.87  11.68  8.29  2.878  2.888  2.866  2.976  3.061  5/2^{+}  0^{+} 
28  15  231.04  232.68  8.25  8.31  24.04  23.78  8.60  9.54  9.31  9.39  2.944  2.989  2.892  3.001  5/2^{+}  1/2^{+}  
29  16  239.85  242.11  8.27  8.35  17.41  25.48  8.81  10.28  8.29  10.52  3.006  3.074  2.920  3.028  5/2^{+}  0^{+}  
30  17  245.81  247.84  8.19  8.26  14.77  27.44  5.96  11.18  8.19  11.67  3.067  3.156  2.947  3.054  5/2^{+}  3/2^{+}  
31  18  254.05  254.99  8.20  8.23  14.20  29.55  8.24  12.18  7.16  12.75  3.123  3.220  2.983  3.089  5/2^{+}  0^{+}  
32  19  259.55  259.21  8.11  8.10  13.74  31.59  5.51  13.15  5.23  13.82  3.178  3.287  3.012  3.116  5/2^{+}  3/2^{+}  
33  20  267.21  264.65  8.10  8.02  13.16  33.67  7.65  14.15  5.43  14.86  3.226  3.342  3.041  3.144  5/2^{+}  0^{+}  
34  21  268.61  267.32  7.90  7.86  9.05  35.44  1.40  15.16  5.28  15.85  3.273  3.402  3.053  3.156  5/2^{+}  7/2^{}  
35  22  273.90  272.55  7.83  7.78  6.70  37.18  5.30  15.95  3.63  16.73  3.315  3.453  3.067  3.169  5/2^{+}  0^{+}  
36  23  274.45  274.45  7.62  7.63  5.84  38.12  0.54  16.59  2.98  17.64  3.363  3.514  3.079  3.182  5/2^{+}  7/2^{}  
37  24  279.08  278.66  7.54  7.53  5.18  40.12  4.64  17.45  2.79  18.40  3.409  3.568  3.093  3.093  5/2^{+}  0^{+}  
38  25  279.36  280.33  7.35  7.38  4.91  40.86  0.28  17.81  2.78  18.87  3.515  3.713  3.101  3.203  5/2^{+}  3/2^{}  
39  26  283.47  7.27  4.39  42.71  4.11  18.77  2.43  19.88  3.506  3.685  3.118  3.219  5/2^{+}  0^{+}  
40  27  283.83  7.10  4.47  43.64  0.36  19.25  2.37  20.41  3.595  3.801  3.126  3.227  5/2^{+}  3/2^{}  
41  28  287.36  7.01  3.89  45.04  3.53  19.94  2.16  21.18  3.606  3.802  3.141  3.241  5/2^{+}  0^{+}  
42  29  287.71  6.85  3.89  45.92  0.35  20.49  2.03  21.75  3.682  3.898  3.148  3.248  5/2^{+}  3/2^{}  
43  30  290.90  6.77  3.54  47.21  3.19  21.03  1.92  22.33  3.706  3.919  3.160  3.260  5/2^{+}  0^{+}  
44  31  291.31  6.62  3.60  48.14  0.41  21.47  1.87  22.83  3.783  4.014  3.165  3.264  5/2^{+}  1/2^{}  
45  32  294.10  6.54  3.20  49.29  2.79  22.07  1.63  23.37  3.806  4.034  3.176  3.275  5/2^{+}  0^{+}  
46  33  294.48  6.40  3.17  50.57  0.37  22.76  1.02  23.99  3.869  4.110  3.179  3.279  5/2^{+}  1/2^{}  
47  34  296.65  6.31  2.55  51.28  2.17  23.06  0.94  24.30  3.907  4.148  3.192  3.291  5/2^{+}  0^{+}  
48  35  295.84  6.16  1.36  51.87  0.81  23.04  0.96  24.31  3.974  4.228  3.192  3.291  5/2^{+}  1/2^{}  
49  36  297.04  6.06  0.39  52.38  1.20  23.49  0.26  25.10  4.003  4.244  3.243  3.340  5/2^{+}  0^{+}  
50  37  296.28  5.93  0.44  0.76  0.28  25.12  4.293  4.605  3.244  3.341  5/2^{+}  1/2^{}  
51  38  296.77  5.82  0.27  0.49  0.02  25.88  4.100  4.342  3.293  3.389  5/2^{+}  0^{+}  
\sigma  2.81  0.085  
Z=14~{}~{}(\rm{Si})  
22  8  135.95  6.18  0.28  1.25  21.08  0.38  3.035  2.660  3.230  3.328  0^{+}  0^{+}  
23  9  150.93  6.56  2.51  14.98  3.06  23.58  0.81  2.991  2.728  3.148  3.248  0^{+}  5/2^{+}  
24  10  170.17  172.01  7.09  7.17  34.22  5.54  19.24  4.66  16.48  2.04  2.966  2.779  3.093  3.195  0^{+}  0^{+}  
25  11  184.22  187.01  7.37  7.48  33.29  8.87  14.05  6.00  15.95  3.33  2.939  2.814  3.033  3.137  0^{+}  5/2^{+}  
26  12  202.45  206.05  7.79  7.92  32.28  12.47  18.23  7.59  15.61  4.65  2.917  2.838  2.983  3.089  0^{+}  0^{+}  
27  13  216.67  219.36  8.02  8.12  32.45  16.56  14.23  9.67  13.45  5.97  2.887  2.846  2.924  3.032  0^{+}  5/2^{+}  
28  14  234.04  236.54  8.36  8.45  31.59  20.48  17.37  11.60  13.13  7.35  2.875  2.862  2.889  2.998  3.122  0^{+}  0^{+} 
29  15  243.44  245.01  8.39  8.45  26.77  21.95  9.40  12.41  13.25  8.61  2.942  2.964  2.919  3.026  3.117  0^{+}  1/2^{+} 
30  16  253.15  255.62  8.44  8.52  19.11  23.58  9.70  13.30  9.27  9.90  3.003  3.049  2.949  3.056  3.133  0^{+}  0^{+} 
31  17  259.99  262.21  8.39  8.46  16.55  25.36  6.85  14.18  9.20  11.08  3.068  3.135  2.985  3.090  0^{+}  3/2^{+}  
32  18  269.34  271.41  8.42  8.48  16.20  27.48  9.35  15.30  8.24  12.13  3.123  3.200  3.020  3.124  0^{+}  0^{+}  
33  19  275.89  275.91  8.36  8.36  15.90  29.49  6.55  16.34  7.20  13.19  3.176  3.266  3.050  3.153  0^{+}  3/2^{+}  
34  20  284.63  283.43  8.37  8.34  15.29  31.58  8.74  17.43  6.47  14.23  3.223  3.320  3.078  3.180  0^{+}  0^{+}  
35  21  287.11  285.90  8.20  8.17  11.21  33.66  2.47  18.50  7.26  15.25  3.263  3.375  3.088  3.190  0^{+}  7/2^{}  
36  22  293.31  292.01  8.15  8.11  8.68  35.36  6.21  19.41  4.25  16.17  3.302  3.424  3.101  3.202  0^{+}  0^{+}  
37  23  294.85  294.28  7.97  7.95  7.74  37.27  1.54  20.40  3.84  17.13  3.344  3.478  3.112  3.213  0^{+}  7/2^{}  
38  24  300.29  299.93  7.90  7.89  6.98  38.66  5.44  21.21  3.59  17.95  3.386  3.529  3.126  3.226  0^{+}  0^{+}  
39  25  301.26  301.51  7.72  7.74  6.41  40.36  0.97  21.90  3.22  18.82  3.429  3.583  3.136  3.236  0^{+}  7/2^{}  
40  26  306.25  306.47  7.66  7.66  5.96  41.54  4.99  22.78  3.12  19.54  3.473  3.635  3.150  3.250  0^{+}  0^{+}  
41  27  307.15  307.85  7.49  7.51  5.89  42.92  0.90  23.32  3.10  20.10  3.548  3.734  3.157  3.257  0^{+}  3/2^{}  
42  28  311.48  7.42  5.24  44.06  4.82  24.12  2.76  20.92  3.564  3.745  3.171  3.271  0^{+}  0^{+}  
43  29  312.46  7.27  5.31  45.24  0.98  24.75  2.65  21.54  3.627  3.825  3.179  3.278  0^{+}  3/2^{}  
44  30  316.17  7.19  4.69  46.30  3.71  25.27  2.44  22.10  3.659  3.859  3.190  3.289  0^{+}  0^{+}  
45  31  317.05  7.05  4.59  47.21  0.88  25.74  2.42  22.61  3.725  3.941  3.194  3.293  0^{+}  1/2^{}  
46  32  320.39  6.97  4.22  48.35  3.34  26.29  2.10  23.13  3.757  3.974  3.206  3.304  0^{+}  0^{+}  
47  33  321.47  6.84  4.42  49.75  1.08  27.00  1.42  23.73  3.810  4.040  3.203  3.301  0^{+}  1/2^{}  
48  34  323.86  6.75  3.47  50.27  2.39  27.21  1.33  24.03  3.857  4.090  3.219  3.317  0^{+}  0^{+}  
49  35  323.28  6.60  1.80  50.50  0.59  27.64  1.40  24.43  3.918  4.160  3.234  3.332  0^{+}  5/2^{}  
50  36  324.95  6.50  1.09  51.40  1.68  27.91  0.63  24.84  3.952  4.185  3.275  3.371  0^{+}  0^{+}  
51  37  324.20  6.36  0.92  0.75  28.47  0.65  24.86  4.238  4.550  3.275  3.372  0^{+}  1/2^{}  
52  38  325.46  6.26  0.51  1.26  28.69  0.35  25.67  4.043  4.276  3.329  3.423  0^{+}  0^{+}  
53  39  324.77  6.13  0.57  0.69  29.42  0.38  25.69  4.300  4.599  3.329  3.424  0^{+}  1/2^{+}  
54  40  325.46  6.03   0.00  0.69  0.17  26.42  4.150  4.389  3.374  3.467  0^{+}  0^{+}  
\sigma  1.84  0.099  
Z=15~{}~{}(\rm{P})  
22  7  108.82  4.95  19.90  1.73  4.105  2.594  4.645  4.713  1/2^{+}  1/2^{}  
23  8  133.70  5.81  0.22  24.88  2.25  21.93  1.01  3.464  2.661  3.824  3.907  1/2^{+}  0^{+}  
24  9  150.46  6.27  41.64  2.82  16.76  0.47  22.53  0.52  3.053  2.720  3.237  3.334  1/2^{+}  5/2^{+}  
25  10  170.98  6.84  37.27  5.46  20.52  0.81  17.82  0.97  3.022  2.773  3.177  3.276  1/2^{+}  0^{+}  
26  11  186.05  7.16  35.60  7.84  15.08  1.83  17.12  2.40  2.997  2.816  3.124  3.225  1/2^{+}  5/2^{+}  
27  12  205.05  206.91  7.59  7.66  34.07  10.18  18.99  2.60  16.56  3.52  2.985  2.854  3.086  3.188  1/2^{+}  0^{+}  
28  13  219.68  221.41  7.85  7.91  33.63  12.68  14.64  3.01  14.58  4.26  2.963  2.872  3.039  3.143  1/2^{+}  5/2^{+}  
29  14  237.73  239.29  8.20  8.25  32.69  15.30  18.05  3.69  14.31  6.98  2.951  2.892  3.005  3.110  1/2^{+}  0^{+}  
30  15  249.13  250.60  8.30  8.35  29.45  18.10  11.40  5.69  13.99  8.05  3.002  2.984  3.020  3.124  1/2^{+}  1/2^{+}  
31  16  260.83  262.92  8.41  8.48  23.10  20.99  11.70  7.69  10.78  8.99  3.047  3.060  3.033  3.137  3.189  1/2^{+}  0^{+} 
32  17  269.07  270.85  8.41  8.46  19.94  23.26  8.24  9.08  10.77  9.84  3.098  3.134  3.057  3.160  1/2^{+}  3/2^{+}  
33  18  279.40  280.96  8.47  8.51  18.57  25.36  10.33  10.06  9.20  10.80  3.146  3.197  3.083  3.185  1/2^{+}  0^{+}  
34  19  286.98  287.24  8.44  8.45  17.91  27.42  7.57  11.08  7.64  11.76  3.192  3.259  3.106  3.208  1/2^{+}  3/2^{+}  
35  20  296.59  295.62  8.47  8.45  17.19  28.46  9.62  11.96  7.36  13.41  3.235  3.311  3.130  3.230  1/2^{+}  0^{+}  
36  21  299.94  299.08  8.33  8.31  12.96  31.33  3.35  12.83  7.21  13.94  3.272  3.365  3.138  3.238  1/2^{+}  7/2^{}  
37  22  307.06  305.90  8.30  8.27  10.47  32.59  7.12  13.75  5.11  15.36  3.308  3.412  3.148  3.248  1/2^{+}  0^{+}  
38  23  309.42  309.64  8.14  8.15  9.48  34.64  2.36  14.57  4.66  16.30  3.346  3.463  3.157  3.257  1/2^{+}  7/2^{}  
39  24  315.70  315.88  8.09  8.10  8.64  36.62  6.28  15.41  4.39  17.12  3.384  3.512  3.169  3.268  1/2^{+}  0^{+}  
40  25  317.45  319.19  7.94  7.98  8.03  38.09  1.75  16.19  3.99  18.01  3.423  3.562  3.178  3.277  1/2^{+}  7/2^{}  
41  26  323.21  324.17  7.88  7.91  7.51  39.74  5.76  16.96  3.86  18.74  3.462  3.611  3.189  3.288  1/2^{+}  0^{+}  
42  27  324.82  326.25  7.73  7.77  7.37  40.99  1.61  17.67  3.42  19.32  3.523  3.692  3.196  3.294  1/2^{+}  3/2^{}  
43  28  329.89  330.65  7.67  7.66  6.68  42.53  5.07  18.41  3.43  20.16  3.545  3.713  3.208  3.306  1/2^{+}  0^{+}  
44  29  331.55  7.54  6.73  43.84  1.66  19.09  3.30  20.87  3.595  3.778  3.214  3.312  1/2^{+}  3/2^{}  
45  30  335.90  7.46  6.01  45.00  4.35  19.73  3.05  21.46  3.631  3.818  3.225  3.322  1/2^{+}  0^{+}  
46  31  337.32  7.33  5.77  46.20  1.42  20.26  2.70  22.12  3.680  3.879  3.230  3.328  1/2^{+}  3/2^{}  
47  32  341.26  7.26  5.36  47.15  3.94  20.87  2.61  22.51  3.721  3.927  3.239  3.337  1/2^{+}  0^{+}  
48  33  343.05  7.15  5.73  48.57  1.79  21.57  1.96  23.01  3.770  3.989  3.237  3.334  1/2^{+}  1/2^{}  
49  34  345.61  7.05  4.35  48.96  2.56  21.75  1.76  23.41  3.818  4.041  3.254  3.351  1/2^{+}  0^{+}  
50  35  345.45  6.91  2.40  49.61  0.16  22.18  1.81  23.89  3.874  4.106  3.268  3.365  1/2^{+}  5/2^{}  
51  36  347.58  6.82  1.97  50.54  2.13  22.63  1.01  24.30  3.906  4.132  3.301  3.397  1/2^{+}  0^{+}  
52  37  347.04  6.67  1.58  50.76  0.55  23.01  0.81  24.76  3.957  4.187  3.323  3.418  1/2^{+}  5/2^{}  
53  38  348.78  6.58  1.20  52.01  1.75  23.32  0.71  25.24  3.993  4.220  3.352  3.446  1/2^{+}  0^{+}  
54  39  348.14  6.45  1.10  52.79  0.64  23.38  0.74  25.25  4.239  4.534  3.352  3.446  1/2^{+}  1/2^{}  
55  40  349.42  6.35  0.63  1.28  23.96  0.09  26.09  4.087  4.316  3.403  3.495  1/2^{+}  0^{+}  
56  41  348.87  6.23  0.73  0.55  0.05  26.13  4.300  4.584  3.405  3.497  1/2^{+}  1/2^{+}  
57  42  348.62  6.12  0.79  0.24  0.33  26.40  4.301  4.576  3.415  3.508  1/2^{+}  0^{+}  
\sigma  1.36  0.052  
Z=16~{}~{}(\rm{S})  
23  7  107.31  4.67  1.51  22.17  2.90  3.371  2.627  3.649  3.735  0^{+}  1/2^{}  
24  8  133.19  5.55  2.76  25.88  0.51  23.07  2.08  3.212  2.677  3.448  3.540  0^{+}  0^{+}  
25  9  150.32  6.01  43.00  0.61  17.13  0.14  23.76  1.18  3.142  2.743  3.346  3.440  0^{+}  5/2^{+}  
26  10  171.70  6.60  38.51  1.53  21.38  0.72  18.61  0.27  3.103  2.796  3.280  3.376  0^{+}  0^{+}  
27  11  187.51  6.94  37.19  3.29  15.81  1.46  17.92  0.60  3.072  2.840  3.223  3.321  0^{+}  5/2^{+}  
28  12  207.36  209.41  7.41  7.48  35.66  4.91  19.85  2.31  17.37  1.52  3.057  2.880  3.183  3.282  0^{+}  0^{+}  
29  13  222.61  224.71  7.68  7.66  35.10  5.93  15.25  2.92  15.75  2.27  3.031  2.899  3.134  3.234  0^{+}  5/2^{+}  
30  14  241.42  243.68  8.05  7.91  34.06  7.38  18.81  3.69  15.54  3.13  3.019  2.923  3.099  3.201  0^{+}  0^{+}  
31  15  254.85  256.74  8.22  8.25  32.24  11.40  13.43  5.71  15.31  4.63  3.054  3.000  3.104  3.205  0^{+}  1/2^{+}  
32  16  269.38  271.78  8.42  8.35  27.96  16.24  14.54  8.55  12.85  6.75  3.077  3.058  3.058  3.198  3.261  0^{+}  0^{+} 
33  17  278.50  280.42  8.44  8.48  23.65  18.51  9.12  9.43  12.50  7.84  3.124  3.131  3.117  3.218  0^{+}  3/2^{+}  
34  18  289.64  291.84  8.52  8.46  20.26  20.30  11.14  10.24  10.03  8.88  3.168  3.194  3.138  3.238  3.285  0^{+}  0^{+} 
35  19  297.92  298.82  8.51  8.51  19.42  22.02  8.27  10.94  9.19  9.79  3.218  3.262  3.166  3.266  0^{+}  3/2^{+}  
36  20  308.55  308.71  8.57  8.45  18.90  23.91  10.63  11.95  8.45  10.82  3.262  3.317  3.190  3.289  3.298  0^{+}  0^{+} 
37  21  312.96  313.02  8.46  8.45  15.04  25.85  4.41  13.02  9.14  11.82  3.297  3.370  3.199  3.297  0^{+}  7/2^{}  
38  22  321.11  321.05  8.45  8.31  12.57  27.80  8.16  14.05  6.12  12.79  3.330  3.416  3.208  3.307  0^{+}  0^{+}  
39  23  324.58  325.43  8.32  8.27  11.63  29.74  3.47  15.16  5.68  13.77  3.364  3.463  3.217  3.315  0^{+}  7/2^{}  
40  24  331.84  333.17  8.30  8.15  10.73  31.55  7.26  16.14  5.35  14.69  3.399  3.508  3.228  3.325  0^{+}  0^{+}  
41  25  334.72  337.42  8.16  8.10  10.13  33.46  2.87  17.27  7.43  15.65  3.432  3.551  3.235  3.333  0^{+}  7/2^{}  
42  26  341.34  344.12  8.13  7.98  9.50  35.09  6.62  18.13  4.72  16.49  3.467  3.596  3.246  3.344  0^{+}  0^{+}  
43  27  343.51  346.74  7.99  8.06  8.79  36.36  2.17  18.69  4.74  17.11  3.523  3.673  3.253  3.350  0^{+}  3/2^{}  
44  28  349.71  351.82  7.95  7.99  8.37  38.23  6.20  19.82  4.10  18.09  3.540  3.688  3.264  3.361  0^{+}  0^{+}  
45  29  352.09  354.68  7.82  7.88  8.58  39.63  2.38  20.54  4.00  18.78  3.586  3.749  3.270  3.366  0^{+}  3/2^{}  
46  30  356.92  7.76  7.20  40.74  4.83  21.02  3.57  19.35  3.624  3.794  3.281  3.377  0^{+}  0^{+}  
47  31  358.97  7.64  6.88  41.92  2.05  21.65  3.18  19.93  3.671  3.855  3.286  3.382  0^{+}  3/2^{}  
48  32  363.22  7.57  6.31  42.83  4.25  21.97  3.07  20.40  3.712  3.903  3.296  3.392  0^{+}  0^{+}  
49  33  365.41  7.47  6.44  43.93  2.18  22.36  2.53  20.87  3.760  3.967  3.291  3.387  0^{+}  1/2^{}  
50  34  368.40  7.37  5.18  44.54  3.00  22.79  2.30  21.32  3.804  4.014  3.315  3.410  0^{+}  0^{+}  
51  35  368.75  7.23  3.34  45.47  0.34  23.30  2.33  21.81  3.856  4.074  3.329  3.424  0^{+}  5/2^{}  
52  36  371.76  7.15  3.35  46.80  3.01  24.17  1.67  22.46  3.891  4.105  3.360  3.454  0^{+}  0^{+}  
53  37  371.88  7.02  3.13  47.85  0.12  24.84  1.49  23.00  3.938  4.156  3.381  3.474  0^{+}  5/2^{}  
54  38  374.38  6.93  2.63  48.92  2.50  25.60  1.35  23.60  3.973  4.188  3.408  3.501  0^{+}  0^{+}  
55  39  374.29  6.81  2.41  50.23  0.09  26.15  0.42  24.14  4.018  4.235  3.432  3.524  0^{+}  5/2^{}  
56  40  376.43  6.72  2.05  50.97  2.14  27.02  0.50  24.68  4.054  4.269  3.456  3.548  0^{+}  0^{+}  
57  41  375.89  6.59  1.60  0.55  27.02  0.50  24.71  4.269  4.547  3.457  3.548  0^{+}  1/2^{+}  
58  42  375.91  6.48  0.52  0.03  27.29  0.18  25.23  4.199  4.444  3.473  3.564  0^{+}  0^{+}  
\sigma  1.93  0.046  
Z=17~{}~{}(\rm{Cl})  
28  11  185.67  6.63  0.40  1.84  18.80  0.32  3.160  2.867  3.336  3.430  3/2^{+}  5/2^{+}  
29  12  206.45  7.12  1.40  20.78  0.91  18.27  1.23  3.134  2.909  3.283  3.379  3/2^{+}  0^{+}  
30  13  222.45  7.42  36.78  2.77  16.00  0.16   