Source author record

B. G. Carlsson

B. G. Carlsson 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

23works
4topics
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

23 published item(s)

preprint2022arXiv

Nuclear spectra from low-energy interactions

A method to describe spectra starting from nuclear density functionals is explored. The idea is based on postulating an effective Hamiltonian that reproduces the stiffness associated with collective modes. The method defines a simple form of such an effective Hamiltonian and a mapping to go from a density functional to the corresponding Hamiltonian. In order to test the method, the Hamiltonian is constrained using a Skyrme functional and solved with the generator-coordinate method to describe low-lying levels and electromagnetic transitions in $^{48,49,50,52}$Cr and $^{24}$Mg.

preprint2020arXiv

Calculated fission-fragment mass yields and average total kinetic energies of heavy and superheavy nuclei

Fission-fragment mass and total-kinetic-energy (TKE) distributions following fission of even-even nuclides in the region $74 \leq Z \leq 126$ and $92 \leq N \leq 230$, comprising 896 nuclides have been calculated using the Brownian shape-motion method. The emphasis is the region of superheavy nuclei. To show compatibility with earlier results the calculations are extended to include earlier studied regions. An island of asymmetric fission is obtained in the superheavy region, $106\leq Z\leq114$ and $162\leq N\leq 176$, where the heavy fragment is found to be close to $^{208}$Pb and the light fragment adjusts accordingly. Most experimentally observed $α$-decay chains of superheavy nuclei with $Z > 113 $ terminate by spontaneous fission in our predicted region of asymmetric fission. In these cases, the pronounced large asymmetry is accompanied by a low TKE value compatible with measurements.

preprint2020arXiv

Correlation studies of fission fragment neutron multiplicities

We calculate neutron multiplicities from fission fragments with specified mass numbers for events having a specified total fragment kinetic energy. The shape evolution from the initial compound nucleus to the scission configurations is obtained with the Metropolis walk method on the five-dimensional potential-energy landscape, calculated with the macroscopic-microscopic method for the three-quadratic-surface shape family. Shape-dependent microscopic level densities are used to guide the random walk, to partition the intrinsic excitation energy between the two proto-fragments at scission, and to determine the spectrum of the neutrons evaporated from the fragments. The contributions to the total excitation energy of the resulting fragments from statistical excitation and shape distortion at scission is studied. Good agreement is obtained with available experimental data on neutron multiplicities in correlation with fission fragments from $^{235}$U(n$_{\rm th}$,f). At higher neutron energies a superlong fission mode appears which affects the dependence of the observables on the total fragment kinetic energy.

preprint2016arXiv

The effect of realistic nuclear charge distributions on isotope shifts and towards the extraction of higher order nuclear radial moments

Background: Atomic spectral lines from different isotopes display a small shift in energy, commonly referred to as the line isotope shift. One of the components of the isotope shift is the field shift, which depends on the extent and the shape of the nuclear charge density distribution. Purpose: To investigate how sensitive field shifts are with respect to variations in the nuclear size and shape and what information of nuclear charge distributions that can be extracted from measured field shifts. Methods: Nuclear properties are obtained from nuclear density functional theory calculations based on the Skyrme-Hartree-Fock-Bogoliubov approach. These results are combined with multiconfiguration Dirac-Hartree-Fock methods to obtain realistic field shifts. Results: Phenomena such as nuclear deformation and variations in the diffuseness of nuclear charge distributions give measurable contributions to the field shifts. Using a novel approach, we demonstrate the possibility to extract new information concerning the nuclear charge densities from the observed field shifts. Conclusions: Combining methods used in atomic and nuclear structure theory gives an improved description of field shifts. Extracting additional nuclear information from measured field shifts is possible in the near future with improved experimental methods.

preprint2015arXiv

Recoil-alpha-fission and recoil-alpha-alpha-fission events observed in the reaction Ca-48 + Am-243

Products of the fusion-evaporation reaction Ca-48 + Am-243 were studied with the TASISpec set-up at the gas-filled separator TASCA at the GSI Helmholtzzentrum für Schwerionenforschung. Amongst the detected thirty correlated alpha-decay chains associated with the production of element Z=115, two recoil-alpha-fission and five recoil-alpha-alpha-fission events were observed. The latter are similar to four such events reported from experiments performed at the Dubna gas-filled separator. Contrary to their interpretation, we propose an alternative view, namely to assign eight of these eleven decay chains of recoil-alpha(-alpha)-fission type to start from the 3n-evaporation channel 115-288. The other three decay chains remain viable candidates for the 2n-evaporation channel 115-289.

preprint2014arXiv

$α$-decay calculations of heavy nuclei using an effective Skyrme interaction

