Density Functional Theory studies of cluster states in nuclei

The framework of nuclear energy density functionals is applied to a study of the formation and evolution of cluster states in nuclei. The relativistic functional DD-ME2 is used in triaxial and reflection-asymmetric relativistic Hartree-Bogoliubov calculations of relatively light $N = Z$ and neutron-rich nuclei. The role of deformation and degeneracy of single-nucleon states in the formation of clusters is analysed, and interesting cluster structures are predicted in excited configurations of Be, C, O, Ne, Mg, Si, S, Ar and Ca $N = Z$ nuclei. Cluster phenomena in neutron-rich nuclei are discussed, and it is shown that in neutron-rich Be and C nuclei cluster states occur as a result of molecular bonding of $alpha$-particles by the excess neutrons, and also that proton covalent bonding can occur in $^{10}$C.

DIRHB -- a relativistic self-consistent mean-field framework for atomic nuclei

The DIRHB package consists of three Fortran computer codes for the calculation of the ground-state properties of even-even atomic nuclei using the framework of relativistic self-consistent mean-field models. Each code corresponds to a particular choice of spatial symmetry: the DIRHBS, DIRHBZ and DIRHBT codes are used to calculate nuclei with spherical symmetry, axially symmetric quadrupole deformation, and triaxial quadrupole shapes, respectively. Reflection symmetry is assumed in all three cases. The latest relativistic nuclear energy density functionals are implemented in the codes, thus enabling efficient and accurate calculations over the entire nuclide chart.

Microscopic description of octupole shape-phase transitions in light actinides and rare-earth nuclei

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.

An efficient method for computing the Thouless-Valatin inertia parameters

Starting from the adiabatic time-dependent Hartree-Fock approximation (ATDHF), we propose an efficient method to calculate the Thouless-Valatin moments of inertia for the nuclear system. The method is based on the rapid convergence of the expansion of the inertia matrix. The accuracy of the proposed method is verified in the rotational case by comparing the results with the exact Thouless-Valatin moments of inertia calculated using the self-consistent cranking model. The proposed method is computationally much more efficient than the full ATDHF calculation, yet it retains a high accuracy of the order of 1%.

Localization and clustering in the nuclear Fermi liquid

Using the framework of nuclear energy density functionals we examine the conditions for single-nucleon localization and formation of cluster structures in finite nuclei. We propose to characterize localization by the ratio of the dispersion of single-nucleon wave functions to the average inter-nucleon distance. This parameter generally increases with mass and describes the gradual transition from a hybrid phase in light nuclei, characterized by the spatial localization of individual nucleon states that leads to the formation of cluster structures, toward the Fermi liquid phase in heavier nuclei. Values of the localization parameter that correspond to a crystal phase cannot occur in finite nuclei. Typical length and energy scales in nuclei allow the formation of liquid drops, clusters, and halo structures.

Relativistic Energy Density Functional Description of Shape Transition in Superheavy Nuclei

Relativistic energy density functionals (REDF) provide a complete and accurate, global description of nuclear structure phenomena. A modern semi-empirical functional, adjusted to the nuclear matter equation of state and to empirical masses of deformed nuclei, is applied to studies of shapes of superheavy nuclei. The theoretical framework is tested in a comparison of calculated masses, quadrupole deformations, and potential energy barriers to available data on actinide isotopes. Self-consistent mean-field calculations predict a variety of spherical, axial and triaxial shapes of long-lived superheavy nuclei, and their alpha-decay energies and half-lives are compared to data. A microscopic, REDF-based, quadrupole collective Hamiltonian model is used to study the effect of explicit treatment of collective correlations in the calculation of Q{alpha} values and half-lives.

Microscopic Analysis of Nuclear Quantum Phase Transitions in the N=90 region

The analysis of shape transitions in Nd isotopes, based on the framework of relativistic energy density functionals and restricted to axially symmetric shapes in Ref. cite{PRL99}, is extended to the region $Z = 60$, 62, 64 with $N approx 90$, and includes both $beta$ and $gamma$ deformations. Collective excitation spectra and transition probabilities are calculated starting from a five-dimensional Hamiltonian for quadrupole vibrational and rotational degrees of freedom, with parameters determined by constrained self-consistent relativistic mean-field calculations for triaxial shapes. The results reproduce available data, and show that there is an abrupt change of structure at N=90 that can be approximately characterized by the X(5) analytic solution at the critical point of the first-order quantum phase transition between spherical and axially deformed shapes.

Relativistic QRPA calculation of muon capture rates

The relativistic proton-neutron quasiparticle random phase approximation (PN-RQRPA) is applied in the calculation of total muon capture rates on a large set of nuclei from $^{12}$C to $^{244}$Pu, for which experimental values are available. The microscopic theoretical framework is based on the Relativistic Hartree-Bogoliubov (RHB) model for the nuclear ground state, and transitions to excited states are calculated using the PN-RQRPA. The calculation is fully consistent, i.e., the same interactions are used both in the RHB equations that determine the quasiparticle basis, and in the matrix equations of the PN-RQRPA. The calculated capture rates are sensitive to the in-medium quenching of the axial-vector coupling constant. By reducing this constant from its free-nucleon value $g_A = 1.262$ by 10% for all multipole transitions, the calculation reproduces the experimental muon capture rates to better than 10% accuracy.

Beyond the relativistic mean-field approximation (III): collective Hamiltonian in five dimensions

The framework of relativistic energy density functionals is extended to include correlations related to restoration of broken symmetries and fluctuations of collective variables. A new implementation is developed for the solution of the eigenvalue problem of a five-dimensional collective Hamiltonian for quadrupole vibrational and rotational degrees of freedom, with parameters determined by constrained self-consistent relativistic mean-field calculations for triaxial shapes. The model is tested in a series of illustrative calculations of potential energy surfaces and the resulting collective excitation spectra and transition probabilities of the chain of even-even gadolinium isotopes.

Relativistic Nuclear Energy Density Functionals: adjusting parameters to binding energies

We study a particular class of relativistic nuclear energy density functionals in which only nucleon degrees of freedom are explicitly used in the construction of effective interaction terms. Short-distance (high-momentum) correlations, as well as intermediate and long-range dynamics, are encoded in the medium (nucleon density) dependence of the strength functionals of an effective interaction Lagrangian. Guided by the density dependence of microscopic nucleon self-energies in nuclear matter, a phenomenological ansatz for the density-dependent coupling functionals is accurately determined in self-consistent mean-field calculations of binding energies of a large set of axially deformed nuclei. The relationship between the nuclear matter volume, surface and symmetry energies, and the corresponding predictions for nuclear masses is analyzed in detail. The resulting best-fit parametrization of the nuclear energy density functional is further tested in calculations of properties of spherical and deformed medium-heavy and heavy nuclei, including binding energies, charge radii, deformation parameters, neutron skin thickness, and excitation energies of giant multipole resonances.