BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

T. Duguet

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

ResearcherUnclaimed source record
nucl-thnucl-excond-mat.str-elcond-mat.supr-conquant-phphysics.chem-phcond-mat.quant-gasphysics.comp-ph

Catalog footprint

What is connected

45works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

45 published item(s)

preprint2022arXiv

A First Glimpse at the Shell Structure beyond $^{54}$Ca: Spectroscopy of $^{55}$K, $^{55}$Ca, and $^{57}$Ca

States in the $N=35$ and 37 isotopes $^{55,57}$Ca have been populated by direct proton-induced nucleon removal reactions from $^{56,58}$Sc and $^{56}$Ca beams at the RIBF. In addition, the $(p,2p)$ quasi-free single-proton removal reaction from $^{56}$Ca was studied. Excited states in $^{55}$K, $^{55}$Ca, and $^{57}$Ca were established for the first time via in-beam $γ$-ray spectroscopy. Results for the proton and neutron removal reactions from $^{56}$Ca to states in $^{55}$K and $^{55}$Ca for the level energies, excited state lifetimes, and exclusive cross sections agree well with state-of-the-art theoretical calculations using different approaches. The observation of a short-lived state in $^{57}$Ca suggests a transition in the calcium isotopic chain from single-particle dominated states at $N=35$ to collective excitations at $N=37$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Angular-momentum projection in coupled-cluster theory: structure of $^{34}$Mg

Single-reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for $^8$Be and $^{20}$Ne using a two-nucleon interaction from chiral effective field theory and for $pf$-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus $^{34}$Mg and find agreement with data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rooting the EDF method into the ab initio framework. PGCM-PT formalism based on MR-IMSRG pre-processed Hamiltonians

Recently, ab initio techniques have been successfully connected to the traditional valence-space shell model. In doing so, they can either explicitly provide ab initio shell-model effective Hamiltonians or constrain the construction of empirical ones. In the present work, the possibility to follow a similar path for the nuclear energy density functional (EDF) method is analyzed. For this connection to be actualized, two theoretical techniques are instrumental: the recently proposed ab initio PGCM-PT many-body formalism and the MR-IMSRG pre-processing of the nuclear Hamiltonian. Based on both formal arguments and numerical results, possible new lines of research are briefly discussed, namely to compute ab initio EDF effective Hamiltonians at low computational cost, to constrain empirical ones or to produce them directly via an effective field theory that remains to be invented.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

ADG: Automated generation and evaluation of many-body diagrams III. Bogoliubov in-medium similarity renormalization group formalism

The goal of the present paper is twofold. First, a novel expansion many-body method applicable to superfluid open-shell nuclei, the so-called Bogoliubov in-medium similarity renormalization group (BIMSRG) theory, is formulated. This generalization of standard single-reference IMSRG theory for closed-shell systems parallels the recent extensions of coupled cluster, self-consistent Green's function or many-body perturbation theory. Within the realm of IMSRG theories, BIMSRG provides an interesting alternative to the already existing multi-reference IMSRG (MR-IMSRG) method applicable to open-shell nuclei. The algebraic equations for low-order approximations, i.e., BIMSRG(1) and BIMSRG(2), can be derived manually without much difficulty. However, such a methodology becomes already impractical and error prone for the derivation of the BIMSRG(3) equations, which are eventually needed to reach high accuracy. Based on a diagrammatic formulation of BIMSRG theory, the second objective of the present paper is thus to describe the third version (v3.0.0) of the ADG code that automatically (1) generates all valid BIMSRG(n) diagrams and (2) evaluates their algebraic expressions in a matter of seconds. This is achieved in such a way that equations can easily be retrieved for both the flow equation and the Magnus expansion formulations of BIMSRG. Expanding on this work, the first future objective is to numerically implement BIMSRG(2) (eventually BIMSRG(3)) equations and perform ab initio calculations of mid-mass open-shell nuclei.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Novel chiral Hamiltonian and observables in light and medium-mass nuclei

A novel parameterisation of a Hamiltonian based on chiral effective field theory is introduced. Specifically, three-nucleon operators at next-to-next-to-leading order are combined with an existing (and successful) two-body interaction containing terms up to next-to-next-to-next-to-leading order. The resulting potential is labelled $N\!N\!$+$3N\text{(lnl)}$. The objective of the present work is to investigate the performance of this new Hamiltonian across light and medium-mass nuclei. Binding energies, nuclear radii and excitation spectra are computed using no-core shell model and self-consistent Green's function approaches. Calculations with $N\!N\!$+$3N\text{(lnl)}$ are compared to two other representative Hamiltonians currently in use, namely NNLO$_{\text{sat}}$ and the older $N\!N\!$+$3N(400)$. Overall, the performance of the novel interaction is very encouraging. In light nuclei, total energies are generally in good agreement with experimental data. Known spectra are also well reproduced with a few notable exceptions. The good description of ground-state energies carries on to heavier nuclei, all the way from oxygen to nickel isotopes. Except for those involving excitation processes across the $N=20$ gap, which is overestimated by the new interaction, spectra are of very good quality, in general superior to those obtained with NNLO$_{\text{sat}}$. Although largely improving on $N\!N\!$+$3N(400)$ results, charge radii calculated with $N\!N\!$+$3N\text{(lnl)}$ still underestimate experimental values, as opposed to the ones computed with NNLO$_{\text{sat}}$ that successfully reproduce available data on nickel. On the whole, the new two- plus three-nucleon Hamiltonian introduced in the present work represents a promising alternative to existing nuclear interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Particle-number projected Bogoliubov coupled cluster theory. Application to the pairing Hamiltonian

