Source author record

J. Dobaczewski

J. Dobaczewski 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

58works
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

58 published item(s)

preprint2022arXiv

Precision measurement of the magnetic octupole moment in $^{45}$Sc as a test for state-of-the-art atomic- and nuclear-structure theory

We report on measurements of the hyperfine $A, B$ and $C$-constants of the $3d4s^2 ~^2D_{5/2}$ and $3d4s^2 ~^2D_{3/2}$ atomic states in $^{45}$Sc. High-precision atomic calculations of the hyperfine fields of these states and second-order corrections are performed, and are used to extract $C_{5/2}=-0.06(6)$ kHz and $C_{3/2}=+0.04(3)$ kHz from the data. These results are one order of magnitude more precise than the available literature. From the combined analysis of both atomic states, we infer the nuclear magnetic octupole moment $Ω= -0.07(53) μ_N b$, including experimental and atomic structure-related uncertainties. With a single valence proton outside of a magic calcium core, scandium is ideally suited to test a variety of nuclear models, and to investigate in-depth the many intriguing nuclear structure phenomena observed within the neighboring isotopes of calcium. We perform nuclear shell-model calculations of $Ω$, and furthermore explore the use of Density Functional Theory for evaluating $Ω$. From this, mutually consistent theoretical values of $Ω$ are obtained, which are in agreement with the experimental value. This confirms atomic structure calculations possess the accuracy and precision required for magnetic octupole moment measurements, and shows that modern nuclear theory is capable of providing meaningful insight into this largely unexplored observable.

preprint2020arXiv

Model nuclear energy density functionals derived from ab initio calculations

We present the first application of a new approach, proposed in [Journal of Physics G: Nuclear and Particle Physics, 43, 04LT01 (2016)] to derive coupling constants of the Skyrme energy density functional (EDF) from ab initio Hamiltonian. By perturbing the ab initio Hamiltonian with several functional generators defining the Skyrme EDF, we create a set of metadata that is then used to constrain the coupling constants of the functional. We use statistical analysis to obtain such an ab initio-equivalent Skyrme EDF. We find that the resulting functional describes properties of atomic nuclei and infinite nuclear matter quite poorly. This may point out to the necessity of building up the ab initio-equivalent functionals from more sophisticated generators. However, we also indicate that the current precision of the ab initio calculations may be insufficient for deriving meaningful nuclear EDFs.

preprint2019arXiv

Regularized pseudopotential for mean-field calculations

We present preliminary results obtained with a finite-range two-body pseudopotential complemented with zero-range spin-orbit and density-dependent terms. After discussing the penalty function used to adjust parameters, we discuss predictions for binding energies of spherical nuclei calculated at the mean-field level, and we compare them with those obtained using the standard Gogny D1S finite-range effective interaction.

preprint2018arXiv

Towards a novel energy density functional for beyond-mean-field calculations with pairing and deformation

We take an additional step towards the optimization of the novel finite-range pseudopotential at constrained Hartree-Fock-Bogolyubov level and implement an optimization procedure within an axial code using harmonic oscillator basis. We perform the optimization using three different numbers of the harmonic oscillator shells. We apply the new parameterizations in the O-Kr part of the nuclear chart and isotopic chain of Sn, and we compare the results with experimental values and those given by a parameterization obtained using a spherical code.

preprint2015arXiv

Ab Initio Derivation of Model Energy Density Functionals

I propose a simple and manageable method that allows for deriving coupling constants of model energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion systems. A proof-of-principle application allows for linking properties of finite nuclei, determined by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme functional. The method does not rely on properties of infinite fermion systems but on the ab initio calculations in finite systems. It also allows for quantifying merits of different model EDFs in describing the ab initio results.

preprint2015arXiv

Approximate restoration of translational and rotational symmetries within the Lipkin method

Background: Nuclear self-consistent mean-field approaches are rooted in the density functional theory and, through the spontaneous symmetry breaking mechanism, allow for including important correlations, while keeping the simplicity of the approach. Because real ground states should have all symmetries of the nuclear Hamiltonian, these methods require subsequent symmetry restoration. Purpose: We implement and study Lipkin method of approximate variation after projection applied to the restoration of the translational or rotational symmetries. Methods: We use Lipkin operators up to quadratic terms in momenta or angular momenta with self-consistently determined values of the Peierls-Yoccoz translational masses or moments of inertia, respectively. Calculations based on Skyrme energy-density functional are performed for heavy, deformed, and paired nuclei. Results: In deformed nuclei, the Peierls-Yoccoz translational masses along three different principal-axes directions of the intrinsic system can be different, which illustrates different widths of the total-momentum distributions. Numerically, the differences are only of the order of a few per cent. For the rotational-symmetry restoration, the Lipkin method compares favorably with the exact angular-momentum projection, which requires much larger computational effort. Conclusions: The Lipkin method of translational and rotational symmetry restoration is a practical low-cost method of determining the corresponding correlation energies. It allows for a simultaneous restoration of several symmetries and can be relatively easily implemented in the standard self-consistent mean-field calculations, including those required for the parameter adjustments.

preprint2015arXiv

Properties of nuclei in the nobelium region studied within the covariant, Skyrme, and Gogny energy density functionals

We calculate properties of the ground and excited states of nuclei in the nobelium region for proton and neutron numbers of 92 <= Z <= 104 and 144 <= N <= 156, respectively. We use three different energy-density-functional (EDF) approaches, based on covariant, Skyrme, and Gogny functionals, each within two different parameter sets. A comparative analysis of the results obtained for odd-even mass staggerings, quasiparticle spectra, and moments of inertia allows us to identify single-particle and shell effects that are characteristic to these different models and to illustrate possible systematic uncertainties related to using the EDF modelling

preprint2015arXiv

Recoil-alpha-fission and recoil-alpha-alpha-fission events observed in the reaction Ca-48 + Am-243

