
We present a systematic study of neutronproton scattering in Nuclear Lattice
Effective Field Theory (NLEFT), in terms of the computationally efficient
radial Hamiltonian method. Our leadingorder (LO) interaction consists of
smeared, local contact terms and static onepion exchange. We show results for
a fully nonperturbative analysis up to nexttonexttoleading order (NNLO),
followed by a perturbative treatment of contributions beyond LO. The latter
analysis anticipates practical Monte Carlo simulations of heavier nuclei. We
explore how our results depend on the lattice spacing a, and estimate sources
of uncertainty in the determination of the lowenergy constants of the
nexttoleadingorder (NLO) twonucleon force. We give results for lattice
spacings ranging from a = 1.97 fm down to a = 0.98 fm, and discuss the effects
of lattice artifacts on the scattering observables. At a = 0.98 fm, lattice
artifacts appear small, and our NNLO results agree well with the Nijmegen
partialwave analysis for Swave and Pwave channels. We expect the peripheral
partial waves to be equally well described once the lattice momenta in the
pionnucleon coupling are taken to coincide with the continuum dispersion
relation, and higherorder (N3LO) contributions are included. We stress that
for centerofmass momenta below 100 MeV, the physics of the twonucleon system
is independent of the lattice spacing.

Nuclear clustering describes the appearance of structures resembling smaller
nuclei such as alpha particles (4He nuclei) within the interior of a larger
nucleus. While clustering is important for several wellknown examples, little
is known about the general nature of clustering in nuclei. In this letter we
present lattice Monte Carlo calculations based on chiral effective field theory
for the ground states of helium, beryllium, carbon, and oxygen isotopes. By
computing modelindependent measures that probe three and fournucleon
correlations at short distances, we determine the effective number of alpha
clusters in any nucleus as well as their shape compared to alpha particles in
vacuum. We also introduce a new computational approach called the pinhole
algorithm, which solves a longstanding deficiency of auxiliaryfield Monte
Carlo simulations in computing density correlations relative to the center of
mass. We use the pinhole algorithm to determine the proton and neutron density
distributions and the geometry of cluster correlations in 12C, 14C, and 16C.
The structural similarities among the carbon isotopes suggest that 14C and 16C
have excitations analogous to the wellknown Hoyle state resonance in 12C.

We develop a multidimensionally constrained relativistic HartreeBogoliubov
(MDCRHB) model in which the pairing correlations are taken into account by
making the Bogoliubov transformation. In this model, the nuclear shape is
assumed to be invariant under the reversion of $x$ and $y$ axes; i.e., the
intrinsic symmetry group is $V_4$ and all shape degrees of freedom
$\beta_{\lambda\mu}$ with even $\mu$ are included selfconsistently. The RHB
equation is solved in an axially deformed harmonic oscillator basis. A
separable pairing force of finite range is adopted in the MDCRHB model. The
potential energy curves of neutronrich eveneven Zr isotopes are calculated
with relativistic functionals DDPC1 and PCPK1 and possible tetrahedral shapes
in the ground and isomeric states are investigated. The ground state shape of
$^{110}$Zr is predicted to be tetrahedral with both functionals and so is that
of $^{112}$Zr with the functional DDPC1. The tetrahedral ground states are
caused by large energy gaps around $Z=40$ and $N=70$ when $\beta_{32}$
deformation is included. Although the inclusion of the $\beta_{30}$ deformation
can also reduce the energy around $\beta_{20}=0$ and lead to minima with
pearlike shapes for nuclei around $^{110}$Zr, these minima are unstable due to
their shallowness.

How do protons and neutrons bind to form nuclei? This is the central question
of ab initio nuclear structure theory. While the answer may seem as simple as
the fact that nuclear forces are attractive, the full story is more complex and
interesting. In this work we present numerical evidence from ab initio lattice
simulations showing that nature is near a quantum phase transition, a
zerotemperature transition driven by quantum fluctuations. Using lattice
effective field theory, we perform Monte Carlo simulations for systems with up
to twenty nucleons. For even and equal numbers of protons and neutrons, we
discover a firstorder transition at zero temperature from a Bosecondensed gas
of alpha particles (4He nuclei) to a nuclear liquid. Whether one has an
alphaparticle gas or nuclear liquid is determined by the strength of the
alphaalpha interactions, and we show that the alphaalpha interactions depend
on the strength and locality of the nucleonnucleon interactions. This insight
should be useful in improving calculations of nuclear structure and important
astrophysical reactions involving alpha capture on nuclei. Our findings also
provide a tool to probe the structure of alpha cluster states such as the Hoyle
state responsible for the production of carbon in red giant stars and point to
a connection between nuclear states and the universal physics of bosons at
large scattering length.

