Showing 101–150 of 255 results for author: Nazarewicz, W

Nuclear Energy Density Optimization: UNEDF2

Authors: M. Kortelainen, J. McDonnell, W. Nazarewicz, E. Olsen, P. -G. Reinhard, J. Sarich, N. Schunck, S. M. Wild, D. Davesne, J. Erler, A. Pastore

Abstract: 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.… ▽ More 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. △ Less

Submitted 30 October, 2014; originally announced October 2014.

Comments: 6 pages, 2 figures. Talk given at the ARIS 2014 conference, Tokyo, Japan, June 1-6, 2014

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

Authors: K. Fossez, N. Michel, W. Nazarewicz, M. Płoszajczak, Y. Jaganathen

Abstract: 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… ▽ More 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. △ Less

Submitted 24 October, 2014; originally announced October 2014.

Comments: 11 pages, 13 figures

Journal ref: Phys. Rev. A 91, 012503 (2015)

Pairing-induced speedup of nuclear spontaneous fission

Authors: Jhilam Sadhukhan, J. Dobaczewski, W. Nazarewicz, J. A. Sheikh, A. Baran

Abstract: 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 fissi… ▽ More 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. △ Less

Submitted 6 October, 2014; originally announced October 2014.

Comments: Submitted to PRC

Journal ref: Phys. Rev. C 90, 061304(R) (2014)

Adaptive Multi-resolution 3D Hartree-Fock-Bogoliubov Solver for Nuclear Structure

Authors: Junchen Pei, George Fann, Robert Harrison, Witold Nazarewicz, Yue Shi, Scott Thornton

Abstract: Complex many-body systems, such as triaxial and reflection-asymmetric nuclei, weakly-bound halo states, cluster configurations, nuclear fragments produced in heavy-ion fusion reactions, cold Fermi gases, and pasta phases in neutron star crust, they are all characterized by large sizes and complex topologies, in which many geometrical symmetries characteristic of ground-state configurations are bro… ▽ More Complex many-body systems, such as triaxial and reflection-asymmetric nuclei, weakly-bound halo states, cluster configurations, nuclear fragments produced in heavy-ion fusion reactions, cold Fermi gases, and pasta phases in neutron star crust, they are all characterized by large sizes and complex topologies, in which many geometrical symmetries characteristic of ground-state configurations are broken. A tool of choice to study such complex forms of matter is an adaptive multi-resolution wavelet analysis. This method has generated much excitement since it provides a common framework linking many diversified methodologies across different fields, including signal processing, data compression, harmonic analysis and operator theory, fractals, and quantum field theory. To describe complex superfluid many-fermion systems, we introduce an adaptive pseudo-spectral method for solving self-consistent equations of nuclear density functional theory in three dimensions, without symmetry restrictions. The new adaptive multi-resolution Hartree-Fock-Bogoliubov (HFB) solver {\madnesshfb} is benchmarked against a two-dimensional coordinate-space solver {\hfbax} based on B-spline technique and three-dimensional solver {\hfodd} based on the harmonic oscillator basis expansion. Several examples are considered, including self-consistent HFB problem for spin-polarized trapped cold fermions and Skyrme-Hartree-Fock (+BCS) problem for triaxial deformed nuclei. The new {\madnesshfb} framework has many attractive features when applied to nuclear and atomic problems involving many-particle superfluid systems. Of particular interest are weakly-bound nuclear configurations close to particle drip lines, strongly elongated and dinuclear configurations such as those present in fission and heavy ion fusion, and exotic pasta phases that appear in the neutron star crust. △ Less

Submitted 14 July, 2014; originally announced July 2014.

Comments: 9 pages, 5 figures, submitted to PRC

Excitation energy dependence of fission in the mercury region

Authors: J. D. McDonnell, W. Nazarewicz, J. A. Sheikh, A. Staszczak, M. Warda

Abstract: 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 e… ▽ More 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. △ Less

Submitted 26 June, 2014; originally announced June 2014.

Comments: 6 pages, 6 figures

Journal ref: Phys. Rev. C 90, 021302 (2014)

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