Products of the fusion-evaporation reaction Ca-48 + Am-243 were studied with the TASISpec set-up at the gas-filled separator TASCA at the GSI Helmholtzzentrum für Schwerionenforschung. Amongst the detected thirty correlated alpha-decay chains associated with the production of element Z=115, two recoil-alpha-fission and five recoil-alpha-alpha-fission events were observed. The latter are similar to four such events reported from experiments performed at the Dubna gas-filled separator. Contrary to their interpretation, we propose an alternative view, namely to assign eight of these eleven decay chains of recoil-alpha(-alpha)-fission type to start from the 3n-evaporation channel 115-288. The other three decay chains remain viable candidates for the 2n-evaporation channel 115-289.

preprint2014arXiv

Error Estimates of Theoretical Models: a Guide

This guide offers suggestions/insights on uncertainty quantification of nuclear structure models. We discuss a simple approach to statistical error estimates, strategies to assess systematic errors, and show how to uncover inter-dependencies by correlation analysis. The basic concepts are illustrated through simple examples. By providing theoretical error bars on predicted quantities and using statistical methods to study correlations between observables, theory can significantly enhance the feedback between experiment and nuclear modeling.

preprint2014arXiv

Isospin-invariant Skyrme energy-density-functional approach with axial symmetry

We develop the isospin-invariant Skyrme-EDF method by considering local densities in all possible isospin channels and proton-neutron (p-n) mixing terms as mandated by the isospin symmetry. The EDF employed has the most general form that depends quadratically on the isoscalar and isovector densities. We test and benchmark the resulting p-n EDF approach, and study the general properties of the new scheme by means of the cranking in the isospin space. We extend the existing axial DFT solver HFBTHO to the case of isospin-invariant EDF approach with all possible p-n mixing terms. Explicit expressions have been derived for all the densities and potentials that appear in the isospin representation. In practical tests, we consider the Skyrme EDF SkM* and, as a first application, concentrate on Hartree-Fock aspects of the problem, i.e., pairing has been disregarded. Calculations have been performed for the (A=78, T~11), (A=40, T~8), and (A=48, T~4) isobaric analog chains. Isospin structure of self-consistent p-n mixed solutions has been investigated with and without the Coulomb interaction, which is the sole source of isospin symmetry breaking in our approach. The extended axial HFBTHO solver has been benchmarked against the symmetry-unrestricted HFODD code for deformed and spherical states. We developed and tested a general isospin-invariant Skyrme-EDF framework. The new approach permits spin-isospin densities that may give rise to, hitherto, unexplored modes in the excitation spectrum. The new formalism has been tested in the Hartree-Fock limit. A systematic comparison between HFODD and HFBTHO results show a maximum deviation of about 10 keV on the total binding energy for deformed nuclei when the Coulomb term is included. Without this term, the results of both solvers agree down to a ~10 eV level.

preprint2014arXiv

Lipkin method of particle-number restoration to higher orders

Background: On the mean-field level, pairing correlations are incorporated through the Bogoliubov-Valatin transformation, whereupon the particle degrees of freedom are replaced by quasiparticles. This approach leads to a spontaneous breaking of the particle-number symmetry and mixing of states with different particle numbers. In order to restore the particle number, various methods have been employed, which are based on projection approaches before or after variation. Approximate variation-after-projection (VAP) schemes, utilizing the Lipkin method, have mostly been used within the Lipkin-Nogami prescription. Purpose: Without recurring to the Lipkin-Nogami prescription, and using instead states rotated in the gauge space, we derive the Lipkin method of particle-number restoration up to sixth order and we test the convergence and accuracy of the obtained expansion. Methods: We perform self-consistent calculations using the higher-order Lipkin method to restore the particle-number symmetry in the framework of superfluid nuclear energy-density functional theory. We also apply the Lipkin method to a schematic exactly solvable two-level pairing model. Results: Calculations performed in open-shell tin and lead isotopes show that the Lipkin method converges at fourth order and satisfactorily reproduces the VAP ground-state energies and energy kernels. Near closed shells, the higher-order Lipkin method cannot be applied because of a non-analytic kink in the ground-state energies in function of the particle number. Conclusions: In open-shell nuclei, the higher-order Lipkin method provides a good approximation to the exact VAP energies. The method is computationally inexpensive, making it particularly suitable, for example, for future optimizations of the nuclear energy-density functionals and simultaneous restoration of different symmetries.

preprint2014arXiv

Nonlocal energy density functionals for low-energy nuclear structure

We introduce a finite-range pseudopotential built as an expansion in derivatives up to next-to-next-to-next-to-leading order (N$^3$LO) and we calculate the corresponding nonlocal energy density functional (EDF). The coupling constants of the nonlocal EDF, for both finite nuclei and infinite nuclear matter, are expressed through the parameters of the pseudopotential. All central, spin-orbit, and tensor terms of the pseudopotential are derived both in the spherical-tensor and Cartesian representation. At next-to-leading order (NLO), we also derive relations between the nonlocal EDF expressed in the spherical-tensor and Cartesian formalism. Finally, a simplified version of the finite-range pseudopotential is considered, which generates the EDF identical to that generated by a local potential.

preprint2014arXiv

Pairing-induced speedup of nuclear spontaneous fission

