Source author record

Witold Nazarewicz

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2023arXiv

Combined theoretical analysis of the parity-violating asymmetry for ${}^{48}$Ca and ${}^{208}Pb$

The recent experimental determination of the parity violating asymmetry $A_{\rm pv}$ in ${}^{48}$Ca and ${}^{208}$Pb at Jefferson Lab is important for our understanding on how neutrons and protons arrange themselves inside the atomic nucleus. To better understand the impact of these measurements, we present a rigorous theoretical investigation of $A_{\rm pv}$ in ${}^{48}$Ca and ${}^{208}$Pb and assess the associated uncertainties. We complement our study by inspecting the static electric dipole polarizability in these nuclei. The analysis is carried out within nuclear energy density functional theory with quantified input. We conclude that the simultaneous accurate description of $A_{\rm pv}$ in ${}^{48}$Ca and ${}^{208}$Pb cannot be achieved by our models that accommodate a pool of global nuclear properties, such as masses and charge radii, throughout the nuclear chart, and describe -- within one standard deviation -- the experimental dipole polarizabilities $α_{\rm D}$ in these nuclei.

preprint2022arXiv

Evidence of Two-Source King Plot Nonlinearity in Spectroscopic Search for New Boson

Optical precision spectroscopy of isotope shifts can be used to test for new forces beyond the Standard Model, and to determine basic properties of atomic nuclei. We measure isotope shifts on the highly forbidden ${}^2S_{1/2} \rightarrow {}^2F_{7/2}$ octupole transition of trapped $^{168,170,172,174,176}$Yb ions. When combined with previous measurements in Yb$^+$ and very recent measurements in Yb, the data reveal a King plot nonlinearity of up to 240$σ$. The trends exhibited by experimental data are explained by nuclear density functional theory calculations with the Fayans functional. We also find, with 4.3$σ$ confidence, that there is a second distinct source of nonlinearity, and discuss its possible origin.

preprint2022arXiv

Information content of the differences in the charge radii of mirror nuclei

Differences in the charge radii of mirror nuclei have been recently suggested to contain information on the slope of the symmetry energy L. To test this hypothesis, we perform statistical correlation analysis using quantified energy density functionals that are consistent with our previous knowledge on global nuclear observables such as binding energies and charge radii. We conclude that the difference in charge radii between a mirror pair, R_mir, is an inferior isovector indicator compared to other observables, such at the neutron skin or electric dipole polarizability. In particular, this quantity correlates poorly with both the neutron skin and L. We demonstrate that R_mir is influenced by pairing correlations in the presence of low-lying proton continuum in the proton-rich mirror-partner nucleus. Considering the large theoretical uncertainties on R_mir, we conclude that the precise data on mirror charge radii cannot provide a stringent constraint on L.

preprint2022arXiv

Nudged elastic band approach to nuclear fission pathways

The nuclear fission process is a dramatic example of the large-amplitude collective motion in which the nucleus undergoes a series of shape changes before splitting into distinct fragments. This motion can be represented by a pathway in the many-dimensional space of collective coordinates. The collective action along the fission pathway determines the spontaneous fission half-lives as well as mass and charge distributions of fission fragments. We study the performance and precision of various methods to determine the minimum action and minimum-energy fission trajectories in the collective space. We apply the nudged elastic band method (NEB), grid-based methods, and Euler Lagrange approach to the collective action minimization in two and three dimensional collective spaces. The performance of various approaches to the fission pathway problem is assessed by studying the collective motion along both analytic energy surfaces and realistic potential energy surfaces obtained with the Hartree-Fock-Bogoliubov theory. The uniqueness and stability of the solutions is studied. The NEB method is capable of efficient determination of the exit points on the outer turning surface that characterize the most probable fission pathway and constitute the key input for fission studies. This method can also be used to accurately compute the critical points (i.e., local minima and saddle points) on the potential energy surface of the fissioning nucleus that determine the static fission path. The NEB method is the tool of choice for finding the least-action and minimum energy fission trajectories. It will be particularly useful in large-scale fission calculation of superheavy nuclei and neutron-rich fissioning nuclei contributing to the astrophysical r-process recycling.

preprint2022arXiv

Theoretical description of fission yields: towards a fast and efficient global model