While coupled cluster theory accurately models weakly correlated quantum systems, it often fails in the presence of strong correlations where the standard mean-field picture is qualitatively incorrect. In many cases, these failures can be largely ameliorated by permitting the mean-field reference to break physical symmetries. Symmetry-broken coupled cluster, e.g. Bogoliubov coupled cluster, theory can indeed provide reasonably accurate energetic predictions, but the broken symmetry can compromise the quality of the resulting wave function and predictions of observables other than the energy. Merging symmetry projection and coupled cluster theory is therefore an appealing way to describe strongly correlated systems. Independently, two different but related formalisms have been recently proposed to achieve this goal. The two formalisms are contrasted in this manuscript, with results tested on the Richardson pairing Hamiltonian. Both formalisms are based on the disentangled cluster representation of the symmetry-rotated coupled cluster wavefunction. However, they differ in the way that the disentangled clusters are solved. One approach sets up angle-dependent coupled cluster equations, while the other involves first-order ordinary differential equations. The latter approach yields energies and occupation probabilities significantly better than those of number-projected BCS and BCS coupled cluster and, when the disentangled clusters are truncated at low excitation levels, has a computational cost not too much larger than that of BCS coupled cluster. The high quality of results presented in this manuscript indicates that symmetry-projected coupled cluster is a promising method that can accurately describe both weakly and strongly correlated finite many-fermion systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab initio-driven nuclear energy density functional method. A proposal for safe/correlated/improvable parametrizations of the off-diagonal EDF kernels

This programmatic paper lays down the possibility to reconcile the necessity to resum many-body correlations into the energy kernel with the fact that safe multi-reference energy density functional (EDF) calculations cannot be achieved whenever the Pauli principle is not strictly enforced, as is for example the case when many-body correlations are parametrized under the form of empirical density dependencies. Our proposal is to exploit a newly developed ab initio many-body formalism to guide the construction of safe, explicitly correlated and systematically improvable parametrizations of the {\it off-diagonal} energy and norm kernels that lie at the heart of the nuclear EDF method. The many-body formalism of interest relies on the concepts of symmetry breaking {\it and} restoration that have made the fortune of the nuclear EDF method and is, as such, amenable to this guidance. After elaborating on our proposal, we briefly outline the project we plan to execute in the years to come.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-observable nature of the nuclear shell structure. Meaning, illustrations and consequences

The concept of single-nucleon shells constitutes a basic pillar of our understanding of nuclear structure. Effective single-particle energies (ESPEs) introduced by French and Baranger represent the most appropriate tool to relate many-body observables to a single-nucleon shell structure. As briefly discussed in [T. Duguet, G. Hagen, Phys. Rev. C {\bf 85}, 034330 (2012)], the dependence of ESPEs on one-nucleon transfer probability matrices makes them purely theoretical quantities that "run" with the non-observable resolution scale $λ$ employed in the calculation. Given that ESPEs provide a way to interpret the many-body problem in terms of simpler theoretical ingredients, the goal is to specify the terms, i.e. the exact sense and conditions, in which this interpretation can be conducted meaningfully. State-of-the-art multi-reference in-medium similarity renormalization group and self-consistent Gorkov Green's function many-body calculations are employed to corroborate the formal analysis. This is done by comparing the behavior of several observables and of non-observable ESPEs (and spectroscopic factors) under (quasi) unitary similarity renormalization group transformations of the Hamiltonian parameterized by the resolution scale $λ$. The non-observable nature of the nuclear shell structure, i.e. the fact that it constitutes an intrinsically theoretical object with no counterpart in the empirical world, must be recognized and assimilated. Eventually, practitioners can refer to nuclear shells and spectroscopic factors in their analyses of nuclear phenomena if, and only if, they use consistent structure and reaction theoretical schemes based on a fixed resolution scale they have agreed on prior to performing their analysis and comparisons.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symmetry broken and restored coupled-cluster theory. II. Global gauge symmetry and particle number