Collective inertia is strongly influenced at the level crossing at which quantum system changes diabatically its microscopic configuration. Pairing correlations tend to make the large-amplitude nuclear collective motion more adiabatic by reducing the effect of those configuration changes. Competition between pairing and level crossing is thus expected to have a profound impact on spontaneous fission lifetimes. To elucidate the role of nucleonic pairing on spontaneous fission, we study the dynamic fission trajectories of $^{264}$Fm and $^{240}$Pu using the state-of-the-art self-consistent framework. We employ the superfluid nuclear density functional theory with the Skyrme energy density functional SkM$^*$ and a density-dependent pairing interaction. Along with shape variables, proton and neutron pairing correlations are taken as collective coordinates. The collective inertia tensor is calculated within the nonperturbative cranking approximation. The fission paths are obtained by using the least action principle in a four-dimensional collective space of shape and pairing coordinates. Pairing correlations are enhanced along the minimum-action fission path. For the symmetric fission of $^{264}$Fm, where the effect of triaxiality on the fission barrier is large, the geometry of fission pathway in the space of shape degrees of freedom is weakly impacted by pairing. This is not the case for $^{240}$Pu where pairing fluctuations restore the axial symmetry of the dynamic fission trajectory. The minimum-action fission path is strongly impacted by nucleonic pairing. In some cases, the dynamical coupling between shape and pairing degrees of freedom can lead to a dramatic departure from the static picture. Consequently, in the dynamical description of nuclear fission, particle-particle correlations should be considered on the same footing as those associated with shape degrees of freedom.

preprint2014arXiv

Rotational properties of nuclei around $^{254}$No investigated using a spectroscopic-quality Skyrme energy density functional

Nuclei in the $Z\approx100$ mass region represent the heaviest systems where detailed spectroscopic information is experimentally available. Although microscopic-macroscopic and self-consistent models have achieved great success in describing the data in this mass region, a fully satisfying precise theoretical description is still missing. By using fine-tuned parametrizations of the energy density functionals, the present work aims at an improved description of the single-particle properties and rotational bands in the nobelium region. Such locally optimized parameterizations may have better properties when extrapolating towards the superheavy region. Skyrme-Hartree-Fock-Bogolyubov and Lipkin-Nogami methods were used to calculate the quasiparticle energies and rotational bands of nuclei in the nobelium region. Starting from the most recent Skyrme parametrization, UNEDF1, the spin-orbit coupling constants and pairing strengths have been tuned, so as to achieve a better agreement with the excitation spectra and odd-even mass differences in $^{251}$Cf and $^{249}$Bk. The quasiparticle properties of $^{251}$Cf and $^{249}$Bk were very well reproduced. At the same time, crucial deformed neutron and proton shell gaps open up at $N=152$ and $Z=100$, respectively. Rotational bands in Fm, No, and Rf isotopes, where experimental data are available, were also fairly well described. To help future improvements towards a more precise description, small deficiencies of the approach were carefully identified. In the $Z\approx100$ mass region, larger spin-orbit strengths than those from global adjustments lead to improved agreement with data. Puzzling effects of particle-number restoration on the calculated moment of inertia, at odds with the experimental behaviour, require further scrutiny.

preprint2014arXiv

Spectroscopic properties of nuclear Skyrme energy density functionals

We address the question of how to improve the agreement between theoretical nuclear single-particle energies (SPEs) and experiment. Empirically, in doubly magic nuclei, the SPEs can be deduced from spectroscopic properties of odd nuclei that have one more, or one less neutron or proton. Theoretically, bare SPEs, before being confronted with experiment, must be corrected for the effects of the particle-vibration-coupling (PVC). In the present work, we determine the PVC corrections in a fully self-consistent way. Then, we adjust the SPEs, with PVC corrections included, to empirical data. In this way, the agreement with experiment, on average, improves; nevertheless, large discrepancies still remain. We conclude that the main source of disagreement is still in the underlying mean fields, and not in including or neglecting the PVC corrections.

preprint2014arXiv

Structure of Superheavy Nuclei Along Element 115 Decay Chains

A recent high-resolution $α$, $X$-ray, and $γ$-ray coincidence-spectroscopy experiment offered first glimpse of excitation schemes of isotopes along $α$-decay chains of $Z=115$. To understand these observations and to make predictions about shell structure of superheavy nuclei below $^{288}115$, we employ two complementary mean-field models: self-consistent Skyrme Energy Density Functional approach and the macroscopic-microscopic Nilsson model. We discuss the spectroscopic information carried by the new data. In particular, candidates for the experimentally observed $E1$ transitions in $^{276}$Mt are proposed. We find that the presence and nature of low-energy $E1$ transitions in well-deformed nuclei around $Z=110, N=168$ strongly depends on the strength of the spin-orbit coupling; hence, it provides an excellent constraint on theoretical models of superheavy nuclei. To clarify competing theoretical scenarios, an experimental search for $E1$ transitions in odd-$A$ systems $^{275,277}$Mt, $^{275}$Hs, and $^{277}$Ds is strongly recommended.

preprint2013arXiv

Isospin mixing within the multi-reference nuclear density functional theory and beyond - selected aspects

The results of systematic calculations of isospin-symmetry-breaking corrections to superallowed beta-decays based on the self-consistent isospin- and angular-momentum-projected nuclear density functional theory (DFT) are reviewed with an emphasis on theoretical uncertainties of the model. Extensions of the formalism towards no core shell model approach with basis cutoff scheme dictated by the self-consistent particle-hole DFT solutions will be also discussed.

preprint2013arXiv

Isospin mixing within the symmetry restored density functional theory and beyond

We present results of systematic calculations of the isospin-symmetry-breaking corrections to the superallowed I=$0+,T=1 --> I=0+,T=1 beta-decays, based on the self-consistent isospin- and angular-momentum-projected nuclear density functional theory (DFT). We discuss theoretical uncertainties of the formalism related to the basis truncation, parametrization of the underlying energy density functional, and ambiguities related to determination of Slater determinants in odd-odd nuclei. A generalization of the double-projected DFT model towards a no core shell-model-like configuration-mixing approach is formulated and implemented. We also discuss new opportunities in charge-symmetry- and charge-independence-breaking studies offered by the newly developed DFT formalism involving proton-neutron mixing in the particle-hole channel.

preprint2013arXiv

Kerman-Onishi conditions in self-consistent tilted-axis-cranking mean-field calculations

