Researcher profile

Witold Nazarewicz

Witold Nazarewicz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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