Source author record

D. Vretenar

D. Vretenar 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

57works
6topics
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

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

preprint2022arXiv

Model for collective motion

Collective motion is a manifestation of emergent phenomena in medium-heavy and heavy nuclei. A relatively large number of constituent nucleons contribute coherently to nuclear excitations (vibrations, rotations) that are characterized by large electromagnetic moments and transition rates. Basic features of collective excitations are reviewed, and a simple model introduced that describes large-amplitude quadrupole and octupole shape dynamics, as well as the dynamics of induced fission. Modern implementations of the collective Hamiltonian model are based on the microscopic framework of energy density functionals, that provide an accurate global description of nuclear ground states and collective excitations. Results of illustrative calculations are discussed in comparison with available data.

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

Low-energy cluster vibrations in N = Z nuclei

Significant transition strength in light $α$-conjugate nuclei at low energy, typically below 10 MeV, has been observed in many experiments. In this work the isoscalar low-energy response of N=Z nuclei is explored using the Finite Amplitude Method (FAM) based on the microscopic framework of nuclear energy density functionals. Depending on the multipolarity of the excitation and the equilibrium deformation of a particular isotope, the low-energy strength functions display prominent peaks that can be attributed to vibration of cluster structures: $α$+$^{12}$C+$α$ and $α$+$^{16}$O in $^{20}$Ne, $^{12}$C+$^{12}$C in $^{24}$Mg, 4$α$+$^{12}$C in $^{28}$Si, etc. Such cluster excitations are favored in light nuclei with large deformation.

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.

preprint2020arXiv

Shape phase transitions in odd-A Zr isotopes

Spectroscopic properties that characterize shape phase transitions in neutron-rich odd-A Zr isotopes are investigated using the framework of nuclear density functional theory and particle-core coupling. The interacting-boson Hamiltonian of the even-even core nuclei, and the single-particle energies and occupation probabilities of the unpaired neutron are completely determined by deformation constrained self-consistent mean-field calculations based on the relativistic Hartree-Bogoliubov model with a choice of a universal energy density functional and pairing interaction. The triaxial $(β,γ)$ deformation energy surfaces for even-even $^{94-102}$Zr indicates transition from triaxial or $γ$-soft ($^{94,96}$Zr) to prolate ($^{98}$Zr), and triaxial ($^{100,102}$Zr) shapes. The corresponding low-energy excitation spectra of the odd-A Zr isotopes are in very good agreement with recent experimental results. Consistent with the structural evolution of the neighboring even-even Zr nuclei, the state-dependent effective deformations and their fluctuations in the odd-A isotopes indicate a pronounced discontinuity around the transitional nucleus $^{99}$Zr.

preprint2016arXiv

Beyond mean-field boson-fermion model for odd-mass nuclei

A novel method for calculating spectroscopic properties of medium-mass and heavy atomic nuclei with an odd number of nucleons is introduced, based on the framework of nuclear energy density functional theory and the particle-core coupling scheme. The deformation energy surface of the even-even core, as well as the spherical single-particle energies and occupation probabilities of the odd particle(s), are obtained in a self-consistent mean-field calculation determined by the choice of the energy density functional and pairing interaction. This method uniquely determines the parameters of the Hamiltonian of the boson core, and only the strength of the particle-core coupling is specifically adjusted to selected data for a particular nucleus. The approach is illustrated in a systematic study of low-energy excitation spectra and transition rates of axially deformed odd-mass Eu isotopes.

preprint2016arXiv

Spin-orbit interaction in relativistic nuclear structure models

Relativistic self-consistent mean-field (SCMF) models naturally account for the coupling of the nucleon spin to its orbital motion, whereas non-relativistic SCMF methods necessitate a phenomenological ansatz for the effective spin-orbit potential. Recent experimental studies aim to explore the isospin properties of the effective spin-orbit interaction in nuclei. SCMF models are very useful in the interpretation of the corresponding data, however standard relativistic mean-field and non-relativistic Hartree-Fock models use effective spin-orbit potentials with different isovector properties, mainly because exchange contributions are not treated explicitly in the former. The impact of exchange terms on the effective spin-orbit potential in relativistic mean-field models is analysed, and it is shown that it leads to an isovector structure similar to the one used in standard non-relativistic Hartree-Fock. Data on the isospin dependence of spin-orbit splittings in spherical nuclei could be used to constrain the isovector-scalar channel of relativistic mean-field models. The reproduction of the empirical kink in the isotope shifts of even Pb nuclei by relativistic effective interactions points to the occurrence of pseudospin symmetry in the single-neutron spectra in these nuclei.