\item[Background] For cranked mean-field calculations with arbitrarily oriented rotational frequency vector $\boldsymbolω$ in the intrinsic frame, one has to employ constraints on average values of the quadrupole-moment tensor, so as to keep the nucleus in the principal-axis reference frame. Kerman and Onishi [Nucl. Phys. A {\bf 361}, 179 (1981)] have shown that the Lagrangian multipliers that correspond to the required constraints are proportional to $\boldsymbolω \times \boldsymbol{J}$, where $\boldsymbol{J}$ is the average angular momentum vector. \item[Purpose] We study the validity and consequences of the Kerman-Onishi conditions in the context of self-consistent tilted-axis-cranking (TAC) mean-field calculations. \item[Methods] We perform a two-dimensional self-consistent calculations (with and without pairing) utilizing the symmetry-unrestricted solver {\sc hfodd}. At each tilting angle, we compare the calculated values of quadrupole-moment-tensor Lagrangian multipliers and $\boldsymbolω \times \boldsymbol{J}$. \item[Results] We show that in self-consistent calculations, the Kerman-Onishi conditions are obeyed with high precision. Small deviations seen in the calculations with pairing can be attributed to the truncation of the quasiparticle spectrum. We also provide results of systematic TAC calculations for triaxial strongly deformed bands in $^{160}$Yb. \item[Conclusions] For non-stationary TAC solutions, Kerman-Onishi conditions link the non-zero values of the angle between rotational-frequency and angular-momentum vectors to the constraints on off-diagonal components of the quadrupole-moment tensor. To stabilize the convergence of self-consistent iterations, such constraints have to be taken into account. Only then one can determine the Routhian surfaces as functions of the tilting angles.

preprint2013arXiv

New density-independent interactions for nuclear structure calculations

We present a new two-body finite-range and momentum-dependent but density-independent effective interaction, which can be interpreted as a regularized zero-range force. We show that no three-body or density-dependent terms are needed for a correct description of saturation properties in infinite matter, that is, on the level of low-energy density functional, the physical three-body effects can be efficiently absorbed in effective two-body terms. The new interaction gives a very satisfying equation of state of nuclear matter and opens up extremely interesting perspectives for the mean-field and beyond-mean-field descriptions of atomic nuclei.

preprint2013arXiv

Polarization corrections to single-particle energies studied within the energy-density-functional and QRPA approaches

Background: Models based on using perturbative polarization corrections and mean-field blocking approximation give conflicting results for masses of odd nuclei. Purpose: Systematically investigate the polarization and mean-field models, implemented within self-consistent approaches that use identical interactions and model spaces, so as to find reasons for the conflicts between them. Methods: For density-dependent interactions and with pairing correlations included, we derive and study links between the mean-field and polarization results obtained for energies of odd nuclei. We also identify and discuss differences between the polarization-correction and full particle-vibration-coupling (PVC) models. Numerical calculations are performed for the mean-field ground-state properties of deformed odd nuclei and then compared to the polarization corrections determined by using the approach that conserves spherical symmetry. Results: We have identified and numerically evaluated self-interaction (SI) energies that are at the origin of different results obtained within the mean-field and polarization-correction approaches. Conclusions: Mean-field energies of odd nuclei are polluted by the SI energies, and this makes them different from those obtained by using polarization-correction methods. A comparison of both approaches allows for the identification and determination of the SI terms, which then can be calculated and removed from the mean-field results, giving the self-interaction-free energies. The simplest deformed mean-field approach that does not break parity symmetry is unable to reproduce full PVC effects.

preprint2013arXiv

Propagation of uncertainties in the Skyrme energy-density-functional model

Parameters of nuclear energy-density-functionals (EDFs) are always derived by an optimization to experimental data. For the minima of appropriately defined penalty functions, a statistical sensitivity analysis provides the uncertainties of the EDF parameters. To quantify theoretical errors of observables given by the model, we studied the propagation of uncertainties within the UNEDF0 Skyrme-EDF approach. We found that typically the standard errors rapidly increase towards neutron rich nuclei. This can be linked to large uncertainties of the isovector coupling constants of the currently used EDFs.

preprint2013arXiv

Spontaneous fission lifetimes from the minimization of self-consistent collective action

The spontaneous fission lifetime of 264Fm has been studied within nuclear density functional theory by minimizing the collective action integral for fission in a two-dimensional quadrupole collective space representing elongation and triaxiality. The collective potential and inertia tensor are obtained self-consistently using the Skyrme energy density functional and density-dependent pairing interaction. The resulting spontaneous fission lifetimes are compared with the static result obtained with the minimum-energy pathway. We show that fission pathways strongly depend on assumptions underlying collective inertia. With the non-perturbative mass parameters, the dynamic fission pathway becomes strongly triaxial and it approaches the static fission valley. On the other hand, when the standard perturbative cranking inertia tensor is used, axial symmetry is restored along the path to fission; an effect that is an artifact of the approximation used.

preprint2012arXiv

Effective theory for low-energy nuclear energy density functionals

We introduce a new class of effective interactions to be used within the energy-density-functional approaches. They are based on regularized zero-range interactions and constitute a consistent application of the effective-theory methodology to low-energy phenomena in nuclei. They allow for defining the order of expansion in terms of the order of derivatives acting on the finite-range potential. Numerical calculations show a rapid convergence of the expansion and independence of results of the regularization scale.

preprint2012arXiv

Giant Monopole Resonances and nuclear incompressibilities studied for the zero-range and separable pairing interactions

Background: Following the 2007 precise measurements of monopole strengths in tin isotopes, there has been a continuous theoretical effort to obtain a precise description of the experimental results. Up to now, there is no satisfactory explanation of why the tin nuclei appear to be significantly softer than 208Pb. Purpose: We determine the influence of finite-range and separable pairing interactions on monopole strength functions in semi-magic nuclei. Methods: We employ self-consistently the Quasiparticle Random Phase Approximation on top of spherical Hartree-Fock-Bogolyubov solutions. We use the Arnoldi method to solve the linear-response problem with pairing. Results: We found that the difference between centroids of Giant Monopole Resonances measured in lead and tin (about 1 MeV) always turns out to be overestimated by about 100%. We also found that the volume incompressibility, obtained by adjusting the liquid-drop expression to microscopic results, is significantly larger than the infinite-matter incompressibility. Conclusions: The zero-range and separable pairing forces cannot induce modifications of monopole strength functions in tin to match experimental data.