Background: For nuclei heavier than $^{208}$Pb $α$ decay is a dominating decay mode, and in the search of new superheavy elements one often observes chains of $α$ decays. Purpose: Explore and test microscopic descriptions of $α$ decay based on theories with effective nuclear interactions. Methods: The nuclear ground states are calculated with the Hartree-Fock-Bogoliubov (HFB) method using the Skyrme interaction. Microscopic $α$-decay formation amplitudes are calculated from the HFB wave functions, and the $R$-matrix formalism is utilized to obtain decay probabilities. Results: Using a large harmonic-oscillator basis we obtain converged $α$-decay widths. A comparison with experiment including all spherical even-even $α$ emitting nuclei shows that the model consistently predicts too small formation amplitudes while relative values are in good agreement with experiment. Conclusions: The method was found to be numerically practical even with a large basis size. The comparison of formation amplitudes suggests that the pairing type correlations included in the HFB approach cannot produce sufficient $α$-particle clustering.

preprint2014arXiv

Alpha-particle formation and decay rates from Skyrme-HFB wave functions

$α$ decay is treated microscopically, where the unstable mother nucleus and residual daughter nucleus are described using HFB wave functions, obtained with the Skyrme effective interaction. From these wave functions the amplitude for forming $α$ particles in the mother nucleus is computed. Two different Skyrme parametrizations with different pairing properties are compared, and we find good agreement with experiment for relative decay rates in both cases. The absolute values of the decay rates are under-estimated.

preprint2014arXiv

Low-lying states in near-magic odd-odd nuclei and the effective interaction

The iterative quasi-particle-random-phase approximation (QRPA) method we previously developed to accurately calculate properties of individual nuclear states is extended so that it can be applied for nuclei with odd numbers of neutrons and protons. The approach is based on the proton-neutron-QRPA (pnQRPA) and uses an iterative non-hermitian Arnoldi diagonalization method where the QRPA matrix does not have to be explicitly calculated and stored. The method is used to calculate excitation energies of proton-neutron multiplets for several nuclei. The influence of a pairing interaction in the $T=0$ channel is studied.

preprint2014arXiv

Spectroscopic properties of nuclear Skyrme energy density functionals

We address the question of how to improve the agreement between theoretical nuclear single-particle energies (SPEs) and experiment. Empirically, in doubly magic nuclei, the SPEs can be deduced from spectroscopic properties of odd nuclei that have one more, or one less neutron or proton. Theoretically, bare SPEs, before being confronted with experiment, must be corrected for the effects of the particle-vibration-coupling (PVC). In the present work, we determine the PVC corrections in a fully self-consistent way. Then, we adjust the SPEs, with PVC corrections included, to empirical data. In this way, the agreement with experiment, on average, improves; nevertheless, large discrepancies still remain. We conclude that the main source of disagreement is still in the underlying mean fields, and not in including or neglecting the PVC corrections.

preprint2014arXiv

Structure of Superheavy Nuclei Along Element 115 Decay Chains

A recent high-resolution $α$, $X$-ray, and $γ$-ray coincidence-spectroscopy experiment offered first glimpse of excitation schemes of isotopes along $α$-decay chains of $Z=115$. To understand these observations and to make predictions about shell structure of superheavy nuclei below $^{288}115$, we employ two complementary mean-field models: self-consistent Skyrme Energy Density Functional approach and the macroscopic-microscopic Nilsson model. We discuss the spectroscopic information carried by the new data. In particular, candidates for the experimentally observed $E1$ transitions in $^{276}$Mt are proposed. We find that the presence and nature of low-energy $E1$ transitions in well-deformed nuclei around $Z=110, N=168$ strongly depends on the strength of the spin-orbit coupling; hence, it provides an excellent constraint on theoretical models of superheavy nuclei. To clarify competing theoretical scenarios, an experimental search for $E1$ transitions in odd-$A$ systems $^{275,277}$Mt, $^{275}$Hs, and $^{277}$Ds is strongly recommended.

preprint2013arXiv

Polarization corrections to single-particle energies studied within the energy-density-functional and QRPA approaches

