Source author record

M. Kortelainen

M. Kortelainen 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
6topics
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)

preprint2021arXiv

Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni

Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).

preprint2019arXiv

Regularized pseudopotential for mean-field calculations

We present preliminary results obtained with a finite-range two-body pseudopotential complemented with zero-range spin-orbit and density-dependent terms. After discussing the penalty function used to adjust parameters, we discuss predictions for binding energies of spherical nuclei calculated at the mean-field level, and we compare them with those obtained using the standard Gogny D1S finite-range effective interaction.

preprint2018arXiv

Towards a novel energy density functional for beyond-mean-field calculations with pairing and deformation

We take an additional step towards the optimization of the novel finite-range pseudopotential at constrained Hartree-Fock-Bogolyubov level and implement an optimization procedure within an axial code using harmonic oscillator basis. We perform the optimization using three different numbers of the harmonic oscillator shells. We apply the new parameterizations in the O-Kr part of the nuclear chart and isotopic chain of Sn, and we compare the results with experimental values and those given by a parameterization obtained using a spherical code.

preprint2016arXiv

Fayans functional for deformed nuclei. Uranium region

Fayans energy density functional (EDF) FaNDF^0 has been applied to the nuclei around uranium region. Ground state characteristics of the Th, U and Pu isotopic chains, up to the two-neutron drip line, are found and compared with predictions from several Skyrme EDFs. The two-neutron drip line is found for FaNDF^0, SLy4 and SkM^* EDFs for a set of elements with even proton number, from Pb up to Fm.

preprint2016arXiv

Inelastic WIMP-nucleus scattering to the first excited state in $^{125}$Te

The direct detection of dark matter constituents, in particular the weakly interacting massive particles (WIMPs), is considered central to particle physics and cosmology. In this paper we study transitions to the excited states, possible in some nuclei, which have sufficiently low lying excited states. Examples considered previously were the first excited states of $^{127}$I and $^{129}$Xe and $^{83}$Kr. Here we examine $^{125}$Te, which offers some advantages and is currently being considered as a target.In all these cases the extra signature of the gamma rays following the de-excitation of these states has definite advantages over the purely nuclear recoil and, in principle, such a signature can be exploited experimentally. A brief discussion of the experimental feasibility is given in the context of the CUORE experiment.

preprint2015arXiv

First applications of Fayans functional to deformed nuclei

First calculations for deformed nuclei with the Fayans functional are carried out for the uranium and lead isotopic chains. The ground state deformations and deformation energies are compared to Skyrme-Hartree-Fock-Bogolyubov results of HFB-17 and HFB-27 functionals. For the uranium isotopic chain, the Fayans functional predictions are rather similar properties compared to HFB-17 and HFB-27. However, there is a disagreement for the lead isotopic chain. Both of the Skyrme HFB functionals predict rather strong deformations for the light Pb isotopes which does not agree with the experimental data on charge radii and magnetic moments of the odd Pb isotopes. On the other hand, the Fayans functional predicts a spherical ground state for all of the lead isotopes, in accordance with the data and the known in literature results obtained with the Gogny D1S force and SLy6 functional as well. The deformation energy curves are calculated and compared to four Skyrme functionals, SLy4, Sly6, SkM* and UNEDF1, for $^{238}$U nucleus and several lead deficient Pb isotopes. In the first case, the Fayans functional result is rather close to SkM* and UNEDF1 which, in particularly the latter one, describe the first and second barriers in $^{238}$U rather well. For the light lead isotopes, the Fayans deformation energy curves are qualitatively close to those of the SLy6 functional.

preprint2015arXiv

Multipole modes in deformed nuclei within the finite amplitude method