Background: A quantitative microscopic understanding of the fission-fragment yield distributions represents a major challenge for nuclear theory as it involves the intricate competition between large-amplitude nuclear collective motion and single-particle nucleonic motion. Purpose: A recently proposed approach to global modeling of fission fragment distributions is extended to account for odd-even staggering in charge yields and for neutron evaporation. Method: Fission trajectories are obtained within the density functional theory framework, allowing for a microscopic determination of the most probable fission prefragment configurations. Mass and charge yields distributions are constructed by means of a statistical approach rooted in a microcanonical ensemble. Result: We show that the proposed hybrid model can reproduce experimental mass and charge fragment yields, including the odd-even staggering, for a wide range of fissioning nuclei. Experimental isotopic yields can be described within a simple neutron evaporation scheme. We also explore fission fragment distributions of exotic neutron-rich and superheavy systems, and compare our predictions with other state-of-the art global calculations. Conclusion: Our study suggests that the microscopic rearrangement of nucleons into fission fragments occurs well before the scission, and that the subsequent dynamics is mainly driven by the thermal excitations and bulk features of the nuclear binding. The proposed simple hybrid approach is well suited for large-scale calculations involving hundreds of fissioning nuclei.

preprint2022arXiv

Three-dimensional Skyrme Hartree-Fock-Bogoliubov solver in coordinate-space representation

The coordinate-space representation of the Hartree-Fock-Bogoliubov theory is the method of choice to study weakly bound nuclei whose properties are affected by the quasiparticle continuum space. To describe such systems, we developed a three-dimensional Skyrme-Hartree-Fock-Bogoliubov solver HFBFFT based on the existing, highly optimized and parallelized Skyrme-Hartree-Fock code Sky3D. The code does not impose any self-consistent spatial symmetries such as mirror inversions or parity. The underlying equations are solved in HFBFFT directly in the canonical basis using the fast Fourier transform. To remedy the problems with pairing collapse, we implemented the soft energy cutoff and pairing annealing. The convergence of HFB solutions was improved by a sub-iteration method. The Hermiticity violation of differential operators brought by Fourier-transform-based differentiation has also been solved. The accuracy and performance of HFBFFT were tested by benchmarking it against other HFB codes, both spherical and deformed, for a set of nuclei, both well-bound and weakly-bound.

preprint2021arXiv

Microscopic origin of reflection-asymmetric nuclear shapes

Background: The presence of nuclear ground states with stable reflection-asymmetric shapes is supported by rich experimental evidence. Theoretical surveys of odd-multipolarity deformations predict the existence of pear-shaped isotopes in several fairly localized regions of the nuclear landscape in the vicinity of near-lying single-particle shells with $Δ\ell=Δj=3$. Purpose: We analyze the role of isoscalar, isovector, neutron-proton, neutron-neutron, and proton-proton multipole interaction energies in inducing the onset of reflection-asymmetric ground-state deformations. Methods: The calculations are performed in the framework of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory using two Skyrme energy density functionals and density-dependent pairing force. Results: We show that reflection-asymmetric ground-state shapes of atomic nuclei are driven by the odd-multipolarity neutron-proton (or isoscalar) part of the nuclear interaction energy. This result is consistent with the particle-vibration picture, in which the main driver of octupole instability is the isoscalar octupole-octupole interaction giving rise to large $E3$ polarizability.

preprint2021arXiv

Universal trend of charge radii of even-even Ca-Zn nuclei

Radii of nuclear charge distributions carry information about the strong and electromagnetic forces acting inside the atomic nucleus. While the global behavior of nuclear charge radii is governed by the bulk properties of nuclear matter, their local trends are affected by quantum motion of proton and neutron nuclear constituents. The measured differential charge radii $δ\langle r^2_c\rangle$ between neutron numbers $N=28$ and $N=40$ exhibit a universal pattern as a function of $n=N-28$ that is independent of the atomic number. Here we analyze this remarkable behavior in even-even nuclei from calcium to zinc using two state-of-the-art theories based on quantified nuclear interactions: the ab-initio coupled cluster theory and nuclear density functional theory. Both theories reproduce the smooth rise of differential charge radii and their weak dependence on the atomic number. By considering a large set of isotopic chains, we show that this trend can be captured by just two parameters: the slope and curvature of ${δ\langle r^2_c\rangle(n)}$. We demonstrate that these parameters show appreciable model dependence, and the statistical analysis indicates that they are not correlated with any single model property, i.e., they are impacted by both bulk nuclear properties as well as shell structure.