Authors: J. A. Sheikh, N. Hinohara, J. Dobaczewski, T. Nakatsukasa, W. Nazarewicz, K. Sato

Abstract: 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… ▽ More 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. △ Less

Submitted 20 May, 2014; v1 submitted 10 March, 2014; originally announced March 2014.

Comments: 13 pages, 10 figures, submitted to Phys. Rev. C

Journal ref: Phys.Rev.C89, 054317 (2014)

Error Estimates of Theoretical Models: a Guide

Authors: J. Dobaczewski, W. Nazarewicz, P. -G. Reinhard

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 19 February, 2014; originally announced February 2014.

Comments: 24 pages, 9 figures, submitted to Journal of Physics G

Report number: INT-PUB-14-004

Journal ref: J. Phys. G: Nucl. Part. Phys. 41 (2014) 074001

Structure of Superheavy Nuclei Along Element 115 Decay Chains

Authors: Yue Shi, D. E. Ward, B. G. Carlsson, J. Dobaczewski, W. Nazarewicz, I. Ragnarsson, D. Rudolph

Abstract: 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 approac… ▽ More 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. △ Less

Submitted 26 June, 2014; v1 submitted 30 January, 2014; originally announced January 2014.

Comments: 9 pages, 8 figures, 6 table; accepted for publication in Physical Review C

Journal ref: Phys. Rev. C 90, 014308 (2014)

Nuclear Theory and Science of the Facility for Rare Isotope Beams

Authors: A. B Balantekin, J. Carlson, D. J. Dean, G. M. Fuller, R. J. Furnstahl, M. Hjorth-Jensen, R. V. F. Janssens, Bao-An Li, W. Nazarewicz, F. M. Nunes, W. E. Ormand, S. Reddy, B. M. Sherrill

Abstract: 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 scien… ▽ More 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. △ Less

Submitted 24 January, 2014; originally announced January 2014.

Comments: 20 pages, 7 figures

Journal ref: Mod. Phys. Lett. A, Vol. 29, No. 11 (2014) 1430010

Nuclear energy density optimization: Shell structure

Authors: M. Kortelainen, J. McDonnell, W. Nazarewicz, E. Olsen, P. -G. Reinhard, J. Sarich, N. Schunck, S. M. Wild, D. Davesne, J. Erler, A. Pastore

Abstract: 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 Har… ▽ More 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. △ Less

Submitted 28 April, 2014; v1 submitted 5 December, 2013; originally announced December 2013.

Comments: 18 pages, 13 figures, 12 tables; resubmitted for publication to Phys. Rev. C after second review by referee

Journal ref: Phys. Rev. C 89, 054314 (2014)

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

Authors: A. T. Kruppa, G. Papadimitriou, W. Nazarewicz, N. Michel

Abstract: 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 exceedingl… ▽ More 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. △ Less

Submitted 29 October, 2013; v1 submitted 28 October, 2013; originally announced October 2013.

Comments: 15 pages, 16 figures

Near-threshold correlations of neutrons

Authors: J. Okołowicz, M. Płoszajczak, W. Nazarewicz

Abstract: The appearance of charged-particle clustering in near-threshold configuration is a phenomenon that can be explained in the Open Quantum System description of the atomic nucleus. In this work we apply the realistic Shell Model Embedded in the Continuum to elucidate the emergence of neutron correlations in near-threshold many-body states coupled to l=1,2 neutron decay channels. Spectral consequences… ▽ More The appearance of charged-particle clustering in near-threshold configuration is a phenomenon that can be explained in the Open Quantum System description of the atomic nucleus. In this work we apply the realistic Shell Model Embedded in the Continuum to elucidate the emergence of neutron correlations in near-threshold many-body states coupled to l=1,2 neutron decay channels. Spectral consequences of such continuum coupling are briefly discussed together with the emergence of complex multi-neutron correlations. △ Less

Submitted 24 October, 2013; originally announced October 2013.

Comments: Invited talk at XXXIII Mazurian Lakes Conference on Physics

Journal ref: Acta Physica Polonica B 45, 331 (2014)