Background: To access selected excited states of nuclei, within the framework of nuclear density functional theory, the quasiparticle random phase approximation (QRPA) is commonly used. Purpose: We present a computationally efficient, fully self-consistent framework to compute the QRPA transition strength function of an arbitrary multipole operator in axially-deformed superfluid nuclei. Methods: The method is based on the finite amplitude method (FAM) QRPA, allowing fast iterative solution of QRPA equations. A numerical implementation of the FAM-QRPA solver module has been carried out for deformed nuclei. Results: The practical feasibility of the deformed FAM module has been demonstrated. In particular, we calculate the quadrupole and octupole strength in a heavy deformed nucleus $^{240}$Pu, without any truncations in the quasiparticle space. To demonstrate the capability to calculate individual QRPA modes, we also compute low-lying negative-parity collective states in $^{154}$Sm. Conclusions: The new FAM implementation enables calculations of the QRPA strength function throughout the nuclear landscape. This will facilitate global surveys of multipole modes and beta decays, and will open new avenues for constraining the nuclear energy density functional.

preprint2015arXiv

Theoretical direct WIMP detection rates for transitions to nuclear excited states

The recent WMAP and Planck data have confirmed that exotic dark matter together with the vacuum energy (cosmological constant) dominate in the flat Universe. Many extensions of the standard model provide dark matter candidates, in particular Weakly Interacting Massive Particles (WIMPs). Thus the direct dark matter detection is central to particle physics and cosmology. Most of the research on this issue has hitherto focused on the detection of the recoiling nucleus. In this paper we study transitions to the excited states, possible in some nuclei, which have sufficiently low lying excited states. Examples considered previously were the first excited states of $^{127}$I and $^{129}$Xe. We examine here $^{83}$Kr, which offers some kinematical advantages and is currently considered as a possible target. We find appreciable branching ratios for the inelastic scattering mediated by the spin cross sections, with an inelastic event rate of $4.4\times 10^{-4}$kg$^{-1}$d$^{-1}$. So, the extra signature of the gamma ray following the de-excitation of these states can, in principle, be exploited experimentally. A brief discussion of the experimental feasibility is given

preprint2014arXiv

Emergent Soft Monopole Modes in Weakly-Bound Deformed Nuclei

Based on the Hartree-Fock-Bogoliubov solutions in large deformed coordinate spaces, the finite amplitude method for quasiparticle random phase approximation (FAM-QRPA) has been implemented, providing a suitable approach to probe collective excitations of weakly-bound nuclei embedded in the continuum. The monopole excitation modes in Magnesium isotopes up to the neutron drip line have been studied with the FAM-QRPA framework on both the coordinate-space and harmonic oscillator basis methods. Enhanced soft monopole strengths and collectivity as a result of weak-binding effects have been unambiguously demonstrated.

preprint2014arXiv

Lipkin method of particle-number restoration to higher orders

Background: On the mean-field level, pairing correlations are incorporated through the Bogoliubov-Valatin transformation, whereupon the particle degrees of freedom are replaced by quasiparticles. This approach leads to a spontaneous breaking of the particle-number symmetry and mixing of states with different particle numbers. In order to restore the particle number, various methods have been employed, which are based on projection approaches before or after variation. Approximate variation-after-projection (VAP) schemes, utilizing the Lipkin method, have mostly been used within the Lipkin-Nogami prescription. Purpose: Without recurring to the Lipkin-Nogami prescription, and using instead states rotated in the gauge space, we derive the Lipkin method of particle-number restoration up to sixth order and we test the convergence and accuracy of the obtained expansion. Methods: We perform self-consistent calculations using the higher-order Lipkin method to restore the particle-number symmetry in the framework of superfluid nuclear energy-density functional theory. We also apply the Lipkin method to a schematic exactly solvable two-level pairing model. Results: Calculations performed in open-shell tin and lead isotopes show that the Lipkin method converges at fourth order and satisfactorily reproduces the VAP ground-state energies and energy kernels. Near closed shells, the higher-order Lipkin method cannot be applied because of a non-analytic kink in the ground-state energies in function of the particle number. Conclusions: In open-shell nuclei, the higher-order Lipkin method provides a good approximation to the exact VAP energies. The method is computationally inexpensive, making it particularly suitable, for example, for future optimizations of the nuclear energy-density functionals and simultaneous restoration of different symmetries.

preprint2014arXiv

Nuclear energy density optimization: Shell structure

