Showing 1–22 of 22 results for author: Sarich, J

Axially deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator basis (III) hfbtho (v3.00): a new version of the program

Authors: R. Navarro Perez, N. Schunck, R. -D. Lasseri, C. Zhang, J. Sarich

Abstract: We describe the new version 3.00 of the code HFBTHO that solves the nuclear Hartree-Fock (HF) or 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 full Gogny force in both particle-hole and particle-particle channels, (ii) the calculation of the nuclear collective in… ▽ More We describe the new version 3.00 of the code HFBTHO that solves the nuclear Hartree-Fock (HF) or 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 full Gogny force in both particle-hole and particle-particle channels, (ii) the calculation of the nuclear collective inertia at the perturbative cranking approximation, (iii) the calculation of fission fragment charge, mass and deformations based on the determination of the neck (iv) the regularization of zero-range pairing forces (v) the calculation of localization functions (vi)MPI interface for large-scale mass table calculations. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: 29 pages, 3 figures, 4 tables; Submitted to Computer Physics Communications

Journal ref: Comput. Phys. Commun. 220, 363 (2017)

xSDK Foundations: Toward an Extreme-scale Scientific Software Development Kit

Authors: Roscoe Bartlett, Irina Demeshko, Todd Gamblin, Glenn Hammond, Michael Heroux, Jeffrey Johnson, Alicia Klinvex, Xiaoye Li, Lois Curfman McInnes, J. David Moulton, Daniel Osei-Kuffuor, Jason Sarich, Barry Smith, Jim Willenbring, Ulrike Meier Yang

Abstract: Extreme-scale computational science increasingly demands multiscale and multiphysics formulations. Combining software developed by independent groups is imperative: no single team has resources for all predictive science and decision support capabilities. Scientific libraries provide high-quality, reusable software components for constructing applications with improved robustness and portability.… ▽ More Extreme-scale computational science increasingly demands multiscale and multiphysics formulations. Combining software developed by independent groups is imperative: no single team has resources for all predictive science and decision support capabilities. Scientific libraries provide high-quality, reusable software components for constructing applications with improved robustness and portability. However, without coordination, many libraries cannot be easily composed. Namespace collisions, inconsistent arguments, lack of third-party software versioning, and additional difficulties make composition costly. The Extreme-scale Scientific Software Development Kit (xSDK) defines community policies to improve code quality and compatibility across independently developed packages (hypre, PETSc, SuperLU, Trilinos, and Alquimia) and provides a foundation for addressing broader issues in software interoperability, performance portability, and sustainability. The xSDK provides turnkey installation of member software and seamless combination of aggregate capabilities, and it marks first steps toward extreme-scale scientific software ecosystems from which future applications can be composed rapidly with assured quality and scalability. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: 14 pages

ACM Class: D.2.0; D.2.2; D.2.11

Uncertainty Quantification for Optical Model Parameters

Authors: A. E. Lovell, F. M. Nunes, J. Sarich, S. M. Wild

Abstract: Although uncertainty quantification has been making its way into nuclear theory, these methods have yet to be explored in the context of reaction theory. For example, it is well known that different parameterizations of the optical potential can result in different cross sections, but these differences have not been systematically studied and quantified. The purpose of this work is to investigate… ▽ More Although uncertainty quantification has been making its way into nuclear theory, these methods have yet to be explored in the context of reaction theory. For example, it is well known that different parameterizations of the optical potential can result in different cross sections, but these differences have not been systematically studied and quantified. The purpose of this work is to investigate the uncertainties in nuclear reactions that result from fitting a given model to elastic-scattering data, as well as to study how these uncertainties propagate to the inelastic and transfer channels. We use statistical methods to determine a best fit and create corresponding 95\% confidence bands. A simple model of the process is fit to elastic-scattering data and used to predict either inelastic or transfer cross sections. In this initial work, we assume that our model is correct, and the only uncertainties come from the variation of the fit parameters. We study a number of reactions involving neutron and deuteron projectiles with energies in the range of 5-25 MeV/u, on targets with mass $A$=12-208. We investigate the correlations between the parameters in the fit. The case of deuterons on $^{12}$C is discussed in detail: the elastic-scattering fit and the prediction of $^{12}$C(d,p)$^{13}$C transfer angular distributions, using both uncorrelated and correlated $χ^2$ minimization functions. The general features for all cases are compiled in a systematic manner to identify trends. Our work shows that, in many cases, the correlated $χ^2$ functions (in comparison to the uncorrelated $χ^2$ functions) provide a more natural parameterization of the process. These correlated functions do, however, produce broader confidence bands. Further optimization may require improvement in the models themselves and/or more information included in the fit. △ Less