We have recently extended many-body perturbation theory and coupled-cluster theory performed on top of a Slater determinant breaking rotational symmetry to allow for the restoration of the angular momentum at any truncation order [T. Duguet, J. Phys. G: Nucl. Part. Phys. 42 (2015) 025107]. Following a similar route, we presently extend Bogoliubov many-body perturbation theory and Bogoliubov coupled cluster theory performed on top of a Bogoliubov reference state breaking global gauge symmetry to allow for the restoration of the particle number at any truncation order. Eventually, formalisms can be merged to handle $SU(2)$ and $U(1)$ symmetries at the same time. Several further extensions of the newly proposed many-body formalisms can be foreseen in the mid-term future. The long-term goal relates to the ab initio description of near-degenerate finite quantum systems with an open-shell character.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Shell structure of potassium isotopes deduced from their magnetic moments

\item[Background] Ground-state spins and magnetic moments are sensitive to the nuclear wave function, thus they are powerful probes to study the nuclear structure of isotopes far from stability. \item[Purpose] Extend our knowledge about the evolution of the $1/2^+$ and $3/2^+$ states for K isotopes beyond the $N = 28$ shell gap. \item[Method] High-resolution collinear laser spectroscopy on bunched atomic beams. \item[Results] From measured hyperfine structure spectra of K isotopes, nuclear spins and magnetic moments of the ground states were obtained for isotopes from $N = 19$ up to $N = 32$. In order to draw conclusions about the composition of the wave functions and the occupation of the levels, the experimental data were compared to shell-model calculations using SDPF-NR and SDPF-U effective interactions. In addition, a detailed discussion about the evolution of the gap between proton $1d_{3/2}$ and $2s_{1/2}$ in the shell model and {\it{ab initio}} framework is also presented. \item[Conclusions] The dominant component of the wave function for the odd-$A$ isotopes up to $^{45}$K is a $π1d_{3/2}^{-1}$ hole. For $^{47,49}$K, the main component originates from a $π2s_{1/2}^{-1}$ hole configuration and it inverts back to the $π1d_{3/2}^{-1}$ in $^{51}$K. For all even-$A$ isotopes, the dominant configuration arises from a $π1d_{3/2}^{-1}$ hole coupled to a neutron in the $ν1f_{7/2}$ or $ν2p_{3/2}$ orbitals. Only for $^{48}$K, a significant amount of mixing with $π2s_{1/2}^{-1} \otimes ν(pf)$ is observed leading to a $I^π=1^{-}$ ground state. For $^{50}$K, the ground-state spin-parity is $0^-$ with leading configuration $π1d_{3/2}^{-1} \otimes ν2p_{3/2}^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Symmetry broken and restored coupled-cluster theory I. Rotational symmetry and angular momentum

We extend coupled-cluster theory performed on top of a Slater determinant breaking rotational symmetry to allow for the exact restoration of the angular momentum at any truncation order. The main objective relates to the description of near-degenerate finite quantum systems with an open-shell character. As such, the newly developed many-body formalism offers a wealth of potential applications and further extensions dedicated to the ab initio description of, e.g., doubly open-shell atomic nuclei and molecule dissociation. The formalism, which encompasses both single-reference coupled cluster theory and projected Hartree-Fock theory as particular cases, permits the computation of usual sets of connected diagrams while consistently incorporating static correlations through the highly non-perturbative restoration of rotational symmetry. Interestingly, the yrast spectroscopy of the system, i.e. the lowest energy associated with each angular momentum, is accessed within a single calculation. A key difficulty presently overcome relates to the necessity to handle generalized energy {\it and} norm kernels for which naturally terminating coupled-cluster expansions could be eventually obtained. The present work focuses on $SU(2)$ but can be extended to any (locally) compact Lie group and to discrete groups, such as most point groups. In particular, the formalism will be soon generalized to $U(1)$ symmetry associated with particle number conservation. This is relevant to Bogoliubov coupled cluster theory that was recently formulated and applied to singly open-shell nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ab-initio self-consistent Gorkov-Green's function calculations of semi-magic nuclei - II. Numerical implementation at second order with a two-nucleon interaction

The newly developed Gorkov-Green's function approach represents a promising path to the ab initio description of medium-mass open-shell nuclei. We discuss the implementation of the method at second order with a two-body interaction, with particular attention to the numerical solution of Gorkov's equation. Different sources of theoretical error and degrees of self-consistency are investigated. We show that Krylov projection techniques with a multi-pivot Lanczos algorithm efficiently handle the growth of poles in the one-body Green's function when Gorkov's equation is solved self-consistently. The end result is a tractable, accurate and gently scaling ab initio scheme applicable to full isotopic chains in the medium-mass region.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Leading chiral three-nucleon forces along isotope chains in the calcium region

