BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

W. Nazarewicz

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

ResearcherUnclaimed source record
nucl-thnucl-exastro-ph.SRphysics.atm-clusphysics.comp-phcond-mat.supr-conquant-phastro-ph.COastro-ph.GAastro-ph.HEcond-mat.quant-gasphysics.atom-phphysics.chem-ph

Catalog footprint

What is connected

92works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

92 published item(s)

preprint2023arXiv

Pushing the Limits of the Periodic Table -- A Review on Atomic Relativistic Electronic Structure Theory and Calculations for the Superheavy Elements

We review the progress in atomic structure theory with a focus on superheavy elements and the aim to predict their ground state configuration and element's placement in the periodic table. To understand the electronic structure and correlations in the regime of large atomic numbers, it is important to correctly solve the Dirac equation in strong Coulomb fields, and also to take into account quantum electrodynamic effects. We specifically focus on the fundamental difficulties encountered when dealing with the many-particle Dirac equation. We further discuss the possibility for future many-electron atomic structure calculations going beyond the critical nuclear charge \(Z_{\rm crit}\approx 170\), where levels such as the \(1s\) shell dive into the negative energy continuum (\(E_{nκ}<-m_ec^2\)). The nature of the resulting Gamow states within a rigged Hilbert space formalism is highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

$β^-{\rm p}$ and $β^-α$ decay of the $^{11}$Be neutron halo ground state

Beta-delayed proton emission from the neutron halo ground state of $^{11}$Be raised much attention due to the unusually high decay rate. It was argued that this may be due to the existence of a resonance just above the proton decay threshold. In this Letter, we use the lenses of real-energy continuum shell model to describe several observables including the Gamow-Teller rates for the $β^-$-delayed $α$ and proton decays, and argue that, within our model, the large $β^-{\rm p}$ branching ratio cannot be reconciled with other data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni

Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nucleon-nucleon correlations in the extreme oxygen isotopes

There has been an upsurge of interest in two-nucleon decays thanks to the studies of nucleon-nucleon correlations. In our previous work, based on a novel time-dependent three-body approach, we demonstrated that the energy and angular correlations of the emitted nucleons can shed light on the structure of nucleonic pairs formed inside the nucleus. In this work, we apply the new framework to study the decay dynamics and properties of some extreme proton-rich and neutron-rich oxygen isotopes, including two-proton ($2p$) decays of $^{11,12}$O and two-neutron ($2n$) decay of $^{26}$O. Here we show that the low-$\ell$ components of $^{11,12}$O wave functions, which are affected by continuum and configuration-interaction effects, strongly impact decay dynamics and asymptotic correlations. In the calculated wave functions of $^{11,12}$O, diproton and cigarlike structures merge together during the tunneling process and the resulting energy- and angular correlations are very consistent with the experimental data. The asymptotic correlations of the $2n$ decay of $^{26}$O dramatically change as the two-neutron decay energy approaches the zero-energy threshold. The small reported value of $Q_{2n}$ suggests that the $2n$ decay of this nucleus can be understood in terms of the universal phase-space limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy Spectrum of Neutron-Rich Helium Isotopes: Complex Made Simple

We demonstrate that the intricate energy spectrum of neutron-rich helium isotopes can be straightforwardly described by taking advantage of the low-energy properties of neutron-neutron interaction and the scale separation that is present in diluted dripline systems. By using arguments based on the halo effective field theory, we carry out a parameter reduction of the complex-energy configuration interaction framework in the $spd$ space, including resonant and scattering states. By adjusting only one parameter, the strength of the spin-singlet central neutron-neutron interaction, we reproduce experimental energies and widths of $^{5-8}$He within tens of keV precision. We predict a parity inversion of narrow resonances in $^{9}$He and show that the ground state of $^{10}$He is an $s$-wave-dominated threshold configuration that could decay through two-neutron emission.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gamow Shell Model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we construct the minimal complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 $\leq$ A $\leq$ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid $^4$He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the sp space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: Our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei $^{10}$Li-$^{10}$N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measurement and microscopic description of odd-even staggering of charge radii of exotic copper isotopes

The mesoscopic nature of the atomic nucleus gives rise to a wide array of macroscopic and microscopic phenomena. The size of the nucleus is a window into this duality: while the charge radii globally scale as $A^{1/3}$, their evolution across isotopic chains reveals unanticipated structural phenomena [1-3]. The most ubiquitous of these is perhaps the Odd-Even Staggering (OES) [4]: isotopes with an odd number of neutrons are usually smaller in size than the trend of their even-neutron neighbours suggests. This OES effect varies with the number of protons and neutrons and poses a significant challenge for nuclear theory [5-7]. Here, we examine this problem with new measurements of the charge radii of short-lived copper isotopes up to the very exotic $^{78}$Cu $(Z=29, N=49)$, produced at only 20 ions/s, using the highly-sensitive Collinear Resonance Ionisation Spectroscopy (CRIS) method at ISOLDE-CERN. Due to the presence of a single proton outside of the closed Z=28 shell, these measurements provide crucial insights into the single-particle proton structure and how this affects the charge radii. We observe an unexpected reduction in the OES for isotopes approaching the $N=50$ shell gap. To describe the data, we applied models based on nuclear Density Functional Theory [2,8] (DFT) and ab-initio Valence-Space In-Medium Similarity Renormalization Group (VS-IMSRG) theory [9,10]. Through these comparisons, we demonstrate a relation between the global behavior of charge radii and the saturation density of nuclear matter, and show that the local charge radii variations, which reflect the many-body polarization effects due to the odd neutron, naturally emerge from the VS-IMSRG calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Beyond the charge radius: the information content of the fourth radial moment

Measurements of atomic transitions in different isotopes offer key information on the nuclear charge radius. The anticipated high-precision experimental techniques, augmented by atomic calculations, will soon enable extraction of the higher-order radial moments of the charge density distribution. To assess the value of such measurements for nuclear structure research, we study the information content of the fourth radial moment $\rf$ by means of nuclear density functional theory and a multiple correlation analysis. We show that $\rf$ can be directly related to the surface thickness of nuclear density, a fundamental property of the atomic nucleus that is difficult to obtain for radioactive systems. Precise knowledge of these radial moments is essential to establish reliable constraints on the existence of new forces from precision isotope shift measurements .

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Convenient location of a near-threshold proton-emitting resonance in $^{11}$B

The presence of cluster-like narrow resonances in the vicinity of reaction/decay thresholds is a ubiquitous phenomenon with profound consequences. We argue that the continuum coupling, present in the open quantum system description of the atomic nucleus, can profoundly impact the nature of near-threshold states. In this Letter, we discuss the structure of the recently observed near-threshold resonance in $^{11}$B, whose very existence explains the puzzling beta-delayed proton emission of the neutron-rich $^{11}$Be.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear charge and neutron radii and nuclear matter: trend analysis