Submitted 27 January, 2017; v1 submitted 4 November, 2016; originally announced November 2016.

Comments: 11 pages, 7 figures, accepted to Phys. Rev. C

Journal ref: Phys. Rev. C 95, 024611 (2017)

One-nucleon transfer reactions and the optical potential

Authors: F. M. Nunes, A. Lovell, A. Ross, L. J. Titus, R. J. Charity, W. H. Dickhoff, M. H. Mahzoon, J. Sarich, S. M. Wild

Abstract: We provide a summary of new developments in the area of direct reaction theory with a particular focus on one-nucleon transfer reactions. We provide a status of the methods available for describing (d,p) reactions. We discuss the effects of nonlocality in the optical potential in transfer reactions. The results of a purely phenomenological potential and the optical potential obtained from the disp… ▽ More We provide a summary of new developments in the area of direct reaction theory with a particular focus on one-nucleon transfer reactions. We provide a status of the methods available for describing (d,p) reactions. We discuss the effects of nonlocality in the optical potential in transfer reactions. The results of a purely phenomenological potential and the optical potential obtained from the dispersive optical model are compared; both point toward the importance of including nonlocality in transfer reactions explicitly. Given the large ambiguities associated with optical potentials, we discuss some new developments toward the quantification of this uncertainty. We conclude with some general comments and a brief account of new advances that are in the pipeline. △ Less

Submitted 15 September, 2015; originally announced September 2015.

Comments: 7 pages, 5 figures, proceedings for the 14th International Conference on Nuclear Reaction Mechanisms, Varenna, June 2015

Uncertainty Quantification and Propagation in Nuclear Density Functional Theory

Authors: N. Schunck, J. D. McDonnell, D. Higdon, J. Sarich, S. M. Wild

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

Submitted 19 March, 2015; originally announced March 2015.

Comments: Proceedings of the "Second International Workshop on Perspectives on Nuclear Data for the Next Decade", 14-17 October 2014, Bruyères-le-Châtel, France

Journal ref: Eur. Phys. J. A (2015) 51: 169

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

Authors: J. D. McDonnell, N. Schunck, D. Higdon, J. Sarich, S. M. Wild, W. Nazarewicz

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

Submitted 14 January, 2015; originally announced January 2015.

Comments: 5 pages, 4 figures, 1 table; Supplemental material with 2 tables; Submitted to Physical Review Letters

Journal ref: Phys. Rev. Lett. 114, 122501 (2015)

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

A Bayesian Approach for Parameter Estimation and Prediction using a Computationally Intensive Model

Authors: Dave Higdon, Jordan D. McDonnell, Nicolas Schunck, Jason Sarich, Stefan M. Wild

Abstract: Bayesian methods have been very successful in quantifying uncertainty in physics-based problems in parameter estimation and prediction. In these cases, physical measurements y are modeled as the best fit of a physics-based model $η(θ)$ where $θ$ denotes the uncertain, best input setting. Hence the statistical model is of the form $y = η(θ) + ε$, where $ε$ accounts for measurement, and possibly oth… ▽ More Bayesian methods have been very successful in quantifying uncertainty in physics-based problems in parameter estimation and prediction. In these cases, physical measurements y are modeled as the best fit of a physics-based model $η(θ)$ where $θ$ denotes the uncertain, best input setting. Hence the statistical model is of the form $y = η(θ) + ε$, where $ε$ accounts for measurement, and possibly other error sources. When non-linearity is present in $η(\cdot)$, the resulting posterior distribution for the unknown parameters in the Bayesian formulation is typically complex and non-standard, requiring computationally demanding computational approaches such as Markov chain Monte Carlo (MCMC) to produce multivariate draws from the posterior. While quite generally applicable, MCMC requires thousands, or even millions of evaluations of the physics model $η(\cdot)$. This is problematic if the model takes hours or days to evaluate. To overcome this computational bottleneck, we present an approach adapted from Bayesian model calibration. This approach combines output from an ensemble of computational model runs with physical measurements, within a statistical formulation, to carry out inference. A key component of this approach is a statistical response surface, or emulator, estimated from the ensemble of model runs. We demonstrate this approach with a case study in estimating parameters for a density functional theory (DFT) model, using experimental mass/binding energy measurements from a collection of atomic nuclei. We also demonstrate how this approach produces uncertainties in predictions for recent mass measurements obtained at Argonne National Laboratory (ANL). △ Less

