• The excitation function for the 7 alpha de-excitation of $^{28}$Si nuclei excited to high excitation energies in the collisions of 35 MeV/nucleon $^{28}$Si with $^{12}$C reveals resonance structures that may indicate the population of high spin toroidal isomers such as those predicted by a number of recent theoretical calculations. This interpretation is supported by extended theoretical analyses.
  • The paradox on the previous interpretation for the nuclear chiral geometry based on the effective angle has been clarified by reexamining the system with the particle-hole configuration $\pi (1h_{11/2})^1 \otimes \nu(1h_{11/2})^{-1}$ and rotor with deformation parameter $\gamma=30^\circ$. It is found that the paradox is caused by the fact that the angular momentum of the rotor is much smaller than those of the proton and the neutron near the bandhead. Hence, it does not support a chiral rotation interpretation near the bandhead. The nuclear chiral geometry based on the effective angle makes sense only when the angular momentum of the rotor becomes comparable with those of the proton and the neutron at the certain spin region.
  • The effective field theory (EFT) for triaxially deformed even-even nuclei is generalized to include the vibrational degrees of freedom. The pertinent Hamiltonian is constructed up to next-to-leading order. The leading order part describes the vibrational motion, and the NLO part couples rotations to vibrations. The applicability of the EFT Hamiltonian is examined through the description of the energy spectra of the ground state bands, $\gamma$-bands, and $K=4$ bands in the $^{108, 110, 112}$Ru isotopes. It is found that by taking into account the vibrational degrees of freedom, the deviations for high-spin states in the $\gamma$-band observed in the EFT with only rotational degrees of freedom disappear. This supports the importance of including vibrational degrees of freedom in the EFT formulation for the collective motion of triaxially deformed nuclei.
  • Microscopic signatures of nuclear ground-state shape phase transitions in odd-mass Eu isotopes are explored starting from excitation spectra and collective wave functions obtained by diagonalization of a core-quasiparticle coupling Hamiltonian based on energy density functionals. As functions of the physical control parameter -- the number of nucleons -- theoretical low-energy spectra, two-neutron separation energies, charge isotope shifts, spectroscopic quadrupole moments, and $E2$ reduced transition matrix elements accurately reproduce available data, and exhibit more pronounced discontinuities at neutron number $N=90$, compared to the adjacent even-even Sm and Gd isotopes. The enhancement of the first-order quantum phase transition in odd-mass systems can be attributed to a shape polarization effect of the unpaired proton which, at the critical neutron number, starts predominantly coupling to Gd core nuclei that are characterized by larger quadrupole deformation and weaker proton pairing correlations compared to the corresponding Sm isotopes.
  • The shell-model-like approach is implemented to treat the cranking many-body Hamiltonian based on the covariant density functional theory including pairing correlations with exact particle number conservation. The self-consistency is achieved by iterating the single-particle occupation probabilities back to the densities and currents. As an example, the rotational structures observed in the neutron-rich nucleus $^{60}$Fe are investigated and analyzed. Without introducing any \emph{ad hoc} parameters, the bandheads, the rotational spectra, and the relations between the angular momentum and rotational frequency for the positive parity band A, and negative parity bands B and C are well reproduced. The essential role of the pairing correlations is revealed. It is found that for band A, the bandcrossing is due to the change of the last two occupied neutrons from the $1f_{5/2}$ signature partners to the $1g_{9/2}$ signature partners. For the two negative parity signature partner bands B and C, the bandcrossings are due to the pseudo-crossing between the $1f_{7/2,~5/2}$ and the $1f_{5/2,~1/2}$ orbitals. Generally speaking, the deformation parameters $\beta$ for bands A, B, and C decrease with rotational frequency. For band A, the deformation jumps from $\beta\sim0.19$ to $\beta\sim0.29$ around the bandcrossing. In comparison with its signature partner band C, band B exhibits appreciable triaxial deformation.
  • The ground-state properties of nuclei with 8 $\le$ Z $\le$ 120 from the proton drip line to the neutron drip line have been investigated using the relativistic continuum Hartree-Bogoliubov (RCHB) theory with the relativistic density functional PC-PK1. 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, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, ground-state spins and parities are tabulated. The extension of the nuclear landscape obtained with RCHB is discussed in detail, in particular for the neutron-rich 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 alpha-decay energies and proton emitters based on the RCHB calculations are investigated.
  • The pairing-plus-quadrupole Hamiltonian is diagonalized in a symmetry restored basis, i.e., the triaxial quasiparticle-states with angular momentum and particle number projections, and applied for chiral doublet bands in 128Cs. The observed energy spectra and electromagnetic transition probabilities are reproduced well without introducing any parameter. The orientation of the angular momentum in the intrinsic frame is investigated by the distributions of its components on the three principle axes as well as those of its tilted angles. The evolution of the chirality with spin is illustrated and the chiral geometry is demonstrated in the angular momentum projected model for the first time.
  • Effective field theory (EFT) is generalized to investigate the rotational motion of triaxially deformed even-even nuclei. A Hamiltonian, called the triaxial rotor model (TRM), is obtained up to next-to-leading order (NLO) within the EFT formalism. Its applicability is examined by comparing with a five-dimensional collective Hamiltonian (5DCH) for the description of the energy spectra of the ground state and $\gamma$ band in Ru isotopes. It is found that by taking into account the NLO corrections, the ground state band in the whole spin region and the $\gamma$ band in the low spin region are well described. The results presented here indicate that it should be possible to further generalize the EFT to triaxial nuclei with odd mass number.
  • A new method to solve the Dirac equation on a 3D lattice is proposed, in which the variational collapse problem is avoided by the inverse Hamiltonian method and the fermion doubling problem is avoided by performing spatial derivatives in momentum space with the help of the discrete Fourier transform, i.e., the spectral method. This method is demonstrated in solving the Dirac equation for a given spherical potential in 3D lattice space. In comparison with the results obtained by the shooting method, the differences in single particle energy are smaller than $10^{-4}$~MeV, and the densities are almost identical, which demonstrates the high accuracy of the present method. The results obtained by applying this method without any modification to solve the Dirac equations for an axial deformed, non-axial deformed, and octupole deformed potential are provided and discussed.
  • Background: The discovery of neutrinoless double-beta ($0\nu\beta\beta$) decay would demonstrate the nature of neutrinos, have profound implications for our understanding of matter-antimatter mystery, and solve the mass hierarchy problem of neutrinos. The calculations for the nuclear matrix elements $M^{0\nu}$ of $0\nu\beta\beta$ decay are crucial for the interpretation of this process. Purpose: We study the effects of relativity and nucleon-nucleon short-range correlations on the nuclear matrix elements $M^{0\nu}$ by assuming the mechanism of exchanging light or heavy neutrinos for the $0\nu\beta\beta$ decay. Methods: The nuclear matrix elements $M^{0\nu}$ are calculated within the framework of covariant density functional theory, where the beyond-mean-field correlations are included in the nuclear wave functions by configuration mixing of both angular-momentum and particle-number projected quadrupole deformed mean-field states. Results: The nuclear matrix elements $M^{0\nu}$ are obtained for ten $0\nu\beta\beta$-decay candidate nuclei. The impact of relativity is illustrated by adopting relativistic or nonrelativistic decay operators. The effects of short-range correlations are evaluated. Conclusions: The effects of relativity and short-range correlations play an important role in the mechanism of exchanging heavy neutrinos though the influences are marginal for light neutrinos. Combining the nuclear matrix elements $M^{0\nu}$ with the observed lower limits on the $0\nu\beta\beta$-decay half-lives, the predicted strongest limits on the effective masses are $|\langle m_\nu\rangle|<0.06~\mathrm{eV}$ for light neutrinos and $|\langle m_{\nu_h}^{-1}\rangle|^{-1}>3.065\times 10^8~\mathrm{GeV}$ for heavy neutrinos.
  • The recently reported wobbling bands in $^{135}$Pr are investigated by the collective Hamiltonian, in which the collective parameters, including the collective potential and the mass parameter, are respectively determined from the tilted axis cranking (TAC) model and the harmonic frozen alignment (HFA) formula. It is shown that the experimental energy spectra of both yrast and wobbling bands are well reproduced by the collective Hamiltonian. It is confirmed that the wobbling mode in $^{135}$Pr changes from transverse to longitudinal with the rotational frequency. The mechanism of this transition is revealed by analyzing the effective moments of inertia of the three principal axes, and the corresponding variation trend of the wobbling frequency is determined by the softness and shapes of the collective potential.
  • Motivated by the recently observed hypernucleus (Kiso event) $^{15}_{\Xi}$C ($^{14}$N$+\Xi^-$), we identify the state of this system theoretically within the framework of the relativistic-mean-field and Skyrme-Hartree-Fock models. The $\Xi N$ interactions are constructed to reproduce the two possibly observed $\Xi^-$ removal energies, $4.38\pm 0.25$ MeV or $1.11\pm 0.25$ MeV. The present result is preferable to be $^{14}{\rm N}({\rm g.s.})+\Xi^-(1p)$, corresponding to the latter value.
  • A new method to calculate spectroscopic properties of deformed nuclei is proposed: configuration interaction on top of projected density functional theory (CI-PDFT). The general concept of this approach is discussed in the framework of covariant density functional theory and its validity is illustrated in an application to the yrast band of the nucleus Cr-54. It is found that the experimentally observed excitation energies for the yrast band in Cr-54 can be well reproduced. In contrast to conventional shell-model calculations, there is no core and only a relatively small number of configurations is sufficient for a satisfying description. No new parameters are necessary, because the effective interaction is derived from an universal density functional given in the literature.
  • A two-dimensional collective Hamiltonian (2DCH) on both azimuth and polar motions in triaxial nuclei is proposed to investigate the chiral and wobbling modes. In the 2DCH, the collective potential and the mass parameters are determined from three-dimensional tilted axis cranking (TAC) calculations. The broken chiral and signature symmetries in the TAC solutions are restored by the 2DCH. The validity of the 2DCH is illustrated with a triaxial rotor ($\gamma=-30^\circ$) coupling to one $h_{11/2}$ proton particle and one $h_{11/2}$ neutron hole. By diagonalizing the 2DCH, the angular momenta and energy spectra are obtained. These results agree with the exact solutions of the particle rotor model (PRM) at high rotational frequencies. However, at low frequencies, the energies given by the 2DCH are larger than those by the PRM due to the underestimation of the mass parameters. In addition, with increasing angular momentum, the transitions from the chiral vibration to chiral rotation and further to longitudinal wobbling motion have been presented in the 2DCH.
  • The self-consistent quasiparticle random-phase approximation (QRPA) approach is formulated in the canonical single-nucleon basis of the relativistic Hatree-Fock-Bogoliubov (RHFB) theory. This approach is applied to study the isobaric analog states (IAS) and Gamov-Teller resonances (GTR) by taking Sn isotopes as examples. It is found that self-consistent treatment of the particle-particle residual interaction is essential to concentrate the IAS in a single peak for open-shell nuclei and the Coulomb exchange term is very important to predict the IAS energies. For the GTR, the isovector pairing can increase the calculated GTR energy, while the isoscalar pairing has an important influence on the low-lying tail of the GT transition. Furthermore, the QRPA approach is employed to predict nuclear $\beta$-decay half-lives. With an isospin-dependent pairing interaction in the isoscalar channel, the RHFB+QRPA approach almost completely reproduces the experimental $\beta$-decay half-lives for nuclei up to the Sn isotopes with half-lives smaller than one second. Large discrepancies are found for the Ni, Zn, and Ge isotopes with neutron number smaller than $50$, as well as the Sn isotopes with neutron number smaller than $82$. The potential reasons for these discrepancies are discussed in detail.
  • We present a beyond mean-field study of clusters and molecular structures in low-spin states of $^{20}$Ne with a multireference relativistic energy density functional, where the dynamical correlation effects of symmetry restoration and quadrupole-octupole shapes fluctuation are taken into account with projections on parity, particle number and angular momentum in the framework of the generator coordinate method. Both the energy spectrum and the electric multipole transition strengths for low-lying parity-doublet bands are better reproduced after taking into account the dynamical octupole vibration effect. Consistent with the finding in previous studies, a rotation-induced dissolution of the $\alpha+^{16}$O molecular structure in $^{20}$Ne is predicted.
  • For the first time, the tilted axis cranking covariant density functional theory with pairing correlations has been formulated and implemented in a fully self-consistent and microscopic way to investigate the evolution of the spin axis and the pairing effects in rotating triaxial nuclei. The measured energy spectrum and transition probabilities for the Nd-135 yrast band are reproduced well without any ad hoc renormalization factors when pairing effects are taken into account. A transition from collective to chiral rotation has been demonstrated. It is found that pairing correlations introduce additional admixtures in the single-particle orbitals, and, thus, influence the structure of tilted axis rotating nuclei by reducing the magnitude of the proton and neutron angular momenta while merging their direction.
  • The self-consistent tilted axis cranking covariant density functional theory based on the point- coupling interaction PC-PK1 is applied to investigate the possible existence of antimagnetic ro- tation in the nucleus Fe-58. The observed data for Band 3 and Band 4 are reproduced well with two assigned configurations. It is found that both bands correspond to a rotation of antimagnetic character, but, due to the presence of considerable deformation, the interplay between antimag- netic rotation and collective motion plays an essential role. In particular for Band 4, collective rotation is dominant in the competition with antimagnetic rotation. Moreover, it is shown that the behavior of the ratios between the dynamic moments of inertia and the B(E2) values reflects the interplay between antimagnetic and collective rotation.
  • The anomalous rod shape in carbon isotopes has been investigated in the framework of the cranking covariant density functional theory, and two mechanisms to stabilize such a novel shape with respect to the bending motion, extreme spin, and isospin, are simultaneously discussed for the first time in a self-consistent and microscopic way. By adding valence neutrons and rotating the system, we have found the mechanism stabilizing the rod shape; i.e., the $\sigma$ orbitals (parallel to the symmetry axis) of the valence neutrons, important for the rod shape, are lowered by the rotation due to the Coriolis term. The spin and isospin effects enhance the stability of the rod-shaped configuration. This provides a strong hint that a rod shape could be realized in nuclei towards extreme spin and isospin.
  • 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 from $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.
  • We report a systematic study of nuclear matrix elements (NMEs) in neutrinoless double-beta decays with a state-of-the-art beyond mean-field covariant density functional theory. The dynamic effects of particle-number and angular-momentum conservations as well as quadrupole shape fluctuations are taken into account with projections and generator coordinate method for both initial and final nuclei. The full relativistic transition operator is adopted to calculate the NMEs. The present systematic studies show that in most of the cases there is a much better agreement with the previous non-relativistic calculation based on the Gogny force than in the case of the nucleus $^{150}$Nd found in Song et al. [Phys. Rev. C 90, 054309 (2014)]. In particular, we find that the total NMEs can be well approximated by the pure axial-vector coupling term with a considerable reduction of the computational effort.
  • Neutrinoless double-$\beta$ ($0\nu\beta\beta$) decay is related to many fundamental concepts in nuclear and particle physics beyond the standard model. Currently there are many experiments searching for this weak process. An accurate knowledge of the nuclear matrix element for the $0\nu\beta\beta$ decay is essential for determining the effective neutrino mass once this process is eventually measured. We report the first full relativistic description of the $0\nu\beta\beta$ decay matrix element based on a state-of-the-art nuclear structure model. We adopt the full relativistic transition operators which are derived with the charge-changing nucleonic currents composed of the vector coupling, axial-vector coupling, pseudoscalar coupling, and weak-magnetism coupling terms. The wave functions for the initial and final nuclei are determined by the multireference covariant density functional theory (MR-CDFT) based on the point-coupling functional PC-PK1. The low-energy spectra and electric quadrupole transitions in ${}^{150}$Nd and its daughter nucleus ${}^{150}$Sm are well reproduced by the MR-CDFT calculations. The $0\nu\beta\beta$ decay matrix elements for both the $0_1^+\rightarrow 0_1^+$ and $0_1^+\rightarrow 0_2^+$ decays of ${}^{150}$Nd are evaluated. The effects of particle number projection, static and dynamic deformations, and the full relativistic structure of the transition operators on the matrix elements are studied in detail. The resulting $0\nu\beta\beta$ decay matrix element for the $0_1^+\rightarrow 0_1^+$ transition is $5.60$, which gives the most optimistic prediction for the next generation of experiments searching for the $0\nu\beta\beta$ decay in ${}^{150}$Nd.
  • Relativistic mean field theory is formulated with the Green's function method in coordinate space to investigate the single-particle bound states and resonant states on the same footing. Taking the density of states for free particle as a reference, the energies and widths of single-particle resonant states are extracted from the density of states without any ambiguity. As an example, the energies and widths for single-neutron resonant states in $^{120}$Sn are compared with those obtained by the scattering phase-shift method, the analytic continuation in the coupling constant approach, the real stabilization method and the complex scaling method. Excellent agreements are found for the energies and widths of single-neutron resonant states.
  • In this talk we address two topics: The first one is an empirical explanation in terms of a new state $h_1$ of the peak in the $K^{*0}\bar{K}^{*0}$ invariant mass distribution close to threshold of this channel in the $J/\psi \to \eta K^{*0}\bar{K}^{*0}$ decay. The second one is a theoretical description of the isospin $I=1$ $Z_c(3900)$ and $Z_c(4020)$ states in terms of molecular states of $D \bar D^*+ cc$ and $D^* \bar D^*$.
  • A study of 4$\alpha$ linear-chain structure in high-lying collective excitation states of $^{16}$O with a covariant density functional theory is presented. The low-spin states are obtained by configuration mixing of particle-number and angular-momentum projected quadrupole deformed mean-field states with generator coordinate method. The high-spin states are determined by cranking calculations. These two calculations are based on the same energy density functional PC-PK1. We have found a rotational band at low-spin with the dominated intrinsic configuration considered to be the one that 4$\alpha$ clusters stay along a common axis. The strongly deformed rod shape also appears in the high-spin region with the angular momentum $13-18\hbar$; however whether the state is pure $4\alpha$ linear chain or not is less obvious than that in the low-spin states.