Radii of charge and neutron distributions are fundamental nuclear properties. They depend on both nuclear interaction parameters related to the equation of state of infinite nuclear matter and on quantal shell effects, which are strongly impacted by the presence of nuclear surface. In this work, by studying the dependence of charge and neutron radii, and neutron skin, on nuclear matter parameters, we assess different mechanisms that drive nuclear sizes. We apply nuclear density functional theory using a family of Skyrme functionals obtained by means of different optimization protocols targeting specific nuclear properties. By performing the Monte-Carlo sampling of reasonable functionals around the optimal parametrization, we study correlations between nuclear matter paramaters and observables characterizing charge and neutron distributions. We demonstrate the existence of the strong converse relation between the nuclear charge radii and the saturation density of symmetric nuclear matter and also between the neutron skins and the slope of the symmetry energy. For functionals optimized to experimental binding energies only, proton and neutron radii are weakly correlated due to canceling trends from different nuclear matter parameters. We show that by requiring that the nuclear functional reproduces the empirical saturation point of symmetric nuclear matter practically fixes the charge (or proton) radii, and vice versa. The neutron skin uncertainty primarily depends on the slope of the symmetry energy. Consequently, imposing a constraint on both $ρ_0$ and $L$ practically determines the nuclear size, modulo small variations due to shell effects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear rotation in the continuum

${\textbf{Background:}}$ Atomic nuclei often exhibit collective rotational-like behavior in highly excited states, well above the particle emission threshold. What determines the existence of collective motion in the continuum region, is not fully understood. ${\textbf{Purpose:}}$ In this work, by studying the collective rotation of the positive-parity deformed configurations of the one-neutron halo nucleus $^{11}$Be, we assess different mechanisms that stabilize collective behavior beyond the limits of particle stability. ${\textbf{Method:}}$ To solve a particle-plus-core problem, we employ a non-adiabatic coupled-channel formalism and the Berggren single-particle ensemble, which explicitly contains bound states, narrow resonances, and the scattering continuum. We study the valence-neutron density in the intrinsic rotor frame to assess the validity of the adiabatic approach as the excitation energy increases. ${\textbf{Results:}}$ We demonstrate that collective rotation of the ground band of $^{11}$Be is stabilized by (i) the fact that the $\ell=0$ one-neutron decay channel is closed, and (ii) the angular momentum alignment, which increases the parentage of high-$\ell$ components at high spins; both effects act in concert to decrease decay widths of ground-state band members. This is not the case for higher-lying states of $^{11}$Be, where the $\ell=0$ neutron-decay channel is open and often dominates. ${\textbf{Conclusion:}}$ We demonstrate that long-lived collective states can exist at high excitation energy in weakly bound neutron drip-line nuclei such as $^{11}$Be.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Resonant spectra of quadrupolar anions

In quadrupole-bound anions, an extra electron is attached at a sufficiently large quadrupole moment of a neutral molecule, which is lacking a permanent dipole moment. The nature of the bound states and low-lying resonances of such anions is of interest for understanding the threshold behavior of open quantum systems in general. In this work, we investigate the properties of quadrupolar anions as extreme halo systems, the formation of rotational bands, and the transition from a subcritical to supercritical electric quadrupole moment. We solve the electron-plus-molecule problem using a non-adiabatic coupled-channel formalism by employing the Berggren ensemble, which explicitly contains bound states, narrow resonances, and the scattering continuum. We demonstrate that binding energies and radii of quadrupolar anions strictly follow the scaling laws for two-body halo systems. Contrary to the case of dipolar anions, ground-state band of quadrupolar anions smoothly extend into the continuum, and many rotational bands could be identified above the detachment threshold. We study the evolution of a bound state of an anion as to dives into the continuum at a critical quadrupole moment and we show that the associated critical exponent is consistent with the second-order phase transition. Everything considered, quadrupolar anions represent a perfect laboratory for the studies of marginally bound open quantum systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Single-particle and collective motion in unbound deformed $ ^{39}\text{Mg} $

Background: Deformed neutron-rich magnesium isotopes constitute a fascinating territory where the interplay between collective rotation and single-particle motion is strongly affected by the neutron continuum. The unbound $fp$-shell nucleus $ ^{39}\text{Mg} $ is an ideal candidate to study this interplay. Purpose: In this work, we predict the properties of low-lying resonant states of $ ^{39}\text{Mg} $, using a suite of realistic theoretical approaches rooted in the open quantum system framework. Method: To describe the spectrum and decay modes of $ ^{39}\text{Mg} $ we use the conventional Shell Model, Gamow Shell Model, Resonating Group Method, Density Matrix Renormalization Group method, and the non-adiabatic Particle-Plus-Rotor model formulated in the Berggren basis. Results: The unbound ground state of $ ^{39}\text{Mg} $ is predicted to be either a $ J^π = {7/2}^- $ state or a $ {3/2}^- $ state. A narrow $ J^π = {7/2}^- $ ground-state candidate exhibits a resonant structure reminiscent of that of its one-neutron halo neighbor $ ^{37}\text{Mg} $, which is dominated by the $ f_{7/2} $ partial wave at short distances and a $ p_{3/2} $ component at large distances. A $ J^π={3/2}^- $ ground-state candidate is favored by the large deformation of the system. It can be associated with the $ {1/2}^- [321] $ Nilsson orbital dominated by the $\ell=1$ wave; hence its predicted width is large. The excited $J^π = {1/2}^-$ and $5/2^-$ states are expected to be broad resonances, while the $ J^π = {9/2}^- $ and $ {11/2}^- $ members of the ground-state rotational band are predicted to have very small neutron decay widths. Conclusion: We demonstrate that the subtle interplay between deformation, shell structure, and continuum coupling can result in a variety of excitations in an unbound nucleus just outside the neutron drip line.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Time-dependent density functional theory with twist-averaged boundary conditions

Time-dependent density functional theory is widely used to describe excitations of many-fermion systems. In its many applications, 3D coordinate-space representation is used, and infinite-domain calculations are limited to a finite volume represented by a box. For finite quantum systems (atoms, molecules, nuclei), the commonly used periodic or reflecting boundary conditions introduce spurious quantization of the continuum states and artificial reflections from boundary; hence, an incorrect treatment of evaporated particles. These artifacts can be practically cured by introducing absorbing boundary conditions (ABC) through an absorbing potential in a certain boundary region sufficiently far from the described system. But also the calculations of infinite matter (crystal electrons, quantum fluids, neutron star crust) suffer artifacts from a finite computational box. In this regime, twist- averaged boundary conditions (TABC) have been used successfully to diminish the finite-volume effects. In this work, we extend TABC to time-dependent framework and apply it to resolve the box artifacts for finite quantum systems using as test case small- and large-amplitude nuclear vibrations. We demonstrate that by using such a method, one can reduce finite volume effects drastically without adding any additional parameters. While they are almost equivalent in the linear regime, TABC and ABC differ in the nonlinear regime in their treatment of evaporated particles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Twist-averaged boundary conditions for nuclear pasta Hartree-Fock calculations