preprint2012arXiv

Isospin-breaking corrections to superallowed Fermi beta-decay in isospin- and angular-momentum-projected nuclear Density Functional Theory

Background: The superallowed beta-decay rates provide stringent constraints on physics beyond the Standard Model of particle physics. To extract crucial information about the electroweak force, small isospin-breaking corrections to the Fermi matrix element of superallowed transitions must be applied. Purpose: We perform systematic calculations of isospin-breaking corrections to superallowed beta-decays and estimate theoretical uncertainties related to the basis truncation, time-odd polarization effects related to the intrinsic symmetry of the underlying Slater determinants, and to the functional parametrization. Methods: We use the self-consistent isospin- and angular-momentum-projected nuclear density functional theory employing two density functionals derived from the density independent Skyrme interaction. Pairing correlations are ignored. Our framework can simultaneously describe various effects that impact matrix elements of the Fermi decay: symmetry breaking, configuration mixing, and long-range Coulomb polarization. Results: The isospin-breaking corrections to the I=0+,T=1 --> I=0+,T=1 pure Fermi transitions are computed for nuclei from A=10 to A=98 and, for the first time, to the Fermi branch of the I,T=1/2 --> I,T=1/2 transitions in mirror nuclei from A=11 to A=49. We carefully analyze various model assumptions impacting theoretical uncertainties of our calculations and provide theoretical error bars on our predictions. Conclusions: The overall agreement with empirical isospin-breaking corrections is very satisfactory. Using computed isospin-breaking corrections we show that the unitarity of the CKM matrix is satisfied with a precision better than 0.1%.

preprint2012arXiv

Precision Mass Measurements beyond $^{132}$Sn: Anomalous behaviour of odd-even staggering of binding energies

Atomic masses of the neutron-rich isotopes $^{121-128}$Cd, $^{129,131}$In, $^{130-135}$Sn, $^{131-136}$Sb, and $^{132-140}$Te have been measured with high precision (10 ppb) using the Penning trap mass spectrometer JYFLTRAP. Among these, the masses of four r-process nuclei $^{135}$Sn, $^{136}$Sb, and $^{139,140}$Te were measured for the first time. The data reveals a strong $N$=82 shell gap at $Z$=50 but indicates the importance of correlations for $Z>50$. An empirical neutron pairing gap expressed as the odd-even staggering of isotopic masses shows a strong quenching across $N$=82 for Sn, with the $Z$-dependence that is unexplainable by the current theoretical models.

preprint2011arXiv

Continuity equation and local gauge invariance for the N3LO nuclear Energy Density Functionals

Background: The next-to-next-to-next-to-leading order (N3LO) nuclear energy density functional extends the standard Skyrme functional with new terms depending on higher-order derivatives of densities, introduced to gain better precision in the nuclear many-body calculations. A thorough study of the transformation properties of the functional with respect to different symmetries is required, as a step preliminary to the adjustment of the coupling constants. Purpose: Determine to which extent the presence of higher-order derivatives in the functional can be compatible with the continuity equation. In particular, to study the relations between the validity of the continuity equation and invariance of the functional under gauge transformations. Methods: Derive conditions for the validity of the continuity equation in the framework of time-dependent density functional theory. The conditions apply separately to the four spin-isospin channels of the one-body density matrix. Results: We obtained four sets of constraints on the coupling constants of the N3LO energy density functional that guarantee the validity of the continuity equation in all spin-isospin channels. In particular, for the scalar-isoscalar channel, the constraints are the same as those resulting from imposing the standard U(1) local-gauge-invariance conditions. Conclusions: Validity of the continuity equation in the four spin-isospin channels is equivalent to the local-gauge invariance of the energy density functional. For vector and isovector channels, such validity requires the invariance of the functional under local rotations in the spin and isospin spaces.

preprint2011arXiv

Effective pseudopotential for energy density functionals with higher order derivatives

We derive a zero-range pseudopotential that includes all possible terms up to sixth order in derivatives. Within the Hartree-Fock approximation, it gives the average energy that corresponds to a quasi-local nuclear Energy Density Functional (EDF) built of derivatives of the one-body density matrix up to sixth order. The direct reference of the EDF to the pseudopotential acts as a constraint that divides the number of independent coupling constants of the EDF by two. This allows, e.g., for expressing the isovector part of the functional in terms of the isoscalar part, or vice versa. We also derive the analogous set of constraints for the coupling constants of the EDF that is restricted by spherical, space-inversion, and time-reversal symmetries.

preprint2011arXiv

Microscopic calculations of isospin-breaking corrections to superallowed beta-decay

The superallowed beta-decay rates that provide stringent constraints on physics beyond the Standard Model of particle physics are affected by nuclear structure effects through isospin-breaking corrections. The self-consistent isospin- and angular-momentum-projected nuclear density functional theory is used for the first time to compute those corrections for a number of Fermi transitions in nuclei from A=10 to A=74. The resulting leading element of the CKM matrix, |V_{ud}|= 0.97447(23), agrees well with the recent result by Towner and Hardy [Phys. Rev. C {\bf 77}, 025501 (2008)].

preprint2011arXiv

Self-consistent tilted-axis-cranking study of triaxial strongly deformed bands in $^{158}$Er at ultrahigh spin