Three-nucleon forces (3NFs), and in particular terms of the Fujita-Miyazawa type, strongly influence the structure of neutron-rich exotic isotopes. Ab-initio calculations have shown that chiral two- and three-nucleon interactions correctly reproduce binding energy systematics and neutron driplines of oxygen and nearby isotopes. Exploiting the novel self-consistent Gorkov-Green's function approach, we present the first investigation of Ar, K, Ca, Sc and Ti isotopic chains. Leading chiral 3N interactions are mandatory to reproduce the trend of binding energies throughout these chains and to obtain a good description of two-neutron separation energies. At the same time, nuclei in this mass region are systematically overbound by about 40 MeV and the $N=20$ magic gap is significantly overestimated. We conclude that ab-initio many-body calculations of mid-mass isotopic chains challenge modern theories of elementary nuclear interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Skyrme functional from a three-body pseudo-potential of second-order in gradients. Formalism for central terms

In one way or the other, all modern parametrizations of the nuclear energy density functional (EDF) do not respect the exchange symmetry associated with Pauli's principle. It has been recently shown that this practice jeopardizes multi-reference (MR) EDF calculations by contaminating the energy with spurious self-interactions that, for example, lead to finite steps or even divergences when plotting it as a function of collective coordinates. As of today, the only viable option to bypass these pathologies is to rely on EDF kernels that enforce Pauli's principle from the outset by strictly and exactly deriving from a genuine, i.e. density-independent, Hamilton operator. We wish to develop the most general Skyrme-like EDF parametrization containing linear, bilinear and trilinear terms in the density matrices with up to two gradients, under the key constraint that it derives strictly from an effective Hamilton operator. The most general three-body Skyrme-like pseudo-potential containing up to two gradient operators is constructed to generate the trilinear part. The present study is limited to central terms. Spin-orbit and tensor will be addressed in a forthcoming paper. (See paper for full abstract)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spurious finite-size instabilities in nuclear energy density functionals

It is known that some well-established parametrizations of the EDF do not always provide converged results for nuclei and a qualitative link between this finding and the appearance of finite-size instabilities of SNM near saturation density when computed within the RPA has been pointed out. We seek for a quantitative and systematic connection between the impossibility to converge self-consistent calculations of nuclei and the occurrence of finite-size instabilities in SNM for the example of scalar-isovector (S=0, T=1) instabilities of the standard Skyrme EDF. We aim to establish a stability criterion based on computationally-friendly RPA calculations of SNM that is independent on the functional form of the EDF and that can be utilized during the adjustment of its coupling constants. Tuning the coupling constant $C^{ρΔρ}_{1}$ of the gradient term that triggers scalar-isovector instabilities of the standard Skyrme EDF, we find that the occurrence of instabilities in finite nuclei depends strongly on the numerical scheme used to solve the self-consistent mean-field equations. The link to instabilities of SNM is made by extracting the lowest density $ρ_{\text{crit}}$ at which a pole appears at zero energy in the RPA response function when employing the critical value of the coupling constant $C^{ρΔρ}_{1}$ extracted in nuclei. Our analysis suggests a two-fold stability criterion to avoid scalar-isovector instabilities: (i) The density $ρ_{\text{min}}$ corresponding to the lowest pole in the RPA response function should be larger than about 1.2 times the saturation density; (ii) one needs to verify that $ρ_{p}(q_{\text{pq}})$ exhibits a distinct global minimum and is not a decreasing function for large transferred momenta.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The nuclear energy density functional formalism

The present document focuses on the theoretical foundations of the nuclear energy density functional (EDF) method. As such, it does not aim at reviewing the status of the field, at covering all possible ramifications of the approach or at presenting recent achievements and applications. The objective is to provide a modern account of the nuclear EDF formalism that is at variance with traditional presentations that rely, at one point or another, on a {\it Hamiltonian-based} picture. The latter is not general enough to encompass what the nuclear EDF method represents as of today. Specifically, the traditional Hamiltonian-based picture does not allow one to grasp the difficulties associated with the fact that currently available parametrizations of the energy kernel $E[g',g]$ at play in the method do not derive from a genuine Hamilton operator, would the latter be effective. The method is formulated from the outset through the most general multi-reference, i.e. beyond mean-field, implementation such that the single-reference, i.e. "mean-field", derives as a particular case. As such, a key point of the presentation provided here is to demonstrate that the multi-reference EDF method can indeed be formulated in a {\it mathematically} meaningful fashion even if $E[g',g]$ does {\it not} derive from a genuine Hamilton operator. In particular, the restoration of symmetries can be entirely formulated without making {\it any} reference to a projected state, i.e. within a genuine EDF framework. However, and as is illustrated in the present document, a mathematically meaningful formulation does not guarantee that the formalism is sound from a {\it physical} standpoint. The price at which the latter can be enforced as well in the future is eventually alluded to.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-nucleon forces in exotic open-shell isotopes

Advances in the self-consistent Green's function approach to finite nuclei are discussed, including the implementation of three-nucleon forces and the extension to the Gorkov formalism. We report results on binding energies in the nitrogen and fluorine isotopic chains, as well as spectral functions of 22O. The application to medium-mass open-shell systems is illustrated by separation energy spectra of two argon isotopes, which are compared to one-neutron removal experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab-initio approach to effective single-particle energies in doubly closed shell nuclei