Background: Nuclear pasta phases, present in the inner crust of neutron stars, are associated with nucleonic matter at sub-saturation densities arranged in regular shapes. Those complex phases, residing in a layer which is approximately 100 m thick, impact many features of neutron stars. Theoretical quantum-mechanical simulations of nuclear pasta are usually carried out in finite 3D boxes assuming periodic boundary conditions (PBC). The resulting solutions are affected by spurious finite-size effects. Purpose: In order to remove spurious finite-size effects, it is convenient to employ twist-averaged boundary conditions (TABC) used in condensed matter, nuclear matter, and lattice QCD applications. In this work, we study the effectiveness of TABC in the context of pasta phases simulations within nuclear density functional theory. Methods: We perform Skyrme-Hartree-Fock calculations in three dimensions by implementing Bloch boundary conditions. The TABC averages are obtained by means of Gauss-Legendre integration over twist angles. Results: We benchmark the TABC for a free nucleonic gas and apply it to simple cases such as the rod and slab phases, as well as to more elaborate P-surface and gyroidal phases. Conclusions: We demonstrate that by applying TABC reliable results can be obtained from calculations performed in relatively small volumes. By studying various contributions to the total energy, we gain insights into pasta phases in mid-density range.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unexpectedly large charge radii of neutron-rich calcium isotopes

Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successfully described by nuclear theory, it is still a challenge to predict their charge radii evolution. Here we present the first measurements of the charge radii of 49,51,52Ca, obtained from laser spectroscopy experiments at ISOLDE, CERN. The experimental results are complemented by state-of-the-art theoretical calculations. The large and unexpected increase of the size of the neutron-rich calcium isotopes beyond N = 28 challenges the doubly-magic nature of 52Ca and opens new intriguing questions on the evolution of nuclear sizes away from stability, which are of importance for our understanding of neutron-rich atomic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate nuclear radii and binding energies from a chiral interaction

With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO$_{\rm sat}$, yield accurate binding energies and radii of nuclei up to $^{40}$Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective $J^π=3^-$ states in $^{16}$O and $^{40}$Ca are described accurately, while spectra for selected $p$- and $sd$-shell nuclei are in reasonable agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Charge, neutron, and weak size of the atomic nucleus

What is the size of the atomic nucleus? This deceivably simple question is difficult to answer. While the electric charge distributions in atomic nuclei were measured accurately already half a century ago, our knowledge of the distribution of neutrons is still deficient. In addition to constraining the size of atomic nuclei, the neutron distribution also impacts the number of nuclei that can exist and the size of neutron stars. We present an ab initio calculation of the neutron distribution of the neutron-rich nucleus $^{48}$Ca. We show that the neutron skin (difference between radii of neutron and proton distributions) is significantly smaller than previously thought. We also make predictions for the electric dipole polarizability and the weak form factor; both quantities are currently targeted by precision measurements. Based on ab initio results for $^{48}$Ca, we provide a constraint on the size of a neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multipole modes in deformed nuclei within the finite amplitude method

Background: To access selected excited states of nuclei, within the framework of nuclear density functional theory, the quasiparticle random phase approximation (QRPA) is commonly used. Purpose: We present a computationally efficient, fully self-consistent framework to compute the QRPA transition strength function of an arbitrary multipole operator in axially-deformed superfluid nuclei. Methods: The method is based on the finite amplitude method (FAM) QRPA, allowing fast iterative solution of QRPA equations. A numerical implementation of the FAM-QRPA solver module has been carried out for deformed nuclei. Results: The practical feasibility of the deformed FAM module has been demonstrated. In particular, we calculate the quadrupole and octupole strength in a heavy deformed nucleus $^{240}$Pu, without any truncations in the quasiparticle space. To demonstrate the capability to calculate individual QRPA modes, we also compute low-lying negative-parity collective states in $^{154}$Sm. Conclusions: The new FAM implementation enables calculations of the QRPA strength function throughout the nuclear landscape. This will facilitate global surveys of multipole modes and beta decays, and will open new avenues for constraining the nuclear energy density functional.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Uncertainty Quantification for Nuclear Density Functional Theory and Information Content of New Measurements

Statistical tools of uncertainty quantification can be used to assess the information content of measured observables with respect to present-day theoretical models; to estimate model errors and thereby improve predictive capability; to extrapolate beyond the regions reached by experiment; and to provide meaningful input to applications and planned measurements. To showcase new opportunities offered by such tools, we make a rigorous analysis of theoretical statistical uncertainties in nuclear density functional theory using Bayesian inference methods. By considering the recent mass measurements from the Canadian Penning Trap at Argonne National Laboratory, we demonstrate how the Bayesian analysis and a direct least-squares optimization, combined with high-performance computing, can be used to assess the information content of the new data with respect to a model based on the Skyrme energy density functional approach. Employing the posterior probability distribution computed with a Gaussian process emulator, we apply the Bayesian framework to propagate theoretical statistical uncertainties in predictions of nuclear masses, two-neutron dripline, and fission barriers. Overall, we find that the new mass measurements do not impose a constraint that is strong enough to lead to significant changes in the model parameters. The example discussed in this study sets the stage for quantifying and maximizing the impact of new measurements with respect to current modeling and guiding future experimental efforts, thus enhancing the experiment-theory cycle in the scientific method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bound and resonance states of the dipolar anion of hydrogen cyanide: competition between threshold effects and rotation in an open quantum system

Bound and resonance states of the dipole-bound anion of hydrogen cyanide HCN$^-$ are studied using a non-adiabatic pseudopotential method and the Berggren expansion technique involving bound states, decaying resonant states, and non-resonant scattering continuum. We devise an algorithm to identify the resonant states in the complex energy plane. To characterize spatial distributions of electronic wave functions, we introduce the body-fixed density and use it to assign families of resonant states into collective rotational bands. We find that the non-adiabatic coupling of electronic motion to molecular rotation results in a transition from the strong-coupling to weak-coupling regime. In the strong coupling limit, the electron moving in a subthreshold, spatially extended halo state follows the rotational motion of the molecule. Above the ionization threshold, electron's motion in a resonance state becomes largely decoupled from molecular rotation. Widths of resonance-band members depend primarily on the electron orbital angular momentum.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Excitation energy dependence of fission in the mercury region

