Source author record

T. Niksic

T. Niksic appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Dynamical synthesis of 4He in the scission phase of nuclear fission

In the exothermic process of fission decay, an atomic nucleus splits into two or more independent fragments. Several aspects of nuclear fission are not properly understood, in particular the formation of the neck between the nascent fragments, and the subsequent mechanism of scission into two or more independent fragments. Using an implementation of time-dependent density functional theory, based on a relativistic energy density functional and including pairing correlations, we analyze the final phase of the process of induced fission of $^{240}$Pu, and show that the time-scale of neck formation coincides with the assembly of two $α$-like clusters (less than 1 zs = 10$^{-21}$ s). Because of its much larger binding energy, the dynamical synthesis of 4He in the neck predominates over other light clusters, e.g., $^3$H and $^6$He. At the instant of scission the neck ruptures exactly between the two $α$-like clusters, which separate because of the Coulomb repulsion and are eventually absorbed by the two emerging fragments. The newly proposed mechanism of light charged clusters formation at scission provides a natural explanation of ternary fission.

preprint2022arXiv

Microscopic analysis of induced nuclear fission dynamics

The dynamics of low-energy induced fission is explored using a consistent microscopic framework that combines the time-dependent generator coordinate method (TDGCM) and time-dependent nuclear density functional theory (TDDFT). While the former presents a fully quantum mechanical approach that describes the entire fission process as an adiabatic evolution of collective degrees of freedom, the latter models the dissipative dynamics of the final stage of fission by propagating the nucleons independently toward scission and beyond. By combining the two methods, based on the same nuclear energy density functional and pairing interaction, we perform an illustrative calculation of the charge distribution of yields and total kinetic energy for induced fission of $^{240}$Pu. For the saddle-to-scission phase a set of initial points for the TDDFT evolution is selected along an iso-energy curve beyond the outer fission barrier on the deformation energy surface, and the TDGCM is used to calculate the probability that the collective wave function reaches these points at different times. Fission observables are computed with both methods and compared with available data. The relative merits of including quantum fluctuations (TDGCM) and the one-body dissipation mechanism (TDDFT) are discussed.

preprint2020arXiv

Coupling of shape and pairing vibrations in a collective Hamiltonian based on nuclear energy density functionals

The quadrupole collective Hamiltonian, based on relativistic energy density functionals, is extended to include a pairing collective coordinate. In addition to quadrupole shape vibrations and rotations, the model describes pairing vibrations and the coupling between shape and pairing degrees of freedom. The parameters of the collective Hamiltonian are determined by constrained self-consistent relativistic mean-field plus Bardeen-Cooper-Schrieffer (RMF+BCS) calculations in the space of intrinsic shape and pairing deformations. The effect of coupling between shape and pairing degrees of freedom is analyzed in a study of low-energy spectra and transition rates of four axially symmetric $N=92$ rare-earth isotones. When compared to results obtained with the standard quadrupole collective Hamiltonian, the inclusion of dynamical pairing increases the moment of inertia, lowers the energies of excited $0^+$ states and reduces the E0-transition strengths, in better agreement with data.

preprint2020arXiv

Microscopic description of the self-conjugate $^{108}$Xe and $^{104}$Te $α$-decay chain

A microscopic calculation of half-lives for the recently observed $^{108}$Xe $\to$ $^{104}$Te $\to$ $^{100}$Sn $α$-decay chain is performed using a self-consistent framework based on energy density functionals. The relativistic density functional DD-PC1 and a separable pairing interaction of finite range are used to compute axially-symmetric deformation energy surfaces of $^{104}$Te and $^{108}$Xe as functions of quadrupole, octupole and hexadecupole collective coordinates. Dynamic least-action paths are determined that trace the $α$-particle emission from the equilibrium deformation to the point of scission. The calculated half-lives: 197 ns for $^{104}$Te and 50 $μ$s for $^{108}$Xe, are compared to recent experimental values of the half-lives of superallowed $α$-decay of $^{104}$Te: $< 18$ ns, and $^{108}$Xe: 58$^{+106}_{-23}$ $μ$s.

preprint2014arXiv

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 $α$-particles by the excess neutrons, and also that proton covalent bonding can occur in $^{10}$C.

