Search

The predictive power of the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with density functional PC-PK1 is demonstrated for superheavy region ($101 \leqslant Z \leqslant 120$) by comparing with available experimental and evaluated data in the AME2020. The DRHBc theory predicts 93 bound nuclei beyond the drip line $N = 258$ in the region of $106 \leqslant Z \leqslant 112$, which form a stability peninsula. The odd-even differences between odd-$N$ and even-$N$ nuclei are remarkable in the stability peninsula; the number of bound odd-$N$ nuclei is less than that of bound even-$N$ nuclei, and the one-neutron separation energy of an odd-$N$ nucleus is smaller than those of its neighboring even-$N$ nuclei due to the blocking effect. The deformation effect is indispensable for the reentrant stability beyond the drip line by significantly affecting the structure of single-particle levels around the Fermi energy. The interplay between deformation and pairing effects affects the position where the odd-$N$ nucleus becomes bound in the stability peninsula. By examining the deformation effect at different orders, it is found that quadrupole deformation makes leading contribution to the appearance of stability peninsula and the effects of hexadecapole and hexacontatetrapole deformations are nonnegligible., Comment: 6 figures

Within the supersymmetry scheme, which includes many-body interactions and a perturbation possessing the SO(5) (or SU(5)) symmetry, the rotational bands of the $A\sim 250$ mass region are studied systematically. A novel modification is introduced, extending the Arima coefficient to the third order. This study is dedicated to the quantitative analysis of evolving trends in intraband $\gamma$-transition energy, kinematic, and the dynamic moment of inertia within the rotational bands of $^{244}$Pu and $^{248}$Cm. The computed outcomes exhibit an exceptional degree of agreement with experimental observations across various conditions. The significance of including a higher-order Arima coefficient is further examined by contrasting it with the previously suggested model. The calculated results demonstrate the significance of both the anti-pairing and pairing effects in the evolution of the dynamic moment of inertia. Additionally, these insights reveal the importance of a newly introduced parameter in accurately depicting complex nuclear behaviors such as back-bending, up-bending, and the downturn in the moment of inertia.

The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $\alpha$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $\alpha$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed., Comment: 394 pages, 17 figures, 2 tables, published in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

For the first time, we conducted a comprehensive analysis of the rotational bands in the \( A\sim250 \) mass region. Utilizing a variety of rotational energy models and formulas, we have extracted free parameters for 36 rotational bands within this mass region, encompassing neutron numbers from N = 148 to 152. A significant enhancement has been made to the vibrational distortion model, leading to an exceptional match with the experimental data on dynamic moment of inertia across most of the rotational bands studied. Through the application of this revised model, we have categorized the dynamic moment of inertia into three distinct groups. The implications of the modified vibrational distortion term and how it evolves with changes in rotational frequency is also investigated. Our systematic investigation brings to light the peculiar behavior of Pu isotopes, especially in ground state bands of \( ^{242,244}\text{Pu}\), within this mass region. We provide a thorough discussion on various aspects such as the role of shape-fluctuation energy, the softness parameter, band-head and average moment of inertia, the prominence of the anti-pairing effect, and the transition from pairing to anti-pairing effects in Pu isotopes. Moreover, these findings are compared with other isotones, offering a comprehensive understanding of their unique characteristics.

Based on the reflection-asymmetric Nilsson potential, the parity-doublet rotational bands in odd-$A$ isotones \element{237}U and \element{239}Pu have been investigated by using the particle-number-conserving (PNC) method in the framework of the cranked shell model (CSM). The experimental kinematic moments of inertia (MOIs) and angular momentum alignments are reproduced very well by the PNC-CSM calculations. The significant differences of rotational properties between \element{237}U and \element{239}Pu are explained with the contribution of nucleons occupying proton octupole-correlation pairs of $\pi ^{2} i_{13/2}f_{7/2}$. The upbendings of moments of inertia of the parity-doublet bands in \element{237}U are due to the interference terms of alignments of protons occupying $\pi f_{7/2}$ (${1/2}$) and the high-$j$ intruder $\pi i_{13/2}$ $(1/2, 3/2)$ orbitals. The splittings between the simplex partner bands of the parity-doublet bands in both \element{237}U and \element{239}Pu result from the contribution of alignment of neutron occupying the $\nu d_{5/2}$ ($1/2$) orbital.

