
The excitation function for the 7 alpha deexcitation 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 particlehole 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 eveneven nuclei is
generalized to include the vibrational degrees of freedom. The pertinent
Hamiltonian is constructed up to nexttoleading 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
highspin 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 groundstate shape phase transitions in
oddmass Eu isotopes are explored starting from excitation spectra and
collective wave functions obtained by diagonalization of a corequasiparticle
coupling Hamiltonian based on energy density functionals. As functions of the
physical control parameter  the number of nucleons  theoretical lowenergy
spectra, twoneutron 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 eveneven Sm and Gd isotopes.
The enhancement of the firstorder quantum phase transition in oddmass 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 shellmodellike approach is implemented to treat the cranking manybody
Hamiltonian based on the covariant density functional theory including pairing
correlations with exact particle number conservation. The selfconsistency is
achieved by iterating the singleparticle occupation probabilities back to the
densities and currents. As an example, the rotational structures observed in
the neutronrich 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 pseudocrossing 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 groundstate 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 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.

The pairingplusquadrupole Hamiltonian is diagonalized in a symmetry
restored basis, i.e., the triaxial quasiparticlestates 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 eveneven nuclei. A Hamiltonian, called the
triaxial rotor model (TRM), is obtained up to nexttoleading order (NLO)
within the EFT formalism. Its applicability is examined by comparing with a
fivedimensional 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, nonaxial deformed, and octupole
deformed potential are provided and discussed.

Background: The discovery of neutrinoless doublebeta ($0\nu\beta\beta$)
decay would demonstrate the nature of neutrinos, have profound implications for
our understanding of matterantimatter 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 nucleonnucleon
shortrange 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
beyondmeanfield correlations are included in the nuclear wave functions by
configuration mixing of both angularmomentum and particlenumber projected
quadrupole deformed meanfield 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 shortrange correlations are evaluated.
Conclusions: The effects of relativity and shortrange 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 halflives, 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 relativisticmeanfield and SkyrmeHartreeFock 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 (CIPDFT). 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 Cr54. It is
found that the experimentally observed excitation energies for the yrast band
in Cr54 can be well reproduced. In contrast to conventional shellmodel
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 twodimensional 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 threedimensional 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 selfconsistent quasiparticle randomphase approximation (QRPA) approach
is formulated in the canonical singlenucleon basis of the relativistic
HatreeFockBogoliubov (RHFB) theory. This approach is applied to study the
isobaric analog states (IAS) and GamovTeller resonances (GTR) by taking Sn
isotopes as examples. It is found that selfconsistent treatment of the
particleparticle residual interaction is essential to concentrate the IAS in a
single peak for openshell 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 lowlying tail of the GT transition. Furthermore,
the QRPA approach is employed to predict nuclear $\beta$decay halflives. With
an isospindependent pairing interaction in the isoscalar channel, the
RHFB+QRPA approach almost completely reproduces the experimental $\beta$decay
halflives for nuclei up to the Sn isotopes with halflives 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 meanfield study of clusters and molecular structures in
lowspin states of $^{20}$Ne with a multireference relativistic energy density
functional, where the dynamical correlation effects of symmetry restoration and
quadrupoleoctupole 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 lowlying paritydoublet bands are better
reproduced after taking into account the dynamical octupole vibration effect.
Consistent with the finding in previous studies, a rotationinduced 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
selfconsistent 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 Nd135 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 singleparticle 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 selfconsistent tilted axis cranking covariant density functional theory
based on the point coupling interaction PCPK1 is applied to investigate the
possible existence of antimagnetic ro tation in the nucleus Fe58. 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
selfconsistent 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 rodshaped
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 eveneven nuclei with proton
numbers ranging from $Z=8$ to $Z=108$ by solving a fivedimensional collective
Hamiltonian, the collective parameters of which are determined from triaxial
relativistic meanfield plus BCS calculation using the PCPK1 force. After
taking into account these beyond meanfield DCEs, the rootmeansquare (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),
DDME2 (2.39 MeV), DDME$\delta$ (2.29 MeV) and DDPC1 (2.01 MeV). Moreover,
the rms deviation for twonucleon separation energies is reduced by $\sim34\%$
in comparison with cranking prescription.

We report a systematic study of nuclear matrix elements (NMEs) in
neutrinoless doublebeta decays with a stateoftheart beyond meanfield
covariant density functional theory. The dynamic effects of particlenumber and
angularmomentum 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 nonrelativistic
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 axialvector 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 stateoftheart nuclear
structure model. We adopt the full relativistic transition operators which are
derived with the chargechanging nucleonic currents composed of the vector
coupling, axialvector coupling, pseudoscalar coupling, and weakmagnetism
coupling terms. The wave functions for the initial and final nuclei are
determined by the multireference covariant density functional theory (MRCDFT)
based on the pointcoupling functional PCPK1. The lowenergy spectra and
electric quadrupole transitions in ${}^{150}$Nd and its daughter nucleus
${}^{150}$Sm are well reproduced by the MRCDFT 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 singleparticle 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 singleparticle resonant
states are extracted from the density of states without any ambiguity. As an
example, the energies and widths for singleneutron resonant states in
$^{120}$Sn are compared with those obtained by the scattering phaseshift
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 singleneutron 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$ linearchain structure in highlying collective
excitation states of $^{16}$O with a covariant density functional theory is
presented. The lowspin states are obtained by configuration mixing of
particlenumber and angularmomentum projected quadrupole deformed meanfield
states with generator coordinate method. The highspin states are determined by
cranking calculations. These two calculations are based on the same energy
density functional PCPK1. We have found a rotational band at lowspin 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 highspin region with the angular momentum $1318\hbar$; however whether
the state is pure $4\alpha$ linear chain or not is less obvious than that in
the lowspin states.