preprint2015arXiv

Modeling nuclear weak-interaction processes with relativistic energy density functionals

Relativistic energy density functionals have become a standard framework for nuclear structure studies of ground-state properties and collective excitations over the entire nuclide chart. We review recent developments in modeling nuclear weak-interaction processes: charge-exchange excitations and the role of isoscalar proton-neutron pairing, charged-current neutrino-nucleus reactions relevant for supernova evolution and neutrino detectors, and calculation of beta-decay rates for r-process nucleosynthesis.

preprint2015arXiv

Spin-orbit coupling rule in bound fermions systems

Spin-orbit coupling characterizes quantum systems such as atoms, nuclei, hypernuclei, quarkonia, etc., and is essential for understanding their spectroscopic properties. Depending on the system, the effect of spin-orbit coupling on shell structure is large in nuclei, small in quarkonia, perturbative in atoms. In the standard non-relativistic reduction of the single-particle Dirac equation, we derive a universal rule for the relative magnitude of the spin-orbit effect that applies to very different quantum systems, regardless of whether the spin-orbit coupling originates from the strong or electromagnetic interaction. It is shown that in nuclei the near equality of the mass of the nucleon and the difference between the large repulsive and attractive potentials explains the fact that spin-orbit splittings are comparable to the energy spacing between major shells. For a specific ratio between the particle mass and the effective potential whose gradient determines the spin-orbit force, we predict the occurrence of giant spin-orbit energy splittings that dominate the single-particle excitation spectrum.

preprint2015arXiv

The neutron skin thickness from the measured electric dipole polarizability in $^{68}$Ni, $^{120}$Sn, and $^{208}$Pb

The information on the symmetry energy and its density dependence is deduced by comparing the available data on the electric dipole polarizability $α_D$ of ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb with the predictions of the Random Phase Approximation, using a representative set of nuclear energy density functionals. The calculated values of $α_D$ are used to validate different correlations involving $α_D$, the symmetry energy at the saturation density $J$, the corresponding slope parameter $L$ and the neutron skin thickness $Δr_{\!np}$, as suggested by the Droplet Model. A subset of models that reproduce simultaneously the measured polarizabilities in ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb are employed to predict the values of the symmetry energy parameters at saturation density and $Δr_{\!np}$. The resulting intervals are: $J\!=\!30 \text{-}35$ MeV, $L\!=\!20 \text{-} 66$ MeV; and the values for $Δr_{\!np}$ in ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb are in the ranges: 0.15\text{-}0.19 fm, 0.12\text{-}0.16 fm, and 0.13\text{-}0.19 fm, respectively. The strong correlation between the electric dipole polarizabilities of two nuclei is instrumental to predict the values of electric dipole polarizabilities in other nuclei.

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.

preprint2014arXiv

Neutron star structure and collective excitations of finite nuclei

We study relationships between properties of collective excitations in finite nuclei and the phase transition density $n_t$ and pressure $P_t$ at the inner edge separating the liquid core and the solid crust of a neutron star. A theoretical framework that includes the thermodynamic method, relativistic nuclear energy density functionals and the quasiparticle random-phase approximation is employed in a self-consistent calculation of $(n_t,P_t)$ and collective excitations in nuclei. The covariance analysis shows that properties of charge-exchange dipole transitions, isovector giant dipole and quadrupole resonances and pygmy dipole transitions are correlated with the core-crust transition density and pressure. A set of relativistic nuclear energy density functionals, characterized by systematic variation of the density dependence of the symmetry energy of nuclear matter, is used to constrain possible values for $(n_t,P_t)$. By comparing the calculated excitation energies of giant resonances, energy weighted pygmy dipole strength, and dipole polarizability with available data, we obtain the weighted average values: $n_t = 0.0955 \pm 0.0007$ fm$^{-3}$ and $P_t = 0.59 \pm 0.05$ MeV fm$^{-3}$.