The aim of this work is to extend the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) based on the point-coupling density functionals to odd-$A$ and odd-odd nuclei and examine its applicability by taking odd-$A$ Nd isotopes as examples. In the DRHBc theory, the densities and potentials with axial deformation are expanded in terms of Legendre polynomials, and the relativistic Hartree-Bogoliubov equations are solved in a Dirac Woods-Saxon basis to include the continuum effects. For an odd-$A$ or odd-odd nucleus, the blocking effect of unpaired nucleon(s) is taken into account with the equal filling approximation. To determine its ground state, an automatic blocking procedure is adopted, in which the orbital with the lowest quasiparticle energy is blocked during the iteration. This procedure is justified by comparing with the results from the orbital-fixed blocking calculations, in which the blocked orbital near the Fermi surface is fixed during the iteration. The ground states for both light and heavy nuclei can be provided by the automatic blocking procedure as the orbital-fixed blocking procedure, but with considerably reduced computational cost. The numerical details for even-even nuclei are found to be valid for odd-$A$ and odd-odd nuclei as well. Taking Nd isotopes including both even-even and odd-$A$ ones as examples, the calculated ground-state properties with PC-PK1 are in good agreement with the available experimental data. This work paves the way to construct the DRHBc mass table including all even-even, odd-$A$ and odd-odd nuclei in the nuclear chart., Comment: 39 pages, 13 figures, and 1 table

We build and train the artificial neural network model (ANN) based on the experimental $\alpha$-decay energy ($Q_{\alpha}$) data. Besides decays between the ground states of parent and daughter nuclei, decays from the ground state of parent nuclei to the excited state of daughter nuclei are also included. By this way, the number of samples are increased dramatically. The results calculated by ANN model reproduce the experimental data with a good accuracy. The root-mean-square (rms) relative to the experiment data is 0.105 MeV. The influence of different input is investigated. It is found that either the shell effect or the pairing effect results in an obvious improvement of the predictive power of ANN model, and the shell effect plays a more important role. The optimal result can be obtained as both the shell and pairing effects are considered simultaneously. Application of ANN model in prediction of the $\alpha$-decay energy shows the neutron magic number at $N=184$, and a possible sub-shell gap around $N=174$ or 176 in the superheavy nuclei region.

Ground-state properties of even-even nuclei with $8\le Z\le120$ from the proton drip line to the neutron drip line have been investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the density functional PC-PK1. With the effects of deformation and continuum included simultaneously, 2583 even-even nuclei are predicted to be bound. The calculated binding energies, two-nucleon separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. The rms deviation from the 637 mass data is 1.518 MeV, providing one of the best microscopic descriptions for nuclear masses. The drip lines obtained from DRHBc calculations are compared with other calculations, including the spherical relativistic continuum Hartree-Bogoliubov (RCHB) and triaxial relativistic Hartree-Bogoliubov (TRHB) calculations with PC-PK1. The deformation and continuum effects on the limits of the nuclear landscape are discussed. Possible peninsulas consisting of bound nuclei beyond the two-neutron drip line are predicted. The systematics of the two-nucleon separation energies, two-nucleon gaps, rms radii, quadrupole deformations, potential energy curves, neutron densities, neutron mean-field potentials, and pairing energies in the DRHBc calculations are also discussed. In addition, the $\alpha$ decay energies extracted are in good agreement with available data., Comment: 217 pages, 15 figures, 2 tables, accepted for publication in Atomic Data and Nuclear Data Tables, data file in the TXT form is available for download under "Ancillary files"