The present work discusses, from an ab-initio standpoint, the definition, the meaning, and the usefulness of effective single-particle energies (ESPEs) in doubly closed shell nuclei. We perform coupled-cluster calculations to quantify to what extent selected closed-shell nuclei in the oxygen and calcium isotopic chains can effectively be mapped onto an effective independent-particle picture. To do so, we revisit in detail the notion of ESPEs in the context of strongly correlated many-nucleon systems and illustrate the necessity to extract ESPEs through the diagonalization of the centroid {\it matrix}, as originally argued by Baranger. For the purpose of illustration, we analyse the impact of correlations on observable one-nucleon separation energies and non-observable ESPEs in selected closed-shell oxygen and calcium isotopes. We then state and illustrate the non-observability of ESPEs. Similarly to spectroscopic factors, ESPEs can indeed be modified by a redefinition of inaccessible quantities while leaving actual observables unchanged. This leads to the absolute necessity to employ consistent structure and reaction models based on the same nuclear Hamiltonian to extract the shell structure in a meaningful fashion from experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab-initio Gorkov-Green's function calculations of open-shell nuclei

We present results from a new ab-initio method that uses the self-consistent Gorkov Green's function theory to address truly open-shell systems. The formalism has been recently worked out up to second order and is implemented here in nuclei for the first time on the basis of realistic nuclear forces. We find good convergence of the results with respect to the basis size in Ca44 and Ni74 and discuss quantities of experimental interest including ground-state energies, pairing gaps and particle addition/removal spectroscopy. These results demonstrate that the Gorkov method is a valid alternative to multireference approaches for tackling degenerate or near degenerate quantum systems. In particular, it increases the number of mid-mass nuclei accessible in an ab-initio fashion from a few tens to a few hundreds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Skyrme pseudo-potential-based EDF parametrisation for spuriousity-free MR-EDF calculations

First exploratory steps towards a pseudo-potential-based Skyrme energy density functional for spuriousity-free multi-reference calculations are presented. A qualitatively acceptable fit can be accomplished by adding simple three- and four-body contact terms to the standard central plus spin-orbit two-body terms. To achieve quantitative predictive power, higher-order terms, e.g. velocity-dependent three-body terms, will be required.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Toward the Ab-initio Description of Medium Mass Nuclei

As ab-initio calculations of atomic nuclei enter the A=40-100 mass range, a great challenge is how to approach the vast majority of open-shell (degenerate) isotopes. We add realistic three-nucleon interactions to the state of the art many-body Green's function theory of closed-shells, and find that physics of neutron driplines is reproduced with very good quality. Further, we introduce the Gorkov formalism to extend ab-initio theory to semi-magic, fully open-shell, isotopes. Proof-of-principle calculations for Ca-44 and Ni-74 confirm that this approach is indeed feasible. Combining these two advances (open-shells and three-nucleon interactions) requires longer, technical, work but it is otherwise within reach.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ab-initio self-consistent Gorkov-Green's function calculations of semi-magic nuclei - I. Formalism at second order with a two-nucleon interaction

An ab-initio calculation scheme for finite nuclei based on self-consistent Green's functions in the Gorkov formalism is developed. It aims at describing properties of doubly-magic and semi-magic nuclei employing state-of-the-art microscopic nuclear interactions and explicitly treating pairing correlations through the breaking of U(1) symmetry associated with particle number conservation. The present paper introduces the formalism, necessary to undertake applications at (self-consistent) second-order using two-nucleon interactions, in a detailed and self-contained fashion. First applications of such a scheme will be reported soon in a forthcoming publication. Future works will extend the present scheme to include three-nucleon interactions and implement more advanced truncation schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An improved density matrix expansion for spin-unsaturated nuclei

A current objective of low-energy nuclear theory is to build non-empirical nuclear energy density functionals (EDFs) from underlying inter-nucleon interactions and many-body perturbation theory (MBPT). The density matrix expansion (DME) of Negele and Vautherin is a convenient method to map highly non-local Hartree-Fock expressions into the form of a quasi-local Skyrme functional with density-dependent couplings. In this work, we assess the accuracy of the DME at reproducing the non-local exchange (Fock) contribution to the energy. In contrast to the scalar part of the density matrix for which the original formulation of Negele and Vautherin is reasonably accurate, we demonstrate the necessity to reformulate the DME for the vector part of the density matrix, which is needed for an accurate description of spin-unsaturated nuclei. Phase-space averaging techniques are shown to yield a significant improvement for the vector part of the density matrix compared to the original formulation of Negele and Vautherin. The key to the improved accuracy is to take into account the anisotropy that characterizes the local-momentum distribution in the surface region of finite Fermi systems. Optimizing separately the DME for the central, tensor and spin-orbit contributions to the Fock energy, one reaches a few-percent accuracy over a representative set of semi-magic nuclei. With such an accuracy at hand, one can envision using the corresponding Skyrme-like energy functional as a microscopically-constrained starting point around which future phenomenological parameterizations can be built and refined.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Breaking and restoring symmetries within the nuclear energy density functional method

