Researcher profile

W. Satula

W. Satula contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
4topics
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

18 published item(s)

preprint2015arXiv

Beta-decay study within multi-reference density functional theory and beyond

Pioneering study of Gamow-Teller (GT) and Fermi matrix elements (MEs) using no-core-configuration-interaction formalism rooted in multi-reference density functional theory is presented. After successful test performed for 6He -> 6Li beta-decay, the model is applied to compute MEs in the sd- and pf-shell T=1/2 mirror nuclei. The calculated GT MEs and the isospin-symmetry-breaking corrections to the Fermi branch are found to be in a very good agreement with shell-model predictions in spite of fundamental differences between these models concerning model space, treatment of correlations or inclusion of a core. This result indirectly supports the two-body current based scenarios behind the quenching of axial-vector coupling constant.

preprint2013arXiv

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

The results of systematic calculations of isospin-symmetry-breaking corrections to superallowed beta-decays based on the self-consistent isospin- and angular-momentum-projected nuclear density functional theory (DFT) are reviewed with an emphasis on theoretical uncertainties of the model. Extensions of the formalism towards no core shell model approach with basis cutoff scheme dictated by the self-consistent particle-hole DFT solutions will be also discussed.

preprint2013arXiv

Isospin mixing within the symmetry restored density functional theory and beyond

We present results of systematic calculations of the isospin-symmetry-breaking corrections to the superallowed I=$0+,T=1 --> I=0+,T=1 beta-decays, based on the self-consistent isospin- and angular-momentum-projected nuclear density functional theory (DFT). We discuss theoretical uncertainties of the formalism related to the basis truncation, parametrization of the underlying energy density functional, and ambiguities related to determination of Slater determinants in odd-odd nuclei. A generalization of the double-projected DFT model towards a no core shell-model-like configuration-mixing approach is formulated and implemented. We also discuss new opportunities in charge-symmetry- and charge-independence-breaking studies offered by the newly developed DFT formalism involving proton-neutron mixing in the particle-hole channel.

preprint2013arXiv

Symmetry energy in nuclear density functional theory

The nuclear symmetry energy represents a response to the neutron-proton asymmetry. In this survey we discuss various aspects of symmetry energy in the framework of nuclear density functional theory, considering both non-relativistic and relativistic self-consistent mean-field realizations side-by-side. Key observables pertaining to bulk nucleonic matter and finite nuclei are reviewed. Constraints on the symmetry energy and correlations between observables and symmetry-energy parameters, using statistical covariance analysis, are investigated. Perspectives for future work are outlined in the context of ongoing experimental efforts.

preprint2012arXiv

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

Background: The superallowed beta-decay rates provide stringent constraints on physics beyond the Standard Model of particle physics. To extract crucial information about the electroweak force, small isospin-breaking corrections to the Fermi matrix element of superallowed transitions must be applied. Purpose: We perform systematic calculations of isospin-breaking corrections to superallowed beta-decays and estimate theoretical uncertainties related to the basis truncation, time-odd polarization effects related to the intrinsic symmetry of the underlying Slater determinants, and to the functional parametrization. Methods: We use the self-consistent isospin- and angular-momentum-projected nuclear density functional theory employing two density functionals derived from the density independent Skyrme interaction. Pairing correlations are ignored. Our framework can simultaneously describe various effects that impact matrix elements of the Fermi decay: symmetry breaking, configuration mixing, and long-range Coulomb polarization. Results: The isospin-breaking corrections to the I=0+,T=1 --> I=0+,T=1 pure Fermi transitions are computed for nuclei from A=10 to A=98 and, for the first time, to the Fermi branch of the I,T=1/2 --> I,T=1/2 transitions in mirror nuclei from A=11 to A=49. We carefully analyze various model assumptions impacting theoretical uncertainties of our calculations and provide theoretical error bars on our predictions. Conclusions: The overall agreement with empirical isospin-breaking corrections is very satisfactory. Using computed isospin-breaking corrections we show that the unitarity of the CKM matrix is satisfied with a precision better than 0.1%.

preprint2011arXiv

Microscopic calculations of isospin mixing in N~Z nuclei and isospin-symmetry-breaking corrections to the superallowed beta-decay