We introduce a general and accurate method for determining lattice phase
shifts and mixing angles, which is applicable to arbitrary, noncubic lattices.
Our method combines angular momentum projection, spherical wall boundaries and
an adjustable auxiliary potential. This allows us to construct radial lattice
wave functions and to determine phase shifts at arbitrary energies. For coupled
partial waves, we use a complexvalued auxiliary potential that breaks
timereversal invariance. We benchmark our method using a system of two
spin1/2 particles interacting through a finiterange potential with a strong
tensor component. We are able to extract phase shifts and mixing angles for all
angular momenta and energies, with precision greater than that of extant
methods. We discuss a wide range of applications from nuclear lattice
simulations to optical lattice experiments.

Studies of fission dynamics, based on nuclear energy density functionals,
have shown that the coupling between shape and pairing degrees of freedom has a
pronounced effect on the nonperturbative collective inertia and, therefore, on
dynamic (leastaction) spontaneous fission paths and halflives.
Collective potentials and nonperturbative cranking collective inertia tensors
are calculated using the multidimensionallyconstrained relativistic meanfield
(MDCRMF) model. Pairing correlations are treated in the BCS approximation
using a separable pairing force of finite range. Pairing fluctuations are
included as a collective variable using a constraint on particlenumber
dispersion. Fission paths are determined with the dynamic programming method by
minimizing the action in multidimensional collective spaces.
The dynamics of spontaneous fission of $^{264}$Fm and $^{250}$Fm are
explored. Fission paths, action integrals and corresponding halflives computed
in the threedimensional collective space of shape and pairing coordinates,
using the relativistic functional DDPC1 and a separable pairing force of
finite range, are compared with results obtained without pairing fluctuations.
Results for $^{264}$Fm are also discussed in relation with those recently
obtained using the HFB model.
The inclusion of pairing correlations in the space of collective coordinates
favors axially symmetric shapes along the dynamic path of the fissioning
system, amplifies pairing as the path traverses the fission barriers,
significantly reduces the action integral and shortens the corresponding SF
halflife.

Recent microscopic studies, based on the theoretical framework of nuclear
energy density functionals, have analyzed dynamic (least action) and static
(minimum energy) fission paths, and it has been shown that in addition to the
important role played by nonaxial and/or octupole collective degrees of
freedom, fission paths crucially depend on the approximations adopted in
calculating the collective inertia.
The dynamics of spontaneous fission of $^{264}$Fm and $^{250}$Fm is explored.
The fission paths, action integrals and the corresponding halflives predicted
by the functionals PCPK1 and DDPC1 are compared and, in the case of
$^{264}$Fm, discussed in relation with recent results obtained using the HFB
model based on the Skyrme functional SkM$^*$ and a density dependent mixed
pairing interaction.
Deformation energy surfaces, collective potentials, and perturbative and
nonperturbative cranking collective inertia tensors are calculated using the
multidimensionallyconstrained relativistic HartreeBogoliubov (MDCRHB) model,
with the energy density functionals PCPK1 and DDPC1. Pairing correlations are
treated in the Bogoliubov approximation using a separable pairing force of
finite range. The leastaction principle is employed to determine dynamic
spontaneous fission paths.

We study the breaking of rotational symmetry on the lattice for irreducible
tensor operators and practical methods for suppressing this breaking. We
illustrate the features of the general problem using an $\alpha$ cluster model
for $^{8}$Be. We focus on the lowest states with nonzero angular momentum and
examine the matrix elements of multipole moment operators. We show that the
physical reduced matrix element is well reproduced by averaging over all
possible orientations of the quantum state, and this is expressed as a sum of
matrix elements weighted by the corresponding ClebschGordan coefficients. For
our $\alpha$ cluster model we find that the effects of rotational symmetry
breaking can be largely eliminated for lattice spacings of $a\leq 1.7$ fm, and
we expect similar improvement for actual lattice Monte Carlo calculations.