We review the notion of symmetry breaking and restoration within the frame of nuclear energy density functional methods. We focus on key differences between wave-function- and energy-functional-based methods. In particular, we point to difficulties to formulate the restoration of symmetries within the energy functional framework. The problems tackled recently in connection with particle-number restoration serve as a baseline to the present discussion. Reaching out to angular-momentum restoration, we identify an exact mathematical property of the energy density $E^{LM}(\vec{R})$ that could be used to constrain energy density functional kernels. Consequently, we suggest possible routes towards a better formulation of symmetry restorations within energy density functional methods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lowest-order contributions of chiral three-nucleon interactions to pairing properties of nuclear ground states

We perform a systematic study of the odd-even mass staggering generated using a pairing interaction computed at first order in low-momentum interactions. Building on previous work including the (nuclear plus Coulomb) two-nucleon interaction only, we focus here on the first-order contribution from chiral three-nucleon forces. We observe a significant repulsive effect from the three-nucleon interaction. The combined contribution from two- and three-nucleon interactions accounts for approximately 70% of the experimental pairing gaps. This leaves room for higher-order contributions to the pairing kernel and the normal self-energy that need to be computed consistently.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Microscopically-constrained Fock energy density functionals from chiral effective field theory. I. Two-nucleon interactions

The density matrix expansion (DME) of Negele and Vautherin is a convenient tool to map finite-range physics associated with vacuum two- and three-nucleon interactions into the form of a Skyme-like energy density functional (EDF) with density-dependent couplings. In this work, we apply the improved formulation of the DME proposed recently in arXiv:0910.4979 by Gebremariam {\it et al.} to the non-local Fock energy obtained from chiral effective field theory (EFT) two-nucleon (NN) interactions at next-to-next-to-leading-order (N$^2$LO). The structure of the chiral interactions is such that each coupling in the DME Fock functional can be decomposed into a cutoff-dependent coupling {\it constant} arising from zero-range contact interactions and a cutoff-independent coupling {\it function} of the density arising from the universal long-range pion exchanges. This motivates a new microscopically-guided Skyrme phenomenology where the density-dependent couplings associated with the underlying pion-exchange interactions are added to standard empirical Skyrme functionals, and the density-independent Skyrme parameters subsequently refit to data. A Mathematica notebook containing the novel density-dependent couplings is provided.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the breaking and restoration of symmetries within the nuclear energy density functional formalism

We review the notion of symmetry breaking and restoration within the frame of nuclear energy density functional methods. We focus on key differences between wave-function- and energy-functional-based methods. In particular, we point to difficulties encountered within the energy functional framework and discuss new potential constraints on the underlying energy density functional that could make the restoration of broken symmetries better formulated within such a formalism. We refer to Ref.~\cite{duguet10a} for details.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Regularization of Multi-Reference Energy Density Functional Calculations

We report the first application of a recently proposed regularization procedure for multi-reference energy density functionals, which removes spurious divergent or non-continuous contributions to the binding energy, to a general configuration mixing. As an example, we present a calculation that corresponds to the particle-number and angular momentum projection of axially symmetric time-reversal invariant quasiparticle vacua of different quadrupole deformation for the nucleus O18. The SIII parameterization of the Skyrme energy functional is used.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Configuration Mixing within the Energy Density Functional Formalism: Removing Spurious Contributions from Non-Diagonal Energy Kernels

Multi-reference calculations along the lines of the Generator Coordinate Method or the restoration of broken symmetries within the nuclear Energy Density Functional (EDF) framework are becoming a standard tool in nuclear structure physics. These calculations rely on the extension of a single-reference energy functional, of the Gogny or the Skyrme types, to non-diagonal energy kernels. There is no rigorous constructive framework for this extension so far. The commonly accepted way proceeds by formal analogy with the expressions obtained when applying the generalized Wick theorem to the non-diagonal matrix element of a Hamilton operator between two product states. It is pointed out that this procedure is ill-defined when extended to EDF calculations as the generalized Wick theorem is taken outside of its range of applicability. In particular, such a procedure is responsible for the appearance of spurious divergences and steps in multi-reference EDF energies, as was recently observed in calculations restoring particle number or angular momentum. In the present work, we give a formal analysis of the origin of this problem for calculations with and without pairing, i.e. constructing the density matrices from either Slater determinants or quasi-particle vacua. We propose a correction to energy kernels that removes the divergences and steps, and which is applicable to calculations based on any symmetry restoration or generator coordinate. The method is formally illustrated for particle number restoration and is specified to configuration mixing calculations based on Slater determinants.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Halo phenomenon in finite many-fermion systems. Atom-positron complexes and large-scale study of atomic nuclei