Background: Models based on using perturbative polarization corrections and mean-field blocking approximation give conflicting results for masses of odd nuclei. Purpose: Systematically investigate the polarization and mean-field models, implemented within self-consistent approaches that use identical interactions and model spaces, so as to find reasons for the conflicts between them. Methods: For density-dependent interactions and with pairing correlations included, we derive and study links between the mean-field and polarization results obtained for energies of odd nuclei. We also identify and discuss differences between the polarization-correction and full particle-vibration-coupling (PVC) models. Numerical calculations are performed for the mean-field ground-state properties of deformed odd nuclei and then compared to the polarization corrections determined by using the approach that conserves spherical symmetry. Results: We have identified and numerically evaluated self-interaction (SI) energies that are at the origin of different results obtained within the mean-field and polarization-correction approaches. Conclusions: Mean-field energies of odd nuclei are polluted by the SI energies, and this makes them different from those obtained by using polarization-correction methods. A comparison of both approaches allows for the identification and determination of the SI terms, which then can be calculated and removed from the mean-field results, giving the self-interaction-free energies. The simplest deformed mean-field approach that does not break parity symmetry is unable to reproduce full PVC effects.

preprint2012arXiv

Collective vibrational states with fast iterative QRPA method

An iterative method we previously proposed to compute nuclear strength functions is developed to allow it to accurately calculate properties of individual nuclear states. The approach is based on the quasi-particle-random-phase approximation (QRPA) and uses an iterative non-hermitian Arnoldi diagonalization method where the QRPA matrix does not have to be explicitly calculated and stored. The method gives substantial advantages over conventional QRPA calculations with regards to the computational cost. The method is used to calculate excitation energies and decay rates of the lowest lying 2+ and 3- states in Pb, Sn, Ni and Ca isotopes using three different Skyrme interactions and a separable gaussian pairing force.

preprint2012arXiv

Fluctuating parts of nuclear ground state correlation energies

Background: Heavy atomic nuclei are often described using the Hartree-Fock-Bogoliubov (HFB) method. In principle, this approach takes into account Pauli effects and pairing correlations while other correlation effects are mimicked through the use of effective density-dependent interactions. Purpose: Investigate the influence of higher order correlation effects on nuclear binding energies using Skyrme's effective interaction. Methods: A cut-off in relative momenta is introduced in order to remove ultraviolet divergences caused by the zero-range character of the interaction. Corrections to binding energies are then calculated using the quasiparticle-random-phase approximation (QRPA) and second order many-body perturbation theory (MBPT2). Result: Contributions to the correlation energies are evaluated for several isotopic chains and an attempt is made to disentangle which parts give rise to fluctuations that may be difficult to incorporate on the HFB level. The dependence of the results on the cut-off is also investigated. Conclusions: The improved interaction allows explicit summations of perturbation series which is useful for the description of some nuclear observables. However, refits of the interaction parameters are needed to obtain more quantitative results.

preprint2012arXiv

Giant Monopole Resonances and nuclear incompressibilities studied for the zero-range and separable pairing interactions

Background: Following the 2007 precise measurements of monopole strengths in tin isotopes, there has been a continuous theoretical effort to obtain a precise description of the experimental results. Up to now, there is no satisfactory explanation of why the tin nuclei appear to be significantly softer than 208Pb. Purpose: We determine the influence of finite-range and separable pairing interactions on monopole strength functions in semi-magic nuclei. Methods: We employ self-consistently the Quasiparticle Random Phase Approximation on top of spherical Hartree-Fock-Bogolyubov solutions. We use the Arnoldi method to solve the linear-response problem with pairing. Results: We found that the difference between centroids of Giant Monopole Resonances measured in lead and tin (about 1 MeV) always turns out to be overestimated by about 100%. We also found that the volume incompressibility, obtained by adjusting the liquid-drop expression to microscopic results, is significantly larger than the infinite-matter incompressibility. Conclusions: The zero-range and separable pairing forces cannot induce modifications of monopole strength functions in tin to match experimental data.

preprint2011arXiv

Continuity equation and local gauge invariance for the N3LO nuclear Energy Density Functionals

Background: The next-to-next-to-next-to-leading order (N3LO) nuclear energy density functional extends the standard Skyrme functional with new terms depending on higher-order derivatives of densities, introduced to gain better precision in the nuclear many-body calculations. A thorough study of the transformation properties of the functional with respect to different symmetries is required, as a step preliminary to the adjustment of the coupling constants. Purpose: Determine to which extent the presence of higher-order derivatives in the functional can be compatible with the continuity equation. In particular, to study the relations between the validity of the continuity equation and invariance of the functional under gauge transformations. Methods: Derive conditions for the validity of the continuity equation in the framework of time-dependent density functional theory. The conditions apply separately to the four spin-isospin channels of the one-body density matrix. Results: We obtained four sets of constraints on the coupling constants of the N3LO energy density functional that guarantee the validity of the continuity equation in all spin-isospin channels. In particular, for the scalar-isoscalar channel, the constraints are the same as those resulting from imposing the standard U(1) local-gauge-invariance conditions. Conclusions: Validity of the continuity equation in the four spin-isospin channels is equivalent to the local-gauge invariance of the energy density functional. For vector and isovector channels, such validity requires the invariance of the functional under local rotations in the spin and isospin spaces.