Isospin mixing within the symmetry restored density functional theory and beyond

Authors: W. Satula, J. Dobaczewski, M. Konieczka, W. Nazarewicz

Abstract: 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 functio… ▽ More 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. △ Less

Submitted 15 October, 2013; originally announced October 2013.

Comments: Invited talk presented by WS at the XXXIII Mazurian Lakes Conference on Physics, Piaski, Poland, 1-7 September 2013

Journal ref: Acta Phys. Pol. B45, 167 (2014)

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

Authors: Jhilam Sadhukhan, K. Mazurek, A. Baran, J. Dobaczewski, W. Nazarewicz, J. A. Sheikh

Abstract: 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 inter… ▽ More 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. △ Less

Submitted 8 October, 2013; originally announced October 2013.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. C 88, 064314 (2013)

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

Authors: A. V. Afanasjev, Yue Shi, W. Nazarewicz

Abstract: 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… ▽ More 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. △ Less

Submitted 12 September, 2013; originally announced September 2013.

Comments: 6 pages, 4 figures

Journal ref: Physical Review C 86, 031304 (2012)

Information content of the weak-charge form factor

Authors: P. -G. Reinhard, J. Piekarewicz, W. Nazarewicz, B. K. Agrawal, N. Paar, X. Rocca-Maza

Abstract: 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 nu… ▽ More 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. △ Less

Submitted 7 August, 2013; originally announced August 2013.

Comments: 10 pages, 4 figures

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

Authors: Yue Shi, C. L. Zhang, J. Dobaczewski, W. Nazarewicz

Abstract: \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 cor… ▽ More \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. △ Less

Submitted 11 September, 2013; v1 submitted 1 August, 2013; originally announced August 2013.

Comments: 7 pages, 4 figures

Journal ref: Physical Review C 88, 034311 (2013)

Symmetry energy in nuclear density functional theory

Authors: W. Nazarewicz, P. -G. Reinhard, W. Satula, D. Vretenar

Abstract: 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… ▽ More 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. △ Less

Submitted 22 July, 2013; originally announced July 2013.

Comments: 14 pages, 8 figures, submitted to the Special EPJA Issue on "Symmetry Energy"

Journal ref: Eur. Phys. J. A50, 20 (2014)

Neutron skin uncertainties of Skyrme energy density functionals

Authors: M. Kortelainen, J. Erler, W. Nazarewicz, N. Birge, Y. Gao, E. Olsen

Abstract: 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 violatin… ▽ More 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. △ Less

Submitted 16 July, 2013; originally announced July 2013.

Comments: 5 pages, 4 figures

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

Authors: W. Satula, J. Dobaczewski, M. Konieczka, W. Nazarewicz

Abstract: 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… ▽ More 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. △ Less

Submitted 5 July, 2013; originally announced July 2013.

Comments: Invited talk presented by WS at the EURISOL User Group Topical Meeting 2013, July 1-3, Krakow, Poland

Low energy collective modes of deformed superfluid nuclei within the finite amplitude method

Authors: Nobuo Hinohara, Markus Kortelainen, Witold Nazarewicz

Abstract: Background: The major challenge for nuclear theory is to describe and predict global properties and collective modes of atomic nuclei. Of particular interest is the response of the nucleus to a time-dependent external field that impacts the low-energy multipole and beta-decay strength. Purpose: We propose a method to compute low-lying collective modes in deformed nuclei within the finite amplitu… ▽ More Background: The major challenge for nuclear theory is to describe and predict global properties and collective modes of atomic nuclei. Of particular interest is the response of the nucleus to a time-dependent external field that impacts the low-energy multipole and beta-decay strength. Purpose: We propose a method to compute low-lying collective modes in deformed nuclei within the finite amplitude method (FAM) based on the quasiparticle random-phase approximation (QRPA). By using the analytic property of the response function, we find the QRPA amplitudes by computing the residua of the FAM amplitudes by means of a contour integration around the QRPA poles in a complex frequency plane. Methods: We use the superfluid nuclear density functional theory with Skyrme energy density functionals, FAM-QRPA approach, and the conventional matrix formulation of the QRPA (MQRPA). Results: We demonstrate that the complex-energy FAM-QRPA method reproduces low-lying collective states obtained within the conventional matrix formulation of the QRPA theory. Illustrative calculations are performed for the isoscalar monopole strength in deformed 24Mg and for low-lying K = 0 quadrupole vibrational modes of deformed Yb and Er isotopes. Conclusions: The proposed FAM-QRPA approach allows one to efficiently calculate low-lying collective modes in spherical and deformed nuclei throughout the entire nuclear landscape, including shape-vibrational excitations, pairing vibrational modes, and beta-decay rates. △ Less