Recently, we have applied for the first time the angular momentum and isospin projected nuclear density functional theory to calculate the isospin-symmetry breaking (ISB) corrections to the superallowed beta-decay. With the calculated set of the ISB corrections we found |V_{ud}|=0.97447(23) for the leading element of the Cabibbo-Kobayashi-Maskawa matrix. This is in nice agreement with both the recent result of Towner and Hardy [Phys. Rev. {\bf C77}, 025501 (2008)] and the central value deduced from the neutron decay. In this work we extend our calculations of the ISB corrections covering all superallowed transitions A,I^π=0^+,T=1,T_z \rightarrow A,I^π=0^+,T=1,T_z+1 with T_z =-1,0 and A ranging from 10 to 74.

preprint2011arXiv

Microscopic calculations of isospin-breaking corrections to superallowed beta-decay

The superallowed beta-decay rates that provide stringent constraints on physics beyond the Standard Model of particle physics are affected by nuclear structure effects through isospin-breaking corrections. The self-consistent isospin- and angular-momentum-projected nuclear density functional theory is used for the first time to compute those corrections for a number of Fermi transitions in nuclei from A=10 to A=74. The resulting leading element of the CKM matrix, |V_{ud}|= 0.97447(23), agrees well with the recent result by Towner and Hardy [Phys. Rev. C {\bf 77}, 025501 (2008)].

preprint2011arXiv

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

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

preprint2010arXiv

Isospin mixing in nuclei around N=Z and the superallowed beta-decay

Theoretical approaches that use one-body densities as dynamical variables, such as Hartree-Fock or the density functional theory (DFT), break isospin symmetry both explicitly, by virtue of charge-dependent interactions, and spontaneously. To restore the spontaneously broken isospin symmetry, we implemented the isospin-projection scheme on top of the Skyrme-DFT approach. This development allows for consistent treatment of isospin mixing in both ground and exited nuclear states. In this study, we apply this method to evaluate the isospin impurities in ground states of even-even and odd-odd N~Z nuclei. By including simultaneous isospin and angular-momentum projection, we compute the isospin-breaking corrections to the 0+ --> 0+ superallowed beta-decay.

preprint2010arXiv

Isospin mixing in the vicinity of the N=Z line

We present the isospin- and angular-momentum-projected nuclear density functional theory (DFT) and its applications to the isospin-breaking corrections to the superallowed beta-decay rates in the vicinity of the N=Z line. A preliminary value obtained for the Cabbibo-Kobayashi-Maskawa matrix element, |V_{ud}|=0.97463(24), agrees well with the recent estimate by Towner and Hardy [Phys. Rev. C{\bf 77}, 025501 (2008)]. We also discuss new opportunities to study the symmetry energy by using the isospin-projected DFT.

preprint2010arXiv

Surface-peaked effective mass in the nuclear energy density functional and its influence on single-particle spectra

Calculations for infinite nuclear matter with realistic nucleon-nucleon interactions suggest that the isoscalar effective mass of a nucleon at the saturation density, m*/m, equals 0.8 +/- 0.1. This result is at variance with empirical data on the level density in finite nuclei, which are consistent with m*/m ~ 1. Ma and Wambach suggested that these two contradicting results may be reconciled within a single theoretical framework by assuming a radial-dependent effective mass, peaked at the nuclear surface. The aim of this exploratory work is to investigate this idea within the density functional theory by using a Skyrme-type local functional enriched with new terms, $τ(\mathbf{\nabla}ρ)^2$ and $τ\frac{dρ}{dr}$, where $τ$ and $ρ$ denote the kinetic and particle densities, respectively. We show that each of these terms can give rise to a surface peak in the effective mass, but of a limited height. We investigate the influence of the radial profile of the effective mass on the spin-orbit splittings and centroids. In particular, we demonstrate that the $τ\frac{dρ}{dr}$ term quenches the 1f5/2-1f7/2 splitting in 40Ca, which is strongly overestimated within conventional Skyrme parametrizations.

preprint2009arXiv

Isospin-symmetry restoration within the nuclear density functional theory: formalism and applications

Isospin symmetry of atomic nuclei is explicitly broken by the charge-dependent interactions, primarily the Coulomb force. Within the nuclear density functional theory, isospin is also broken spontaneously. We propose a projection scheme rooted in a mean field theory, that allows the consistent treatment of isospin breaking in both ground and exited nuclear states. We demonstrate that this scheme is essentially free from spurious divergences plaguing particle-number and angular-momentum restoration approaches. Applications of the new technique include excited high-spin states in medium-mass N=Z nuclei, such as superdeformed bands and many-particle-many-hole terminating states.

preprint2009arXiv

The nuclear energy density functionals with modified radial dependence of the isoscalar effective mass