Nuclear density functional theory is the only microscopical theory that can be applied throughout the entire nuclear landscape. Its key ingredient is the energy density functional. In this work, we propose a new parameterization UNEDF2 of the Skyrme energy density functional. The functional optimization is carried out using the POUNDerS optimization algorithm within the framework of the Skyrme Hartree-Fock-Bogoliubov theory. Compared to the previous parameterization UNEDF1, restrictions on the tensor term of the energy density have been lifted, yielding a very general form of the energy density functional up to second order in derivatives of the one-body density matrix. In order to impose constraints on all the parameters of the functional, selected data on single-particle splittings in spherical doubly-magic nuclei have been included into the experimental dataset. The agreement with both bulk and spectroscopic nuclear properties achieved by the resulting UNEDF2 parameterization is comparable with UNEDF1. While there is a small improvement on single-particle spectra and binding energies of closed shell nuclei, the reproduction of fission barriers and fission isomer excitation energies has degraded. As compared to previous UNEDF parameterizations, the parameter confidence interval for UNEDF2 is narrower. In particular, our results overlap well with those obtained in previous systematic studies of the spin-orbit and tensor terms. UNEDF2 can be viewed as an all-around Skyrme EDF that performs reasonably well for both global nuclear properties and shell structure. However, after adding new data aiming to better constrain the nuclear functional, its quality has improved only marginally. These results suggest that the standard Skyrme energy density has reached its limits and significant changes to the form of the functional are needed.

preprint2014arXiv

Nuclear Energy Density Optimization: UNEDF2

The parameters of the UNEDF2 nuclear energy density functional (EDF) model were obtained in an optimization to experimental data consisting of nuclear binding energies, proton radii, odd-even mass staggering data, fission-isomer excitation energies, and single particle energies. In addition to parameter optimization, sensitivity analysis was done to obtain parameter uncertainties and correlations. The resulting UNEDF2 is an all-around EDF. However, the sensitivity analysis also demonstrated that the limits of current Skyrme-like EDFs have been reached and that novel approaches are called for.

preprint2013arXiv

Axially deformed solution of the Skyrme-Hartree-Fock-Bogolyubov equations using the transformed harmonic oscillator basis (II) HFBTHO v2.00d: a new version of the program

We describe the new version 2.00d of the code HFBTHO that solves the nuclear Skyrme Hartree-Fock (HF) or Skyrme Hartree-Fock-Bogolyubov (HFB) problem by using the cylindrical transformed deformed harmonic-oscillator basis. In the new version, we have implemented the following features: (i) the modified Broyden method for non-linear problems, (ii) optional breaking of reflection symmetry, (iii) calculation of axial multipole moments, (iv) finite temperature formalism for the HFB method, (v) linear constraint method based on the approximation of the Random Phase Approximation (RPA) matrix for multi-constraint calculations, (vi) blocking of quasi-particles in the Equal Filling Approximation (EFA), (vii) framework for generalized energy density with arbitrary density-dependences, and (viii) shared memory parallelism via OpenMP pragmas.

preprint2013arXiv

Neutron skin uncertainties of Skyrme energy density functionals

Background: Neutron-skin thickness is an excellent indicator of isovector properties of atomic nuclei. As such, it correlates strongly with observables in finite nuclei that depend on neutron-to-proton imbalance and the nuclear symmetry energy that characterizes the equation of state of neutron-rich matter. A rich worldwide experimental program involving studies with rare isotopes, parity violating electron scattering, and astronomical observations is devoted to pinning down the isovector sector of nuclear models. Purpose: We assess the theoretical systematic and statistical uncertainties of neutron-skin thickness and relate them to the equation of state of nuclear matter, and in particular to nuclear symmetry energy parameters. Methods: We use the nuclear superfluid Density Functional Theory with several Skyrme energy density functionals and density dependent pairing. To evaluate statistical errors and their budget, we employ the statistical covariance technique. Results: We find that the errors on neutron skin increase with neutron excess. Statistical errors due to uncertain coupling constants of the density functional are found to be larger than systematic errors, the latter not exceeding 0.06 fm in most neutron-rich nuclei across the nuclear landscape. The single major source of uncertainty is the poorly determined slope L of the symmetry energy that parametrizes its density dependence. Conclusions: To provide essential constraints on the symmetry energy of the nuclear energy density functional, next-generation measurements of neutron skins are required to deliver precision better than 0.06 fm.