Submitted 15 April, 2013; originally announced April 2013.

Comments: 9 pages, 2 figures, submitted to Phys. Rev. C

Journal ref: Phys. Rev. C 87, 064309 (2013)

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

Authors: Scott Bogner, Aurel Bulgac, Joseph A. Carlson, Jonathan Engel, George Fann, Richard J. Furnstahl, Stefano Gandolfi, Gaute Hagen, Mihai Horoi, Calvin W. Johnson, Markus Kortelainen, Ewing Lusk, Pieter Maris, Hai Ah Nam, Petr Navratil, Witold Nazarewicz, Esmond G. Ng, Gustavo P. A. Nobre, Erich Ormand, Thomas Papenbrock, Junchen Pei, Steven C. Pieper, Sofia Quaglioni, Kenneth J. Roche, Jason Sarich , et al. (6 additional authors not shown)

Abstract: The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists can lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science… ▽ More The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists can lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science results to illustrate this interplay. △ Less

Submitted 12 April, 2013; originally announced April 2013.

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

Authors: A. Ekström, G. Baardsen, C. Forssén, G. Hagen, M. Hjorth-Jensen, G. R. Jansen, R. Machleidt, W. Nazarewicz, T. Papenbrock, J. Sarich, S. M. Wild

Abstract: 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 us… ▽ More 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. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 110, 192502 (2013)

Bound states of dipolar molecules studied with the Berggren expansion method

Authors: K. Fossez, N. Michel, W. Nazarewicz, M. Płoszajczak

Abstract: 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-b… ▽ More 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. △ Less

Submitted 8 March, 2013; originally announced March 2013.

Comments: 12 pages, 10 figures

Journal ref: Phys. Rev. A 87, 042515 (2013)

The landscape of two-proton radioactivity

Authors: E. Olsen, M. Pfützner, N. Birge, M. Brown, W. Nazarewicz, A. Perhac

Abstract: 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 a… ▽ More 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. △ Less

Submitted 16 September, 2013; v1 submitted 5 March, 2013; originally announced March 2013.

Comments: 1 page, 2 encapsulated-postscript figures; this work has been accepted for publication in Physical Review Letters

Reflection-asymmetric nuclear deformations within the Density Functional Theory

Authors: E. Olsen, J. Erler, W. Nazarewicz, M. Stoitsov

Abstract: 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. Illus… ▽ More 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. △ Less

Submitted 4 March, 2013; originally announced March 2013.

Comments: 7 pages, 5 encapsulated-postscript figures

Journal ref: Journal of Physics: Conference Series 402 (2012) 012034

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

Authors: J. D. McDonnell, W. Nazarewicz, J. A. Sheikh

Abstract: 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 interp… ▽ More 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. △ Less

Submitted 23 May, 2013; v1 submitted 5 February, 2013; originally announced February 2013.

Comments: 9 pages, 10 figures; published in Physical Review C

Journal ref: Phys. Rev. C 87, 054327 (2013)

Microscopic Description of Nuclear Fission: Fission Barrier Heights of Even-Even Actinides

Authors: J. McDonnell, N. Schunck, W. Nazarewicz

Abstract: We evaluate the performance of modern nuclear energy density functionals for predicting inner and outer fission barrier heights and energies of fission isomers of even-even actinides. For isomer energies and outer barrier heights, we find that the self-consistent theory at the HFB level is capable of providing quantitative agreement with empirical data. We evaluate the performance of modern nuclear energy density functionals for predicting inner and outer fission barrier heights and energies of fission isomers of even-even actinides. For isomer energies and outer barrier heights, we find that the self-consistent theory at the HFB level is capable of providing quantitative agreement with empirical data. △ Less

