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

N. Schunck

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

ResearcherUnclaimed source record
nucl-thphysics.comp-phnucl-ex

Catalog footprint

What is connected

28works
3topics
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

28 published item(s)

preprint2022arXiv

Axially-deformed solution of the Skyrme-Hartree-Fock-Bogoliubov equations using the transformed harmonic oscillator basis (IV) hfbtho (v4.0): A new version of the program

We describe the new version 4.0 of the code hfbtho that solves the nuclear Hartree-Fock-Bogoliubov problem by using the deformed harmonic oscillator basis in cylindrical coordinates. In the new version, we have implemented the restoration of rotational, particle number, and reflection symmetry. The restoration of rotational symmetry does not require using bases closed under rotation. Furthermore, the SeaLL1 functional was added and the calculation of the Coulomb potential was improved. Finally, the code has been refactored to facilitate maintenance and future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two-body weak currents in heavy nuclei

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and $β$-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a density-matrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these unexpected effects are related to changes in the total nucleon density as the number of neutrons increases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Calibration of Energy Density Functionals with Deformed Nuclei

Nuclear density functional theory is the prevalent theoretical framework for accurately describing nuclear properties at the scale of the entire chart of nuclides. Given an energy functional and a many-body scheme (e.g., single- or multireference level), the predictive power of the theory depends strongly on how the parameters of the energy functionals have been calibrated with experimental data. Expanded algorithms and computing power have enabled recent optimization protocols to include data in deformed nuclei in order to optimize the coupling constants of the energy functional. The primary motivation of this work is to test the robustness of such protocols with respect to some of the technical and numerical details of the underlying calculations, especially when the calibration explores a large parameter space. To this end, we quantify the effect of these uncertainties on both the optimization and statistical emulation of composite objective functions. We also emphasize that Bayesian calibration can provide better estimates of the theoretical errors used to define objective functions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fission in a microscopic framework: from basic science to support for applications

Recent developments, both in theoretical modeling and computational power, have allowed us to make progress on a goal not fully achieved yet in nuclear theory: a microscopic theory of nuclear fission. Even if the complete microscopic description remains a computationally demanding task, the information that can be provided by current calculations can be extremely useful to guide and constrain more phenomenological approaches, which are simpler to implement. First, a microscopic model that describes the real-time dynamics of the fissioning system can justify or rule out some of the approximations. Second, the microscopic approach can be used to obtain trends, e.g., with increasing excitation energy of the fissioning system, or even to compute observables that cannot be otherwise calculated in phenomenological approaches or that can be hindered by the limitations of the method. We briefly present in this contribution the time-dependent superfluid local density approximation (TDSLDA) approach to nuclear fission, approach that has become a very successful theoretical model in many areas of many-body research. The TDSLDA incorporates the effects of the continuum, the dynamics of the pairing field, and the numerical solution is implemented with controlled approximations and negligible numerical corrections. The main part of the current contribution will be dedicated to discussing the method, and recent results concerning the fission dynamics. In addition, we present results on the excitation energy sharing between the fragments, which are in agreement with a qualitative conclusions extracted from a limited number of experimental measurements of properties of prompt neutrons.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Microscopically based energy density functionals for nuclei using the density matrix expansion

While ab initio many-body techniques have been able to successfully describe the properties of light and intermediate mass nuclei based on chiral effective field theory interactions, neutron-rich nuclei still remain out of reach for these methods. Conversely, energy density functional approaches can be used to calculate properties of heavy nuclei but rely mostly on phenomenological interactions. A usable form of the nuclear energy density functional that is rooted in the modern theory of nuclear forces was presented recently. The first component of this new set of functionals corresponds to the direct part (Hartree term) of the expectation value of local chiral potentials on a Slater determinant. The exchange term, which is a functional of the non-local density, is transformed into a local functional by applying the density matrix expansion. In order to reduce the computational cost due to the direct implementation of non-separable, local interactions in the Hartree term, we use an approximation to represent the regularized Yukawa functions in terms of a sum of (separable) Gaussian functions. These proceedings analyze the accuracy of such an approximation in terms of the number of Gaussian functions and look for an optimal value that gives an acceptable level of accuracy while maintaining the computational memory requirements in a many-body calculation as low as possible.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Propagation of Statistical Uncertainties of Skyrme Mass Models to Simulations of $r$-Process Nucleosynthesis

