We present a newly enhanced version of the Monte Carlo Shell Model method by incorporating the conjugate gradient method and energy-variance extrapolation. This new method enables us to perform large-scale shell-model calculations that the direct dia
gonalization method cannot reach. This new generation framework of the MCSM provides us with a powerful tool to perform most-advanced large-scale shell-model calculations on current massively parallel computers such as the K computer. We discuss the validity of this method in ab initio calculations of light nuclei, and propose a new method to describe the intrinsic wave function in terms of the shell-model picture. We also apply this new MCSM to the study of neutron-rich Cr and Ni isotopes using the conventional shell-model calculations with an inert 40Ca core and discuss how the magicity of N = 28, 40, 50 remains or is broken.
We discuss a variational calculation for nuclear shell-model calculations and propose a new procedure for the energy-variance extrapolation (EVE) method using a sequence of the approximated wave functions obtained by the variational calculation. The
wave functions are described as linear combinations of the parity, angular-momentum projected Slater determinants, the energy of which is minimized by the conjugate gradient method obeying the variational principle. The EVE generally works well using the wave functions, but we found some difficult cases where the EVE gives a poor estimation. We discuss the origin of the poor estimation concerning shape coexistence. We found that the appropriate reordering of the Slater determinants allows us to overcome this difficulty and to reduce the uncertainty of the extrapolation.