The analysis method proposed in Ref. \cite{rotival07a} is applied to characterize halo properties in finite many-fermion systems. First, the versatility of the method is highlighted by applying it to light and medium-mass nuclei as well as to atom-positron and ion-positronium complexes. Second, the dependence of nuclear halo properties on the characteristics of the energy density functional used in self-consistent Hartree-Fock-Bogoliubov calculations is studied. It is found that (a) the low-density behavior of the pairing functional and the regularization/renormalization scheme must be chosen coherently and with care to provide meaningful predictions, (b) the impact of pairing correlations on halo properties is significant and is the result of two competing effects, (c) the detailed characteristics of the pairing functional has however only little importance, (d) halo properties depend significantly on any ingredient of the energy density functional that influences the location of single-particle levels; i.e. the effective mass, the tensor terms and the saturation density of nuclear matter. The latter dependencies give insights to how experimental data on medium-mass drip-line nuclei can be used in the distant future to constrain some characteristics of the nuclear energy density functional. Last but not least, large scale predictions of halos among all spherical even-even nuclei are performed using specific sets of particle-hole and particle-particle energy functionals. It is shown that halos in the ground state of medium-mass nuclei will only be found at the very limit of neutron stability and for a limited number of elements.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

New analysis method of the halo phenomenon in finite many-fermion systems. First applications to medium-mass atomic nuclei

A new analysis method to investigate halos in finite many-fermion systems is designed, as existing characterization methods are proven to be incomplete/inaccurate. A decomposition of the internal wave-function of the {$N$-body} system in terms of overlap functions allows a model-independent analysis of medium-range and asymptotic properties of the internal one-body density. The existence of a spatially decorrelated region in the density profile is related to the existence of three typical energy scales in the excitation spectrum of the {$(N-1)$-body} system. A series of model-independent measures, taking the internal density as the only input, are introduced. The new measures allow a quantification of the potential halo in terms of the average number of fermions participating to it and of its impact on the system extension. Those new "halo factors" are validated through simulations and applied to results obtained through energy density functional calculations of medium-mass nuclei. Performing spherical Hartree-Fock-Bogoliubov calculations with state-of-the-art Skyrme plus pairing functionals, a collective halo is predicted in drip-line Cr isotopes, whereas no such effect is seen in Sn isotopes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-empirical nuclear energy functionals, pairing gaps and odd-even mass differences

First, we briefly outline some aspects of the starting project to design non-empirical energy functionals based on low-momentum vacuum interactions and many-body perturbation theory. Second, we present results obtained within an approximation of such a scheme where the pairing part of the energy density functional is constructed at first order in the nuclear plus Coulomb two-body interaction. We discuss in detail the physics of the odd-even mass staggering and the necessity to compute actual odd-even mass differences to analyze it meaningfully.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Non-empirical pairing energy density functional. First order in the nuclear plus Coulomb two-body interaction

We perform systematic calculations of pairing gaps in semi-magic nuclei across the nuclear chart using the Energy Density Functional method and a {\it non-empirical} pairing functional derived, without further approximation, at lowest order in the two-nucleon vacuum interaction, including the Coulomb force. The correlated single-particle motion is accounted for by the SLy4 semi-empirical functional. Rather unexpectedly, both neutron and proton pairing gaps thus generated are systematically close to experimental data. Such a result further suggests that missing effects, i.e. higher partial-waves of the NN interaction, the NNN interaction and the coupling to collective fluctuations, provide an overall contribution that is sub-leading as for generating pairing gaps in nuclei. We find that including the Coulomb interaction is essential as it reduces proton pairing gaps by up to 40%.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Particle-Number Restoration within the Energy Density Functional Formalism

We give a detailed analysis of the origin of spurious divergences and finite steps that have been recently identified in particle-number restoration calculations within the nuclear energy density functional framework. We isolate two distinct levels of spurious contributions to the energy. The first one is encoded in the definition of the basic energy density functional itself whereas the second one relates to the canonical procedure followed to extend the use of the energy density functional to multi-reference calculations. The first level of spuriosity relates to the long-known self-interaction problem and to the newly discussed self-pairing interaction process which might appear when describing paired systems with energy functional methods using auxiliary reference states of Bogoliubov or BCS type. A minimal correction to the second level of spuriosity to the multi-reference nuclear energy density functional proposed in [D. Lacroix, T. Duguet, M. Bender, arXiv:0809.2041] is shown to remove completely the anomalies encountered in particle-number restored calculations. In particular, it restores sum-rules over (positive) particle numbers that are to be fulfilled by the particle-number-restored formalism. The correction is found to be on the order of several hundreds of keVs up to about 1 MeV in realistic calculations, which is small compared to the total binding energy, but often accounts for a substantial percentage of the energy gain from particle-number restoration and is on the same energy scale as the excitations one addresses with multi-reference energy density functional methods.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Particle-Number Restoration within the Energy Density Functional formalism: Nonviability of terms depending on noninteger powers of the density matrices