preprint2014arXiv

Pairing transitions in the finite-temperature relativistic Hartree-Bogoliubov theory

We formulate the finite-temperature relativistic Hartree-Bogoliubov theory for spherical nuclei based on a point-coupling functional, with the Gogny or separable pairing force. Using the functional PC-PK1, the framework is applied to the study of pairing transitions in Ca, Ni, Sn, and Pb isotopic chains. The separable pairing force reproduces the gaps calculated with the Gogny force not only at zero temperature, but also at finite temperatures. By performing a systematic calculation of the even-even Ca, Ni, Sn, and Pb isotopes, it is found that the critical temperature for a pairing transition generally follows the rule $T_c =0.6 Δ_n(0)$, where $Δ_n(0)$ is the neutron pairing gap at zero temperature. This rule is further verified by adjusting the pairing gap at zero temperature with a strength parameter.

preprint2013arXiv

$β$-decay half-lives of neutron-rich nuclei and matter flow in the $r$-process

The $β$-decay half-lives of neutron-rich nuclei with $20 \leqslant Z \leqslant 50$ are systematically investigated using the newly developed fully self-consistent proton-neutron quasiparticle random phase approximation (QRPA), based on the spherical relativistic Hartree-Fock-Bogoliubov (RHFB) framework. Available data are reproduced by including an isospin-dependent proton-neutron pairing interaction in the isoscalar channel of the RHFB+QRPA model. With the calculated $β$-decay half-lives of neutron-rich nuclei a remarkable speeding up of $r$-matter flow is predicted. This leads to enhanced $r$-process abundances of elements with $A \gtrsim 140$, an important result for the understanding of the origin of heavy elements in the universe.

preprint2013arXiv

Anti-analog giant dipole resonances and the neutron skin of nuclei

We examine a method to determine the neutron-skin thickness of nuclei using data on the charge-exchange anti-analog giant dipole resonance (AGDR). Calculations performed using the relativistic proton-neutron quasiparticle random-phase approximation (RQRPA) reproduce the isotopic trend of the excitation energies of the AGDR, as well as that of the spin-flip giant dipole resonances (IVSGDR), in comparison to available data for the even-even isotopes 112-124Sn. It is shown that the excitation energies of the AGDR, obtained using a set of density-dependent effective interactions which span a range of the symmetry energy at saturation density, supplemented with the experimental values, provide a stringent constraint on the value of the neutron-skin thickness. For 124Sn, in particular, we determine the value DRpn=0.21(5) fm. The result of the present study shows that a measurement of the excitation energy of the AGDR in (p,n) reactions using rare-isotope beams in inverse kinematics, provides a valuable method for the determination of the neutron-skin thickness in exotic nuclei.

preprint2013arXiv

Electric Dipole Polarizability in ${}^{208}$Pb: insights from the Droplet Model

We study the electric dipole polarizability $α_D$ in ${}^{208}$Pb based on the predictions of a large and representative set of relativistic and non-relativistic nuclear mean field models. We adopt the droplet model as a guide to better understand the correlations between $α_D$ and other isovector observables. Insights from the droplet model suggest that the product of $α_D$ and the nuclear symmetry energy at saturation density $J$ is much better correlated with the neutron skin thickness $Δr_{np}$ of ${}^{208}$Pb than the polarizability alone. Correlations of $α_D J$ with $Δr_{np}$ and with the symmetry energy slope parameter $L$ suggest that $α_D J$ is a strong isovector indicator. Hence, we explore the possibility of constraining the isovector sector of thenuclear energy density functional by comparing our theoretical predictions against measurements of both $α_D$ and the parity-violating asymmetry in ${}^{208}$Pb. We find that the recent experimental determination of $α_D$ in ${}^{208}$Pb in combination with the range for the symmetry energy at saturation density $J=[31\pm (2)_{\rm est.}]$\,MeV suggests $Δr_{np}({}^{208}{\rm Pb}) = 0.165 \pm (0.009)_{\rm exp.} \pm (0.013)_{\rm theo.} \pm (0.021)_{\rm est.} {\rm fm}$ and $L= 43 \pm(6)_{\rm exp.} \pm (8)_{\rm theo.}\pm(12)_{\rm est.}$ MeV.