Based on the relativistic calculations of the nuclear masses in the transfermium region from No $(Z=102)$ to Ds $(Z=110)$ by the deformed relativistic Hartree-Bogoliubov theory in continuum, the possible existence of the bound nuclei beyond the neutron drip lines is studied. The two-neutron and multi-neutron emission bound nuclei beyond the primary neutron drip line of $N=258$ are predicted in $Z=106,108$ and $110$ isotopes. Detailed microscopic mechanism investigation reveals that nuclear deformation plays a vital role in the existence of the bound nuclei beyond the drip line. Furthermore, not only the quadrupole deformation $\beta_{2}$, but also the higher orders of deformation are indispensible in the reliable description of the phenomenon of the reentrant binding.

We study the chiral phase transition inside a rotating cylinder within the framework of the Namb--Jona-Lasinio model. A spectral boundary condition is imposed to avoid faster than light. We investigate how the geometry of the cylinder and rotation influence the chiral phase transition at finite temperature and chemical potential. The inhomogeneous effects caused by the finite size and rotation are also taken into account. It is found that finite size will reduce the chiral transition temperature and raises the chiral transition chemical potential, while the rotation reduces both the chiral transition temperature and chemical potential. In addition, we discuss the implications of our results in heavy-ion collisions and equation of states of neutron star.

Within the framework of Dyson--Schwinger equations of QCD, we study the effect of finite volume on the chiral phase transition in a sphere with the MIT boundary condition. We find that the chiral quark condensate $\langle\bar{\psi} \psi\rangle$ and pseudotransition temperature $T_{pc}$ of the crossover decreases as the volume decreases, until there is no chiral crossover transition at last. We find that the system for $R = \infty $\ fm is indistinguishable from $R=10$ fm and there is a significant decrease in $T_{pc}$ with $R$ as $R<4$ fm. When $R<1.5$ fm, there is no chiral transition in the system., Comment: 5 pages, 2 figures

The particle-number conserving (PNC) method in the framework of cranked shell model (CSM) is developed to deal with the reflection-asymmetric nuclear system by applying the $S_x$ symmetry. Based on an octupole-deformed Nilsson potential, the alternating-parity bands in \element{236,238}{U} and \element{238,240}{Pu} are investigated. The experimental kinematic moments of inertia (MoI) and the angular momentum alignments of all studied bands are reproduced well in the PNC-CSM calculations. The striking difference of rotational behaviors between U and Pu isotopes can be linked to the strength of octupole correlations. The upbendings of the alternating-parity bands in\element{236,238}{U} are due to the alignments of pairs of nucleons occupying $\nu g_{9/2}$, $\pi f_{7/2}$ orbitals and $\nu j_{15/2}$, $\pi i_{13/2}$ high-$j$ intruder orbitals. Particularly, the interference terms of nucleon occupying the octupole-correlation pairs of $\nu^2 j_{15/2} g_{9/2}$ and of $\pi^2 i_{13/2} f_{7/2}$ give a very important contribution to the suddenly gained alignments.

We study the chirally imbalanced hot and dense strongly interacting matter by means of the Dyson-Schwinger equations (DSEs). The chiral phase diagram is studied in the presence of chiral chemical potential $\mu_5$. The chiral quark condensate $\langle \bar{\psi} \psi \rangle$ is obtained with the Cornwall-Jackiw-Tomboulis (CJT) effective action in concert with the Rainbow truncation. Catalysis effect of dynamical chiral symmetry breaking (DCSB) by $\mu_5$ is observed. We examine with two popular gluon models and consistency is found within the DSE approach, as well as in comparison with lattice QCD. The CEP location $(\mu_E,T_E)$ shifts toward larger $T_E$ but constant $\mu_E$ as $\mu_5$ increases. A technique is then introduced to compute the chiral charge density $n_5$ from the fully dressed quark propagator. We find the $n_5$ generally increases with temperature $T$, quark number chemical potential $\mu$ and $\mu_5$. Since the chiral magnetic effect (CME) is typically investigated with peripheral collisions, we also investigate the finite size effect on $n_5$ and find an increase in $n_5$ with smaller system size.