We discuss the origin of pathological behaviors that have been recently identified in particle-number-restoration calculations performed within the nuclear energy density functional framework. A regularization method that removes the problematic terms from the multi-reference energy density functional and which applies (i) to any symmetry restoration- and/or generator-coordinate-method-based configuration mixing calculation and (ii) to energy density functionals depending only on integer powers of the density matrices, was proposed in [D. Lacroix, T. Duguet, M. Bender, arXiv:0809.2041] and implemented for particle-number restoration calculations in [M. Bender, T. Duguet, D. Lacroix, arXiv:0809.2045]. In the present paper, we address the viability of non-integer powers of the density matrices in the nuclear energy density functional. Our discussion builds upon the analysis already carried out in [J. Dobaczewski \emph{et al.}, Phys. Rev. C \textbf{76}, 054315 (2007)]. First, we propose to reduce the pathological nature of terms depending on a non-integer power of the density matrices by regularizing the fraction that relates to the integer part of the exponent using the method proposed in [D. Lacroix, T. Duguet, M. Bender, arXiv:0809.2041]. Then, we discuss the spurious features brought about by the remaining fractional power. Finally, we conclude that non-integer powers of the density matrices are not viable and should be avoided in the first place when constructing nuclear energy density functionals that are eventually meant to be used in multi-reference calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Symbolic computation of the Hartree-Fock energy from a chiral EFT three-nucleon interaction at N$^2$LO

We present the first of a two-part Mathematica notebook collection that implements a symbolic approach for the application of the density matrix expansion (DME) to the Hartree-Fock (HF) energy from a chiral effective field theory (EFT) three-nucleon interaction at N$^2$LO. The final output from the notebooks is a Skyrme-like energy density functional that provides a quasi-local approximation to the nonlocal HF energy. In this paper, we discuss the derivation of the HF energy and its simplification in terms of the scalar/vector-isoscalar/isovector parts of the one-body density matrix. Furthermore, a set of steps is described and illustrated on how to extend the approach to other three-nucleon interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Symbolic integration of a product of two spherical bessel functions with an additional exponential and polynomial factor

We present a mathematica package that performs the symbolic calculation of integrals of the form \int^{\infty}_0 e^{-x/u} x^n j_ν (x) j_μ (x) dx where $j_ν (x)$ and $j_μ (x)$ denote spherical Bessel functions of integer orders, with $ν\ge 0$ and $μ\ge 0$. With the real parameter $u>0$ and the integer $n$, convergence of the integral requires that $n+ν+μ\ge 0$. The package provides analytical result for the integral in its most simplified form. The novel symbolic method employed enables the calculation of a large number of integrals of the above form in a fraction of the time required for conventional numerical and Mathematica based brute-force methods. We test the accuracy of such analytical expressions by comparing the results with their numerical counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The tensor part of the Skyrme energy density functional. I. Spherical nuclei

We perform a systematic study of the impact of the J^2 tensor term in the Skyrme energy functional on properties of spherical nuclei. In the Skyrme energy functional, the tensor terms originate both from zero-range central and tensor forces. We build a set of 36 parameterizations, which covers a wide range of the parameter space of the isoscalar and isovector tensor term coupling constants, with a fit protocol very similar to that of the successful SLy parameterizations. We analyze the impact of the tensor terms on a large variety of observables in spherical mean-field calculations, such as the spin-orbit splittings and single-particle spectra of doubly-magic nuclei, the evolution of spin-orbit splittings along chains of semi-magic nuclei, mass residuals of spherical nuclei, and known anomalies of charge radii. Our main conclusion is that the currently used central and spin-orbit parts of the Skyrme energy density functional are not flexible enough to allow for the presence of large tensor terms.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Excitation Energies and Spins of the Yrast Superdeformed Band in ^{191}Hg

The excitation energies and spins of the levels in the yrast superdeformed band of $^{191}$Hg have been determined from two single-step $γ$ transitions and the quasi-continuum spectrum connecting the superdeformed and normal-deformed states. The results are compared with those from theoretical mean-field calculations with different interactions. A discussion of pairing in superdeformed states is also included.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Rotational properties of $^{252, 253, 254}$No. Influence of pairing correlations

Rotational bands of $^{252, 253, 254}$No and the fission barriers of $^{254}$No at spin 0$\hbar$ and 20$\hbar$ are calculated with the Hartree-Fock-Bogolyubov theory and the Lipkin-Nogami approximate particle number projection. The SLy4 Skyrme force is used in the particle-hole channel. A zero-range force with and without density-dependence is used in the particle-particle channel. The experimental ground state deformation (${\cal{Q}}_{20}$ = 32.8 b) is reproduced as well as the increase of the dynamical moment of inertia with frequency both for $^{252}$No and $^{254}$No. The rotational band of $^{253}$No is also calculated. Fission barriers of $^{254}$No at spin 0$\hbar$ and 20$\hbar$ show the robustness of shell-corrections against rotation in these heavy nuclei.

0 citations0 reviewsNo code linkedOpen access