preprint2013arXiv

Giant Quadrupole Resonances in 208Pb, the nuclear symmetry energy and the neutron skin thickness

Recent improvements in the experimental determination of properties of the Isovector Giant Quadrupole Resonance (IVGQR), as demonstrated in the A=208 mass region, may be instrumental for characterizing the isovector channel of the effective nuclear interaction. We analyze properties of the IVGQR in 208Pb, using both macroscopic and microscopic approaches. The microscopic method is based on families of non-relativistic and covariant Energy Density Functionals (EDF), characterized by a systematic variation of isoscalar and isovector properties of the corresponding nuclear matter equations of state. The macroscopic approach yields an explicit dependence of the nuclear symmetry energy at some subsaturation density, for instance S(ρ=0.1 fm^{-3}), or the neutron skin thickness Δr_{np} of a heavy nucleus, on the excitation energies of isoscalar and isovector GQRs. Using available data it is found that S(ρ=0.1 fm{}^{-3})=23.3 +/- 0.6 MeV. Results obtained with the microscopic framework confirm the correlation of the Δr_{np} to the isoscalar and isovector GQR energies, as predicted by the macroscopic model. By exploiting this correlation together with the experimental values for the isoscalar and isovector GQR energies, we estimate Δr_{np} = 0.14 +/- 0.03 fm for 208Pb, and the slope parameter of the symmetry energy: L = 37 +/- 18 MeV.

preprint2013arXiv

Incompressibility of finite fermionic systems: stable and exotic atomic nuclei

The incompressibility of finite fermionic systems is investigated using analytical approaches and microscopic models. The incompressibility of a system is directly linked to the zero-point kinetic energy of constituent fermions, and this is a universal feature of fermionic systems. In the case of atomic nuclei, this implies a constant value of the incompressibility in medium-heavy and heavy nuclei. The evolution of nuclear incompressibility along Sn and Pb isotopic chains is analyzed using global microscopic models, based on both non-relativistic and relativistic energy functionals. The result is an almost constant incompressibility in stable nuclei and systems not far from stability, and a steep decrease in nuclei with pronounced neutron excess, caused by the emergence of a soft monopole mode in neutron-rich nuclei.

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.

preprint2013arXiv

Measurement of the 92,93,94,100Mo(g,n) reactions by Coulomb Dissociation

The Coulomb Dissociation (CD) cross sections of the stable isotopes 92,94,100Mo and of the unstable isotope 93Mo were measured at the LAND/R3B setup at GSI Helmholtzzentrum für Schwerionenforschung in Darmstadt, Germany. Experimental data on these isotopes may help to explain the problem of the underproduction of 92,94Mo and 96,98Ru in the models of p-process nucleosynthesis. The CD cross sections obtained for the stable Mo isotopes are in good agreement with experiments performed with real photons, thus validating the method of Coulomb Dissociation. The result for the reaction 93Mo(g,n) is especially important since the corresponding cross section has not been measured before. A preliminary integral Coulomb Dissociation cross section of the 94Mo(g,n) reaction is presented. Further analysis will complete the experimental database for the (g,n) production chain of the p-isotopes of molybdenum.

preprint2013arXiv

Microscopic analysis of the octupole phase transition in Th isotopes

A shape phase transition between stable octupole deformation and octupole vibrations in Th nuclei is analyzed in a microscopic framework based on nuclear density functional theory. The relativistic functional DD-PC1 is used to calculate axially-symmetric quadrupole-octupole constrained energy surfaces. Observables related to order parameters are computed using an interacting-boson Hamiltonian, with parameters determined by mapping the microscopic energy surfaces to the expectation value of the Hamiltonian in the boson condensate. The systematics of constrained energy surfaces and low-energy excitation spectra point to the occurrence of a phase transition between octupole-deformed shapes and shapes characterized by octupole-soft potentials.