Background: Recent experiments on beta-delayed fission reported an asymmetric mass yield in the neutron-deficient nucleus 180Hg. Earlier experiments in the mass region A=190-200 close to the beta-stability line, using the (p,f) and (α,f) reactions, observed a more symmetric distribution of fission fragments. While the beta-delayed fission of 180Hg can be associated with relatively low excitation energy, this is not the case for light-ion reactions, which result in warm compound nuclei. Purpose: To elucidate the roles of proton and neutron numbers and excitation energy in determining symmetric and asymmetric fission yields, we compute and analyze the isentropic potential energy surfaces of 174,180,198Hg and 196,210Po. Methods: We use the finite-temperature superfluid nuclear density functional theory, for excitation energies up to E*=30MeV and zero angular momentum. For our theoretical framework, we consider the Skyrme energy density functional SkM* and a density-dependent pairing interaction. Results: For 174,180Hg, we predict fission pathways consistent with asymmetric fission at low excitation energies, with the symmetric fission pathway opening very gradually as excitation energy is increased. For 198Hg and 196Po, we expect the nearly-symmetric fission channel to dominate. 210Po shows a preference for a slightly asymmetric pathway at low energies, and a preference for a symmetric pathway at high energies. Conclusions: Our self-consistent theory suggests that excitation energy weakly affects the fission pattern of the nuclei considered. The transition from the asymmetric fission in the proton-rich nuclei to a more symmetric fission in the heavier isotopes is governed by the shell structure of pre-scission configurations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear energy density optimization: Shell structure

Nuclear density functional theory is the only microscopical theory that can be applied throughout the entire nuclear landscape. Its key ingredient is the energy density functional. In this work, we propose a new parameterization UNEDF2 of the Skyrme energy density functional. The functional optimization is carried out using the POUNDerS optimization algorithm within the framework of the Skyrme Hartree-Fock-Bogoliubov theory. Compared to the previous parameterization UNEDF1, restrictions on the tensor term of the energy density have been lifted, yielding a very general form of the energy density functional up to second order in derivatives of the one-body density matrix. In order to impose constraints on all the parameters of the functional, selected data on single-particle splittings in spherical doubly-magic nuclei have been included into the experimental dataset. The agreement with both bulk and spectroscopic nuclear properties achieved by the resulting UNEDF2 parameterization is comparable with UNEDF1. While there is a small improvement on single-particle spectra and binding energies of closed shell nuclei, the reproduction of fission barriers and fission isomer excitation energies has degraded. As compared to previous UNEDF parameterizations, the parameter confidence interval for UNEDF2 is narrower. In particular, our results overlap well with those obtained in previous systematic studies of the spin-orbit and tensor terms. UNEDF2 can be viewed as an all-around Skyrme EDF that performs reasonably well for both global nuclear properties and shell structure. However, after adding new data aiming to better constrain the nuclear functional, its quality has improved only marginally. These results suggest that the standard Skyrme energy density has reached its limits and significant changes to the form of the functional are needed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Energy Density Optimization: UNEDF2

The parameters of the UNEDF2 nuclear energy density functional (EDF) model were obtained in an optimization to experimental data consisting of nuclear binding energies, proton radii, odd-even mass staggering data, fission-isomer excitation energies, and single particle energies. In addition to parameter optimization, sensitivity analysis was done to obtain parameter uncertainties and correlations. The resulting UNEDF2 is an all-around EDF. However, the sensitivity analysis also demonstrated that the limits of current Skyrme-like EDFs have been reached and that novel approaches are called for.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Theory and Science of the Facility for Rare Isotope Beams

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An optimized chiral nucleon-nucleon interaction at next-to-next-to-leading order

We optimize the nucleon-nucleon interaction from chiral effective field theory at next-to-next- to-leading order. The resulting new chiral force NNLOopt yields χ^2 \approx 1 per degree of freedom for laboratory energies below approximately 125 MeV. In the A = 3, 4 nucleon systems, the contributions of three-nucleon forces are smaller than for previous parametrizations of chiral interactions. We use NNLOopt to study properties of key nuclei and neutron matter, and demonstrate that many aspects of nuclear structure can be understood in terms of this nucleon-nucleon interaction, without explicitly invoking three-nucleon forces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bound states of dipolar molecules studied with the Berggren expansion method

Bound states of dipole-bound negative anions are studied by using a non-adiabatic pseudopotential method and the Berggren expansion involving bound states, decaying resonant states, and non-resonant scattering continuum. The method is benchmarked by using the traditional technique of direct integration of coupled channel equations. A good agreement between the two methods has been found for well-bound states. For weakly-bound subthreshold states with binding energies comparable with rotational energies of the anion, the direct integration approach breaks down and the Berggren expansion method becomes the tool of choice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Description of 158Er at ultrahigh spin in nuclear density functional theory

Rotational bands in $^{158}$Er at ultra-high spin have been studied in the framework of relativistic and non-relativistic nuclear density functional theories. Consistent results are obtained across the theoretical models used but some puzzles remain when confronted with experiment. Namely, the many-body configurations which provide good description of experimental transition quadrupole moments and dynamic moments of inertia require substantial increase of the spins of observed bands as compared with experimental estimates, which are still subject to large uncertainties. If, however, the theoretical spins assignments turned out to be correct, the experimental band 1 in $^{158}$Er would be the highest-spin structure ever observed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Information content of the weak-charge form factor

Parity-violating electron scattering provides a model-independent determination of the nuclear weak-charge form factor that has widespread implications across such diverse areas as fundamental symmetries, nuclear structure, heavy-ion collisions, and neutron-star structure. We assess the impact of precise measurements of the weak-charge form factor of ${}^{48}$Ca and ${}^{208}$Pb on a variety of nuclear observables, such as the neutron skin and the electric-dipole polarizability. We use the nuclear Density Functional Theory with several accurately calibrated non-relativistic and relativistic energy density functionals. To assess the degree of correlation between nuclear observables and to explore systematic and statistical uncertainties on theoretical predictions, we employ the chi-square statistical covariance technique. We find a strong correlation between the weak-charge form factor and the neutron radius, that allows for an accurate determination of the neutron skin of neutron-rich nuclei. We determine the optimal range of the momentum transfer $q$ that maximizes the information content of the measured weak-charge form factor and quantify the uncertainties associated with the strange quark contribution. Moreover, we confirm the role of the electric-dipole polarizability as a strong isovector indicator. Accurate measurements of the weak-charge form factor of ${}^{48}$Ca and ${}^{208}$Pb will have a profound impact on many aspects of nuclear theory and hadronic measurements of neutron skins of exotic nuclei at radioactive-beam facilities.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutron skin uncertainties of Skyrme energy density functionals