preprint2020arXiv

Beyond the proton drip line: Bayesian analysis of proton-emitting nuclei

The limits of the nuclear landscape are determined by nuclear binding energies. Beyond the proton drip lines, where the separation energy becomes negative, there is not enough binding energy to prevent protons from escaping the nucleus. Predicting properties of unstable nuclear states in the vast territory of proton emitters poses an appreciable challenge for nuclear theory as it often involves far extrapolations. In addition, significant discrepancies between nuclear models in the proton-rich territory call for quantified predictions. With the help of Bayesian methodology, we mix a family of nuclear mass models corrected with statistical emulators trained on the experimental mass measurements, in the proton-rich region of the nuclear chart. Separation energies were computed within nuclear density functional theory using several Skyrme and Gogny energy density functionals. We also considered mass predictions based on two models used in astrophysical studies. Quantified predictions were obtained for each model using Bayesian Gaussian processes trained on separation-energy residuals and combined via Bayesian model averaging. We obtained a good agreement between averaged predictions of statistically corrected models and experiment. In particular, we quantified model results for one- and two-proton separation energies and derived probabilities of proton emission. This information enabled us to produce a quantified landscape of proton-rich nuclei. The most promising candidates for two-proton decay studies have been identified. The methodology used in this work has broad applications to model-based extrapolations of various nuclear observables. It also provides a reliable uncertainty quantification of theoretical predictions.

preprint2020arXiv

Efficient method for estimation of fission fragment yields of r-process nuclei

$\textbf{Background}$ More than half of all the elements heavier than iron are made by the rapid neutron capture process (or r process). For very neutron-rich astrophysical conditions, such at those found in the tidal ejecta of neutron stars, nuclear fission determines the r-process endpoint, and the fission fragment yields shape the final abundances of $110\le A \le 170$ nuclei. The knowledge of fission fragment yields of hundreds of nuclei inhabiting very neutron-rich regions of the nuclear landscape is thus crucial for the modeling of heavy-element nucleosynthesis. $\textbf{Purpose}$ In this study, we propose a model for the fast calculation of fission fragment yields based on the concept of shell-stabilized prefragments defined with help of the nucleonic localization functions. $\textbf{Methods}$ To generate realistic potential energy surfaces and nucleonic localizations, we apply Skyrme Density Functional Theory. The distribution of the neck nucleons among the two prefragments is obtained by means of a statistical model. $\textbf{Results}$ We benchmark the method by studying the fission yields of $^{178}$Pt, $^{240}$Pu, $^{254}$Cf, and $^{254,256,258}$Fm and show that it satisfactorily explains the experimental data. We then make predictions for $^{254}$Pu and $^{290}$Fm as two representative cases of fissioning nuclei that are expected to significantly contribute during the r-process nucleosynthesis occurring in neutron star mergers. $\textbf{Conclusions}$ The proposed framework provides an efficient alternative to microscopic approaches based on the evolution of the system in a space of collective coordinates all the way to scission. It can be used to carry out global calculations of fission fragment distributions across the r-process region.

preprint2020arXiv

Electron and Nucleon Localization Functions of Oganesson: Approaching the Thomas-Fermi Limit

Fermion localization functions are used to discuss electronic and nucleonic shell structure effects in the superheavy element oganesson, the heaviest element discovered to date. Spin-orbit splitting in the $7p$ electronic shell becomes so large ($\sim$ 10 eV) that Og is expected to show uniform-gas-like behavior in the valence region with a rather large dipole polarizability compared to the lighter rare gas elements. The nucleon localization in Og is also predicted to undergo a transition to the Thomas-Fermi gas behavior in the valence region. This effect, particularly strong for neutrons, is due to the high density of single-particle orbitals.

preprint2020arXiv

Landscape of pear-shaped even-even nuclei