preprint2013arXiv

Neutron-skin thickness of $^{208}$Pb, and symmetry-energy constraints from the study of the anti-analog giant dipole resonance

The $^{208}$Pb($p$,$nγ\bar p$) $^{207}$Pb reaction at a beam energy of 30 MeV has been used to excite the anti-analog of the giant dipole resonance (AGDR) and to measure its $γ$-decay to the isobaric analog state in coincidence with proton decay of IAS. The energy of the transition has also been calculated with the self-consistent relativistic random-phase approximation (RRPA), and found to be linearly correlated to the predicted value of the neutron-skin thickness ($ΔR_{pn}$). By comparing the theoretical results with the measured transition energy, the value of 0.190 $\pm$ 0.028 fm has been determined for $ΔR_{pn}$ of $^{208}$Pb, in agreement with previous experimental results. The AGDR excitation energy has also been used to calculate the symmetry energy at saturation ($J=32.7 \pm 0.6$ MeV) and the slope of the symmetry energy ($L=49.7 \pm 4.4$ MeV), resulting in more stringent constraints than most of the previous studies.

preprint2013arXiv

Neutron-skin thickness of 208Pb from the energy of the anti-analog giant dipole resonance

The energy of the charge-exchange Anti-analog Giant Dipole Resonance (AGDR) has been calculated for the 208Pb isotope using the state-of-the-art fully self-consistent relativistic proton-neutron quasiparticle random-phase approximation based on the Relativistic Hartree-Bogoliubov model. It is shown that the AGDR centroid energy is very sensitively related to the corresponding neutron-skin thickness. The neutron-skin thickness of 208Pb has been determined very precisely by comparing the theoretical results with the available experimental data on E(AGDR). The result DR(pn)= 0.161(42) agrees nicely with the previous experimental results.

preprint2013arXiv

Nuclear Symmetry Energy: constraints from Giant Quadrupole Resonances and Parity Violating Electron Scattering

Experimental and theoretical efforts are being devoted to the study of observables that can shed light on the properties of the nuclear symmetry energy. We present our new results on the excitation energy [X. Roca-Maza et al., Phys. Rev. C 87, 034301 (2013)] and polarizability of the Isovector Giant Quadrupole Resonance (IVGQR), which has been the object of new experimental investigation [S. S. Henshaw et al., Phys. Rev. Lett. 107, 222501 (2011)]. We also present our theoretical analysis on the parity violating asymmetry at the kinematics of the Lead Radius Experiment [S. Abrahamyan et al. (PREx Collaboration), Phys. Rev. Lett. 108, 112502 (2012)] and highlight its relation with the density dependence of the symmetry energy [X. Roca-Maza et al., Phys. Rev. Lett. 106, 252501 (2011)].

preprint2013arXiv

Simultaneous quadrupole and octupole shape phase transitions in Thorium

The evolution of quadrupole and octupole shapes in Th isotopes is studied in the framework of nuclear Density Functional Theory. Constrained energy maps and observables calculated with microscopic collective Hamiltonians indicate the occurrence of a simultaneous quantum shape phase transition between spherical and quadrupole-deformed prolate shapes, and between non-octupole and octupole-deformed shapes, as functions of the neutron number. The nucleus $^{224}$Th is closest to the critical point of a double phase transition. A microscopic mechanism of this phenomenon is discussed in terms of the evolution of single-nucleon orbitals with deformation.

preprint2013arXiv

Symmetry energy in nuclear density functional theory

The nuclear symmetry energy represents a response to the neutron-proton asymmetry. In this survey we discuss various aspects of symmetry energy in the framework of nuclear density functional theory, considering both non-relativistic and relativistic self-consistent mean-field realizations side-by-side. Key observables pertaining to bulk nucleonic matter and finite nuclei are reviewed. Constraints on the symmetry energy and correlations between observables and symmetry-energy parameters, using statistical covariance analysis, are investigated. Perspectives for future work are outlined in the context of ongoing experimental efforts.

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