Background: Neutron-skin thickness is an excellent indicator of isovector properties of atomic nuclei. As such, it correlates strongly with observables in finite nuclei that depend on neutron-to-proton imbalance and the nuclear symmetry energy that characterizes the equation of state of neutron-rich matter. A rich worldwide experimental program involving studies with rare isotopes, parity violating electron scattering, and astronomical observations is devoted to pinning down the isovector sector of nuclear models. Purpose: We assess the theoretical systematic and statistical uncertainties of neutron-skin thickness and relate them to the equation of state of nuclear matter, and in particular to nuclear symmetry energy parameters. Methods: We use the nuclear superfluid Density Functional Theory with several Skyrme energy density functionals and density dependent pairing. To evaluate statistical errors and their budget, we employ the statistical covariance technique. Results: We find that the errors on neutron skin increase with neutron excess. Statistical errors due to uncertain coupling constants of the density functional are found to be larger than systematic errors, the latter not exceeding 0.06 fm in most neutron-rich nuclei across the nuclear landscape. The single major source of uncertainty is the poorly determined slope L of the symmetry energy that parametrizes its density dependence. Conclusions: To provide essential constraints on the symmetry energy of the nuclear energy density functional, next-generation measurements of neutron skins are required to deliver precision better than 0.06 fm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nuclear three-body problem in the complex energy plane: Complex-Scaling-Slater method

The physics of open quantum systems is an interdisciplinary area of research. The nuclear "openness" manifests itself through the presence of the many-body continuum representing various decay, scattering, and reaction channels. As the radioactive nuclear beam experimentation extends the known nuclear landscape towards the particle drip lines, the coupling to the continuum space becomes exceedingly more important. Of particular interest are weakly bound and unbound nuclear states appearing around particle thresholds. Theories of such nuclei must take into account their open quantum nature. To describe open quantum systems, we introduce a Complex Scaling (CS) approach in the Slater basis. We benchmark it with the complex-energy Gamow Shell Model (GSM) by studying energies and wave functions of the bound and unbound states of the two-neutron halo nucleus 6He viewed as an $α$+ n + n cluster system. In the CS approach, we use the Slater basis, which exhibits the correct asymptotic behavior at large distances. To extract particle densities from the back-rotated CS solutions, we apply the Tikhonov regularization procedure, which minimizes the ultraviolet numerical noise. While standard applications of the inverse complex transformation to the complex-rotated solution provide unstable results, the stabilization method fully reproduces the GSM benchmark. We also propose a method to determine the smoothing parameter of the Tikhonov regularization. The combined suite of CS-Slater and GSM techniques has many attractive features when applied to nuclear problems involving weakly-bound and unbound states. While both methods can describe energies, total widths, and wave functions of nuclear states, the CS-Slater method, if it can be applied, can provide an additional information about partial energy widths associated with individual thresholds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reflection-asymmetric nuclear deformations within the Density Functional Theory

Within the nuclear density functional theory (DFT) we study the effect of reflection-asymmetric shapes on ground-state binding energies and binding energy differences. To this end, we developed the new DFT solver AxialHFB that uses an approximate second-order gradient to solve the Hartree-Fock-Bogoliubov equations of superconducting DFT with the quasi-local Skyrme energy density functionals. Illustrative calculations are carried out for even-even isotopes of radium and thorium.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The landscape of two-proton radioactivity

Ground-state two-proton (2p) radioactivity is a decay mode found in isotopes of elements with even atomic numbers located beyond the two-proton drip line. So far, this exotic process has been experimentally observed in a few light and medium-mass nuclides with Z less than or equal to 30. In this study, using state-of-the-art nuclear density functional theory, we globally analyze 2p radioactivity and for the first time identify 2p decay candidates in elements heavier than strontium. We predict a few cases where the competition between 2p emission and alpha decay may be observed. In nuclei above lead, the alpha decay mode is found to be dominating and no measurable candidates for the 2p radioactivity are expected.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Third minima in thorium and uranium isotopes in a self-consistent theory

Background: Deep third minima have been predicted in some non-self-consistent models to impact fission pathways of thorium and uranium isotopes. These predictions have guided the interpretation of resonances seen experimentally. On the other hand, self-consistent calculations consistently predict very shallow potential-energy surfaces in the third minimum region. Purpose: We investigate the interpretation of third-minimum configurations in terms of dimolecular states. We study the isentropic potential-energy surfaces of selected even-even thorium and uranium isotopes at several excitation energies. In order to understand the driving effects behind the presence of third minima, we study the interplay between pairing and shell effects. Methods: We use the finite-temperature superfluid nuclear density functional theory. We consider a traditional functional, SkM*, and a recent functional, UNEDF1, optimized for fission studies. Results: We predict very shallow or no third minima in the potential-energy surfaces of 232Th and 232U. In Th and U isotopes with N=136 and 138, the third minima are deeper. We show that the reflection-asymmetric configurations around the third minimum can be associated with dimolecular states involving the spherical doubly magic 132Sn and a lighter deformed Zr or Mo fragment. The potential-energy surfaces for 228,232Th and 232U at several excitation energies are presented. Conclusions: We show that the neutron shell effect that governs the existence of the dimolecular states around the third minimum is consistent with the spherical-to-deformed shape transition in the Zr and Mo isotopes around N=58. We demonstrate that the thermal reduction of pairing and enhancement of shell effects at small excitation energies help to develop deeper third minima. At large excitation energies, shell effects are washed out and third minima disappear altogether.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Alpha Decay in the Complex Energy Shell Model

Alpha emission from a nucleus is a fundamental decay process in which the alpha particle formed inside the nucleus tunnels out through the potential barrier. We describe alpha decay of $^{212}$Po and $^{104}$Te by means of the configuration interaction approach. To compute the preformation factor and penetrability, we use the complex-energy shell model with a separable T=1 interaction. The single-particle space is expanded in a Woods-Saxon basis that consists of bound and unbound resonant states. Special attention is paid to the treatment of the norm kernel appearing in the definition of the formation amplitude that guarantees the normalization of the channel function. Without explicitly considering the alpha-cluster component in the wave function of the parent nucleus, we reproduce the experimental alpha-decay width of $^{212}$Po and predict an upper limit of T_{1/2}=5.5x10^{-7} sec for the half-life of $^{104}$Te. The complex-energy shell model in a large valence configuration space is capable of providing a microscopic description of the alpha decay of heavy nuclei having two valence protons and two valence neutrons outside the doubly magic core. The inclusion of proton-neutron interaction between the valence nucleons is likely to shorten the predicted half-live of $^{104}$Te.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic normalization coefficients and continuum coupling in mirror nuclei

Background: An asymptotic normalization coefficient (ANC) characterizes the asymptotic form of a one-nucleon overlap integral required for description of nucleon-removal reactions. Purpose: We investigate the impact of the particle continuum on proton and neutron ANCs for mirror systems from $p$- and $sd$-shell regions. Method: We use the real-energy and complex-energy continuum shell model approaches. Results: We studied the general structure of the single-particle ANCs as a function of the binding energy and orbital angular momentum. We computed ANCs in mirror nuclei for different physical situations, including capture reactions to weakly-bound and unbound states. Conclusions: We demonstrated that the single-particle ANCs exhibit generic behavior that is different for charged and neutral particles. We verified the previously proposed relation [N.K. Timofeyuk, R.C. Johnson, and A.M. Mukhamedzhanov, Phys. Rev. Lett. 91, 232501 (2003); Phys. Rev. Lett. 97, 069904(E) (2006); N. K. Timofeyuk and P. Descouvemont, Phys. Rev. C 72, 064324 (2005)] between proton and neutron mirror ANCs. We find minor modifications if the spectroscopic strength is either localized in a single state or broadly distributed. For cases when several states couple strongly to the decay channel, these modifications may reach 30%.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coordinate-Space Hartree-Fock-Bogoliubov Solvers for Superfluid Fermi Systems in Large Boxes