preprint2013arXiv

Propagation of uncertainties in the Skyrme energy-density-functional model

Parameters of nuclear energy-density-functionals (EDFs) are always derived by an optimization to experimental data. For the minima of appropriately defined penalty functions, a statistical sensitivity analysis provides the uncertainties of the EDF parameters. To quantify theoretical errors of observables given by the model, we studied the propagation of uncertainties within the UNEDF0 Skyrme-EDF approach. We found that typically the standard errors rapidly increase towards neutron rich nuclei. This can be linked to large uncertainties of the isovector coupling constants of the currently used EDFs.

preprint2012arXiv

Nuclear energy density optimization: Large deformations

A new Skyrme-like energy density suitable for studies of strongly elongated nuclei has been determined in the framework of the Hartree-Fock-Bogoliubov theory using the recently developed model-based, derivative-free optimization algorithm POUNDerS. A sensitivity analysis at the optimal solution has revealed the importance of states at large deformations in driving the parameterization of the functional. The good agreement with experimental data on masses and separation energies, achieved with the previous parameterization UNEDF0, is largely preserved. In addition, the new energy density UNEDF1 gives a much improved description of the fission barriers in 240Pu and neighboring nuclei.

preprint2012arXiv

Precision Mass Measurements beyond $^{132}$Sn: Anomalous behaviour of odd-even staggering of binding energies

Atomic masses of the neutron-rich isotopes $^{121-128}$Cd, $^{129,131}$In, $^{130-135}$Sn, $^{131-136}$Sb, and $^{132-140}$Te have been measured with high precision (10 ppb) using the Penning trap mass spectrometer JYFLTRAP. Among these, the masses of four r-process nuclei $^{135}$Sn, $^{136}$Sb, and $^{139,140}$Te were measured for the first time. The data reveals a strong $N$=82 shell gap at $Z$=50 but indicates the importance of correlations for $Z>50$. An empirical neutron pairing gap expressed as the odd-even staggering of isotopic masses shows a strong quenching across $N$=82 for Sn, with the $Z$-dependence that is unexplainable by the current theoretical models.

preprint2012arXiv

UNEDF: Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

The demands of cutting-edge science are driving the need for larger and faster computing resources. With the rapidly growing scale of computing systems and the prospect of technologically disruptive architectures to meet these needs, scientists face the challenge of effectively using complex computational resources to advance scientific discovery. Multidisciplinary collaborating networks of researchers with diverse scientific backgrounds are needed to address these complex challenges. The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper describes UNEDF and identifies attributes that classify it as a successful computational collaboration. We illustrate significant milestones accomplished by UNEDF through integrative solutions using the most reliable theoretical approaches, most advanced algorithms, and leadership-class computational resources.

preprint2011arXiv

Computing Heavy Elements

Reliable calculations of the structure of heavy elements are crucial to address fundamental science questions such as the origin of the elements in the universe. Applications relevant for energy production, medicine, or national security also rely on theoretical predictions of basic properties of atomic nuclei. Heavy elements are best described within the nuclear density functional theory (DFT) and its various extensions. While relatively mature, DFT has never been implemented in its full power, as it relies on a very large number (~ 10^9-10^12) of expensive calculations (~ day). The advent of leadership-class computers, as well as dedicated large-scale collaborative efforts such as the SciDAC 2 UNEDF project, have dramatically changed the field. This article gives an overview of the various computational challenges related to the nuclear DFT, as well as some of the recent achievements.

preprint2011arXiv

Multipole strength function of deformed superfluid nuclei made easy

We present an efficient method for calculating strength functions using the finite amplitude method (FAM) for deformed superfluid heavy nuclei within the framework of the nuclear density functional theory. We demonstrate that FAM reproduces strength functions obtained with the fully self-consistent quasi-particle random-phase approximation (QRPA) at a fraction of computational cost. As a demonstration, we compute the isoscalar and isovector monopole strength for strongly deformed configurations in $^{240}$Pu by considering huge quasi-particle QRPA spaces. Our approach to FAM, based on Broyden's iterative procedure, opens the possibility for large-scale calculations of strength distributions in well-bound and weakly bound nuclei across the nuclear landscape.