The ground-state bands (GSBs) in the even-even hafnium isotopes $^{170-184}$Hf are investigated by using the cranked shell model (CSM) with pairing correlations treated by the particle-number conserving (PNC) method. The experimental kinematic moments of inertia are reproduced very well by theoretical calculations. The second upbending of the GSB at high frequency $\hbar\omega\approx0.5$ MeV observed (predicted) in $^{172}$Hf ($^{170,174-178}$Hf) attributes to the sudden alignments of the proton high-$j$ orbitals $\pi1i_{13/2}$ $(1/2^{+}[660])$, $\pi1h_{9/2}$ $(1/2^{-}[541])$ and orbital $\pi1h_{11/2}$ $(7/2^{-}[523])$. The first upbendings of GSBs at low frequency $\hbar\omega=0.2-0.3$ MeV in $^{170-178}$Hf, which locate below the deformed neutron shell $N=108$, attribute to the alignment of the neutron orbital $\nu1i_{13/2}$. For the heavier even-even isotopes $^{180-184}$Hf, compared to the lighter isotopes, the first band-crossing is delayed to the high frequency due to the existence of the deformed shells $N=108,116$. The upbendings of GSBs in $^{180-184}$Hf are predicted to occur at $\hbar\omega\approx0.5$MeV, which come from the sharp raise of the simultaneous alignments of both proton $\pi1i_{13/2}$, $\pi1h_{9/2}$ and neutron $\nu2g_{9/2}$ orbitals. The pairing correlation plays a very important role in the rotational properties of GSBs in even-even isotopes $^{180-184}$Hf. Its effects on upbendings and band-crossing frequencies are investigated., Comment: 14pages,12figures

The multi-particle states and rotational properties of two-particle bands in $^{254}$No are investigated by the cranked shell model (CSM) with pairing correlations treated by a particle-number conserving (PNC) method. For the first time, the rotational bands on top of two-particle $K^{\pi}=3^+,8^-$ and $10^+$ states and the pairing reduction are studied theoretically in $^{254}$No. The experimental excitation energies and moments of inertia for the multi-particle state are reproduced well by the calculation. Better agreement with the data are achieved by including the high-order deformation $\varepsilon_{6}$ which leads to enlarged $Z=100$ and $N=152$ deformed shell gaps. The rise of the $J^{(1)}$ in these two-particle bands compared with the ground-state band is attributed to the pairing reduction due to the Pauli blocking effects.

The newly observed isomer and ground-state band in the odd-Z neutron-rich rare-earth nucleus $^{163}$Eu are investigated by using the cranked shell model (CSM) with pairing treated by the particle-number conserving (PNC) method. This is the first time detailed theoretical investigations are performed of the observed $964(1)$ keV isomer and ground-state rotational band in $^{163}$Eu. The experimental data are reproduced very well by the theoretical results. The configuration of the $964(1)$ keV isomer is assigned as the three-particle state $\frac{13}{2}^{-}(\nu\frac{7}{2}^{+}[633]\otimes\nu\frac{1}{2}^{-}[521]\otimes\pi\frac{5}{2}^{+}[413]$). More low-lying multi-particle states are predicted in $^{163}$Eu. Due to its significant effect on the nuclear mean field, the high-order $\varepsilon_{6}$ deformation plays an important role in the energy and configuration assignment of the multi-particle states. Compared to its neighboring even-even nuclei $^{162}$Sm and $^{164}$Gd, there is a $10\%\sim15\%$ increase of $J^{(1)}$ of the one-particle ground-state band in $^{163}$Eu. This is explained by the pairing reduction due to the blocking of the nucleon on the proton $\pi\frac{5}{2}^{+}$[413] orbital in $^{163}$Eu.