Stimulated by recent experimental discoveries, triaxial strongly deformed (TSD) states in $^{158}$Er at ultrahigh spins have been studied by means of the Skyrme-Hartree-Fock model and the tilted-axis-cranking method. Restricting the rotational axis to one of the principal axes -- as done in previous cranking calculations -- two well-defined TSD minima in the total Routhian surface are found for a given configuration: one with positive and another with negative triaxial deformation $γ$. By allowing the rotational axis to change direction, the higher-energy minimum is shown to be a saddle point. This resolves the long-standing question of the physical interpretation of the two triaxial minima at a very similar quadrupole shape obtained in the principal axis cranking approach. Several TSD configurations have been predicted, including a highly deformed band expected to cross lesser elongated TSD bands at the highest spins. Its transitional quadrupole moment $Q_t \approx 10.5$\,eb is close to the measured value of $\sim$11\,eb; hence, it is a candidate for the structure observed in experiment.

preprint2011arXiv

Solution of the Skyrme-Hartree-Fock-Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis. (VII) HFODD (v2.49t): a new version of the program

We describe the new version (v2.49t) of the code HFODD which solves the nuclear Skyrme Hartree-Fock (HF) or Skyrme Hartree-Fock-Bogolyubov (HFB) problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we have implemented the following physics features: (i) the isospin mixing and projection, (ii) the finite temperature formalism for the HFB and HF+BCS methods, (iii) the Lipkin translational energy correction method, (iv) the calculation of the shell correction. A number of specific numerical methods have also been implemented in order to deal with large-scale multi-constraint calculations and hardware limitations: (i) the two-basis method for the HFB method, (ii) the Augmented Lagrangian Method (ALM) for multi-constraint calculations, (iii) the linear constraint method based on the approximation of the RPA matrix for multi-constraint calculations, (iv) an interface with the axial and parity-conserving Skyrme-HFB code HFBTHO, (v) the mixing of the HF or HFB matrix elements instead of the HF fields. Special care has been paid to using the code on massively parallel leadership class computers. For this purpose, the following features are now available with this version: (i) the Message Passing Interface (MPI) framework, (ii) scalable input data routines, (iii) multi-threading via OpenMP pragmas, (iv) parallel diagonalization of the HFB matrix in the simplex breaking case using the ScaLAPACK library. Finally, several little significant errors of the previous published version were corrected.

preprint2010arXiv

Convergence of density-matrix expansions for nuclear interactions

We extend density-matrix expansions in nuclei to higher orders in derivatives of densities and test their convergence properties. The expansions allow for converting the interaction energies characteristic to finite- and short-range nuclear effective forces into quasi-local density functionals. We also propose a new type of expansion that has excellent convergence properties when benchmarked against the binding energies obtained for the Gogny interaction.

preprint2010arXiv

Current Developments in Nuclear Density Functional Methods

Density functional theory (DFT) became a universal approach to compute ground-state and excited configurations of many-electron systems held together by an external one-body potential in condensed-matter, atomic, and molecular physics. At present, the DFT strategy is also intensely studied and applied in the area of nuclear structure. The nuclear DFT, a natural extension of the self-consistent mean-field theory, is a tool of choice for computations of ground-state properties and low-lying excitations of medium-mass and heavy nuclei. Over the past thirty-odd years, a lot of experience was accumulated in implementing, adjusting, and using the density-functional methods in nuclei. This research direction is still extremely actively pursued. In particular, current developments concentrate on (i) attempts to improve the performance and precision delivered by the nuclear density-functional methods, (ii) derivations of density functionals from first principles rooted in the low-energy chromodynamics and effective theories, and (iii) including effects of low-energy correlations and symmetry restoration. In this study, we present an overview of recent results and achievements gained in nuclear density-functional methods.

preprint2010arXiv

Isospin mixing in nuclei around N=Z and the superallowed beta-decay

Theoretical approaches that use one-body densities as dynamical variables, such as Hartree-Fock or the density functional theory (DFT), break isospin symmetry both explicitly, by virtue of charge-dependent interactions, and spontaneously. To restore the spontaneously broken isospin symmetry, we implemented the isospin-projection scheme on top of the Skyrme-DFT approach. This development allows for consistent treatment of isospin mixing in both ground and exited nuclear states. In this study, we apply this method to evaluate the isospin impurities in ground states of even-even and odd-odd N~Z nuclei. By including simultaneous isospin and angular-momentum projection, we compute the isospin-breaking corrections to the 0+ --> 0+ superallowed beta-decay.

preprint2010arXiv

Isospin mixing in the vicinity of the N=Z line

We present the isospin- and angular-momentum-projected nuclear density functional theory (DFT) and its applications to the isospin-breaking corrections to the superallowed beta-decay rates in the vicinity of the N=Z line. A preliminary value obtained for the Cabbibo-Kobayashi-Maskawa matrix element, |V_{ud}|=0.97463(24), agrees well with the recent estimate by Towner and Hardy [Phys. Rev. C{\bf 77}, 025501 (2008)]. We also discuss new opportunities to study the symmetry energy by using the isospin-projected DFT.

preprint2010arXiv

One-quasiparticle States in the Nuclear Energy Density Functional Theory

We study one-quasiproton excitations in the rare-earth region in the framework of the nuclear Density Functional Theory in the Skyrme-Hartree-Fock-Bogoliubov variant. The blocking prescription is implemented exactly, with the time-odd mean field fully taken into account. The equal filling approximation is compared with the exact blocking procedure. We show that both procedures are strictly equivalent when the time-odd channel is neglected, and discuss how nuclear alignment properties affect the time-odd fields. The impact of time-odd fields on calculated one-quasiproton bandhead energies is found to be rather small, of the order of 100-200 keV; hence, the equal filling approximation is sufficiently precise for most practical applications. The triaxial polarization of the core induced by the odd particle is studied. We also briefly discuss the occurrence of finite-size spin instabilities that are present in calculations for odd-mass nuclei when certain Skyrme functionals are employed.

preprint2010arXiv

The Negele-Vautherin density matrix expansion applied to the Gogny force