Electric dipole polarizability and the neutron skin

The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.

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

Fragmentation of spin-dipole strength in $^{90}$Zr and $^{208}$Pb

An extension of time-dependent covariant density functional theory that includes particle-vibration coupling is applied to the charge-exchange channel. Spin-dipole excitation spectra are calculated an compared to available data for $^{90}$Zr and $^{208}$Pb. A significant fragmentation is found for all three angular-momentum components of the spin-dipole strength as a result of particle-vibration coupling, as well as a shift of a portion of the strength to higher energy. A high-energy tail is formed in the strength distribution that linearly decreases with energy. Using a model-independent sum rule,the corresponding neutron skin thickness is estimated and shown to be consistent with values obtained at the mean-field level.

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

Low-energy isovector and isoscalar dipole response in neutron-rich nuclei

The self-consistent random phase approximation (RPA), based on the framework of relativistic energy density functionals, is employed in the study of isovector and isoscalar dipole response in $^{68}$Ni, $^{132}$Sn, and $^{208}$Pb. The evolution of pygmy dipole states (PDS) in the region of low excitation energies is analyzed as a function of the density-dependence of the symmetry energy for a set of relativistic effective interactions. The occurrence of PDS is predicted in the response to both the isovector and isoscalar dipole operators, and its strength is enhanced with the increase of the symmetry energy at saturation and the slope of the symmetry energy. In both channels the PDS exhausts a relatively small fraction of the energy-weighted sum rule but a much larger percentage of the inverse energy-weighted sum rule. For the isovector dipole operator the reduced transition probability $B(E1)$ of the PDS is generally small because of pronounced cancellation of neutron and proton partial contributions. The isoscalar reduced transition amplitude is predominantly determined by neutron particle-hole configurations, most of which add coherently, and this results in a collective response of the PDS to the isoscalar dipole operator.

preprint2012arXiv

Neutron-skin thickness from the study of the anti-analog giant dipole resonance

The gamma-decay of the anti-analog of the giant dipole resonance (AGDR) has been measured to the isobaric analog state excited in the p(124Sn,n) reaction at a beam energy of 600 MeV/nucleon. The energy of the transition was also calculated with state-of-the-art self-consistent random-phase approximation (RPA) and turned out to be very sensitive to the neutron-skin thickness (\DeltaR_(pn)). By comparing the theoretical results with the measured one, the \DeltaR_(pn) value for 124Sn was deduced to be 0.175 \pm 0.048 fm, which agrees well with the previous results. The energy of the AGDR measured previously for ^(208)Pb was also used to determine the \DeltaR_(pn) for ^(208)Pb. In this way a very precise \DeltaR_(pn) = 0.181 \pm 0.031 neutron-skin thickness has been obtained for 208Pb. The present method offers new possibilities for measuring the neutron-skin thicknesses of very exotic isotopes.

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.

preprint2012arXiv

Role of momentum transfer in the quenching of Gamow-Teller strength

The relativistic proton-neutron quasiparticle random phase approximation (pn-RQRPA) is applied in the calculation of the L=0 strength in charge-exchange reactions on $^{48}$Ca, $^{90}$Zr, $^{208}$Pb and nuclei in the Sn isotopic chain. The microscopic theoretical framework is based on the relativistic Hartree-Bogoliubov (RHB) model for the nuclear ground state. The calculation is fully self-consistent, i.e. the same interaction is used both in the RHB equations that determine the quasiparticle basis, and in the matrix equations of the pn-RQRPA. The inclusion of the higher-order terms that include the effect of finite momentum transfer, primarily the isovector spin monopole (IVSM) term, in the transition operator shifts a portion of the strength to the high-energy region above the Gamow-Teller (GT) resonance. The total strength is slightly enhanced in nuclei with small neutron-to-proton ratio but remains unchanged with increasing neutron excess. Based on the strength obtained using the full L=0 transition operator in the pn-RQRPA calculation, we have estimated the impact of the IVSM on the strength measured in the charge-exchagne reactions on $^{90}$Zr and found that the data are consistent with the Ikeda sum rule.