Submitted 31 January, 2013; originally announced January 2013.

Comments: 8 pages, 6 figures, 1 table; Proceedings of the 5th International Conference on "Fission and properties of neutron-rich nuclei" (ICFN5), Sanibel Island, Nov. 4-10, 2012

Journal ref: Fission and Properties of Neutron-Rich Nuclei, World Scientific (2013); p. 597

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

Authors: C. Forssen, G. Hagen, M. Hjorth-Jensen, W. Nazarewicz, J. Rotureau

Abstract: 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-nucle… ▽ More 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. △ Less

Submitted 27 December, 2012; originally announced December 2012.

Comments: Contribution to proceedings of Nobel Symposium 152: Physics with radioactive beams, June 2012, Gothenburg, Sweden

Journal ref: Physica Scripta T152, 014022 (2013)

Energy density functional for nuclei and neutron stars

Authors: J. Erler, C. J. Horowitz, W. Nazarewicz, M. Rafalski, P. -G. Reinhard

Abstract: 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 f… ▽ More 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. △ Less

Submitted 27 November, 2012; originally announced November 2012.

Comments: 11 pages, 8 figures

Journal ref: Phys. Rev. C 87, 044320 (2013)

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

Authors: P. -G. Reinhard, W. Nazarewicz

Abstract: 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 densi… ▽ More 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. △ Less

Submitted 7 November, 2012; originally announced November 2012.

Comments: 9 paghes, 13 figures

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

Authors: X. Roca-Maza, B. K. Agrawal, G. Colo', W. Nazarewicz, N. Paar, J. Piekarewicz, P. -G. Reinhard, D. Vretenar

Abstract: 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 dependenc… ▽ More 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). △ Less

Submitted 31 October, 2012; originally announced October 2012.

Comments: Proceedings of NSD12, Opatija, Croatia, 9-13 July 2012

Journal ref: AIP Conf. Proc. 1491, 204 (2012)

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

Authors: W. Satula, J. Dobaczewski, W. Nazarewicz, T. R. Werner

Abstract: 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-d… ▽ More 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%. △ Less

Submitted 25 October, 2012; v1 submitted 3 October, 2012; originally announced October 2012.

Comments: 15 pages; 14 figures; 5 tables; submitted to Phys. Rev. C

Journal ref: Phys. Rev. C86, 054316 (2012)

Alpha Decay in the Complex Energy Shell Model

Authors: R. Id Betan, W. Nazarewicz

Abstract: 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-… ▽ More 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. △ Less

Submitted 7 August, 2012; originally announced August 2012.

Comments: 14 two-column pages, 13 figures, 9 tables

Journal ref: Phys. Rev. C 86, 034338 (2012)

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

Authors: A. Staszczak, A. Baran, W. Nazarewicz

Abstract: 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 coll… ▽ More 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. △ Less

Submitted 6 August, 2012; originally announced August 2012.

Comments: 6 pages, 4 figures

Toward understanding the microscopic origin of nuclear clustering

Authors: J. Okołowicz, W. Nazarewicz, M. Płoszajczak

Abstract: 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)… ▽ More 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. △ Less

Submitted 26 July, 2012; originally announced July 2012.

Comments: Special Issue of Fortschritte der Physik - Progress of Physics: Physics of Non-Hermitian Systems (10 pages, 10 figures)

Journal ref: Fortschr. Phys. 61, 66 (2013)

Hartree-Fock-Bogoliubov solution of the pairing Hamiltonian in finite nuclei

Authors: J. Dobaczewski, W. Nazarewicz