preprint2011arXiv

Effective pseudopotential for energy density functionals with higher order derivatives

We derive a zero-range pseudopotential that includes all possible terms up to sixth order in derivatives. Within the Hartree-Fock approximation, it gives the average energy that corresponds to a quasi-local nuclear Energy Density Functional (EDF) built of derivatives of the one-body density matrix up to sixth order. The direct reference of the EDF to the pseudopotential acts as a constraint that divides the number of independent coupling constants of the EDF by two. This allows, e.g., for expressing the isovector part of the functional in terms of the isoscalar part, or vice versa. We also derive the analogous set of constraints for the coupling constants of the EDF that is restricted by spherical, space-inversion, and time-reversal symmetries.

preprint2011arXiv

Evidence for Shape Co-existence at medium spin in 76Rb

Four previously known rotational bands in 76Rb have been extended to moderate spins using the Gammasphere and Microball gamma ray and charged particle detector arrays and the 40Ca(40Ca,3pn) reaction at a beam energy of 165 MeV. The properties of two of the negative-parity bands can only readily be interpreted in terms of the highly successful Cranked Nilsson-Strutinsky model calculations if they have the same configuration in terms of the number of g9/2 particles, but they result from different nuclear shapes (one near-oblate and the other near-prolate). These data appear to constitute a unique example of shape co-existing structures at medium spins.

preprint2010arXiv

Convergence of density-matrix expansions for nuclear interactions

We extend density-matrix expansions in nuclei to higher orders in derivatives of densities and test their convergence properties. The expansions allow for converting the interaction energies characteristic to finite- and short-range nuclear effective forces into quasi-local density functionals. We also propose a new type of expansion that has excellent convergence properties when benchmarked against the binding energies obtained for the Gogny interaction.

preprint2010arXiv

The Negele-Vautherin density matrix expansion applied to the Gogny force

We use the Negele-Vautherin density matrix expansion to derive a quasi-local density functional for the description of systems of fermions interacting with short-ranged interactions composed of arbitrary finite-range central, spin-orbit, and tensor components. Terms that are absent in the original Negele-Vautherin approach owing to the angle averaging of the density matrix are fixed by employing a gauge invariance condition. We obtain the Kohn-Sham interaction energies in all spin-isospin channels, including the exchange terms, expressed as functions of the local densities and their derivatives up to second (next to leading) order. We illustrate the method by determining the coupling constants of the Skyrme functional or Skyrme force that correspond to the finite-range Gogny central force. The resulting self-consistent solutions reproduce the Gogny-force binding energies and radii within the precision of 1-2%.

preprint2009arXiv

Linear response strength functions with iterative Arnoldi diagonalization

We report on an implementation of a new method to calculate RPA strength functions with iterative non-hermitian Arnoldi diagonalization method, which does not explicitly calculate and store the RPA matrix. We discuss the treatment of spurious modes, numerical stability, and how the method scales as the used model space is enlarged. We perform the particle-hole RPA benchmark calculations for double magic nucleus 132Sn and compare the resulting electromagnetic strength functions against those obtained within the standard RPA.

preprint2009arXiv

Solution of self-consistent equations for the N3LO nuclear energy density functional in spherical symmetry. The program HOSPHE (v1.00)

We present solution of self-consistent equations for the N3LO nuclear energy density functional. We derive general expressions for the mean fields expressed as differential operators depending on densities and for the densities expressed in terms of derivatives of wave functions. These expressions are then specified to the case of spherical symmetry. We also present the computer program HOSPHE (v1.00), which solves the self-consistent equations by using the expansion of single-particle wave functions on the spherical harmonic oscillator basis.

preprint2008arXiv

Local nuclear energy density functional at next-to-next-to-next-to-leading order

We construct nuclear energy density functionals in terms of derivatives of densities up to sixth, next-to-next-to-next-to-leading order (N3LO). A phenomenological functional built in this way conforms to the ideas of the density matrix expansion and is rooted in the expansions characteristic to effective theories. It builds on the standard functionals related to the contact and Skyrme forces, which constitute the zero-order (LO) and second-order (NLO) expansions, respectively. At N3LO, the full functional with density-independent coupling constants, and with the isospin degree of freedom taken into account, contains 376 terms, while the functionals restricted by the Galilean and gauge symmetries contain 100 and 42 terms, respectively. For functionals additionally restricted by the spherical, space-inversion, and time-reversal symmetries, the corresponding numbers of terms are equal to 100, 60, and 22, respectively.