We use the Negele-Vautherin density matrix expansion to derive a quasi-local density functional for the description of systems of fermions interacting with short-ranged interactions composed of arbitrary finite-range central, spin-orbit, and tensor components. Terms that are absent in the original Negele-Vautherin approach owing to the angle averaging of the density matrix are fixed by employing a gauge invariance condition. We obtain the Kohn-Sham interaction energies in all spin-isospin channels, including the exchange terms, expressed as functions of the local densities and their derivatives up to second (next to leading) order. We illustrate the method by determining the coupling constants of the Skyrme functional or Skyrme force that correspond to the finite-range Gogny central force. The resulting self-consistent solutions reproduce the Gogny-force binding energies and radii within the precision of 1-2%.

preprint2010arXiv

Title: Quadrupole collective inertia in nuclear fission: cranking approximation

Collective mass tensor derived from the cranking approximation to the adiabatic time-dependent Hartree-Fock-Bogoliubov (ATDHFB) approach is compared with that obtained in the Gaussian Overlap Approximation (GOA) to the generator coordinate method. Illustrative calculations are carried out for one-dimensional quadrupole fission pathways in 256Fm. It is shown that the collective mass exhibits strong variations with the quadrupole collective coordinate. These variations are related to the changes in the intrinsic shell structure. The differences between collective inertia obtained in cranking and perturbative cranking approximations to ATDHFB, and within GOA, are discussed.

preprint2009arXiv

Isospin-symmetry restoration within the nuclear density functional theory: formalism and applications

Isospin symmetry of atomic nuclei is explicitly broken by the charge-dependent interactions, primarily the Coulomb force. Within the nuclear density functional theory, isospin is also broken spontaneously. We propose a projection scheme rooted in a mean field theory, that allows the consistent treatment of isospin breaking in both ground and exited nuclear states. We demonstrate that this scheme is essentially free from spurious divergences plaguing particle-number and angular-momentum restoration approaches. Applications of the new technique include excited high-spin states in medium-mass N=Z nuclei, such as superdeformed bands and many-particle-many-hole terminating states.

preprint2009arXiv

Linear response strength functions with iterative Arnoldi diagonalization

We report on an implementation of a new method to calculate RPA strength functions with iterative non-hermitian Arnoldi diagonalization method, which does not explicitly calculate and store the RPA matrix. We discuss the treatment of spurious modes, numerical stability, and how the method scales as the used model space is enlarged. We perform the particle-hole RPA benchmark calculations for double magic nucleus 132Sn and compare the resulting electromagnetic strength functions against those obtained within the standard RPA.

preprint2009arXiv

Self-consistent symmetries in the proton-neutron Hartree-Fock-Bogoliubov approach

Symmetry properties of densities and mean fields appearing in the nuclear Density Functional Theory with pairing are studied. We consider energy functionals that depend only on local densities and their derivatives. The most important self-consistent symmetries are discussed: spherical, axial, space-inversion, and mirror symmetries. In each case, the consequences of breaking or conserving the time-reversal and/or proton-neutron symmetries are discussed and summarized in a tabulated form, useful in practical applications.

preprint2009arXiv

Solution of self-consistent equations for the N3LO nuclear energy density functional in spherical symmetry. The program HOSPHE (v1.00)

We present solution of self-consistent equations for the N3LO nuclear energy density functional. We derive general expressions for the mean fields expressed as differential operators depending on densities and for the densities expressed in terms of derivatives of wave functions. These expressions are then specified to the case of spherical symmetry. We also present the computer program HOSPHE (v1.00), which solves the self-consistent equations by using the expansion of single-particle wave functions on the spherical harmonic oscillator basis.

preprint2008arXiv

Local nuclear energy density functional at next-to-next-to-next-to-leading order

We construct nuclear energy density functionals in terms of derivatives of densities up to sixth, next-to-next-to-next-to-leading order (N3LO). A phenomenological functional built in this way conforms to the ideas of the density matrix expansion and is rooted in the expansions characteristic to effective theories. It builds on the standard functionals related to the contact and Skyrme forces, which constitute the zero-order (LO) and second-order (NLO) expansions, respectively. At N3LO, the full functional with density-independent coupling constants, and with the isospin degree of freedom taken into account, contains 376 terms, while the functionals restricted by the Galilean and gauge symmetries contain 100 and 42 terms, respectively. For functionals additionally restricted by the spherical, space-inversion, and time-reversal symmetries, the corresponding numbers of terms are equal to 100, 60, and 22, respectively.

preprint2008arXiv

Shell-structure fingerprints of tensor interaction

We address consequences of strong tensor and weak spin-orbit terms in the local energy density functional, resulting from fits to the $f_{5/2} - f_{7/2}$ splittings in $^{40}$Ca, $^{48}$Ca, and $^{56}$Ni. In this study, we focus on nuclear binding energies. In particular, we show that the tensor contribution to the binding energies exhibits interesting topological features closely resembling that of the shell-correction. We demonstrate that in the extreme single-particle scenario at spherical shape, the tensor contribution shows tensorial magic numbers equal to $N(Z)$=14, 32, 56, and 90, and that this structure is smeared out due to configuration mixing caused by pairing correlations and migration of proton/neutron sub-shells with neutron/proton shell filling. Based on a specific Skyrme-type functional SLy4$_T$, we show that the proton tensorial magic numbers shift with increasing neutron excess to $Z$=14, 28, and 50.

preprint2006arXiv

Large-Scale Self-Consistent Nuclear Mass Calculations

The program of systematic large-scale self-consistent nuclear mass calculations that is based on the nuclear density functional theory represents a rich scientific agenda that is closely aligned with the main research directions in modern nuclear structure and astrophysics, especially the radioactive nuclear beam physics. The quest for the microscopic understanding of the phenomenon of nuclear binding represents, in fact, a number of fundamental and crucial questions of the quantum many-body problem, including the proper treatment of correlations and dynamics in the presence of symmetry breaking. Recent advances and open problems in the field of nuclear mass calculations are presented and discussed.