Submitted 17 September, 2014; v1 submitted 10 July, 2014; originally announced July 2014.

Comments: 19 pages, 10 figures, 1 table; Invited paper for the Journal of Physics G: Nuclear and Particle Physics focus section entitled "Enhancing the interaction between nuclear experiment and theory through information and statistics". In press

Journal ref: J. Phys. G: Nucl. Part. Phys. 42, 034009 (2015)

Derivative-free optimization for parameter estimation in computational nuclear physics

Authors: Stefan M. Wild, Jason Sarich, Nicolas Schunck

Abstract: We consider optimization problems that arise when estimating a set of unknown parameters from experimental data, particularly in the context of nuclear density functional theory. We examine the cost of not having derivatives of these functionals with respect to the parameters. We show that the POUNDERS code for local derivative-free optimization obtains consistent solutions on a variety of computa… ▽ More We consider optimization problems that arise when estimating a set of unknown parameters from experimental data, particularly in the context of nuclear density functional theory. We examine the cost of not having derivatives of these functionals with respect to the parameters. We show that the POUNDERS code for local derivative-free optimization obtains consistent solutions on a variety of computationally expensive energy density functional calibration problems. We also provide a primer on the operation of the POUNDERS software in the Toolkit for Advanced Optimization. △ Less

Submitted 17 September, 2014; v1 submitted 20 June, 2014; originally announced June 2014.

Comments: 17 pages, 4 figures, 5 tables; Invited paper for the Journal of Physics G: Nuclear and Particle Physics focus section entitled "Enhancing the interaction between nuclear experiment and theory through information and statistics". In press

Journal ref: J. Phys. G: Nucl. Part. Phys. 42, 034031 (2015)

Error Analysis in Nuclear Density Functional Theory

Authors: Nicolas Schunck, Jordan D. McDonnell, Jason Sarich, Stefan M. Wild, Dave Higdon

Abstract: Nuclear density functional theory (DFT) is the only microscopic, global approach to the structure of atomic nuclei. It is used in numerous applications, from determining the limits of stability to gaining a deep understanding of the formation of elements in the universe or the mechanisms that power stars and reactors. The predictive power of the theory depends on the amount of physics embedded in… ▽ More Nuclear density functional theory (DFT) is the only microscopic, global approach to the structure of atomic nuclei. It is used in numerous applications, from determining the limits of stability to gaining a deep understanding of the formation of elements in the universe or the mechanisms that power stars and reactors. The predictive power of the theory depends on the amount of physics embedded in the energy density functional as well as on efficient ways to determine a small number of free parameters and solve the DFT equations. In this article, we discuss the various sources of uncertainties and errors encountered in DFT and possible methods to quantify these uncertainties in a rigorous manner. △ Less

Submitted 11 July, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Comments: 18 pages, 3 figures, 4 tables; Invited paper for the Journal of Physics G: Nuclear and Particle Physics focus section entitled "Enhancing the interaction between nuclear experiment and theory through information and statistics"; Revised version after comments by the referees: Figure 1 and Table 4 have been corrected

Journal ref: J. Phys. G: Nucl. Part. Phys. 42, 034024 (2015)

Quantification of Uncertainties in Nuclear Density Functional theory

Authors: N. Schunck, J. D. McDonnell, D. Higdon, J. Sarich, S. Wild

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

Submitted 17 September, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Comments: 4 pages, 2 figure; Proceedings of the International Workshop on Nuclear Data Covariances, April 28 - May 1, 2014, Santa Fe, New Mexico, USA - Second revised version after referee's comments

Journal ref: Nucl. Data Sheets 123, 115 (2015)

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)

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)

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

