• We present a systematic study of neutron-proton scattering in Nuclear Lattice Effective Field Theory (NLEFT), in terms of the computationally efficient radial Hamiltonian method. Our leading-order (LO) interaction consists of smeared, local contact terms and static one-pion exchange. We show results for a fully non-perturbative analysis up to next-to-next-to-leading 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 low-energy constants of the next-to-leading-order (NLO) two-nucleon 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 partial-wave analysis for S-wave and P-wave channels. We expect the peripheral partial waves to be equally well described once the lattice momenta in the pion-nucleon coupling are taken to coincide with the continuum dispersion relation, and higher-order (N3LO) contributions are included. We stress that for center-of-mass momenta below 100 MeV, the physics of the two-nucleon 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 well-known 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 model-independent measures that probe three- and four-nucleon 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 long-standing deficiency of auxiliary-field 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 well-known Hoyle state resonance in 12C.
  • We develop a multidimensionally constrained relativistic Hartree-Bogoliubov (MDC-RHB) 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 self-consistently. The RHB equation is solved in an axially deformed harmonic oscillator basis. A separable pairing force of finite range is adopted in the MDC-RHB model. The potential energy curves of neutron-rich even-even Zr isotopes are calculated with relativistic functionals DD-PC1 and PC-PK1 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 DD-PC1. 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 pear-like 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 zero-temperature 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 first-order transition at zero temperature from a Bose-condensed gas of alpha particles (4He nuclei) to a nuclear liquid. Whether one has an alpha-particle gas or nuclear liquid is determined by the strength of the alpha-alpha interactions, and we show that the alpha-alpha interactions depend on the strength and locality of the nucleon-nucleon 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, non-cubic 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 complex-valued auxiliary potential that breaks time-reversal invariance. We benchmark our method using a system of two spin-1/2 particles interacting through a finite-range 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 (least-action) spontaneous fission paths and half-lives. Collective potentials and nonperturbative cranking collective inertia tensors are calculated using the multidimensionally-constrained relativistic mean-field (MDC-RMF) 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 particle-number 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 half-lives computed in the three-dimensional collective space of shape and pairing coordinates, using the relativistic functional DD-PC1 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 half-life.
  • 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 half-lives predicted by the functionals PC-PK1 and DD-PC1 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 multidimensionally-constrained relativistic Hartree-Bogoliubov (MDC-RHB) model, with the energy density functionals PC-PK1 and DD-PC1. Pairing correlations are treated in the Bogoliubov approximation using a separable pairing force of finite range. The least-action 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 non-zero 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 Clebsch-Gordan 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 two-humped barrier structure. At large deformations beyond the second barrier the occurrence of a third one was predicted by Mic-Mac model calculations in the 1970s, but contradictory results were later reported. In this paper, triple-humped barriers in actinide nuclei are investigated with covariant density functional theory (CDFT). Calculations are performed using the multidimensionally-constrained relativistic mean field (MDC-RMF) model, with functionals PC-PK1 and DD-ME2. Pairing correlations are treated in the BCS approximation with a separable pairing force of finite range. Two-dimensional 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 one-dimensional 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. DD-ME2 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 PC-PK1 a third barrier is found only in $^{226,228,230}$Th. Single-nucleon 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 DD-ME2 functional, as compared to the functional PC-PK1. 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 finite-volume 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 non-zero 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 multidimensionally-constrained relativistic mean field (MDC-RMF) models which can accommodate various shape degree of freedom. %The hyperon-nucleon $\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 low-lying quadrupole and octupole collective states is presented, based on the microscopic energy density functional framework. By mapping the deformation constrained self-consistent axially symmetric mean-field 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 DD-PC1. The resulting IBM Hamiltonian is used to calculate excitation spectra and transition rates for the positive- and negative-parity 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 low-lying negative-parity 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 (MDC-RMF) 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 self-consistently. 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 point-coupling nucleon interactions combined with the nonlinear or density-dependent couplings. The pairing effects are taken into account with the BCS approach. The 1-, 2-, and 3-D 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 even-even 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 macroscopic-microscopic and MDC-RMF models, possible discontinuities on PES's from self-consistent approaches, and the restoration of broken symmetries. MDC-RMF models give reasonably good description of fission barriers of even-even 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 multi-dimensional constrained relativistic mean field (MDC-RMF) models. Recently we developed multi-dimensional constrained covariant density functional theories (MDC-CDFT) 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 point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. In MDC-RMF models, the pairing correlations are treated with the BCS method. With MDC-RMF models, the potential energy surfaces of even-even actinide nuclei were investigated and the effect of triaxiality on the fission barriers in these nuclei was discussed. The non-axial reflection-asymmetric $\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 multi-dimensional constrained covariant density functional theories (MDC-CDFT) 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 point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. For the pp channel, either the BCS approach or the Bogoliubov transformation is implemented. The MDC-CDFTs with the BCS approach for the pairing (in the following labelled as MDC-RMF models with RMF standing for "relativistic mean field") have been applied to investigate multi-dimensional potential energy surfaces and the non-axial octupole $Y_{32}$-correlations in N=150 isotones. In this contribution we present briefly the formalism of MDC-RMF 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.
  • Multi-dimensional 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 point-coupling nucleon interactions combined with the non-linear or density-dependent 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 non-axial reflection-asymmetric \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 square-well and Woods-Saxon potentials.
  • The non-axial reflection-asymmetric $\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 density-dependent point coupling covariant density functional theory with the parameter set DD-PC1 in the particle-hole channel, it is found that, for the ground states of $^{248}$Cf and $^{250}$Fm, the non-axial 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 non-axial 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 non-axial 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 square-well and Woods-Saxon 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 multi-dimensional 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 Skyrme-Hartree-Fock 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.