preprint2014arXiv

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.

preprint2014arXiv

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 $β_{2}$-$β_{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 $β_{3}$-soft potentials.

preprint2013arXiv

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.

preprint2012arXiv

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%.

preprint2012arXiv

Energy Density Functional analysis of shape evolution in N=28 isotones

The structure of low-energy collective states in proton-deficient N=28 isotones is analyzed using structure models based on the relativistic energy density functional DD-PC1. The relativistic Hartree-Bogoliubov model for triaxial nuclei is used to calculate binding energy maps in the $β$-$γ$ plane. The evolution of neutron and proton single-particle levels with quadrupole deformation, and the occurrence of gaps around the Fermi surface, provide a simple microscopic interpretation of the onset of deformation and shape coexistence. Starting from self-consistent constrained energy surfaces calculated with the functional DD-PC1, a collective Hamiltonian for quadrupole vibrations and rotations is employed in the analysis of excitation spectra and transition rates of $^{46}$Ar, $^{44}$S, and $^{42}$Si. The results are compared to available data, and previous studies based either on the mean-field approach or large-scale shell-model calculations. The present study is particularly focused on $^{44}$S, for which data have recently been reported that indicate pronounced shape coexistence.

preprint2012arXiv

How atomic nuclei cluster

Nucleonic matter displays a quantum liquid structure, but in some cases finite nuclei behave like molecules composed of clusters of protons and neutrons. Clustering is a recurrent feature in light nuclei, from beryllium to nickel. For instance, in $^{12}$C the Hoyle state, crucial for stellar nucleosynthesis, can be described as a nuclear molecule consisting of three alpha-particles. The mechanism of cluster formation, however, has not yet been fully understood. We show that the origin of clustering can be traced back to the depth of the confining nuclear potential. By employing the theoretical framework of energy density functionals that encompasses both cluster and quantum liquid-drop aspects of nuclei, it is shown that the depth of the potential determines the energy spacings between single-nucleon orbitals, the localization of the corresponding wave functions and, therefore, the degree of nucleonic density clustering. Relativistic functionals, in particular, are characterized by deep single-nucleon potentials. When compared to non-relativistic functionals that yield similar ground-state properties (binding energy, deformation, radii), they predict the occurrence of much more pronounced cluster structures. More generally, clustering is considered as a transitional phenomenon between crystalline and quantum liquid phases of fermionic systems.

preprint2012arXiv

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α values and half-lives.

preprint2012arXiv

Robust Regularity in γ-Soft Nuclei and its Microscopic Realization

γ-softness in atomic nuclei is investigated in the framework of energy density functionals. By mapping constrained microscopic energy surfaces for a set of representative non-axial medium-heavy and heavy nuclei to a Hamiltonian of the proton-neutron interacting boson model (IBM-2) containing up to three-body interactions, low-lying collective spectra and transition rates are calculated. Observables are analyzed that distinguish between the two limiting geometrical pictures of non-axial nuclei: the rigid-triaxial rotor and the γ-unstable rotor. It is shown that neither of these pictures is realized in actual nuclei, and that a microscopic description leads to results that are almost exactly in between the two geometrical limits. This finding points to the optimal choice of the IBM Hamiltonian for γ-soft nuclei.

preprint2011arXiv

Modern applications of covariant density functional theory

Modern applications of Covariant Density Functional Theory (CDFT) are discussed. First we show a systematic investigation of fission barriers in actinide nuclei within constraint relativistic mean field theory allowing for triaxial deformations. In the second part we discuss a microscopic theory of quantum phase transitions (QPT) based on the relativistic generator coordinate method.

preprint2011arXiv

Quadrupole Collective Dynamics from Energy Density Functionals: Collective Hamiltonian and the Interacting Boson Model

Microscopic energy density functionals (EDF) have become a standard tool for nuclear structure calculations, providing an accurate global description of nuclear ground states and collective excitations. For spectroscopic applications this framework has to be extended to account for collective correlations related to restoration of symmetries broken by the static mean field, and for fluctuations of collective variables. In this work we compare two approaches to five-dimensional quadrupole dynamics: the collective Hamiltonian for quadrupole vibrations and rotations, and the Interacting Boson Model. The two models are compared in a study of the evolution of non-axial shapes in Pt isotopes. Starting from the binding energy surfaces of $^{192,194,196}$Pt, calculated with a microscopic energy density functional, we analyze the resulting low-energy collective spectra obtained from the collective Hamiltonian, and the corresponding IBM-2 Hamiltonian. The calculated excitation spectra and transition probabilities for the ground-state bands and the $γ$-vibration bands are compared to the corresponding sequences of experimental states.

preprint2010arXiv

Constraints on the inner edge of neutron star crusts from relativistic nuclear energy density functionals

The transition density $n_t$ and pressure $P_t$ at the inner edge between the liquid core and the solid crust of a neutron star are analyzed using the thermodynamical method and the framework of relativistic nuclear energy density functionals. Starting from a functional that has been carefully adjusted to experimental binding energies of finite nuclei, and varying the density dependence of the corresponding symmetry energy within the limits determined by isovector properties of finite nuclei, we estimate the constraints on the core-crust transition density and pressure of neutron stars: $0.086 \ {\rm fm}^{-3} \leq n_t < 0.090 \ {\rm fm}^{-3}$ and $0.3\ {\rm MeV \ fm}^{-3} < P_t \leq 0.76 \ {\rm MeV \ fm}^{-3}$.

preprint2010arXiv

Microscopic Description of Spherical to Gamma-Soft Shape Transitions in Ba and Xe Nuclei

The rapid transition between spherical and $γ$-soft shapes in Ba and Xe nuclei in the mass region $A \geq 130$ is analyzed using excitation spectra and collective wave functions obtained by diagonalization of 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 the characteristic evolution of excitation spectra and E2 transition probabilities and, in general, a good agreement with available data is obtained. The calculated spectra display fingerprints of a second-order shape phase transition that can approximately be described by analytic solutions corresponding to the E(5) dynamical symmetry.

preprint2009arXiv

3D Relativistic Hartree-Bogoliubov model with a separable pairing interaction

A recently introduced separable pairing force for relativistic Hartree-Bogoliubov (RHB) calculations, adjusted in nuclear matter to the pairing gap of the Gogny force, is employed in the 3D RHB model for triaxial shapes. The pairing force is separable in momentum space but, when transformed to coordinate space in calculations of finite nuclei, it is no longer separable because of translational invariance. The corresponding pairing matrix elements are represented as a sum of a finite number of separable terms in the basis of a 3D harmonic oscillator. The 3D RHB model is applied to the calculation of binding energy surfaces and pairing energy maps for a sequence of even-A Sm isotopes.

preprint2002arXiv

Relativistic Hartree-Bogoliubov model with density-dependent meson-nucleon couplings

The relativistic Hartree-Bogoliubov (RHB) model is extended to include density dependent meson-nucleon couplings. The effective Lagrangian is characterized by a phenomenological density dependence for the $σ$, $ω$ and $ρ$ meson-nucleon vertex functions, adjusted to properties of nuclear matter and finite nuclei. Pairing correlations are described by the pairing part of the finite range Gogny interaction. The new density-dependent effective interaction DD-ME1 is tested in the analysis of the equations of state for symmetric and asymmetric nuclear matter, and of ground-state properties of the Sn and Pb isotopic chains. Results of self-consistent RHB calculations are compared with experimental data, and with results previously obtained in the RHB model with non-linear self-interactions, as well as in the density dependent relativistic hadron field (DDRH) model. Parity-violating elastic electron scattering on Pb and Sn nuclei is calculated using a relativistic optical model with inclusion of Coulomb distortion effects, and the resulting asymmetry parameters are related to the neutron ground-state density distributions.

preprint2002arXiv

Toroidal dipole resonances in the relativistic random phase approximation

The isoscalar toroidal dipole strength distributions in spherical nuclei are calculated in the framework of a fully consistent relativistic random phase approximation. It is suggested that the recently observed "low-lying component of the isoscalar dipole mode" might in fact correspond to the toroidal giant dipole resonance. Although predicted by several theoretical models, the existence of toroidal resonances has not yet been confirmed in experiment. The strong mixing between the toroidal resonance and the dipole compression mode might help to explain the large discrepancy between theory and experiment on the position of isoscalar giant dipole resonances.