The self-consistent Hartree-Fock-Bogoliubov problem in large boxes can be solved accurately in the coordinate space with the recently developed solvers HFB-AX (2D) and MADNESS-HFB (3D). This is essential for the description of superfluid Fermi systems with complicated topologies and significant spatial extend, such as fissioning nuclei, weakly-bound nuclei, nuclear matter in the neutron star rust, and ultracold Fermi atoms in elongated traps. The HFB-AX solver based on B-spline techniques uses a hybrid MPI and OpenMP programming model for parallel computation for distributed parallel computation, within a node multi-threaded LAPACK and BLAS libraries are used to further enable parallel calculations of large eigensystems. The MADNESS-HFB solver uses a novel multi-resolution analysis based adaptive pseudo-spectral techniques to enable fully parallel 3D calculations of very large systems. In this work we present benchmark results for HFB-AX and MADNESS-HFB on ultracold trapped fermions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energy density functional for nuclei and neutron stars

We aim to develop a nuclear energy density functional that can be simultaneously applied to finite nuclei and neutron stars. We use the self-consistent nuclear density functional theory (DFT) with Skyrme energy density functionals and covariance analysis to assess correlations between observables for finite nuclei and neutron stars. In a first step two energy functionals -- a high density energy functional giving reasonable neutron properties, and a low density functional fitted to nuclear properties -- are matched. In a second step, we optimize a new functional using exactly the same protocol as in earlier studies pertaining to nuclei but now including neutron star data. This allows direct comparisons of performance of the new functional relative to the standard one. The new functional TOV-min yields results for nuclear bulk properties (energy, r.m.s. radius, diffraction radius, surface thickness) that are of the same quality as those obtained with the established Skyrme functionals, including SV-min. When comparing SV-min and TOV-min, isoscalar nuclear matter indicators vary slightly while isovector properties are changed considerably. We discuss neutron skins, dipole polarizability, separation energies of the heaviest elements, and proton and neutron drip lines. We confirm a correlation between the neutron skin of $^{208}$Pb and the neutron star radius. We demonstrate that standard energy density functionals optimized to nuclear data do not carry information on the expected maximum neutron star mass, and that predictions can only be made within an extremely broad uncertainty band. For atomic nuclei, the new functional TOV-min performs at least as well as the standard nuclear functionals, but it also reproduces expected neutron star data within assumed error bands.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fission modes of mercury isotopes

Background: Recent experiments on beta-delayed fission in the mercury-lead region and the discovery of asym- metric fission in 180 Hg [1] have stimulated theoretical interest in the mechanism of fission in heavy nuclei. Purpose: We study fission modes and fusion valleys in 180 Hg and 198 Hg to reveal the role of shell effects in pre-scission region and explain the experimentally observed fragment mass asymmetry and its variation with A. Methods: We use the self-consistent nuclear density functional theory employing Skyrme and Gogny energy density functionals. Results: The potential energy surfaces in multi-dimensional space of collective coordinates, including elongation, triaxiality, reflection-asymmetry, and necking, are calculated for 180 Hg and 198 Hg. The asymmetric fission valleys - well separated from fusion valleys associated with nearly spherical fragments - are found in in both cases. The density distributions at scission configurations are studied and related to the experimentally observed mass splits. Conclusions: The energy density functionals SkM\ast and D1S give a very consistent description of the fission process in 180 Hg and 198 Hg. We predict a transition from asymmetric fission in 180 Hg towards more symmetric distribution of fission fragments in 198 Hg. For 180 Hg, both models yield 100 Ru/80 Kr as the most probable split. For 198 Hg, the most likely split is 108 Ru/90 Kr in HFB-D1S and 110 Ru/88 Kr in HFB-SkM\ast.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Information content of the low-energy electric dipole strength: correlation analysis

We study the information content carried by the electric dipole strength with respect to isovector and isoscalar indicators characterizing bulk nuclear matter and finite nuclei. To separate isoscalar and isovector modes, and low-energy strength and giant resonances, we analyze the E1 strength as a function of excitation energy $E$ and momentum transfer $q$. We use the self-consistent nuclear density functional theory with Skyrme energy density functionals, augmented by the random phase, to compute the E1 strength, and covariance analysis to assess correlations between observables. Calculations are performed for spherical, doubly-magic nuclei $^{208}$Pb and $^{132}$Sn. We demonstrate that E1 transition densities in the low-energy region below the giant dipole resonance exhibit appreciable state dependence and multi-nodal structures, which are fingerprints of weak collectivity. The correlation between the accumulated low-energy strength and symmetry energy is weak, and dramatically depends on the energy cutoff assumed. On the other hand, a strong correlation is predicted between isovector indicators and the accumulated isovector strength at $E$ around 20 MeV and momentum transfer $q\sim 0.65$ fm$^{-1}$. Momentum- and coordinate-space pattern of the low-energy dipole modes indicate a strong fragmentation into individual particle-hole excitations. The global measure of low-energy dipole strength poorly correlates with the nuclear symmetry energy and other isovector characteristics. Consequently, our results do not support the suggestion that there exists a collective "pygmy dipole resonance," which is a strong indicator of nuclear isovector properties.

0 citations0 reviewsNo code linkedOpen access
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%.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Living on the edge of stability, the limits of the nuclear landscape

A first-principles description of nuclear systems along the drip lines presents a substantial theoretical and computational challenge. In this paper, we discuss the nuclear theory roadmap, some of the key theoretical approaches, and present selected results with a focus on long isotopic chains. An important conclusion, which consistently emerges from these theoretical analyses, is that three-nucleon forces are crucial for both global nuclear properties and detailed nuclear structure, and that many-body correlations due to the coupling to the particle continuum are essential as one approaches particle drip lines. In the quest for a comprehensive nuclear theory, high performance computing plays a key role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nuclear energy density optimization: Large deformations

A new Skyrme-like energy density suitable for studies of strongly elongated nuclei has been determined in the framework of the Hartree-Fock-Bogoliubov theory using the recently developed model-based, derivative-free optimization algorithm POUNDerS. A sensitivity analysis at the optimal solution has revealed the importance of states at large deformations in driving the parameterization of the functional. The good agreement with experimental data on masses and separation energies, achieved with the previous parameterization UNEDF0, is largely preserved. In addition, the new energy density UNEDF1 gives a much improved description of the fission barriers in 240Pu and neighboring nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the origin of nuclear clustering

