BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

R. J. Furnstahl

R. J. Furnstahl contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
nucl-thhep-phnucl-exhep-latphysics.data-anastro-ph.HEastro-ph.COastro-ph.GAastro-ph.SRcond-mat.mtrl-scicond-mat.otherphysics.atom-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
35works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

35 published item(s)

preprint2022arXiv

Model reduction methods for nuclear emulators

The field of model order reduction (MOR) is growing in importance due to its ability to extract the key insights from complex simulations while discarding computationally burdensome and superfluous information. We provide an overview of MOR methods for the creation of fast & accurate emulators of memory- and compute-intensive nuclear systems. As an example, we describe how "eigenvector continuation" is a special case of a much more general and well-studied MOR formalism for parameterized systems. We continue with an introduction to the Ritz and Galerkin projection methods that underpin many such emulators, while pointing to the relevant MOR theory and its successful applications along the way. We believe that this will open the door to broader applications in nuclear physics and facilitate communication with practitioners in other fields.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Renormalization group evolution of optical potentials: explorations using a toy model

To take full advantage of experimental facilities such as FRIB for applications to nuclear astrophysics, nuclear structure, and explorations of neutrinos and fundamental symmetries, we need a better understanding of the interplay of reaction and structure theory. The renormalization group (RG) is the natural tool for maintaining a consistent treatment of reaction and structure. Here we make a first study of RG for optical potentials, which are important ingredients for direct reactions. To simplify the analysis, we use a pedagogical one-dimensional model and evolve toward low RG resolution using the similarity RG (or SRG). We show how SRG decoupling at low resolution carries over to the optical potential and enhances perturbative approximations, and how induced SRG nonlocality compares to the nonlocality of the optical potential. We discuss the results in the larger context of consistent SRG evolution of operators and wave functions in the analysis of direct reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables

We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($χ$EFT) that are provided by bound-state observables in the $A=3$ and $A=4$ sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the $χ$EFT expansion at next-to-next-to-leading order. A consistent solution for the ${}^3$H binding energy, the ${}^4$He binding energy and radius, and the ${}^3$H $β$-decay rate can only be obtained if $χ$EFT truncation errors are included in the analysis. All of these except the $β$-decay rate give essentially degenerate constraints on the 3NF low-energy constants, so it is crucial for estimating these parameters. We use eigenvector continuation for fast and accurate emulation of No-Core Shell Model calculations of the considered few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the hyperparameters that quantify the truncation error. We find a $χ$EFT expansion parameter of $Q=0.33 \pm 0.06$ for these observables.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The quasi-deuteron model at low RG resolution

The quasi-deuteron model introduced by Levinger is used to explain cross sections for knocking out high-momentum protons in photo-absorption on nuclei. This is within a framework we characterize as exhibiting high renormalization group (RG) resolution. Assuming a one-body reaction operator, the nuclear wave function must include two-body short-range correlations (SRCs) with deuteron-like quantum numbers. In a previous paper, we showed that SRC physics can be naturally accounted for at low RG resolution. Here we describe the quasi-deuteron model at low RG resolution and determine the Levinger constant, which is proportional to the ratio of nuclear photo-absorption to that for photo-disintegration of a deuteron. We extract the Levinger constant based on the ratio of momentum distributions at high relative momentum. We compute momentum distributions evolved under similarity RG (SRG) transformations, where the SRC physics is shifted into the operator as a universal two-body term. The short-range nature of this operator motivates using local-density approximations with uncorrelated wave functions in evaluating nuclear matrix elements, which greatly simplifies the analysis. The operator must be consistently matched to the RG scale and scheme of the interaction for a reliable extraction. We apply SRG transformations to different nucleon-nucleon (NN) interactions and use the deuteron wave functions and Weinberg eigenvalues to determine approximate matching scales. We predict the Levinger constant for several NN interactions and a wide range of nuclei comparing to experimental extractions. The predictions at low RG resolution are in good agreement with experiment when starting with a hard NN interaction and the initial operator. Similar agreement is found using soft NN interactions when the additional two-body operator induced by evolution from hard to soft is included.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. Here we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Efficient emulators for scattering using eigenvector continuation

Eigenvector continuation EC has been shown to accurately and efficiently reproduce ground states for targeted sets of Hamiltonian parameters. It uses as variational basis vectors the corresponding ground-state eigensolutions from selected other sets of parameters. Here we extend the EC approach to scattering using the Kohn variational principle. We first test it using a model for S-wave nucleon-nucleon scattering and then demonstrate that it also works to give accurate predictions for non-local potentials, charged-particle scattering, complex optical potentials, and higher partial waves. These proofs-of-principle validate EC as an accurate emulator for applying Bayesian inference to parameter estimation constrained by scattering observables. The efficiency of such emulators is because the accuracy is achieved with a small number of variational basis elements and the central computations are just linear algebra calculations in the space spanned by this basis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fast emulation of quantum three-body scattering

We develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from $10^{-13}$ to $10^{-4}$ depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

How well do we know the neutron-matter equation of state at the densities inside neutron stars? A Bayesian approach with correlated uncertainties

We introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory ($χ$EFT). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as $χ$EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the $χ$EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density the predicted symmetry energy and its slope are consistent with experimental constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory ($χ$EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the $χ$EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. We then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium $χ$EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Operator evolution from the similarity renormalization group and the Magnus expansion

The Magnus expansion is an efficient alternative to solving similarity renormalization group (SRG) flow equations with high-order, memory-intensive ordinary differential equation solvers. The numerical simplifications it offers for operator evolution are particularly valuable for in-medium SRG calculations, though challenges remain for difficult problems involving intruder states. Here we test the Magnus approach in an analogous but more accessible situation, which is the free-space SRG treatment of the spurious bound-states arising from a leading-order chiral effective field theory (EFT) potential with very high cutoffs. We show that the Magnus expansion passes these tests and then use the investigations as a springboard to address various aspects of operator evolution that have renewed relevance in the context of the scale and scheme dependence of nuclear processes. These aspects include SRG operator flow with band- versus block-diagonal generators, universality for chiral EFT Hamiltonians and associated operators with different regularization schemes, and the impact of factorization arising from scale separation. Implications for short-range correlations physics and the possibilities for reconciling high- and low-resolution treatments of nuclear structure and reactions are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Turning the nuclear energy density functional method into a proper effective field theory: reflections

Nuclear energy density functionals (EDFs) have a long history of success in reproducing properties of nuclei across the table of the nuclides. They capture quantitatively the emergent features of bound nuclei, such as nuclear saturation and pairing, yet greater accuracy and improved uncertainty quantification are actively sought. Implementations of phenomenological EDFs are suggestive of effective field theory (EFT) formulations and there are hints of an underlying power counting. Multiple paths are possible in trying to turn the nuclear EDF method into a proper EFT. I comment on the current situation and speculate on how to proceed using an effective action formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Infrared extrapolations for atomic nuclei

Harmonic oscillator model-space truncations introduce systematic errors to the calculation of binding energies and other observables. We identify the relevant infrared scaling variable and give values for this nucleus-dependent quantity. We consider isotopes of oxygen computed with the coupled-cluster method from chiral nucleon-nucleon interactions at next-to-next-to-leading order and show that the infrared component of the error is sufficiently understood to permit controlled extrapolations. By employing oscillator spaces with relatively large frequencies, well above the energy minimum, the ultraviolet corrections can be suppressed while infrared extrapolations over tens of MeVs are accurate for ground-state energies. However, robust uncertainty quantification for extrapolated quantities that fully accounts for systematic errors is not yet developed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Theory and Science of the Facility for Rare Isotope Beams

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Systematic expansion for infrared oscillator basis extrapolations

Recent work has demonstrated that the infrared effects of harmonic oscillator basis truncations are well approximated by imposing a partial-wave Dirichlet boundary condition at a properly identified radius L. This led to formulas for extrapolating the corresponding energy E_L and other observables to infinite L and thus infinite basis size. Here we reconsider the energy for a two-body system with a Dirichlet boundary condition at L to identify and test a consistent and systematic expansion for E_L that depends only on observables. We also generalize the energy extrapolation formula to nonzero angular momentum, and apply it to the deuteron. Formulas given previously for extrapolating the radius are derived in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ultraviolet extrapolations in finite oscillator bases

The use of finite harmonic oscillator spaces in many-body calculations introduces both infrared (IR) and ultraviolet (UV) errors. The IR effects are well approximated by imposing a hard-wall boundary condition at a properly identified radius L_eff. We show that duality of the oscillator implies that the UV effects are equally well described by imposing a sharp momentum cutoff at a momentum Lambda_eff complementary to L_eff. By considering two-body systems with separable potentials, we show that the UV energy corrections depend on details of the potential, in contrast to the IR energy corrections, which depend only on the S-matrix. An adaptation of the separable treatment to more general interactions is developed and applied to model potentials as well as to the deuteron with realistic potentials. The previous success with a simple phenomenological form for the UV error is also explained. Possibilities for controlled extrapolations for A > 2 based on scaling arguments are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Universality in Similarity Renormalization Group Evolved Potential Matrix Elements and T-Matrix Equivalence

We examine how the universality of two-nucleon interactions evolved using similarity renormalization group (SRG) transformations correlates with T-matrix equivalence, with the ultimate goal of gaining insight into universality for three-nucleon forces. With sufficient running of the SRG flow equations, the low-energy matrix elements of different realistic potentials evolve to a universal form. Because these potentials are fit to low-energy data, they are (approximately) phase equivalent only up to a certain energy, and we find universality in evolved potentials up to the corresponding momentum. More generally we find universality in local energy regions, reflecting a local decoupling by the SRG. The further requirements for universality in evolved potential matrix elements are explored using two simple alternative potentials. We see evidence that in addition to predicting the same observables, common long-range potentials (i.e., explicit pion physics) is required for universality in the potential matrix elements after SRG flow. In agreement with observations made previously for Vlowk evolution, regions of universal potential matrix elements are restricted to where half-on-shell T-matrix equivalence holds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-resolution probes of low-resolution nuclei

Renormalization group (RG) methods used to soften Hamiltonians allow large-scale computational resources to be used to greater advantage in calculations of nuclear structure and reactions. These RG transformations lower the effective resolution of the nuclei, which raises questions about how to calculate and interpret high-momentum transfer probes of nuclear structure. Such experiments are conventionally explained in terms of short-range correlations, but these disappear with the evolution to low-momentum scales. We highlight the important issues and prospects in the context of recent developments in RG technology, with guidance from the analogous extraction of parton distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutron matter based on consistently evolved chiral three-nucleon interactions

We present the first results for the neutron matter equation of state (EOS) using nucleon-nucleon and three-nucleon chiral effective field theory interactions that are consistently evolved in the framework of the Similarity Renormalization Group (SRG). The dependence of the EOS on the SRG resolution scale is greatly reduced when induced three-nucleon forces (3NF) are included and the residual variation, which in part is from missing induced four-body interactions, is comparable to estimated many-body perturbation theory truncation errors. The relative growth with decreasing resolution of the 3NF contributions to the energy per neutron is of natural size, but it accelerates at the lowest resolutions where strong renormalization of the long-range 3NF matrix elements is also observed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New applications of renormalization group methods in nuclear physics

We review recent developments in the use of renormalization group (RG) methods in low-energy nuclear physics. These advances include enhanced RG technology, particularly for three-nucleon forces, which greatly extends the reach and accuracy of microscopic calculations. We discuss new results for the nucleonic equation of state with applications to astrophysical systems such as neutron stars, new calculations of the structure and reactions of finite nuclei, and new explorations of correlations in nuclear systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

P-shell nuclei using Similarity Renormalization Group evolved three-nucleon interactions

The Similarity Renormalization Group (SRG) is used to soften interactions for ab initio nuclear structure calculations by decoupling low- and high-energy Hamiltonian matrix elements. The substantial contribution of both initial and SRG-induced three-nucleon forces requires their consistent evolution in a three-particle basis space before applying them to larger nuclei. While in principle the evolved Hamiltonians are unitarily equivalent, in practice the need for basis truncation introduces deviations, which must be monitored. Here we present benchmark no-core full configuration calculations with SRG-evolved interactions in p-shell nuclei over a wide range of softening. These calculations are used to assess convergence properties, extrapolation techniques, and the dependence of energies, including four-body contributions, on the SRG resolution scale.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universal properties of infrared oscillator basis extrapolations

Recent work has shown that a finite harmonic oscillator basis in nuclear many-body calculations effectively imposes a hard-wall boundary condition in coordinate space, motivating infrared extrapolation formulas for the energy and other observables. Here we further refine these formulas by studying two-body models and the deuteron. We accurately determine the box size as a function of the model space parameters, and compute scattering phase shifts in the harmonic oscillator basis. We show that the energy shift can be well approximated in terms of the asymptotic normalization coefficient and the bound-state momentum, discuss higher-order corrections for weakly bound systems, and illustrate this universal property using unitarily equivalent calculations of the deuteron.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local Projections of Low-Momentum Potentials

Nuclear interactions evolved via renormalization group methods to lower resolution become increasingly non-local (off-diagonal in coordinate space) as they are softened. This inhibits both the development of intuition about the interactions and their use with some methods for solving the quantum many-body problem. By applying "local projections", a softened interaction can be reduced to a local effective interaction plus a non-local residual interaction. At the two-body level, a local projection after similarity renormalization group (SRG) evolution manifests the elimination of short-range repulsive cores and the flow toward universal low-momentum interactions. The SRG residual interaction is found to be relatively weak at low energy, which motivates a perturbative treatment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Renormalization Group in Nuclear Physics

Modern techniques of the renormalization group (RG) combined with effective field theory (EFT) methods are revolutionizing nuclear many-body physics. In these lectures we will explore the motivation for RG in low-energy nuclear systems and its implementation in systems ranging from the deuteron to neutron stars, both formally and in practice. Flow equation approaches applied to Hamiltonians both in free space and in the medium will be emphasized. This is a conceptually simple technique to transform interactions to more perturbative and universal forms. An unavoidable complication for nuclear systems from both the EFT and flow equation perspective is the need to treat many-body forces and operators, so we will consider these aspects in some detail. We'll finish with a survey of current developments and open problems in nuclear RG.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Similarity Renormalization Group with Novel Generators

The choice of generator in the Similarity Renormalization Group (SRG) flow equation determines the evolution pattern of the Hamiltonian. The kinetic energy has been used in the generator for most prior applications to nuclear interactions, and other options have been largely unexplored. Here we show how variations of this standard choice can allow the evolution to proceed more efficiently without losing its advantages.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Decoupling of Spurious Deep Bound States with the Similarity Renormalization Group

The Similarity Renormalization Group (SRG) is a continuous series of unitary transformations that can be implemented as a flow equation. When the relative kinetic energy ($\Trel$) is used in the SRG generator, nuclear structure calculations have shown greatly improved convergence with basis size because of the decoupling of high-energy and low-energy physics. However this generator can sometimes be problematic. A test case is provided by a study of initial interactions from chiral effective field theories with large cutoffs, which can lead to spurious deep bound states. We would like the SRG to decouple these from the physical shallow bound states. However, with $\Trel$ the high- and low-energy bound states are not decoupled in the usual sense. Replacing $\Trel$ by the momentum-space diagonal of the Hamiltonian ($H_d$) in the SRG generator does produce decoupling, such that the shallow states are in the low-momentum region and the deep bound states are at higher momentum. The flow toward universal low-momentum interactions is also restored.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved nuclear matter calculations from chiral low-momentum interactions

We present new nuclear matter calculations based on low-momentum interactions derived from chiral effective field theory potentials. The current calculations use an improved treatment of the three-nucleon force contribution that includes a corrected combinatorial factor beyond Hartree-Fock that was omitted in previous nuclear matter calculations. We find realistic saturation properties using parameters fit only to few-body data, but with larger uncertainty estimates from cutoff dependence and the 3NF parametrization than in previous calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Operator Evolution via the Similarity Renormalization Group I: The Deuteron

Similarity Renormalization Group (SRG) flow equations can be used to unitarily soften nuclear Hamiltonians by decoupling high-energy intermediate state contributions to low-energy observables while maintaining the natural hierarchy of many-body forces. Analogous flow equations can be used to consistently evolve operators so that observables are unchanged if no approximations are made. The question in practice is whether the advantages of a softer Hamiltonian and less correlated wave functions might be offset by complications in approximating and applying other operators. Here we examine the properties of SRG-evolved operators, focusing in this paper on applications to the deuteron but leading toward methods for few-body systems. We find the advantageous features generally carry over to other operators with additional simplifications in some cases from factorization of the unitary transformation operator.

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolving Nuclear Many-Body Forces with the Similarity Renormalization Group

In recent years, the Similarity Renormalization Group has provided a powerful and versatile means to soften interactions for ab initio nuclear calculations. The substantial contribution of both induced and initial three-body forces to the nuclear interaction has required the consistent evolution of free-space Hamiltonians in the three-particle space. We present the most recent progress on this work, extending the calculational capability to the p-shell nuclei and showing that the hierarchy of induced many-body forces is consistent with previous estimates. Calculations over a range of the flow parameter for 6Li, including fully evolved NN+3N interactions, show moderate contributions due to induced four-body forces and display the same improved convergence properties as in lighter nuclei. A systematic analysis provides further evidence that the hierarchy of many-body forces is preserved.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

From low-momentum interactions to nuclear structure

We present an overview of low-momentum two-nucleon and many-body interactions and their use in calculations of nuclei and infinite matter. The softening of phenomenological and effective field theory (EFT) potentials by renormalization group (RG) transformations that decouple low and high momenta leads to greatly enhanced convergence in few- and many-body systems while maintaining a decreasing hierarchy of many-body forces. This review surveys the RG-based technology and results, discusses the connections to chiral EFT, and clarifies various misconceptions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How should one formulate, extract, and interpret `non-observables' for nuclei?

Nuclear observables such as binding energies and cross sections can be directly measured. Other physically useful quantities, such as spectroscopic factors, are related to measured quantities by a convolution whose decomposition is not unique. Can a framework for these nuclear structure `non-observables' be formulated systematically so that they can be extracted from experiment with known uncertainties and calculated with consistent theory? Parton distribution functions in hadrons serve as an illustrative example of how this can be done. A systematic framework is also needed to address questions of interpretation, such as whether short-range correlations are important for nuclear structure.

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

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.

0 citations0 reviewsNo code linkedOpen access