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

R. J. Furnstahl

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

ResearcherUnclaimed source record
nucl-thhep-phnucl-exphysics.data-anhep-latastro-ph.HEastro-ph.COastro-ph.GAastro-ph.SRcond-matcond-mat.mtrl-scicond-mat.otherphysics.atom-phphysics.comp-ph

Catalog footprint

What is connected

47works
14topics
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

47 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
preprint2016arXiv

Bayesian parameter estimation for effective field theories

We present procedures based on Bayesian statistics for estimating, from data, the parameters of effective field theories (EFTs). The extraction of low-energy constants (LECs) is guided by theoretical expectations in a quantifiable way through the specification of Bayesian priors. A prior for natural-sized LECs reduces the possibility of overfitting, and leads to a consistent accounting of different sources of uncertainty. A set of diagnostic tools are developed that analyze the fit and ensure that the priors do not bias the EFT parameter estimation. The procedures are illustrated using representative model problems, including the extraction of LECs for the nucleon mass expansion in SU(2) chiral perturbation theory from synthetic lattice data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Regulator Artifacts in Uniform Matter for Chiral Interactions

Regulator functions applied to two- and three-nucleon forces are a necessary ingredient in many-body calculations based on chiral effective field theory interactions. These interactions have been developed recently with a variety of different cutoff forms, including regulating both the momentum transfer (local) and the relative momentum (nonlocal). While in principle any regulator that suppresses high momentum modes can be employed, in practice artifacts are inevitable in current power counting schemes. Artifacts from particular regulators may cause significant distortions of the physics or may affect many-body convergence rates, so understanding their nature is important. Here we characterize the differences between cutoff effects using uniform matter at Hartree-Fock and second-order in the interaction as a testbed. This provides a clean laboratory to isolate phase-space effects of various regulators on both two- and three-nucleon interactions. We test the normal-ordering approximation for three-nucleon forces in nuclear matter and find that the relative size of the residual 3N contributions is sensitive to the employed regularization scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Few-nucleon systems with state-of-the-art chiral nucleon-nucleon forces

We apply improved nucleon-nucleon potentials up to fifth order in chiral effective field theory, along with a new analysis of the theoretical truncation errors, to study nucleon-deuteron (Nd) scattering and selected low-energy observables in 3H, 4He, and 6Li. Calculations beyond second order differ from experiment well outside the range of quantified uncertainties, providing truly unambiguous evidence for missing three-nucleon forces within the employed framework. The sizes of the required three-nucleon force contributions agree well with expectations based on Weinberg's power counting. We identify the energy range in elastic Nd scattering best suited to study three-nucleon force effects and estimate the achievable accuracy of theoretical predictions for various observables.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantifying truncation errors in effective field theory

Bayesian procedures designed to quantify truncation errors in perturbative calculations of quantum chromodynamics observables are adapted to expansions in effective field theory (EFT). In the Bayesian approach, such truncation errors are derived from degree-of-belief (DOB) intervals for EFT predictions. Computation of these intervals requires specification of prior probability distributions ("priors") for the expansion coefficients. By encoding expectations about the naturalness of these coefficients, this framework provides a statistical interpretation of the standard EFT procedure where truncation errors are estimated using the order-by-order convergence of the expansion. It also permits exploration of the ways in which such error bars are, and are not, sensitive to assumptions about EFT-coefficient naturalness. We first demonstrate the calculation of Bayesian probability distributions for the EFT truncation error in some representative examples, and then focus on the application of chiral EFT to neutron-proton scattering. Epelbaum, Krebs, and Meißner recently articulated explicit rules for estimating truncation errors in such EFT calculations of few-nucleon-system properties. We find that their basic procedure emerges generically from one class of naturalness priors considered, and that all such priors result in consistent quantitative predictions for 68% DOB intervals. We then explore several methods by which the convergence properties of the EFT for a set of observables may be used to check the statistical consistency of the EFT expansion parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A recipe for EFT uncertainty quantification in nuclear physics

The application of effective field theory (EFT) methods to nuclear systems provides the opportunity to rigorously estimate the uncertainties originating in the nuclear Hamiltonian. Yet this is just one source of uncertainty in the observables predicted by calculations based on nuclear EFTs. We discuss the goals of uncertainty quantification in such calculations and outline a recipe to obtain statistically meaningful error bars for their predictions. We argue that the different sources of theory error can be accounted for within a Bayesian framework, as we illustrate using a toy model.

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
preprint1999arXiv

Parameter Counting in Relativistic Mean-Field Models

Power counting is applied to relativistic mean-field energy functionals to estimate contributions to the energy from individual terms. New estimates for isovector, tensor, and gradient terms in finite nuclei are shown to be consistent with direct, high-quality fits. The estimates establish a hierarchy of model parameters and identify how many parameters are well constrained by bulk nuclear observables. We conclude that four (possibly five) isoscalar, non-gradient parameters, one gradient parameter, and one isovector parameter are well determined by the usual bulk nuclear observables.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Pion-Nucleus Scattering at Medium Energies with Densities from Chiral Effective Field Theories

Recently developed chiral effective field theory models provide excellent descriptions of the bulk characteristics of finite nuclei, but have not been tested with other observables. In this work, densities from both relativistic point-coupling models and mean-field meson models are used in the analysis of meson-nucleus scattering at medium energies. Elastic scattering observables for 790 MeV/$c$ $π^{\pm}$ on $^{208}$Pb are calculated in a relativistic impulse approximation, using the Kemmer-Duffin-Petiau formalism to calculate the $π^{\pm}$ nucleus optical potential.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

QCD Sum Rules and Applications to Nuclear Physics

Applications of QCD sum-rule methods to the physics of nuclei are reviewed, with an emphasis on calculations of baryon self-energies in infinite nuclear matter. The sum-rule approach relates spectral properties of hadrons propagating in the finite-density medium, such as optical potentials for quasinucleons, to matrix elements of QCD composite operators (condensates). The vacuum formalism for QCD sum rules is generalized to finite density, and the strategy and implementation of the approach is discussed. Predictions for baryon self-energies are compared to those suggested by relativistic nuclear physics phenomenology. Sum rules for vector mesons in dense nuclear matter are also considered.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

The Gluon Condensate and Running Coupling of QCD

An expression for the photon condensate in quantum electrodynamics is presented and generalized to deduce a simple relation between the gluon condensate and the running coupling constant of quantum chromodynamics (QCD). Ambiguities in defining the condensates are discussed. The values of the gluon condensate from some Ansätze for the running coupling in the literature are compared with the value determined from QCD sum rules.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Vacuum Contributions in a Chiral Effective Lagrangian for Nuclei

A relativistic hadronic model for nuclear matter and finite nuclei, which incorporates nonlinear chiral symmetry and broken scale invariance, is presented and applied at the one-baryon-loop level to finite nuclei. The model contains an effective light scalar field that is responsible for the mid-range nucleon--nucleon attraction and which has anomalous scaling behavior. One-loop vacuum contributions in this background scalar field at finite density are constrained by low-energy theorems that reflect the broken scale invariance of quantum chromodynamics. A mean-field energy functional for nuclear matter and nuclei is derived that contains small powers of the fields and their derivatives, and the validity of this truncation is discussed. Good fits to the bulk properties of finite nuclei and single-particle spectra are obtained.

0 citations0 reviewsNo code linkedOpen access