Inspired by the newly discovered isomeric states in the rare-earth neutron-rich nuclei, high-$K$ isomeric states in neutron-rich samarium and gadolinium isotopes are investigated within the framework of the cranked shell model (CSM) with pairing correlation treated by a particle-number-conserving (PNC) method. The experimental multi-particle state energies and moments of inertia are reproduced quite well by the PNC-CSM calculations. A remarkable effect from the high-order deformation $\varepsilon_{6}$ is demonstrated. Based on the occupation probabilities, the configurations are assigned to the observed high-$K$ isomeric states. The lower $5^-$ isomeric state in $^{158}$Sm is preferred as the two-proton state with configuration $\pi\frac{5}{2}^{+}[413]\otimes\pi\frac{5}{2}^{-}[532]$. More low-lying two-particle states are predicted. The systematics of the electronic quadrupole transition probabilities, $B(E2)$ values along the neodymium, samarium, gadolinium and dysprosium isotopes and $N=96,98,100,102$ isotones chains is investigated to reveal the midshell collectivities.

Experimentally observed superdeformed (SD) rotational bands in $^{36}$Ar and $^{40}$Ar are studied by the cranked shell model (CSM) with the paring correlations treated by a particle-number-conserving (PNC) method. This is the first time the PNC-CSM calculations are performed on the light nuclear mass region around $A=40$. The experimental kinematic moments of inertia $J^{(1)}$ versus rotational frequency are reproduced well. The backbending of the SD band at frequency around $\hbar\omega=1.5$ MeV in $^{36}$Ar is attributed to the sharp rise of the simultaneous alignments of the neutron and proton $1d_{5/2}[202]5/2$ pairs and $1f_{7/2}[321]3/2$ pairs, which is the consequence of the band crossing between the $1d_{5/2}[202]5/2$ and $1f_{7/2}[321]3/2$ configuration states. The gentle upbending at the low frequency of the SD band in $^{40}$Ar is mainly effected by the alignments of the neutron $1f_{7/2}[321]3/2$ pairs and proton $1d_{5/2}[202]5/2$ pairs. The PNC-CSM calculations show that besides the diagonal parts, the off-diagonal parts of the alignments play an important role in the rotational behavior of the SD bands.

Experimentally observed ground state band based on the $1/2^{-}[521]$ Nilsson state and the first exited band based on the $7/2^{-}[514]$ Nilsson state in the odd-$Z$ nucleus $^{255}$Lr are studied by the cranked shell model (CSM) with the paring correlations treated by the particle-number-conserving (PNC) method. This is the first time the detailed theoretical investigations being performed on these rotational bands. Both the experimental kinematic and dynamic moment of inertia ($\mathcal{J}^{(1)}$ and $\mathcal{J}^{(2)}$) versus rotational frequency are reproduced quite well by the PNC-CSM calculations. By comparing the theoretical kinematic moment of inertia $\mathcal{J}^{(1)}$ with the experimental ones extracted from different spin assignments, the spin $17/2^{-}\rightarrow13/2^{-}$ is assigned to the lowest-lying $196.6(5)$ keV transition of the $1/2^{-}[521]$ band, and $15/2^{-}\rightarrow11/2^{-}$ to the $189(1)$ keV transition of the $7/2^{-}[514]$ band, respectively. The proton $N=7$ major shell is included in the calculations. The intruder of the high$-j$ low$-\Omega$ orbitals $1j_{15/2}$ $ (1/2^{-}[770])$ at the high spin leads to the band-crossing at $\hbar\omega\approx0.20$ ($\hbar\omega\approx0.25$) MeV for the $7/2^{-}[514]$ $\alpha=-1/2$ ($\alpha=+1/2$) band, and at $\hbar\omega\approx0.175$ MeV for the $1/2^{-}[521]$ $\alpha=-1/2$ band, respectively. Further investigations show that the band-crossing frequencies are quadrupole deformation dependent., Comment: arXiv admin note: text overlap with arXiv:1208.1156 by other authors