The phenomenon of reflection-asymmetric nuclear shapes is relevant to nuclear stability, nuclear spectroscopy, nuclear decays and fission, and the search for new physics beyond the standard model. Global surveys of ground-state octupole deformation, performed with a limited number of models, suggest that the number of pear-shaped isotopes is fairly limited across the nuclear landscape. We carry out global analysis of ground-state octupole deformations for particle-bound even-even nuclei with $Z \leq 110$ and $N \leq 210$ using nuclear density functional theory (DFT) with several non-relativistic and covariant energy density functionals. In this way, we can identify the best candidates for reflection-asymmetric shapes. The calculations are performed in the frameworks of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory and relativistic Hartree-Bogoliubov theory using DFT solvers employing harmonic oscillator basis expansion. We consider five Skyrme and four covariant energy density functionals. We predict several regions of ground-state octupole deformation. In addition to the "traditional" regions of neutron-deficient actinide nuclei around $^{224}$Ra and neutron-rich lanthanides around $^{146}$Ba, we identified vast regions of reflecion-asymmetric shapes in very neutron-rich nuclei around $^{200}$Gd and $^{288}$Pu, as well as in several nuclei around $^{112}$Ba. Our analysis suggests several promising candidates with stable ground-state octupole deformation, primarily in the neutron-deficient actinide region, that can be reached experimentally. Detailed comparison between Skyrme and covariant models is performed.

preprint2020arXiv

Neutron drip line in the Ca region from Bayesian model averaging

The region of heavy calcium isotopes forms the frontier of experimental and theoretical nuclear structure research where the basic concepts of nuclear physics are put to stringent test. The recent discovery of the extremely neutron-rich nuclei around $^{60}$Ca [Tarasov, 2018] and the experimental determination of masses for $^{55-57}$Ca (Michimasa, 2018] provide unique information about the binding energy surface in this region. To assess the impact of these experimental discoveries on the nuclear landscape's extent, we use global mass models and statistical machine learning to make predictions, with quantified levels of certainty, for bound nuclides between Si and Ti. Using a Bayesian model averaging analysis based on Gaussian-process-based extrapolations we introduce the posterior probability $p_{ex}$ for each nucleus to be bound to neutron emission. We find that extrapolations for drip-line locations, at which the nuclear binding ends, are consistent across the global mass models used, in spite of significant variations between their raw predictions. In particular, considering the current experimental information and current global mass models, we predict that $^{68}$Ca has an average posterior probability ${p_{ex}\approx76}$% to be bound to two-neutron emission while the nucleus $^{61}$Ca is likely to decay by emitting a neutron (${p_{ex}\approx 46}$ %).

preprint2020arXiv

Quantified limits of the nuclear landscape

The chart of the nuclides is limited by particle drip lines beyond which nuclear stability to proton or neutron emission is lost. Predicting the range of particle-bound isotopes poses an appreciable challenge for nuclear theory as it involves extreme extrapolations of nuclear masses beyond the regions where experimental information is available. Still, quantified extrapolations are crucial for a variety of applications, including the modeling of stellar nucleosynthesis. We use microscopic nuclear mass models and Bayesian methodology to provide quantified predictions of proton and neutron separation energies as well as Bayesian probabilities of existence throughout the nuclear landscape all the way to the particle drip lines. We apply nuclear density functional theory with several energy density functionals. To account for uncertainties, Bayesian Gaussian processes are trained on the separation-energy residuals for each individual model, and the resulting predictions are combined via Bayesian model averaging. This framework allows to account for systematic and statistical uncertainties and propagate them to extrapolative predictions. We characterize the drip-line regions where the probability that the nucleus is particle-bound decreases from $1$ to $0$. In these regions, we provide quantified predictions for one- and two-nucleon separation energies. According to our Bayesian model averaging analysis, 7759 nuclei with $Z\leq 119$ have a probability of existence $\geq 0.5$. The extrapolations obtained in this study will be put through stringent tests when new experimental information on exotic nuclei becomes available. In this respect, the quantified landscape of nuclear existence obtained in this study should be viewed as a dynamical prediction that will be fine-tuned when new experimental information and improved global mass models become available.

preprint2020arXiv

Report from the A.I. For Nuclear Physics Workshop