Abstract: We present an overview of the Hartree-Fock-Bogoliubov (HFB) theory of nucleonic superfluidity for finite nuclei. After introducing basic concepts related to pairing correlations, we show how the correlated pairs are incorporated into the HFB wave function. Thereafter, we present derivation and structure of the HFB equations within the superfluid nuclear density functional formalism and discuss sev… ▽ More We present an overview of the Hartree-Fock-Bogoliubov (HFB) theory of nucleonic superfluidity for finite nuclei. After introducing basic concepts related to pairing correlations, we show how the correlated pairs are incorporated into the HFB wave function. Thereafter, we present derivation and structure of the HFB equations within the superfluid nuclear density functional formalism and discuss several aspects of the theory, including the unitarity of the Bogoliubov transformation in truncated single-particle and quasiparticle spaces, form of the pairing functional, structure of the HFB continuum, regularization and renormalization of pairing fields, and treatment of pairing in systems with odd particle numbers. △ Less

Submitted 1 July, 2012; v1 submitted 12 June, 2012; originally announced June 2012.

Comments: Contributed chapter in "50 Years of Nuclear BCS", edited by R. A. Broglia and V. Zelevinsky, references corrected

Fission modes of mercury isotopes

Authors: M. Warda, A. Staszczak, W. Nazarewicz

Abstract: 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… ▽ More 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. △ Less

Submitted 17 July, 2012; v1 submitted 25 May, 2012; originally announced May 2012.

Comments: 6 pages, 5 figures, to be published in Physical Review C

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

Authors: H. Nam, M. Stoitsov, W. Nazarewicz, A. Bulgac, G. Hagen, M. Kortelainen, P. Maris, J. C. Pei, K. J. Roche, N. Schunck, I. Thompson, J. P. Vary, S. M. Wild

Abstract: 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 resear… ▽ More 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. △ Less

Submitted 1 May, 2012; originally announced May 2012.

Comments: 12 pages, 15 figures, Proceedings for Conference on Computational Physics (CCP2011)

Journal ref: J. Phys.: Conf. Ser. 402 12033 (2012)

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

Authors: J. C. Pei, G. I. Fann, R. J. Harrison, W. Nazarewicz, J. Hill, D. Galindo, J. Jia

Abstract: 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,… ▽ More 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. △ Less

Submitted 23 April, 2012; originally announced April 2012.

Comments: Conference on Computational Physics (CCP 2011) Proceeding

Asymptotic normalization coefficients and continuum coupling in mirror nuclei

Authors: J. Okołowicz, N. Michel, W. Nazarewicz, M. Płoszajczak

Abstract: 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 approa… ▽ More 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%. △ Less

Submitted 17 March, 2012; originally announced March 2012.

Comments: 18 pages, 17 figures, 9 tables

On the origin of nuclear clustering

Authors: J. Okolowicz, M. Ploszajczak, W. Nazarewicz

Abstract: 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 ^… ▽ More 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. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Comments: Proceedings "Frontier Issues in Physics of Exotic Nuclei" (YKIS2011/DCEN2011), 14 pages, 9 figures

Journal ref: Prog. Theor. Phys. Supplement 196 (2012) 230

Electric dipole polarizability and the neutron skin

Authors: J. Piekarewicz, B. K. Agrawal, G. Colo, W. Nazarewicz, N. Paar, P. -G. Reinhard, X. Roca-Maza, D. Vretenar

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2012; originally announced January 2012.

Comments: Manuscript contains 5 pages, 2 figures, and 1 table. Submitted to Physical Review Letters

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

Authors: Yue Shi, J. Dobaczewski, S. Frauendorf, W. Nazarewicz, J. C. Pei, F. R. Xu, N. Nikolov

Abstract: 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… ▽ More 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. △ Less

Submitted 6 December, 2011; originally announced December 2011.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 108, 092501 (2012)

Nuclear energy density optimization: Large deformations

Authors: M. Kortelainen, J. McDonnell, W. Nazarewicz, P. -G. Reinhard, J. Sarich, N. Schunck, M. V. Stoitsov, S. M. Wild

Abstract: 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 funct… ▽ More 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. △ Less

Submitted 20 January, 2012; v1 submitted 18 November, 2011; originally announced November 2011.

Comments: 16 pages, 11 figures, accepted for publication in Phys. Rev. C

Advancing Nuclear Physics Through TOPS Solvers and Tools