preprint2004arXiv

Self-Consistent Description of Multipole Strength in Exotic Nuclei I: Method

We use the canonical Hartree-Fock-Bogoliubov basis to implement a completely self-consistent quasiparticle-random-phase approximation with arbitrary Skyrme energy density functionals and density-dependent pairing functionals. The point of the approach is to accurately describe multipole strength functions in spherical even-even nuclei, including weakly-bound drip-line systems. We describe the method and carefully test its accuracy, particularly in handling spurious modes. To illustrate our approach, we calculate isoscalar and isovector monopole, dipole, and quadrupole strength functions in several Sn isotopes, both in the stable region and at the drip lines.

preprint2003arXiv

Systematic study of deformed nuclei at the drip lines and beyond

An improved prescription for choosing a transformed harmonic oscillator (THO) basis for use in configuration-space Hartree-Fock-Bogoliubov (HFB) calculations is presented. The new HFB+THO framework that follows accurately reproduces the results of coordinate-space HFB calculations for spherical nuclei, including those that are weakly bound. Furthermore, it is fully automated, facilitating its use in systematic investigations of large sets of nuclei throughout the periodic table. As a first application, we have carried out calculations using the Skyrme Force SLy4 and volume pairing, with exact particle number projection following application of the Lipkin-Nogami prescription. Calculations were performed for all even-even nuclei from the proton drip line to the neutron drip line having proton numbers Z=2,4,...,108 and neutron numbers N=2,4,...,188. We focus on nuclei near the neutron drip line and find that there exist numerous particle-bound even-even nuclei (i.e., nuclei with negative Fermi energies) that have at the same time negative two-neutron separation energies. This phenomenon, which was earlier noted for light nuclei, is attributed to bound shape isomers beyond the drip line.

preprint2001arXiv

Pairing Interaction and Self-Consistent Densities in Neutron-Rich Nuclei

Particle and pairing densities in spherical even-even neutron-rich nuclei are studied within the Skyrme-Hartree-Fock-Bogoliubov approach with the density-dependent pairing interaction. The influence of the density dependence of the pairing interaction on asymptotic properties of nucleonic distributions are analyzed. It is demonstrated that the size of the neutron halo dramatically depends on the behavior of the pairing interaction at low density.

preprint1999arXiv

Quadrupole deformations of neutron-drip-line nuclei studied within the Skyrme Hartree-Fock-Bogolyubov approach

We introduce a local-scaling point transformation to allow for modifying the asymptotic properties of the deformed three-dimensional Cartesian harmonic oscillator wave functions. The resulting single-particle bases are very well suited for solving the Hartree-Fock-Bogoliubov equations for deformed drip-line nuclei. We then present results of self-consistent calculations performed for the Mg isotopes and for light nuclei located near the two-neutron drip line. The results suggest that for all even-even elements with $Z$=10--18 the most weakly-bound nucleus has an oblate ground-state shape.

preprint1999arXiv

Structure of nuclei at extreme values of the isospin

Physics of nuclei at extreme values of the isospin is at the focus of present-day nuclear science. Experimentally, thanks to existing and emerging radioactive-ion-beam facilities, we are on the verge of invading the territory of extreme N/Z ratios in an unprecedented way. Theoretically, nuclear exotica represent a formidable challenge for the nuclear many-body theories and their power to predict nuclear properties far from stability. Going to the limits of the nuclear binding is also important for an improvement of our description of normal nuclei from the neighborhood of the beta stability valley. In the present talk, we review several aspects of the present-day mean-field theoretical studies of weakly bound nuclei.

preprint1995arXiv

Time-odd components in the mean field of rotating superdeformed nuclei

Rotation-induced time-odd components in the nuclear mean field are analyzed using the Hartree-Fock cranking approach with effective interactions SIII, SkM*, and SkP. Identical dynamical moments ${{\cal J}^{(2)}}$ are obtained for pairs of superdeformed bands $^{151}$Tb(2)--$^{152}$Dy(1) and $^{150}$Gd(2)--$^{151}$Tb(1). The corresponding relative alignments strongly depend on which time-odd mean-field terms are taken into account in the Hartree-Fock equations.

preprint1994arXiv

Hexadecapole Interaction and the Delta I=4 Staggering Effect in Rotational Bands

A role of the multipole interaction in the description of the $\triangle I$=4 staggering phenomenon is investigated in a model consisting of a single-$j$ shell filled by identical nucleons. Exact diagonalization of the quadrupole-plus-hexadecapole Hamiltonian shows that the hexadecapole-hexadecapole interaction can produce a $\triangle I$=4 periodicity in the yrast sequence.

preprint1994arXiv

SDG fermion-pair algebraic SO(12) and Sp(10) models and their boson realizations

It is shown how the boson mapping formalism may be applied as a useful many-body tool to solve a fermion problem. This is done in the context of generalized Ginocchio models for which we introduce S-, D-, and G-pairs of fermions and subsequently construct the sdg-boson realizations of the generalized Dyson type. The constructed SO(12) and Sp(10) fermion models are solved beyond the explicit symmetry limits. Phase transitions to rotational structures are obtained, also in situations where there is no underlying SU(3) symmetry.

preprint1993arXiv

Superdeformed rotational bands in the Mercury region; A Cranked Skyrme-Hartree-Fock-Bogoliubov study

A study of rotational properties of the ground superdeformed bands in \Hg{0}, \Hg{2}, \Hg{4}, and \Pb{4} is presented. We use the cranked Hartree-Fock-Bogoliubov method with the {\skm} parametrization of the Skyrme force in the particle-hole channel and a seniority interaction in the pairing channel. An approximate particle number projection is performed by means of the Lipkin-Nogami prescription. We analyze the proton and neutron quasiparticle routhians in connection with the present information on about thirty presently observed superdeformed bands in nuclei close neighbours of \Hg{2}.