This report is an outcome of the workshop "AI for Nuclear Physics" held at Thomas Jefferson National Accelerator Facility on March 4-6, 2020. The workshop brought together 184 scientists to explore opportunities for Nuclear Physics in the area of Artificial Intelligence. The workshop consisted of plenary talks, as well as six working groups. The report includes the workshop deliberations and additional contributions to describe prospects for using AI across Nuclear Physics research.

preprint2020arXiv

Statistical aspects of nuclear mass models

We study the information content of nuclear masses from the perspective of global models of nuclear binding energies. To this end, we employ a number of statistical methods and diagnostic tools, including Bayesian calibration, Bayesian model averaging, chi-square correlation analysis, principal component analysis, and empirical coverage probability. Using a Bayesian framework, we investigate the structure of the 4-parameter Liquid Drop Model by considering discrepant mass domains for calibration. We then use the chi-square correlation framework to analyze the 14-parameter Skyrme energy density functional calibrated using homogeneous and heterogeneous datasets. We show that a quite dramatic parameter reduction can be achieved in both cases. The advantage of Bayesian model averaging for improving uncertainty quantification is demonstrated. The statistical approaches used are pedagogically described; in this context this work can serve as a guide for future applications.

preprint2016arXiv

Challenges in Nuclear Structure Theory

The goal of nuclear structure theory is to build a comprehensive microscopic framework in which properties of nuclei and extended nuclear matter, and nuclear reactions and decays can all be consistently described. Due to novel theoretical concepts, breakthroughs in the experimentation with rare isotopes, increased exchange of ideas across different research areas, and the progress in computer technologies and numerical algorithms, nuclear theorists have been quite successful in solving various bits and pieces of the nuclear many-body puzzle and the prospects are exciting. This article contains a brief, personal perspective on the status of the field.

preprint2016arXiv

Clustering and pasta phases in nuclear density functional theory

Nuclear density functional theory (DFT) is the tool of choice in describing properties of complex nuclei and intricate phases of bulk nucleonic matter. It is a microscopic approach based on an energy density functional representing the nuclear interaction. An attractive feature of nuclear DFT is that it can be applied to both finite nuclei and pasta phases appearing in the inner crust of neutron stars. While nuclear pasta clusters in a neutron star can be easily characterized through their density distributions, the level of clustering of nucleons in a nucleus can often be difficult to assess. To this end, we use the concept of nucleonic localization. We demonstrate that the localization measure provides us with fingerprints of clusters in light and heavy nuclei, including fissioning systems. Furthermore we investigate the rod-like pasta phase using twist-averaged boundary conditions, which enable calculations in finite volumes accessible by state of the art DFT solvers.

preprint2016arXiv

Impact of nuclear mass uncertainties on the $r$-process

Nuclear masses play a fundamental role in understanding how the heaviest elements in the Universe are created in the $r$-process. We predict $r$-process nucleosynthesis yields using neutron capture and photodissociation rates that are based on nuclear density functional theory. Using six Skyrme energy density functionals based on different optimization protocols, we determine for the first time systematic uncertainty bands -- related to mass modeling -- for $r$-process abundances in realistic astrophysical scenarios. We find that features of the underlying microphysics make an imprint on abundances especially in the vicinity of neutron shell closures: abundance peaks and troughs are reflected in trends of neutron separation energy. Further advances in nuclear theory and experiments, when linked to observations, will help in the understanding of astrophysical conditions in extreme $r$-process sites.

preprint2016arXiv

Pairing Nambu-Goldstone Modes within Nuclear Density Functional Theory

We show that the Nambu-Goldstone formalism of the broken gauge symmetry in the presence of the $T=1$ pairing condensate offers a quantitative description of the binding-energy differences of open-shell superfluid nuclei. We conclude that the pairing-rotational moments of inertia are excellent pairing indicators, which are free from ambiguities attributed to odd-mass systems. We offer a new, unified interpretation of the binding-energy differences traditionally viewed in the shell model picture as signatures of the valence nucleon properties. We present the first systematic analysis of the off-diagonal pairing-rotational moments of inertia and demonstrate the mixing of the neutron and proton pairing-rotational modes in the ground states of even-even nuclei. Finally, we discuss the importance of mass measurements of neutron-rich nuclei for constraining the pairing energy density functional.