Authors: E Ng, J Sarich, S M Wild, T Munson, H Aktulga, C Yang, P Maris, J P Vary, N Schunck, M G Bertolli, M Kortelainen, W Nazarewicz, T Papenbrock, M V Stoitsov

Abstract: At the heart of many scientific applications is the solution of algebraic systems, such as linear systems of equations, eigenvalue problems, and optimization problems, to name a few. TOPS, which stands for Towards Optimal Petascale Simulations, is a SciDAC applied math center focused on the development of solvers for tackling these algebraic systems, as well as the deployment of such technologies… ▽ More At the heart of many scientific applications is the solution of algebraic systems, such as linear systems of equations, eigenvalue problems, and optimization problems, to name a few. TOPS, which stands for Towards Optimal Petascale Simulations, is a SciDAC applied math center focused on the development of solvers for tackling these algebraic systems, as well as the deployment of such technologies in large-scale scientific applications of interest to the U.S. Department of Energy. In this paper, we highlight some of the solver technologies we have developed in optimization and matrix computations. We also describe some accomplishments achieved using these technologies in UNEDF, a SciDAC application project on nuclear physics. △ Less

Submitted 8 October, 2011; originally announced October 2011.

Comments: SciDAC 2011 Conference, July 10-14, 2011, Denver, CO; 5 pages, 2 tables, 2 figures

FRIB Theory Users Group Report: Joint ATLAS-HRIBF-NSCL-FRIB SuperUsers Meeting 18-20 August 2011

Authors: A. Baha Balantekin, Richard H. Cyburt, W. C. Haxton, Witek Nazarewicz, Filomena Nunes, Thomas Papenbrock, Scott Pratt, James Vary

Abstract: The FRIB (Facility for Rare Isotope Beams) Theory Users Group participated in the Joint ATLAS-HRIBF-NSCL-FRIB SuperUsers Meeting, hosted by Michigan State University August 18-20, 2011. Prior to the meeting a survey of the FRIB Theory Users Group was conducted to assess the health of the low-energy nuclear theory community and to identify perceived areas of need, in anticipation of the new demands… ▽ More The FRIB (Facility for Rare Isotope Beams) Theory Users Group participated in the Joint ATLAS-HRIBF-NSCL-FRIB SuperUsers Meeting, hosted by Michigan State University August 18-20, 2011. Prior to the meeting a survey of the FRIB Theory Users Group was conducted to assess the health of the low-energy nuclear theory community and to identify perceived areas of need, in anticipation of the new demands on theory that will accompany FRIB. Meeting discussions focused on survey results and on possible responses. These discussions are summarized here. △ Less

Submitted 27 September, 2011; originally announced September 2011.

Comments: 14 pages

Report number: UCB-NPAT-11-013,NR-LBNL-11-018

Charge radii and neutron correlations in helium halo nuclei

Authors: G. Papadimitriou, A. T. Kruppa, N. Michel, W. Nazarewicz, M. Płoszajczak, J. Rotureau

Abstract: 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… ▽ More 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. △ Less

Submitted 24 September, 2011; v1 submitted 1 September, 2011; originally announced September 2011.

Comments: 5 pages, 5 figures

Computing Heavy Elements

Authors: N. Schunck, A. Baran, M. Kortelainen, J. McDonnell, J. Moré, W. Nazarewicz, J. Pei, J. Sarich, J. Sheikh, A. Staszczak, M. Stoitsov, S. M. Wild

Abstract: 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) an… ▽ More 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. △ Less

Submitted 25 July, 2011; originally announced July 2011.

Comments: Proceeding of the Invited Talk given at the SciDAC 2011 conference, Jul. 10-15, 2011, Denver, CO

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

Authors: M. Stoitsov, H. Nam, W. Nazarewicz, A. Bulgac, G. Hagen, M. Kortelainen, J. C. Pei, K. J. Roche, N. Schunck, I. Thompson, J. P. Vary, S. M. Wild

Abstract: 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 lea… ▽ More 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. △ Less

Submitted 25 July, 2011; originally announced July 2011.

Comments: 5 pages, 8 figures