Calculations for infinite nuclear matter with realistic nucleon-nucleon interactions suggest that the isoscalar effective mass (IEM) of a nucleon at the saturation density equals $m^*/m\sim 0.8\pm 0.1$, at variance with empirical data on the nuclear level density in finite nuclei which are consistent with $m^*/m\approx 1$. This contradicting results might be reconciled by enriching the radial dependence of IEM. In this work four new terms are introduced into the Skyrme-force inspired local energy-density functional: $τ(\nablaρ)^2$, $τ\frac{dρ}{dr}$, $τ^2$ and $τΔρ$. The aim is to investigate how they influence the radial dependence of IEM and, in turn, the single-particle spectra.

preprint2008arXiv

Shell-structure fingerprints of tensor interaction

We address consequences of strong tensor and weak spin-orbit terms in the local energy density functional, resulting from fits to the $f_{5/2} - f_{7/2}$ splittings in $^{40}$Ca, $^{48}$Ca, and $^{56}$Ni. In this study, we focus on nuclear binding energies. In particular, we show that the tensor contribution to the binding energies exhibits interesting topological features closely resembling that of the shell-correction. We demonstrate that in the extreme single-particle scenario at spherical shape, the tensor contribution shows tensorial magic numbers equal to $N(Z)$=14, 32, 56, and 90, and that this structure is smeared out due to configuration mixing caused by pairing correlations and migration of proton/neutron sub-shells with neutron/proton shell filling. Based on a specific Skyrme-type functional SLy4$_T$, we show that the proton tensorial magic numbers shift with increasing neutron excess to $Z$=14, 28, and 50.

preprint2007arXiv

Critical Survey of Isoscalar and Isovector Contributions to the Spin Orbit Potential in Relativistic Mean Field Theory

The spin-orbit (SO) interaction, emerging naturally from the Relativistic Mean Field (RMF) theory is examined critically in the light of the recently measured excitation energy differences between the terminating states built on two different configurations for nuclei belonging to the lower pf shell. The calculations are carried out using the cranked RMF framework. To further probe the iso-vector dependence of the spin-orbit potential, the energy spacing between the g_{7/2} and h_{11/2} states in the Sb-chain is compared to experiment. It is found that the calculation at the quantitative level deviates strongly from the experiment. In particular the balance of the iso-scalar and iso-vector strengths of the effective one body SO potential indicates that additional terms like tensor couplings may be needed to account for the experimental data.

preprint2007arXiv

Exact solution of the spin-isospin proton-neutron pairing Hamiltonian

The exact solution of proton-neutron isoscalar-isovector (T=0,1) pairing Hamiltonian with non-degenerate single-particle orbits and equal pairing strengths (g_{T=1}= g_{T=0}) is presented for the first time. The Hamiltonian is a particular case of a family of integrable SO(8) Richardson-Gaudin (RG) models. The exact solution of the T=0,1 pairing Hamiltonian is reduced to a problem of 4 sets of coupled non linear equations that determine the spectral parameters of the complete set of eigenstates. The microscopic structure of individual eigenstates is analyzed in terms of evolution of the spectral parameters in the complex plane for system of A=80 nucleons. The spectroscopic trends of the exact solutions are discussed in terms of generalized rotations in isospace.

preprint2002arXiv

Comments on the nuclear symmetry energy

According to standard textbooks, the nuclear symmetry energy originates from the {\it kinetic} energy and the {\it interaction} itself. We argue that this view requires certain modifications. We ascribe the physical origin of the {\it kinetic} term to the discreteness of fermionic levels of, in principle arbitrary binary fermionic systems, and relate its mean value directly to the average level density. Physically it connects this part also to the isoscalar part of the interaction which, at least in self-bound systems like atomic nuclei, decides upon the spatial dimensions of the system. For the general case of binary fermionic systems possible external confining potentials as well as specific boundary conditions will contribute to this part. The reliability of this concept is tested using self-consistent Skyrme Hartree-Fock calculations.

preprint2002arXiv

Mass Number Dependence of Nuclear Pairing

Large scale Hartree-Fock-Bogoliubov (HFB) calculations with the finite-range Gogny force D1S have been performed in order to extract the corresponding theoretical average mass dependence of the nuclear gap values. Good agreement with experimental data from the three-point filter $Δ^{(3)}$(N) with N odd has been found for both the neutron and proton gaps. The study confirms earlier findings [W. Satula, J. Dobaczewski, and W. Nazarewicz, Phys. Rev. Lett. 81 3599 (1998)] that the mass dependence of the gap is much weaker than the so far accepted 12/$\sqrtA$ MeV law.