Clustering is one of the most complex phenomena known to the structure of atomic nuclei. A comprehensive description of this ubiquitous phenomenon goes beyond standard shell model and cluster model frameworks. We argue that clustering is a consequence of an openness of the nuclear many-body system. To illustrate this point, we study the near-threshold behavior of exceptional points in ^16 Ne and ^24 S.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spontaneous fission modes and lifetimes of super-heavy elements in the nuclear density functional theory

Lifetimes of super-heavy (SH) nuclei are primarily governed by alpha decay and spontaneous fission (SF). Here we study the competing decay modes of even-even SH isotopes with 108 <= Z <= 126 and 148 <= N <= 188 using the state-of-the-art self-consistent nuclear density functional theory framework capable of describing the competition between nuclear attraction and electrostatic repulsion. The collective mass tensor of the fissioning superfluid nucleus is computed by means of the cranking approximation to the adiabatic time-dependent Hartree-Fock-Bogoliubov approach. Along the path to fission, our calculations allow for the simultaneous breaking of axial and space inversion symmetries; this may result in lowering SF lifetimes by more than seven orders of magnitude in some cases. We predict two competing SF modes: reflection-symmetric and reflection-asymmetric.The shortest-lived SH isotopes decay by SF; they are expected to lie in a narrow corridor formed by $^{280}$Hs, $^{284}$Fl, and $^{284}_{118}$Uuo that separates the regions of SH nuclei synthesized in "cold fusion" and "hot fusion" reactions. The region of long-lived SH nuclei is expected to be centered on $^{294}$Ds with a total half-life of ?1.5 days.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Toward understanding the microscopic origin of nuclear clustering

Open Quantum System (OQS) description of a many-body system involves interaction of Shell Model (SM) states through the particle continuum. In realistic nuclear applications, this interaction may lead to collective phenomena in the ensemble of SM states. We claim that the nuclear clustering is an emergent, near-threshold phenomenon, which cannot be elucidated within the Closed Quantum System (CQS) framework. We approach this problem by investigating the near-threshold behavior of Exceptional Points (EPs) in the realistic Continuum Shell Model (CSM). The consequences for the alpha-clustering phenomenon are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

UNEDF: Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

The demands of cutting-edge science are driving the need for larger and faster computing resources. With the rapidly growing scale of computing systems and the prospect of technologically disruptive architectures to meet these needs, scientists face the challenge of effectively using complex computational resources to advance scientific discovery. Multidisciplinary collaborating networks of researchers with diverse scientific backgrounds are needed to address these complex challenges. The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper describes UNEDF and identifies attributes that classify it as a successful computational collaboration. We illustrate significant milestones accomplished by UNEDF through integrative solutions using the most reliable theoretical approaches, most advanced algorithms, and leadership-class computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charge radii and neutron correlations in helium halo nuclei

Within the complex-energy configuration interaction framework, we study correlations of valence neutrons to explain the behavior of charge radii in the neutron halo nuclei $^{6,8}$He. We find that the experimentally observed decrease of the charge radius between $^6$He and $^8$He is caused by a subtle interplay between three effects: dineutron correlations, a spin-orbit contribution to the charge radius, and a core swelling effect. We demonstrate that two-neutron angular correlations in the $2^+_1$ resonance of $^6$He differ markedly from the ground-state correlations in $^{6,8}$He.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computing Heavy Elements

Reliable calculations of the structure of heavy elements are crucial to address fundamental science questions such as the origin of the elements in the universe. Applications relevant for energy production, medicine, or national security also rely on theoretical predictions of basic properties of atomic nuclei. Heavy elements are best described within the nuclear density functional theory (DFT) and its various extensions. While relatively mature, DFT has never been implemented in its full power, as it relies on a very large number (~ 10^9-10^12) of expensive calculations (~ day). The advent of leadership-class computers, as well as dedicated large-scale collaborative efforts such as the SciDAC 2 UNEDF project, have dramatically changed the field. This article gives an overview of the various computational challenges related to the nuclear DFT, as well as some of the recent achievements.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Isospin mixing and the continuum coupling in weakly bound nuclei

The isospin breaking effects due to the Coulomb interaction in weakly-bound nuclei are studied using the Gamow Shell Model, a complex-energy configuration interaction approach which simultaneously takes into account many-body correlations between valence nucleons and continuum effects. We investigate the near-threshold behavior of one-nucleon spectroscopic factors and the structure of wave functions along an isomultiplet. Illustrative calculations are carried out for the T=1 isobaric triplet. By using a shell-model Hamiltonian consisting of an isoscalar nuclear interaction and the Coulomb term, we demonstrate that for weakly bound or unbound systems the structure of isobaric analog states varies within the isotriplet and impacts the energy dependence of spectroscopic factors. We discuss the partial dynamical isospin symmetry present in isospin-stretched systems, in spite of the Coulomb interaction that gives rise to large mirror symmetry breaking effects.

0 citations0 reviewsNo code linkedOpen access
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)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multipole strength function of deformed superfluid nuclei made easy

We present an efficient method for calculating strength functions using the finite amplitude method (FAM) for deformed superfluid heavy nuclei within the framework of the nuclear density functional theory. We demonstrate that FAM reproduces strength functions obtained with the fully self-consistent quasi-particle random-phase approximation (QRPA) at a fraction of computational cost. As a demonstration, we compute the isoscalar and isovector monopole strength for strongly deformed configurations in $^{240}$Pu by considering huge quasi-particle QRPA spaces. Our approach to FAM, based on Broyden's iterative procedure, opens the possibility for large-scale calculations of strength distributions in well-bound and weakly bound nuclei across the nuclear landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quasi-particle continuum and resonances in the Hartree-Fock-Bogoliubov theory

The quasi-particle energy spectrum of the Hartree-Fock-Bogoliubov (HFB) equations contains discrete bound states, resonances, and non-resonant continuum states. We study the structure of the unbound quasi-particle spectrum of weakly bound nuclei within several methods that do not rely on imposing scattering or outgoing boundary conditions. Various approximations are examined to estimate resonance widths. It is shown that the stabilization method works well for all HFB resonances except for very narrow ones. The Thomas-Fermi approximation to the non-resonant continuum has been shown to be very effective, especially for coordinate-space HFB calculations in large boxes that involve huge amounts of discretized quasi-particle continuum states.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

UNEDF: Advanced Scientific Computing Transforms the Low-Energy Nuclear Many-Body Problem

The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper illustrates significant milestones accomplished by UNEDF through integration of the theoretical approaches, advanced numerical algorithms, and leadership class computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Augmented Lagrangian Method for Constrained Nuclear Density Functional Theory