Uncertainties in nuclear models have a major impact on simulations that aim at understanding the origin of heavy elements in the universe through the rapid neutron capture process ($r$ process) of nucleosynthesis. Within the framework of the nuclear density functional theory, we use results of Bayesian statistical analysis to propagate uncertainties in the parameters of energy density functionals to the predicted $r$-process abundance pattern, by way not only of the nuclear masses but also through the influence of the masses on $β$-decay and neutron capture rates. We additionally make the first identifications of specific parameters of Skyrme-like energy density functionals which are correlated with particular aspects of the $r$-process abundance pattern. While previous studies have explored the reduction in the abundance pattern uncertainties due to anticipated new measurements of neutron-rich nuclei, here we point out that an even larger reduction will occur when these new measurements are used to reduce the uncertainty of model predictions of masses, which are then propagated through to the abundance pattern. We make a quantitative prediction for how large this reduction will be.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fission fragment charge and mass distributions in 239Pu(n,f) in the adiabatic nuclear energy density functional theory

Accurate knowledge of fission fragment yields is an essential ingredient of numerous applications ranging from the formation of elements in the r-process to fuel cycle optimization for nuclear energy. The need for a predictive theory applicable where no data is available is an incentive to develop a fully microscopic approach to fission dynamics. In this work, we calculate the pre-neutron emission charge and mass distributions of the fission fragments formed in the neutron-induced fission of 239Pu using a microscopic method based on nuclear energy density functional (EDF) method, where large amplitude collective motion is treated adiabatically using the time dependent generator coordinate method (TDGCM) under the Gaussian overlap approximation (GOA). Fission fragment distributions are extracted from the flux of the collective wave packet through the scission line. We find that the main characteristics of the fission charge and mass distributions can be well reproduced by existing energy functionals even in two-dimensional collective spaces. Theory and experiment agree typically within 2 mass units for the position of the asymmetric peak. As expected, calculations are sensitive to the structure of the initial state and the prescription for the collective inertia. We emphasize that results are also sensitive to the continuity of the collective landscape near scission. Our analysis confirms that the adiabatic approximation provides an effective scheme to compute fission fragment yields. It also suggests that, at least in the framework of nuclear DFT, three-dimensional collective spaces may be a prerequisite to reach 10% accuracy in predicting pre-neutron emission fission fragment yields.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Description of Induced Nuclear Fission with Skyrme Energy Functionals: II. Finite Temperature Effects

Understanding the mechanisms of induced nuclear fission for a broad range of neutron energies could help resolve fundamental science issues, such as the formation of elements in the universe, but could have also a large impact on societal applications in energy production or nuclear waste management. The goal of this paper is to set up the foundations of a microscopic theory to study the static aspects of induced fission as a function of the excitation energy of the incident neutron, from thermal to fast neutrons. To account for the high excitation energy of the compound nucleus, we employ a statistical approach based on finite-temperature nuclear density functional theory with Skyrme energy densities, which we benchmark on the 239 Pu(n,f) reaction. We compute the evolution of the least-energy fission pathway across multidimensional potential energy surfaces with up to five collective variables as a function of the nuclear temperature, and predict the evolution of both the inner and outer fission barriers as a function of the excitation energy of the compound nucleus. We show that the coupling to the continuum induced by the finite temperature is negligible in the range of neutron energies relevant for many applications of neutron-induced fission. We prove that the concept of quantum localization introduced recently can be extended to T > 0, and we apply the method to study the interaction energy and total kinetic energy of fission fragments as a function of the temperature for the most probable fission. While large uncertainties in theoretical modeling remain, we conclude that finite-temperature nuclear density functional may provide a useful framework to obtain accurate predictions of fission fragment properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Uncertainty Quantification and Propagation in Nuclear Density Functional Theory

Nuclear density functional theory (DFT) is one of the main theoretical tools used to study the properties of heavy and superheavy elements, or to describe the structure of nuclei far from stability. While on-going efforts seek to better root nuclear DFT in the theory of nuclear forces [see Duguet et al., this issue], energy functionals remain semi-phenomenological constructions that depend on a set of parameters adjusted to experimental data in finite nuclei. In this paper, we review recent efforts to quantify the related uncertainties, and propagate them to model predictions. In particular, we cover the topics of parameter estimation for inverse problems, statistical analysis of model uncertainties and Bayesian inference methods. Illustrative examples are taken from the literature.

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