preprint2011arXiv

Testing the density matrix expansion against ab initio calculations of trapped neutron drops

Microscopic input to a universal nuclear energy density functional can be provided through the density matrix expansion (DME), which has recently been revived and improved. Several DME implementation strategies are tested for neutron drop systems in harmonic traps by comparing to Hartree-Fock (HF) and ab initio no-core full configuration (NCFC) calculations with a model interaction (Minnesota potential). The new DME with exact treatment of Hartree contributions is found to best reproduce HF results and supplementing the functional with fit Skyrme-like contact terms shows systematic improvement toward the full NCFC results.

preprint2011arXiv

UNEDF: Advanced Scientific Computing Transforms the Low-Energy Nuclear Many-Body Problem

The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper illustrates significant milestones accomplished by UNEDF through integration of the theoretical approaches, advanced numerical algorithms, and leadership class computational resources.

preprint2010arXiv

Instabilities in the Nuclear Energy Density Functional

In the field of Energy Density Functionals (EDF) used in nuclear structure and dynamics, one of the unsolved issues is the stability of the functional. Numerical issues aside, some EDFs are unstable with respect to particular perturbations of the nuclear ground-state density. The aim of this contribution is to raise questions about the origin and nature of these instabilities, the techniques used to diagnose and prevent them, and the domain of density functions in which one should expect a nuclear EDF to be stable.

preprint2010arXiv

Microscopically-based energy density functionals for nuclei using the density matrix expansion: Implementation and pre-optimization

In a recent series of papers, Gebremariam, Bogner, and Duguet derived a microscopically based nuclear energy density functional by applying the Density Matrix Expansion (DME) to the Hartree-Fock energy obtained from chiral effective field theory (EFT) two- and three-nucleon interactions. Due to the structure of the chiral interactions, each coupling in the DME functional is given as the sum of a coupling constant arising from zero-range contact interactions and a coupling function of the density arising from the finite-range pion exchanges. Since the contact contributions have essentially the same structure as those entering empirical Skyrme functionals, a microscopically guided Skyrme phenomenology has been suggested in which the contact terms in the DME functional are released for optimization to finite-density observables to capture short-range correlation energy contributions from beyond Hartree-Fock. The present paper is the first attempt to assess the ability of the newly suggested DME functional, which has a much richer set of density dependencies than traditional Skyrme functionals, to generate sensible and stable results for nuclear applications. The results of the first proof-of-principle calculations are given, and numerous practical issues related to the implementation of the new functional in existing Skyrme codes are discussed. Using a restricted singular value decomposition (SVD) optimization procedure, it is found that the new DME functional gives numerically stable results and exhibits a small but systematic reduction of our test $χ^2$ function compared to standard Skyrme functionals, thus justifying its suitability for future global optimizations and large-scale calculations.

preprint2010arXiv

Natural Units For Nuclear Energy Density Functional Theory

Naive dimensional analysis based on chiral effective theory, when adapted to nuclear energy density functionals, prescribes natural units and a hierarchy of contributions that could be used to constrain fits of generalized functionals. By applying these units, a large sample of Skyrme parametrizations is examined for naturalness, which is signaled by dimensionless coupling constants of order one. The bulk of the parameters are found to be natural, with an underlying scale consistent with other determinations. Significant deviations from unity are associated with deficiencies in the corresponding terms of particular functionals or with an incomplete optimization procedure.

preprint2010arXiv

Nuclear Energy Density Optimization

We carry out state-of-the-art optimization of a nuclear energy density of Skyrme type in the framework of the Hartree-Fock-Bogoliubov (HFB) theory. The particle-hole and particle-particle channels are optimized simultaneously, and the experimental data set includes both spherical and deformed nuclei. The new model-based, derivative-free optimization algorithm used in this work has been found to be significantly better than standard optimization methods in terms of reliability, speed, accuracy, and precision. The resulting parameter set UNEDFpre results in good agreement with experimental masses, radii, and deformations and seems to be free of finite-size instabilities. An estimate of the reliability of the obtained parameterization is given, based on standard statistical methods. We discuss new physics insights offered by the advanced covariance analysis.

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

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.