A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the scalar-tensor subset of alternative gravity models. We focus on effects of the slope $L$ of nuclear symmetry energy at saturation density and the high-density behavior of nuclear symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of nuclear symmetry energy leads to large changes in the binding energy of neutron stars. The difference in predictions using the GR and the scalar-tensor theory appears only for massive neutron stars, and even then is significantly smaller than the difference resulting from variations in the symmetry energy., Comment: To appear in the Proceedings of the 13th international symposium on Origin of Matter and Evolution of Galaxies (OMEG2015) , June 24-27, 2015, Beijing, China

There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperons and new phases, the most uncertain part of the nucleonic EOS is currently the density dependence of nuclear symmetry energy. To provide information that may help break the EOS-gravity degeneracy, we investigate effects of symmetry energy within its uncertain range determined by terrestrial nuclear laboratory experiments on the gravitational binding energy and spacetime curvature of neutron stars within GR and the scalar-tensor (ST) theory of gravity. In particular, we focus on effects of the following parameters characterizing the EOS of neutron-rich nucleonic matter: (1) the incompressibility $K_0$ of symmetric nuclear matter, (2) the slope $L$ of symmetry energy at saturation density and (3) the high-density behavior of symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of symmetry energy leads to large changes in both the binding energy and curvature of neutron stars while effects of varying the more constrained $K_0$ are negligibly small. The difference in predictions using the GR and the ST theory appears only for massive neutron stars, and is significantly smaller than the differences resulting from variations in the symmetry energy. We conclude that within the ST theory of gravity, the EOS-gravity degeneracy has been broken by the recent relativistic pulsar measurements, and that measurements of neutron star properties sensitive to the compactness constrain mainly the density dependence of the symmetry energy., Comment: Note a title change in the second version. An additional section has been added. 17 pages, 8 figures (3 additional figures), 1 table

The structure of nuclei with $Z\sim100$ is investigated systematically by the Cranked Shell Model (CSM) with pairing correlations treated by a Particle-Number Conserving (PNC) method. In the PNC method, the particle number is conserved and the Pauli blocking effects are taken into account exactly. By fitting the experimental single-particle spectra in these nuclei, a new set of Nilsson parameters ($\kappa$ and $\mu$) is proposed. The experimental kinematic moments of inertia and the band-head energies are reproduced quite well by the PNC-CSM calculations. The band crossing, the effects of high-$j$ intruder orbitals and deformation are discussed in detail., Comment: To appear in the Proceedings of the International Nuclear Physics Conference (INPC2013), June 2-7, 2013, Florence, Italy

The rotational bands in nuclei with $Z \approx 100$ are investigated systematically by using a cranked shell model (CSM) with the pairing correlations treated by a particle-number conserving (PNC) method, in which the blocking effects are taken into account exactly. By fitting the experimental single-particle spectra in these nuclei, a new set of Nilsson parameters ($\kappa$ and $\mu$) and deformation parameters ($\varepsilon_2$ and $\varepsilon_4$) are proposed. The experimental kinematic moments of inertia for the rotational bands in even-even, odd-$A$ and odd-odd nuclei, and the bandhead energies of the 1-quasiparticle bands in odd-$A$ nuclei, are reproduced quite well by the PNC-CSM calculations. By analyzing the $\omega$-dependence of the occupation probability of each cranked Nilsson orbital near the Fermi surface and the contributions of valence orbitals in each major shell to the angular momentum alignment, the upbending mechanism in this region is understood clearly., Comment: 21 pages, 24 figures, extended version of arXiv: 1101.3607 (Phys. Rev. C83, 011304R); added refs.; added Fig. 4 and discussions; Phys. Rev. C, in press