Description of Induced Nuclear Fission with Skyrme Energy Functionals: I. Static Potential Energy Surfaces and Fission Fragment Properties

Eighty years after its experimental discovery, a microscopic description of induced nuclear fission based solely on the interactions between neutrons and protons and quantum many-body methods still poses formidable challenges. The goal of this paper is to contribute to the development of a predictive microscopic framework for the accurate calculation of static properties of fission fragments for hot fission and thermal or slow neutrons. To this end, we focus on the 239Pu(n,f) reaction and employ nuclear density functional theory with Skyrme energy densities. Potential energy surfaces are computed at the Hartree-Fock-Bogoliubov approximation with up to five collective variables. We find that the triaxial degree of freedom plays an important role, both near the fission barrier and at scission. The impact of the parameterization of the Skyrme energy density on deformation properties from the ground-state up to scission is also quantified. We introduce a general template for the detailed description of fission fragment properties. It is based on the careful analysis of the scission point, using both advanced topological methods and recently proposed quantum many-body techniques. We conclude that an accurate prediction of fission fragment properties at low incident neutron energies, although technologically demanding, should be within the reach of current nuclear density functional theory.

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

Quantification of Uncertainties in Nuclear Density Functional theory

Reliable predictions of nuclear properties are needed as much to answer fundamental science questions as in applications such as reactor physics or data evaluation. Nuclear density functional theory is currently the only microscopic, global approach to nuclear structure that is applicable throughout the nuclear chart. In the past few years, a lot of effort has been devoted to setting up a general methodology to assess theoretical uncertainties in nuclear DFT calculations. In this paper, we summarize some of the recent progress in this direction. Most of the new material discussed here will be be published in separate articles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Axially deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator basis (II) HFBTHO v2.00d: a new version of the program

We describe the new version 2.00d of the code HFBTHO that solves the nuclear Skyrme Hartree-Fock (HF) or Skyrme Hartree-Fock-Bogolyubov (HFB) problem by using the cylindrical transformed deformed harmonic-oscillator basis. In the new version, we have implemented the following features: (i) the modified Broyden method for non-linear problems, (ii) optional breaking of reflection symmetry, (iii) calculation of axial multipole moments, (iv) finite temperature formalism for the HFB method, (v) linear constraint method based on the approximation of the Random Phase Approximation (RPA) matrix for multi-constraint calculations, (vi) blocking of quasi-particles in the Equal Filling Approximation (EFA), (vii) framework for generalized energy density with arbitrary density-dependences, and (viii) shared memory parallelism via OpenMP pragmas.

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

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

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

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

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

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

Microscopically-based energy density functionals for nuclei using the density matrix expansion: Implementation and pre-optimization

In a recent series of papers, Gebremariam, Bogner, and Duguet derived a microscopically based nuclear energy density functional by applying the Density Matrix Expansion (DME) to the Hartree-Fock energy obtained from chiral effective field theory (EFT) two- and three-nucleon interactions. Due to the structure of the chiral interactions, each coupling in the DME functional is given as the sum of a coupling constant arising from zero-range contact interactions and a coupling function of the density arising from the finite-range pion exchanges. Since the contact contributions have essentially the same structure as those entering empirical Skyrme functionals, a microscopically guided Skyrme phenomenology has been suggested in which the contact terms in the DME functional are released for optimization to finite-density observables to capture short-range correlation energy contributions from beyond Hartree-Fock. The present paper is the first attempt to assess the ability of the newly suggested DME functional, which has a much richer set of density dependencies than traditional Skyrme functionals, to generate sensible and stable results for nuclear applications. The results of the first proof-of-principle calculations are given, and numerous practical issues related to the implementation of the new functional in existing Skyrme codes are discussed. Using a restricted singular value decomposition (SVD) optimization procedure, it is found that the new DME functional gives numerically stable results and exhibits a small but systematic reduction of our test $χ^2$ function compared to standard Skyrme functionals, thus justifying its suitability for future global optimizations and large-scale calculations.

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

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