Authors: M. V. Stoitsov, N. Schunck, M. Kortelainen, N. Michel, H. Nam, E. Olsen, J. Sarich, S. Wild

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

Submitted 20 January, 2013; v1 submitted 5 October, 2012; originally announced October 2012.

Comments: 29 pages, 5 figures, 5 tables accepted for publication in Computer Physics Communications; Program will be hosted on CPC Program Library. This version includes a much expanded section on the treatment of the Coulomb potential following the referee's comments, and an updated Acknowledgements section

Journal ref: Comput. Phys. Commun. 184, 1592-1604 (2013)

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

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

Nuclear Energy Density Optimization

Authors: M. Kortelainen, T. Lesinski, J. Moré, W. Nazarewicz, J. Sarich, N. Schunck, M. V. Stoitsov, S. Wild

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

Submitted 27 May, 2010; originally announced May 2010.

Comments: 19 pages, 11 figures

Journal ref: Phys.Rev.C82:024313,2010

One-quasiparticle States in the Nuclear Energy Density Functional Theory

Authors: N. Schunck, J. Dobaczewski, J. McDonnell, J. More, W. Nazarewicz, J. Sarich, M. V. Stoitsov

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

Submitted 28 January, 2010; v1 submitted 12 October, 2009; originally announced October 2009.

Comments: Accepted for publication in Phys. Rev. C after third report from the referee. One reference fixed

Journal ref: Phys.Rev.C81:024316,2010

HFODD (v2.40h) User's Guide

Authors: J. Dobaczewski, B. G. Carlsson, J. Dudek, J. Engel, P. Olbratowski, P. Powałowski, M. Sadziak, J. Sarich, W. Satuła, N. Schunck, A. Staszczak, M. Stoitsov, M. Zalewski, H. Zduńczuk

Abstract: We describe the input data and installation procedures of the code HFODD (v2.40h). The present write-up contains complete and comprehensive information that has originally been given in six independent publications. It is enhanced by the subject index and indexes of variables, input-data keywords, subroutines, and files that are used in this user guide. We describe the input data and installation procedures of the code HFODD (v2.40h). The present write-up contains complete and comprehensive information that has originally been given in six independent publications. It is enhanced by the subject index and indexes of variables, input-data keywords, subroutines, and files that are used in this user guide. △ Less

Submitted 20 September, 2009; originally announced September 2009.

Comments: Comments: 97 LaTeX pages, 4 EPS figures, on-line-only publication, see also http://www.fuw.edu.pl/~dobaczew/hfodd/hfodd.html

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

Authors: J. Dobaczewski, W. Satula, B. G. Carlsson, J. Engel, P. Olbratowski, P. Powalowski, M. Sadziak, J. Sarich, N. Schunck, A. Staszczak, M. Stoitsov, M. Zalewski, H. Zdunczuk

Abstract: We describe the new version (v2.38j) of the code HFODD which solves the nuclear Skyrme-Hartree-Fock or Skyrme-Hartree-Fock-Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we have implemented: (i) projection on good angular momentum (for the Hartree-Fock states), (ii) calculation of the GCM kernels, (iii) calculation of matrix elements of the Yuka… ▽ More We describe the new version (v2.38j) of the code HFODD which solves the nuclear Skyrme-Hartree-Fock or Skyrme-Hartree-Fock-Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we have implemented: (i) projection on good angular momentum (for the Hartree-Fock states), (ii) calculation of the GCM kernels, (iii) calculation of matrix elements of the Yukawa interaction, (iv) the BCS solutions for state-dependent pairing gaps, (v) the HFB solutions for broken simplex symmetry, (vi) calculation of Bohr deformation parameters, (vii) constraints on the Schiff moments and scalar multipole moments, (viii) the D2h transformations and rotations of wave functions, (ix) quasiparticle blocking for the HFB solutions in odd and odd-odd nuclei, (x) the Broyden method to accelerate the convergence, (xi) the Lipkin-Nogami method to treat pairing correlations, (xii) the exact Coulomb exchange term, (xiii) several utility options, and we have corrected two insignificant errors. △ Less

Submitted 5 March, 2009; originally announced March 2009.

Comments: 45 LaTeX pages, 4 figures, submitted to Computer Physics Communications

Journal ref: Comput. Phys. Commun. 180 (2009) 2361