preprint2012arXiv

Sensitivity of the electric dipole polarizability to the neutron skin thickness in ${}^{208}$Pb

The static dipole polarizability, $α_{\rm D}$, in ${}^{208}$Pb has been recently measured with high-resolution via proton inelastic scattering at the Research Center for Nuclear Physics (RCNP). This observable is thought to be intimately connected with the neutron skin thickness, $r_{\rm skin}$, of the same nucleus and, more fundamentally, it is believed to be associated with the density dependence of the nuclear symmetry energy. The impact of $r_{\rm skin}$ on $α_{\rm D}$ in ${}^{208}$Pb is investigated and discussed on the basis of a large and representative set of relativistic and non-relativistic nuclear energy density functionals (EDF).

preprint2011arXiv

Low-energy monopole strength in exotic Nickel isotopes

Low-energy strength is predicted for the isoscalar monopole response of neutron-rich Ni isotopes, in calculations performed using the microscopic Skyrme HF+RPA and relativistic RHB+RQRPA models. Both models, although based on different energy density functionals, predict the occurrence of pronounced monopole states in the energy region between 10 MeV and 15 MeV, well separated from the isoscalar GMR. The analysis of transition densities and corresponding particle-hole configurations shows that these states represent almost pure neutron single hole-particle excitations. Even though their location is not modified with respect to the corresponding unperturbed states, their (Q)RPA strength is considerably enhanced by the residual interaction. The theoretical analysis predicts the gradual enhancement of low-energy monopole strength with neutron excess.

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.

preprint2011arXiv

Relativistic Hartree-Fock-Bogoliubov model for deformed nuclei

The Relativistic Hartree-Fock-Bogoliubov model for axially deformed nuclei (RHFBz) is introduced. The model is based on an effective Lagrangian with density-dependent meson-nucleon couplings in the particle-hole channel, and the pairing part of the Gogny force is used in the pairing channel. The RHFBz quasiparticle equations are solved by expansion in the basis of a deformed harmonic oscillator. Illustrative RHFBz calculations are performed for Carbon, Neon and Magnesium isotopes. The effect of the explicitly including the pion field is investigated for binding energies, deformation parameters, and charge radii.

preprint2011arXiv

Uncertainties in modeling low-energy neutrino induced reactions on iron group nuclei

Charged-current neutrino-nucleus cross sections for 54,56Fe and 58,60Ni are calculated and compared using frameworks based on relativistic and Skyrme energy density functionals, and the shell model. The current theoretical uncertainties in modeling neutrino-nucleus cross sections are assessed in relation to the predicted Gamow-Teller transition strength and available data, multipole decomposition of the cross sections, and cross sections averaged over the Michel flux and Fermi-Dirac distribution. Employing different microscopic approaches and models, the DAR neutrino-56Fe cross section and its theoretical uncertainty are estimated: <sigma>_th=(258+-57) 10^{-42} cm^2, in very good agreement with the experimental value: <sigma>_exp=(256+-108+-43) 10^{-42} cm^2.

preprint2010arXiv

Configuration mixing of angular-momentum projected triaxial relativistic mean-field wave functions

The framework of relativistic energy density functionals is extended to include correlations related to the restoration of broken symmetries and to fluctuations of collective variables. The generator coordinate method is used to perform configuration mixing of angular-momentum projected wave functions, generated by constrained self-consistent relativistic mean-field calculations for triaxial shapes. The effects of triaxial deformation and of $K$-mixing is illustrated in a study of spectroscopic properties of low-spin states in $^{24}$Mg.

preprint2010arXiv

Configuration mixing of angular-momentum projected triaxial relativistic mean-field wave functions. II. Microscopic analysis of low-lying states in magnesium isotopes