preprint2015arXiv

Complex-energy approach to sum rules within nuclear density functional theory

The linear response of the nucleus to an external field contains unique information about the effective interaction, correlations, and properties of its excited states. To characterize the response, it is useful to use its energy-weighted moments, or sum rules. By comparing computed sum rules with experimental values, the information content of the response can be utilized in the optimization process of the nuclear Hamiltonian or EDF. But the additional information comes at a price: compared to the ground state, computation of excited states is more demanding. To establish an efficient framework to compute sum rules of the response that is adaptable to the optimization of the nuclear EDF and large-scale surveys of collective strength, we have developed a new technique within the complex-energy FAM based on the QRPA. To compute sum rules, we carry out contour integration of the response function in the complex-energy plane. We benchmark our results against the conventional matrix formulation of the QRPA theory, the Thouless theorem for the energy-weighted sum rule, and the dielectric theorem for the inverse energy-weighted sum rule. We demonstrate that calculated sum-rule values agree with those obtained from the matrix formulation of the QRPA. We also discuss the applicability of both the Thouless theorem about the energy-weighted sum rule and the dielectric theorem for the inverse energy-weighted sum rule to nuclear density functional theory in cases when the EDF is not based on a Hamiltonian. The proposed sum-rule technique based on the complex-energy FAM is a tool of choice when optimizing effective interactions or energy functionals. The method is very efficient and well-adaptable to parallel computing. The FAM formulation is especially useful when standard theorems based on commutation relations involving the nuclear Hamiltonian and external field cannot be used.

preprint2015arXiv

Isospin effects in N~Z nuclei in extended Density Functional Theory

This paper overviews various phenomena related to the concept of isospin symmetry. The focus is on N~Z nuclei, which are excellent laboratories of isospin physics. The theoretical framework applied is nuclear Density Functional Theory and its isospin- and angular-momentum projected extensions, as well as symmetry-projected multi-reference models. The topics covered include: isospin impurities, superallowed beta decays, beta-transitions in mirror nuclei, isospin breaking hadronic interactions, mirror and triplet binding energy differences, and isoscalar pairing.

preprint2015arXiv

Microscopic modeling of mass and charge distributions in the spontaneous fission of 240Pu

In this letter, we outline a methodology to calculate microscopically mass and charge distributions of spontaneous fission yields. We combine the multi-dimensional minimization of collective action for fission with stochastic Langevin dynamics to track the relevant fission paths from the ground-state configuration up to scission. The nuclear potential energy and collective inertia governing the tunneling motion are obtained with nuclear density functional theory in the collective space of shape deformations and pairing. We obtain a quantitative agreement with experimental data and find that both the charge and mass distributions in the spontaneous fission of 240Pu are sensitive both to the dissipation in collective motion and to adiabatic characteristics.

preprint2014arXiv

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

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.

preprint2013arXiv

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

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.

preprint2013arXiv

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

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.

preprint2010arXiv

Orbital dependent nucleonic pairing in the lightest known isotopes of tin

By studying the 109Xe-->105Te-->101Sn superallowed alpha-decay chain, we observe low-lying states in 101Sn, the one-neutron system outside doubly magic 100Sn. We find that the spins of the ground state (J = 7=2) and first excited state (J = 5=2) in 101Sn are reversed with respect to the traditional level ordering postulated for 103Sn and the heavier tin isotopes. Through simple arguments and state-of-the-art shell model calculations we explain this unexpected switch in terms of a transition from the single-particle regime to the collective mode in which orbital-dependent pairing correlations, dominate.

preprint2002arXiv

Asymptotic Behavior of the Wave Packet Propagation through a Barrier: the Green's Function Approach Revisited

To model the decay of a quasibound state we use the modified two-potential approach introduced by Gurvitz and Kalbermann. This method has proved itself useful in the past for calculating the decay width and the energy shift of an isolated quasistationary state. We follow the same approach in order to propagate the wave-packet in time with the ultimate goal of extracting the momentum-distribution of emitted particles. The advantage of the method is that it provides the time-dependent wave function in a simple semi-analytic form. We intend to apply this method to the modeling of metastable states for which no direct integration of the time-dependent Schroedinger equation is available today.