Potential energy surfaces (PES's) of actinide nuclei are characterized by a
twohumped barrier structure. At large deformations beyond the second barrier
the occurrence of a third one was predicted by MicMac model calculations in
the 1970s, but contradictory results were later reported. In this paper,
triplehumped barriers in actinide nuclei are investigated with covariant
density functional theory (CDFT). Calculations are performed using the
multidimensionallyconstrained relativistic mean field (MDCRMF) model, with
functionals PCPK1 and DDME2. Pairing correlations are treated in the BCS
approximation with a separable pairing force of finite range. Twodimensional
PES's of $^{226,228,230,232}$Th and $^{232,234,236,238}$U are mapped and the
third minima on these surfaces are located. Then onedimensional potential
energy curves along the fission path are analyzed in detail and the energies of
the second barrier, the third minimum, and the third barrier are determined.
DDME2 predicts the occurrence of a third barrier in all Th nuclei and
$^{238}$U. The third minima in $^{230,232}$Th are very shallow, whereas those
in $^{226,228}$Th and $^{238}$U are quite prominent. With PCPK1 a third
barrier is found only in $^{226,228,230}$Th. Singlenucleon levels around the
Fermi surface are analyzed in $^{226}$Th, and it is found that the formation of
the third minimum is mainly due to the $Z=90$ proton energy gap at $\beta_{20}
\approx 1.5$ and $\beta_{30} \approx 0.7$. The possible occurrence of a third
barrier in actinide nuclei depends on the effective interaction used in
multidimensional CDFT calculations. More pronounced minima are predicted by the
DDME2 functional, as compared to the functional PCPK1. The depth of the third
well in Th isotopes decreases with increasing neutron number. The origin of the
third minimum is due to the proton $Z=90$ shell gap at relevant deformations.

We explore the breaking of rotational symmetry on the lattice for bound state
energies and practical methods for suppressing this breaking. We demonstrate
the general problems associated with lattice discretization errors and
finitevolume errors using an $\alpha$ cluster model for $^8$Be and $^{12}$C.
We consider the two and three $\alpha$particle systems and focus on the lowest
states with nonzero angular momentum which split into multiplets corresponding
to different irreducible representations of the cubic group. We examine the
dependence of such splittings on the lattice spacing and box size. We find that
lattice spacing errors are closely related to the commensurability of the
lattice with the intrinsic length scales of the system. We also show that
rotational symmetry breaking effects can be significantly reduced by using
improved lattice actions, and that the physical energy levels are accurately
reproduced by the weighted average of a given spin multiplets.

We study the superdeformed (SD) states and corresponding SD hypernuclei of Ar
isotopes with the multidimensionallyconstrained relativistic mean field
(MDCRMF) models which can accommodate various shape degree of freedom. %The
hyperonnucleon $\sigma$$\Lambda$ and $\omega$$\Lambda $ interactions are
introduced %in the same covariant framework. We found that the density profiles
of SD states in Ar isotopes show a strong localization with a ring structure
near the surface, while the central part of the density is dilute showing a
hole structure. This localization of SD density induces an appreciable
deformation in the hyperon wave function and results in a large overlap between
the core and the hyperon in the SD hypernuclei of Ar isotopes. Then the
$\Lambda$ separation energy of SD state becomes larger than that of normally
deformed or spherical ground state. This feature is different from that found
in other nuclei such as $^{32}$S, $^{56}$Ni, and $^{60}$Zn in which the
$\Lambda$ separation energy of larger deformed state is smaller. In this
context, the measurement of the $\Lambda$ separation energy may provide an
important information on the localization of the density profile of SD states.

A systematic analysis of lowlying quadrupole and octupole collective states
is presented, based on the microscopic energy density functional framework. By
mapping the deformation constrained selfconsistent axially symmetric
meanfield energy surfaces onto the equivalent Hamiltonian of the $sdf$
interacting boson model (IBM), that is, onto the energy expectation value in
the boson condensate state, the Hamiltonian parameters are determined. The
study is based on the global relativistic energy density functional DDPC1. The
resulting IBM Hamiltonian is used to calculate excitation spectra and
transition rates for the positive and negativeparity collective states in
four isotopic chains characteristic for two regions of octupole deformation and
collectivity: Th, Ra, Sm and Ba. Consistent with the empirical trend, the
microscopic calculation based on the systematics of $\beta_{2}$$\beta_{3}$
energy maps, the resulting lowlying negativeparity bands and transition rates
show evidence of a shape transition between stable octupole deformation and
octupole vibrations characteristic for $\beta_{3}$soft potentials.

By breaking both the axial and the spatial reflection symmetries, we develop
multidimensionally constrained relativistic mean field (MDCRMF) models. The
nuclear shape is assumed to be invariant under the reversion of $x$ and $y$
axes, i.e., the intrinsic symmetry group is $V_{4}$ and all shape degrees of
freedom $\beta_{\lambda\mu}$ with even $\mu$, such as $\beta_{20}$,
$\beta_{22}$, $\beta_{30}$, $\beta_{32}$, $\beta_{40}$, $\dots$, are included
selfconsistently. The Dirac wave functions are expanded in an axially deformed
harmonic oscillator basis. The RMF functional can be one of the following four
forms: the meson exchange or pointcoupling nucleon interactions combined with
the nonlinear or densitydependent couplings. The pairing effects are taken
into account with the BCS approach. The 1, 2, and 3D potential energy
surfaces (PES's) of $^{240}$Pu are illustrated for numerical checks and for the
study of the effect of the triaxiality on the fission barriers. Potential
energy curves of eveneven actinide nuclei around the first and second fission
barriers are studied systematically. Besides the first ones, the second fission
barriers in these nuclei are also lowered considerably by the triaxial
deformation. This lowering effect is independent of the effective interactions
used in the RMF functionals. Further discussions are made about different
predictions on the effect of the triaxiality between the
macroscopicmicroscopic and MDCRMF models, possible discontinuities on PES's
from selfconsistent approaches, and the restoration of broken symmetries.
MDCRMF models give reasonably good description of fission barriers of
eveneven actinide nuclei. It is important to include both the nonaxial and the
reflection asymmetric shapes simultaneously for the study of PES's and fission
barriers of actinide nuclei and of those in unknown mass regions such as, e.g.,
superheavy nuclei.

In this contribution we present some results of potential energy surfaces of
actinide and transfermium nuclei from multidimensional constrained
relativistic mean field (MDCRMF) models. Recently we developed
multidimensional constrained covariant density functional theories (MDCCDFT)
in which all shape degrees of freedom $\beta_{\lambda\mu}$ with even $\mu$ are
allowed and the functional can be one of the following four forms: the meson
exchange or pointcoupling nucleon interactions combined with the nonlinear or
densitydependent couplings. In MDCRMF models, the pairing correlations are
treated with the BCS method. With MDCRMF models, the potential energy surfaces
of eveneven actinide nuclei were investigated and the effect of triaxiality on
the fission barriers in these nuclei was discussed. The nonaxial
reflectionasymmetric $\beta_{32}$ shape in some transfermium nuclei with
$N=150$, namely $^{246}$Cm, $^{248}$Cf, $^{250}$Fm, and $^{252}$No were also
studied.

The pseudospin symmetry (PSS) is a relativistic dynamical symmetry directly
connected with the small component of the nucleon Dirac wave function. Much
effort has been made to study this symmetry in bound states. Recently, a
rigorous justification of the PSS in single particle resonant states was
achieved by examining the asymptotic behaviors of the radial Dirac wave
functions: The PSS in single particle resonant states in nuclei is conserved
exactly when the attractive scalar and repulsive vector potentials have the
same magnitude but opposite sign. Several issues related to the exact
conservation and breaking mechanism of the PSS in single particle resonances
were investigated by employing spherical square well potentials in which the
PSS breaking part can be well isolated in the Jost function. A threshold effect
in the energy splitting and an anomaly in the width splitting of pseudospin
partners were found when the depth of the square well potential varies from
zero to a finite value.

The pseudospin symmetry (PSS) has been studied extensively for bound states.
Recently we justified rigorously that the PSS in single particle resonant
states is exactly conserved when the attractive scalar and repulsive vector
potentials of the Dirac Hamiltonian have the same magnitude but opposite sign
[PRL 109, 072501 (2012)]. To understand more deeply the PSS, we focus on
several issues related to the exact conservation and breaking mechanism of the
PSS in single particle resonances. In particular, we are interested in how the
energy and width splittings of PS partners depend on the depth of the scalar
and vector potentials. We investigate the asymptotic behaviors of radial Dirac
wave functions. Spherical square well potentials are employed in which the PSS
breaking part in the Jost function can be well isolated. By examining the zeros
of Jost functions corresponding to small components of the radial Dirac wave
functions, general properties of the PSS are analyzed. By examining the Jost
function, the occurrence of intruder orbitals is explained and it is possible
to trace continuously the PSS partners from the PSS limit to the case with a
finite potential depth. The dependence of the PSS in resonances as well as in
bound states on the potential depth is investigated systematically. We find a
threshold effect in the energy splitting and an anomaly in the width splitting
of pseudospin partners when the depth of the single particle potential varies
from zero to a finite value. The conservation and the breaking of the PSS in
resonant states and bound states share some similar properties. The appearance
of intruder states can be explained by examining the zeros of Jost functions.
Origins of the threshold effect in the energy splitting and the anomaly in the
width splitting of PS partners, together with many other problems, are still
open and should be further investigated.

We have developed multidimensional constrained covariant density functional
theories (MDCCDFT) for finite nuclei in which the shape degrees of freedom
\beta_{\lambda\mu} with even \mu, e.g., \beta_{20}, \beta_{22}, \beta_{30},
\beta_{32}, \beta_{40}, etc., can be described simultaneously. The functional
can be one of the following four forms: the meson exchange or pointcoupling
nucleon interactions combined with the nonlinear or densitydependent
couplings. For the pp channel, either the BCS approach or the Bogoliubov
transformation is implemented. The MDCCDFTs with the BCS approach for the
pairing (in the following labelled as MDCRMF models with RMF standing for
"relativistic mean field") have been applied to investigate multidimensional
potential energy surfaces and the nonaxial octupole $Y_{32}$correlations in
N=150 isotones. In this contribution we present briefly the formalism of
MDCRMF models and some results from these models. The potential energy
surfaces with and without triaxial deformations are compared and it is found
that the triaxiality plays an important role upon the second fission barriers
of actinide nuclei. In the study of Y_{32}correlations in N=150 isotones, it
is found that, for 248Cf and 250Fm, \beta_{32} > 0.03 and the energy is lowered
by the \beta_{32} distortion by more than 300 keV; while for 246Cm and 252No,
the pocket with respect to \beta_{32} is quite shallow.

Multidimensional constrained covariant density functional theories were
developed recently. In these theories, all shape degrees of freedom
\beta_{\lambda\mu} deformations with even \mu are allowed, e.g., \beta_{20},
\beta_{22}, \beta_{30}, \beta_{32}, \beta_{40}, \beta_{42}, \beta_{44}, and so
on and the CDFT functional can be one of the following four forms: the meson
exchange or pointcoupling nucleon interactions combined with the nonlinear or
densitydependent couplings. In this contribution, some applications of these
theories are presented. The potential energy surfaces of actinide nuclei in the
(\beta_{20}, \beta_{22}, \beta_{30}) deformation space are investigated. It is
found that besides the octupole deformation, the triaxiality also plays an
important role upon the second fission barriers. The nonaxial
reflectionasymmetric \beta_{32} shape in some transfermium nuclei with N =
150, namely 246Cm, 248Cf, 250Fm, and 252No are studied.

In this contribution we present some results on the study of pseudospin
symmetry (PSS) in single particle resonant states. The PSS is a relativistic
dynamical symmetry connected with the small component of the nucleon Dirac wave
function. Many efforts have been made to study this symmetry in bound states.
We recently gave a rigorous justification of the PSS in single particle
resonant states by examining the zeros of Jost functions corresponding to the
small components of the radial Dirac wave functions and phase shifts of
continuum states [1, Phys. Rev. Lett. 109 (2012) 072501]. We have shown that
the PSS in single particle resonant states in nuclei is conserved when the
attractive scalar and repulsive vector potentials have the same magnitude but
opposite sign. Examples of exact conservation and breaking of this symmetry in
single particle resonances are given for spherical squarewell and WoodsSaxon
potentials.

The nonaxial reflectionasymmetric $\beta_{32}$ shape in some transfermium
nuclei with N=150, namely $^{246}$Cm, $^{248}$Cf, $^{250}$Fm, and $^{252}$No
are investigated with multidimensional constrained covariant density functional
theories. By using the densitydependent point coupling covariant density
functional theory with the parameter set DDPC1 in the particlehole channel,
it is found that, for the ground states of $^{248}$Cf and $^{250}$Fm, the
nonaxial octupole deformation parameter $\beta_{32} > 0.03$ and the energy
gain due to the $\beta_{32}$ distortion is larger than 300 keV. In $^{246}$Cm
and $^{252}$No, shallow $\beta_{32}$ minima are found. The occurrence of the
nonaxial octupole $\beta_{32}$ correlations is mainly from a pair of neutron
orbitals $[734]9/2$ ($\nu j_{15/2}$) and $[622]5/2$ ($\nu g_{9/2}$) which are
close to the neutron Fermi surface and a pair of proton orbitals $[521]3/2$
($\pi f_{7/2}$) and $[633]7/2$ ($\pi i_{13/2}$) which are close to the proton
Fermi surface. The dependence of the nonaxial octupole effects on the form of
energy density functional and on the parameter set is also studied.

The pseudospin symmetry is a relativistic dynamical symmetry connected with
the small component of the Dirac spinor. The origin of pseudospin symmetry in
single particle bound states in atomic nuclei has been revealed and studied
extensively. By examining the zeros of Jost functions corresponding to the
small components of Dirac wave functions and phase shifts of continuum states,
we show that the pseudospin symmetry in single particle resonant states in
nuclei is conserved when the attractive scalar and repulsive vector potentials
have the same magnitude but opposite sign. The exact conservation and the
breaking of pseudospin symmetry are illustrated for single particle resonances
in spherical squarewell and WoodsSaxon potentials.

For the first time the potential energy surfaces of actinide nuclei in the
$(\beta_{20}, \beta_{22}, \beta_{30})$ deformation space are obtained from a
multidimensional constrained covariant density functional theory. With this
newly developed theory we are able to explore the importance of the triaxial
and octupole shapes simultaneously along the whole fission path. It is found
that besides the octupole deformation, the triaxiality also plays an important
role upon the second fission barriers. The outer barrier as well as the inner
barrier are lowered by the triaxial deformation compared with axially symmetric
results. This lowering effect for the reflection asymmetric outer barrier is
0.5 $\sim$ 1 MeV, accounting for $10 \sim 20%$ of the barrier height. With the
inclusion of the triaxial deformation, a good agreement with the data for the
outer barriers of actinide nuclei is achieved.

The shapes of light normal nuclei and $\Lambda$ hypernuclei are investigated
in the $(\beta, \gamma)$ deformation plane by using a newly developed
constrained relativistic mean field (RMF) model. As examples, the results of
some C, Mg, and Si nuclei are presented and discussed in details. We found that
for normal nuclei the present RMF calculations and previous SkyrmeHartreeFock
models predict similar trends of the shape evolution with the neutron number
increasing. But some quantitative aspects from these two approaches, such as
the depth of the minimum and the softness in the $\gamma$ direction, differ a
lot for several nuclei. For $\Lambda$ hypernuclei, in most cases, the addition
of a $\Lambda$ hyperon alters slightly the location of the ground state minimum
towards the direction of smaller $\beta$ and softer $\gamma$ in the potential
energy surface $E \sim (\beta, \gamma)$. There are three exceptions, namely,
$^{13}_\Lambda$C, $^{23}_\Lambda$C, and $^{31}_\Lambda$Si in which the
polarization effect of the additional $\Lambda$ is so strong that the shapes of
these three hypernuclei are drastically different from their corresponding core
nuclei.