The recently developed structure model that uses the generator coordinate method to perform configuration mixing of angular-momentum projected wave functions, generated by constrained self-consistent relativistic mean-field calculations for triaxial shapes (3DAMP+GCM), is applied in a systematic study of ground states and low-energy collective states in the even-even magnesium isotopes $^{20-40}$Mg. Results obtained using a relativistic point-coupling nucleon-nucleon effective interaction in the particle-hole channel, and a density-independent $δ$-interaction in the pairing channel, are compared to data and with previous axial 1DAMP+GCM calculations, both with a relativistic density functional and the non-relativistic Gogny force. The effects of the inclusion of triaxial degrees of freedom on the low-energy spectra and E2 transitions of magnesium isotopes are examined.

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.

preprint2009arXiv

Calculation of stellar electron-capture cross sections on nuclei based on microscopic Skyrme functionals

A fully self-consistent microscopic framework for evaluation of nuclear weak-interaction rates at finite temperature is introduced, based on Skyrme functionals. The single-nucleon basis and the corresponding thermal occupation factors of the initial nuclear state are determined in the finite-temperature Skyrme Hartree-Fock model, and charge-exchange transitions to excited states are computed using the finite-temperature RPA. Effective interactions are implemented self-consistently: both the finite-temperature single-nucleon Hartree-Fock equations and the matrix equations of RPA are based on the same Skyrme energy density functional. Using a representative set of Skyrme functionals, the model is applied in the calculation of stellar electron-capture cross sections for selected nuclei in the iron mass group and for neutron-rich Ge isotopes.

preprint2009arXiv

Low-energy multipole response in nuclei at finite temperature

The multipole response of nuclei at temperatures T=0-2 MeV is studied using a self-consistent finite-temperature RPA (random phase approximation) based on relativistic energy density functionals. Illustrative calculations are performed for the isoscalar monopole and isovector dipole modes and, in particular, the evolution of low-energy excitations with temperature is analyzed, including the modification of pygmy structures. Both for the monopole and dipole modes, in the temperature range T=1-2 MeV additional transition strength appears at low energies because of thermal unblocking of single-particle orbitals close to the Fermi level. A concentration of dipole strength around 10 MeV excitation energy is predicted in $^{60,62}$Ni, where no low-energy excitations occur at zero temperature. The principal effect of finite temperature on low-energy strength that is already present at zero temperature, e.g. in $^{68}$Ni and $^{132}$Sn, is the spreading of this structure to even lower energy and the appearance of states that correspond to thermally unblocked transitions.

preprint2007arXiv

Inclusive charged-current neutrino-nucleus reactions calculated with the relativistic quasiparticle random phase approximation

Inclusive neutrino-nucleus cross sections are calculated using a consistent relativistic mean-field theoretical framework. The weak lepton-hadron interaction is expressed in the standard current-current form, the nuclear ground state is described with the relativistic Hartree-Bogoliubov model, and the relevant transitions to excited nuclear states are calculated in the relativistic quasiparticle random phase approximation. Illustrative test calculations are performed for charged-current neutrino reactions on $^{12}$C, $^{16}$O, $^{56}$Fe, and $^{208}$Pb, and results compared with previous studies and available data. Using the experimental neutrino fluxes, the averaged cross sections are evaluated for nuclei of interest for neutrino detectors. We analyze the total neutrino-nucleus cross sections, and the evolution of the contribution of the different multipole excitations as a function of neutrino energy. The cross sections for reactions of supernova neutrinos on $^{16}$O and $^{208}$Pb target nuclei are analyzed as functions of the temperature and chemical potential.

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.

preprint1998arXiv

Relativistic Hartree-Bogoliubov description of the deformed N = 28 region

Ground-state properties of neutron-rich $N \approx 28$ nuclei are described in the framework of Relativistic Hartree Bogoliubov (RHB) theory. The model uses the NL3 effective interaction in the mean-field Lagrangian, and describes pairing correlations by the pairing part of the finite range Gogny interaction D1S. Two-neutron separation energies and ground-state quadrupole deformations that result from fully self-consistent RHB solutions, are compared with available experimental data. The model predicts a strong suppression of the spherical N=28 shell gap for neutron-rich nuclei: the $1f7/2 \to fp$ core breaking results in deformed ground states. Shape coexistence is expected for neutron-rich Si, S and Ar isotopes.