We report the first global study of dynamic correlation energies (DCEs) associated with rotational motion and quadrupole shape vibrational motion in a covariant energy density functional (CEDF) for 575 even-even nuclei with proton numbers ranging fro
m $Z=8$ to $Z=108$ by solving a five-dimensional collective Hamiltonian, the collective parameters of which are determined from triaxial relativistic mean-field plus BCS calculation using the PC-PK1 force. After taking into account these beyond mean-field DCEs, the root-mean-square (rms) deviation with respect to nuclear masses is reduced significantly down to 1.14 MeV, which is smaller than those of other successful CEDFs: NL3* (2.96 MeV), DD-ME2 (2.39 MeV), DD-ME$delta$ (2.29 MeV) and DD-PC1 (2.01 MeV). Moreover, the rms deviation for two-nucleon separation energies is reduced by $sim34%$ in comparison with cranking prescription.
Starting from the adiabatic time-dependent Hartree-Fock approximation (ATDHF), we propose an efficient method to calculate the Thouless-Valatin moments of inertia for the nuclear system. The method is based on the rapid convergence of the expansion o
f the inertia matrix. The accuracy of the proposed method is verified in the rotational case by comparing the results with the exact Thouless-Valatin moments of inertia calculated using the self-consistent cranking model. The proposed method is computationally much more efficient than the full ATDHF calculation, yet it retains a high accuracy of the order of 1%.
The structure of low-energy collective states in proton-deficient N=28 isotones is analyzed using structure models based on the relativistic energy density functional DD-PC1. The relativistic Hartree-Bogoliubov model for triaxial nuclei is used to ca
lculate binding energy maps in the $beta$-$gamma$ plane. The evolution of neutron and proton single-particle levels with quadrupole deformation, and the occurrence of gaps around the Fermi surface, provide a simple microscopic interpretation of the onset of deformation and shape coexistence. Starting from self-consistent constrained energy surfaces calculated with the functional DD-PC1, a collective Hamiltonian for quadrupole vibrations and rotations is employed in the analysis of excitation spectra and transition rates of $^{46}$Ar, $^{44}$S, and $^{42}$Si. The results are compared to available data, and previous studies based either on the mean-field approach or large-scale shell-model calculations. The present study is particularly focused on $^{44}$S, for which data have recently been reported that indicate pronounced shape coexistence.
The shape evolution and shape coexistence phenomena in neutron-rich nuclei at $Napprox60$, including Kr, Sr, Zr, and Mo isotopes, are studied in the covariant density functional theory (DFT) with the new parameter set PC-PK1. Pairing correlations are
treated using the BCS approximation with a separable pairing force. Sharp rising in the charge radii of Sr and Zr isotopes at N=60 is observed and shown to be related to the rapid changing in nuclear shapes. The shape evolution is moderate in neighboring Kr and Mo isotopes. Similar as the results of previous Hartree-Fock-Bogogliubov (HFB) calculations with the Gogny force, triaxiality is observed in Mo isotopes and shown to be essential to reproduce quantitatively the corresponding charge radii. In addition, the coexistence of prolate and oblate shapes is found in both $^{98}$Sr and $^{100}$Zr. The observed oblate and prolate minima are related to the low single-particle energy level density around the Fermi surfaces of neutron and proton respectively. Furthermore, the 5-dimensional (5D) collective Hamiltonian determined by the calculations of the PC-PK1 energy functional is solved for $^{98}$Sr and $^{100}$Zr. The resultant excitation energy of $0^+_2$ state and E0 transition strength $rho^2(E0;0^+_2rightarrow0^+_1)$ are in rather good agreement with the data. It is found that the lower barrier height separating the two competing minima along the $gamma$ deformation in $^{100}$Zr gives rise to the larger $rho^2(E0;0^+_2rightarrow0^+_1)$ than that in $^{98}$Sr.
The analysis of shape transitions in Nd isotopes, based on the framework of relativistic energy density functionals and restricted to axially symmetric shapes in Ref. cite{PRL99}, is extended to the region $Z = 60$, 62, 64 with $N approx 90$, and inc
ludes both $beta$ and $gamma$ deformations. Collective excitation spectra and transition probabilities are calculated starting from a five-dimensional Hamiltonian for quadrupole vibrational and rotational degrees of freedom, with parameters determined by constrained self-consistent relativistic mean-field calculations for triaxial shapes. The results reproduce available data, and show that there is an abrupt change of structure at N=90 that can be approximately characterized by the X(5) analytic solution at the critical point of the first-order quantum phase transition between spherical and axially deformed shapes.