The augmented Lagrangiam method (ALM), widely used in quantum chemistry constrained optimization problems, is applied in the context of the nuclear Density Functional Theory (DFT) in the self-consistent constrained Skyrme Hartree-Fock-Bogoliubov (CHFB) variant. The ALM allows precise calculations of multidimensional energy surfaces in the space of collective coordinates that are needed to, e.g., determine fission pathways and saddle points; it improves accuracy of computed derivatives with respect to collective variables that are used to determine collective inertia; and is well adapted to supercomputer applications.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Competition between Normal Superfluidity and Larkin-Ovchinnikov Phases of Polarized Fermi Gases in Elongated Traps

By applying the recently proposed antisymmetric superfluid local density approximation (ASLDA) to strongly interacting polarized atomic gases at unitarity in very elongated traps, we find families of Larkin-Ovchinnikov (LO) type of solutions with prominent transversal oscillation of pairing potential. These LO states coexist with a superfluid state having a smooth pairing potential. We suggest that the LO phase could be accessible experimentally by increasing adiabatically the trap aspect ratio. We show that the local asymmetry effects contained in ASLDA do not support a deformed superfluid core predicted by previous Bogoliubov-de Gennes treatments.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Natural Units For Nuclear Energy Density Functional Theory

Naive dimensional analysis based on chiral effective theory, when adapted to nuclear energy density functionals, prescribes natural units and a hierarchy of contributions that could be used to constrain fits of generalized functionals. By applying these units, a large sample of Skyrme parametrizations is examined for naturalness, which is signaled by dimensionless coupling constants of order one. The bulk of the parameters are found to be natural, with an underlying scale consistent with other determinations. Significant deviations from unity are associated with deficiencies in the corresponding terms of particular functionals or with an incomplete optimization procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nuclear Energy Density Optimization

We carry out state-of-the-art optimization of a nuclear energy density of Skyrme type in the framework of the Hartree-Fock-Bogoliubov (HFB) theory. The particle-hole and particle-particle channels are optimized simultaneously, and the experimental data set includes both spherical and deformed nuclei. The new model-based, derivative-free optimization algorithm used in this work has been found to be significantly better than standard optimization methods in terms of reliability, speed, accuracy, and precision. The resulting parameter set UNEDFpre results in good agreement with experimental masses, radii, and deformations and seems to be free of finite-size instabilities. An estimate of the reliability of the obtained parameterization is given, based on standard statistical methods. We discuss new physics insights offered by the advanced covariance analysis.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Open problems in theory of nuclear open quantum systems

Is there a connection between the branch point singularity at the particle emission threshold and the appearance of cluster states which reveal the structure of a corresponding reaction channel? Which nuclear states are most impacted by the coupling to the scattering continuum? What should be the most important steps in developing the theory that will truly unify nuclear structure and nuclear reactions? The common denominator of these questions is the continuum shell-model approach to bound and unbound nuclear states, nuclear decays, and reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pairing Reentrance Phenomenon in Heated Rotating Nuclei in the Shell Model Monte Carlo Approach

Rotational motion of heated 72-Ge is studied within the microscopic Shell Model Monte Carlo approach. We investigate the the angular momentum alignment and nuclear pairing correlations associated with J-pi Cooper pairs as a function of the rotational frequency and temperature. The reentrance of pairing correlations with temperature is predicted at high rotational frequencies. It manifests itself through the anomalous behavior of specific heat and level density.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Surface Symmetry Energy of Nuclear Energy Density Functionals

We study the bulk deformation properties of the Skyrme nuclear energy density functionals. Following simple arguments based on the leptodermous expansion and liquid drop model, we apply the nuclear density functional theory to assess the role of the surface symmetry energy in nuclei. To this end, we validate the commonly used functional parametrizations against the data on excitation energies of superdeformed band-heads in Hg and Pb isotopes, and fission isomers in actinide nuclei. After subtracting shell effects, the results of our self-consistent calculations are consistent with macroscopic arguments and indicate that experimental data on strongly deformed configurations in neutron-rich nuclei are essential for optimizing future nuclear energy density functionals. The resulting survey provides a useful benchmark for further theoretical improvements. Unlike in nuclei close to the stability valley, whose macroscopic deformability hangs on the balance of surface and Coulomb terms, the deformability of neutron-rich nuclei strongly depends on the surface-symmetry energy; hence, its proper determination is crucial for the stability of deformed phases of the neutron- rich matter and description of fission rates for r-process nucleosynthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The information content of a new observable: the case of the nuclear neutron skin

We address two questions pertaining to the uniqueness and usefulness of a new observable: (i) Considering the current theoretical knowledge, what novel information does new measurement bring in? (ii) How can new data reduce uncertainties of current theoretical models? We illustrate these points by studying the radius of the neutron distribution of a heavy nucleus, a quantity related to the equation of state for neutron matter that determines properties of nuclei and neutron stars. By systematically varying parameters of two theoretical models and studying the resulting confidence ellipsoid, we quantify the relationships between the neutron skin and various properties of finite nuclei and infinite nuclear matter. Using the covariance analysis, we identify observables and pseudo-observables that correlate, and do not correlate, with the neutron skin. By adding the information on the neutron radius to the pool of observables determining the energy functional, we show how precise experimental determination of the neutron radius in $^{208}$Pb would reduce theoretical uncertainties on the neutron matter equation of state.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Odd-even mass differences from self-consistent mean-field theory

We survey odd-even nuclear binding energy staggering using density functional theory with several treatments of the pairing interaction including the BCS, Hartree-Fock-Bogoliubov, and the Hartree-Fock-Bogoliubov with the Lipkin-Nogami approximation. We calculate the second difference of binding energies and compare with 443 measured neutron energy differences in isotope chains and 418 measured proton energy differences in isotone chains. The particle-hole part of the energy functional is taken as the SLy4 Skyrme parametrization and the pairing part of the functional is based on a contact interaction with possible density dependence. An important feature of the data, reproduced by the theory, is the sharp gap quenching at magic numbers. With the strength of the interaction as a free parameter, the theory can reproduce the data to an rms accuracy of about 0.25 MeV. This is slightly better than a single-parameter phenomenological description but slightly poorer than the usual two-parameter phenomenological form C/A^alpha . The following conclusions can be made about the performance of common parametrization of the pairing interaction: (i) there is a weak preference for a surface-peaked neutron-neutron pairing, which might be attributable to many-body effects; (ii) a larger strength is required in the proton pairing channel than in the neutron pairing channel; (iii) pairing strengths adjusted to the well-known spherical isotope chains are too weak to give a good overall fit to the mass differences.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Gamow Shell-Model Description of Weakly Bound and Unbound Nuclear States

Recently, the shell model in the complex k-plane (the so-called Gamow Shell Model) has been formulated using a complex Berggren ensemble representing bound single-particle states, single-particle resonances, and non-resonant continuum states. In this framework, we shall discuss binding energies and energy spectra of neutron-rich helium and lithium isotopes. The single-particle basis used is that of the Hartree-Fock potential generated self-consistently by the finite-range residual interaction.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access