Source author record

Steven C. Pieper

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

24 published item(s)

preprint2016arXiv

Electromagnetic response of 12C: a first-principles calculation

The longitudinal and transverse electromagnetic response functions of $^{12}$C are computed in a "first-principles" Green's function Monte Carlo calculation, based on realistic two- and three-nucleon interactions and associated one- and two-body currents. We find excellent agreement between theory and experiment and, in particular, no evidence for the quenching of measured versus calculated longitudinal response. This is further corroborated by a re-analysis of the Coulomb sum rule, in which the contributions from the low-lying $J^π\,$=$\, 2^+$, $0^+_2$ (Hoyle), and $4^+$ states in $^{12}$C are accounted for explicitly in evaluating the total inelastic strength.

preprint2016arXiv

Local chiral potentials and the structure of light nuclei

We present fully local versions of the minimally non-local nucleon-nucleon potentials constructed in a previous paper [M.\ Piarulli {\it et al.}, Phys.\ Rev.\ C {\bf 91}, 024003 (2015)], and use them in hypersperical-harmonics and quantum Monte Carlo calculations of ground and excited states of $^3$H, $^3$He, $^4$He, $^6$He, and $^6$Li nuclei. The long-range part of these local potentials includes one- and two-pion exchange contributions without and with $Δ$-isobars in the intermediate states up to order $Q^3$ ($Q$ denotes generically the low momentum scale) in the chiral expansion, while the short-range part consists of contact interactions up to order $Q^4$. The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database in two different ranges of laboratory energies, either 0--125 MeV or 0--200 MeV, and to the deuteron binding energy and $nn$ singlet scattering length. Fits to these data are performed for three models characterized by long- and short-range cutoffs, $R_{\rm L}$ and $R_{\rm S}$ respectively, ranging from $(R_{\rm L},R_{\rm S})=(1.2,0.8)$ fm down to $(0.8,0.6)$ fm. The long-range (short-range) cutoff regularizes the one- and two-pion exchange (contact) part of the potential.

preprint2015arXiv

Ab initio calculation of the electromagnetic and neutral-weak response functions of 4He and 12C

Precise measurement of neutrino oscillations, and hence the determination of their masses demands a quantitative understanding of neutrino-nucleus interactions. To this aim, two-body meson-exchange currents have to be accounted for along within realistic models of nuclear dynamics. We summarize our progresses towards the construction of a consistent framework, based on quantum Monte Carlo methods and on the spectral function approach, that can be exploited to accurately describe neutrino interactions with atomic nuclei over the broad kinematical region covered by neutrino experiments.

preprint2015arXiv

Electromagnetic and neutral-weak response functions of 4He and 12C

Ab initio calculations of the quasi-elastic electromagnetic and neutral-weak response functions of 4He and 12C are carried out for the first time. They are based on a realistic approach to nuclear dynamics, in which the strong interactions are described by two- and three-nucleon potentials and the electroweak interactions with external fields include one- and two-body terms. The Green's function Monte Carlo method is used to calculate directly the Laplace transforms of the response functions, and maximum-entropy techniques are employed to invert the resulting imaginary-time correlation functions with associated statistical errors. The theoretical results, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasi-elastic to the dip region and beyond. These findings challenge the conventional picture of quasi-elastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

preprint2015arXiv

Quantum Monte Carlo methods for nuclear physics

Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

preprint2014arXiv

Neutral weak current two-body contributions in inclusive scattering from $^{12}$C

An ab initio calculation of the sum rules of the neutral weak response functions in $^{12}$C is reported, based on a realistic Hamiltonian, including two- and three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (space-like) neutral current are largest and that a significant portion ($\simeq 30$ %) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasi-elastic scattering data on nuclei.

preprint2014arXiv

Quantum Monte Carlo calculations of electromagnetic transitions in 8Be with meson-exchange currents derived from chiral effective field theory

We report quantum Monte Carlo calculations of electromagnetic transitions in 8Be. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon potentials are used to generate the ground state and nine excited states, with energies that are in excellent agreement with experiment. A dozen M1 and eight E2 transition matrix elements between these states are then evaluated. The E2 matrix elements are computed only in impulse approximation, with those transitions from broad resonant states requiring special treatment. The M1 matrix elements include two-body meson-exchange currents derived from chiral effective field theory, which typically contribute 20--30% of the total expectation value. Many of the transitions are between isospin-mixed states; the calculations are performed for isospin-pure states and then combined with empirical mixing coefficients to compare to experiment. Alternate mixings are also explored. In general, we find that transitions between states that have the same dominant spatial symmetry are in reasonable agreement with experiment, but those transitions between different spatial symmetries are often underpredicted.

preprint2013arXiv

Charge form factor and sum rules of electromagnetic response functions in Carbon-12

An "ab initio" calculation of the Carbon-12 elastic form factor, and sum rules of longitudinal and transverse response functions measured in inclusive (e,e') scattering, is reported, based on realistic nuclear potentials and electromagnetic currents. The longitudinal elastic form factor and sum rule are found to be in satisfactory agreement with available experimental data. A direct comparison between theory and experiment is difficult for the transverse sum rule. However, it is shown that the calculated one has large contributions from two-body currents, indicating that these mechanisms lead to a significant enhancement of the quasi-elastic transverse response. This fact may have implications for the anomaly observed in recent neutrino quasi-elastic charge-changing scattering data off Carbon-12.

preprint2013arXiv

Charge-symmetry breaking forces and isospin mixing in 8Be

We report Green's function Monte Carlo calculations of isospin-mixing (IM) matrix elements for the 2+, 1+, and 3+ T=0,1 pairs of states at 16--19 MeV excitation in 8Be. The realistic Argonne v18 (AV18) two-nucleon and Illinois-7 three-nucleon potentials are used to generate the nuclear wave functions. Contributions from the full electromagnetic interaction and strong class III charge-symmetry-breaking (CSB) components of the AV18 potential are evaluated. We also examine two theoretically more complete CSB potentials based on rho-omega mixing, tuned to give the same neutron-neutron scattering length as AV18. The contribution of these different CSB potentials to the 3H-3He, 7Li-7Be, and 8Li-8B isovector energy differences is evaluated and reasonable agreement with experiment is obtained. Finally, for the 8Be IM calculation we add the small class IV CSB terms coming from one-photon, one-pion, and one-rho exchange, as well as rho-omega mixing. The expectation values of the three CSB models vary by up to 20% in the isovector energy differences, but only by 10% or less in the IM matrix element. The total matrix element gives 85--90% of the experimental IM value of -145 keV for the 2+ doublet, with about two thirds coming from the Coulomb interaction. We also report the IM matrix element to the first 2+ state at 3 MeV excitation, which is the final state for various tests of the Standard Model for beta-decay.

preprint2013arXiv

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists can lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science results to illustrate this interplay.

preprint2013arXiv

GFMC calculations of electromagnetic moments and M1 transitions in $A\leq 9$ nuclei

We present recent Green's function Monte Carlo calculations of magnetic moments and M1 transitions in $A \leq 9$ nuclei, which include corrections arising from two-body meson-exchange electromagnetic currents. Two-body effects provide significant corrections to the calculated observables, bringing them in excellent agreement with the experimental data. In particular, we find that two-body corrections are especially large in the A=9, $T=3/2$ systems, in which they account for up to $\sim 20%$ ($\sim 40 %$) of the total predicted value for the $^9$Li ($^9$C) magnetic moment.

preprint2013arXiv

Laser Probing of Neutron-Rich Nuclei in Light Atoms

The neutron-rich 6He and 8He isotopes exhibit an exotic nuclear structure that consists of a tightly bound 4He-like core with additional neutrons orbiting at a relatively large distance, forming a halo. Recent experimental efforts have succeeded in laser trapping and cooling these short-lived, rare helium atoms, and have measured the atomic isotope shifts along the 4He-6He-8He chain by performing laser spectroscopy on individual trapped atoms. Meanwhile, the few-electron atomic structure theory, including relativistic and QED corrections, has reached a comparable degree of accuracy in the calculation of the isotope shifts. In parallel efforts, also by measuring atomic isotope shifts, the nuclear charge radii of lithium and beryllium isotopes have been studied. The techniques employed were resonance ionization spectroscopy on neutral, thermal lithium atoms and collinear laser spectroscopy on beryllium ions. Combining advances in both atomic theory and laser spectroscopy, the charge radii of these light halo nuclei have now been determined for the first time independent of nuclear structure models. The results are compared with the values predicted by a number of nuclear structure calculations, and are used to guide our understanding of the nuclear forces in the extremely neutron-rich environment.

preprint2013arXiv

Nucleon and nucleon-pair momentum distributions in A <= 12 nuclei

We report variational Monte Carlo calculations of single-nucleon momentum distributions for A <= 12 nuclei and nucleon-pair and nucleon-cluster momentum distributions for A <= 8. The wave functions have been generated for a Hamiltonian containing the Argonne v18 two-nucleon and Urbana X three-nucleon potentials. The single-nucleon and nucleon-pair momentum distributions exhibit universal features attributable to the one-pion-exchange tensor interaction. The single-nucleon distributions are broken down into proton and neutron components and spin-up and spin-down components where appropriate. The nucleon-pair momentum distributions are given separately for pp and pn pairs. The nucleon-cluster momentum distributions include dp in 3He, tp and dd in 4He, alpha-d in 6Li, alpha-t in 7Li, and alpha-alpha in 8Be. Detailed tables are provided on-line for download.

preprint2013arXiv

Properties of trapped neutrons interacting with realistic nuclear Hamiltonians

We calculate properties of neutron drops in external potentials using both quantum Monte Carlo and no-core full configuration techniques. The properties of the external wells are varied to examine different density profiles. We compare neutron drop results given by a selection of nuclear Hamiltonians, including realistic two-body interactions as well as several three-body forces. We compute a range of properties for the neutron drops: ground-state energies, spin-orbit splittings, excitation energies, radial densities and rms radii. We compare the equations of state for neutron matter for several of these Hamiltonians. Our results can be used as benchmarks to test other many-body techniques, and to constrain properties of energy-density functionals.

preprint2013arXiv

Quantum Monte Carlo calculations of electromagnetic moments and transitions in A <= 9 nuclei with meson-exchange currents derived from chiral effective field theory

Quantum Monte Carlo calculations of electromagnetic moments and transitions are reported for A <= 9 nuclei. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon potentials are used to generate the nuclear wave functions. Contributions of two-body meson-exchange current (MEC) operators are included for magnetic moments and M1 transitions. The MEC operators have been derived in both a standard nuclear physics approach and a chiral effective field theory formulation with pions and nucleons including up to one-loop corrections. The two-body MEC contributions provide significant corrections and lead to very good agreement with experiment. Their effect is particularly pronounced in the A=9, T=3/2 systems, in which they provide up to 20% (40 %) of the total predicted value for the 9Li (9C) magnetic moment.

preprint2012arXiv

On the lifetime of the 2+ state in 10C

The lifetime of the J=2+ state in 10C was measured using the Doppler Shift Attenuation Method following the inverse kinematics p(10B,n)10C reaction at 95 MeV. The 2+ state, at 3354 keV, has tau = 219\pm(7)stat \pm(10)sys fs corresponding to a B(E2) # of 8.8(3) e2fm4. This measurement,combined with that recently determined for 10Be (9.2(3) e2fm4), provides a unique challenge to abinitio calculations, testing the structure of these states, including the isospin symmetry of the wave functions. Quantum Monte Carlo calculations using realistic two- and three-nucleon Hamiltonians that reproduce the 10Be B(E2) value generally predict a larger 10C B(E2) probability but with considerable sensitivity to the admixture of different spatial symmetry components in the wave functions, and to the three-nucleon potential used.

preprint2012arXiv

Quantum Monte Carlo calculations of magnetic moments and M1 transitions in $A \le 9$ nuclei

We present Quantum Monte Carlo calculations of magnetic moments and M1 transitions in $A\le 9$ nuclei which take into account contributions of two-body electromagnetic currents. The Hamiltonian utilized to generate the nuclear wave functions includes the realistic Argonne-{\it v}$_{18}$ two-nucleon and the Illinois-7 three-nucleon interactions. The nuclear two-body electromagnetic currents are derived from a pionful chiral effective field theory including up to one-loop corrections. These currents involve unknown Low Energy Constants which have been fixed so as to reproduce a number of experimental data for the two- and three-nucleon systems, such as $np$ phase shifts and deuteron, triton, and $^3$He magnetic moments. This preliminary study shows that two-body contributions provide significant corrections which are crucial to bring the theory in agreement with the experimental data in both magnetic moments and M1 transitions.

preprint2011arXiv

Cold neutrons trapped in external fields

The properties of inhomogeneous neutron matter are crucial to the physics of neutron-rich nuclei and the crust of neutron stars. Advances in computational techniques now allow us to accurately determine the binding energies and densities of many neutrons interacting via realistic microscopic interactions and confined in external fields. We perform calculations for different external fields and across several shells to place important constraints on inhomogeneous neutron matter, and hence the large isospin limit of the nuclear energy density functionals that are used to predict properties of heavy nuclei and neutron star crusts. We find important differences between microscopic calculations and current density functionals; in particular the isovector gradient terms are significantly more repulsive than in traditional models, and the spin-orbit and pairing forces are comparatively weaker.

preprint2011arXiv

Quantum Monte Carlo calculations of spectroscopic overlaps in $A \leq 7$ nuclei

We present Green's function Monte Carlo calculations of spectroscopic overlaps for $A \leq 7$ nuclei. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon interactions are used to generate the nuclear states. The overlap matrix elements are extrapolated from mixed estimates between variational Monte Carlo and Green's function Monte Carlo wave functions. The overlap functions are used to obtain spectroscopic factors and asymptotic normalization coefficients, and they can serve as an input for low-energy reaction calculations.

preprint2009arXiv

Precise electromagnetic tests of ab-initio calculations of light nuclei: States in $^{10}$Be

In order to test {\it ab-initio} calculations of light nuclei, we have remeasured lifetimes in $^{10}$Be using the Doppler Shift Attenuation Method (DSAM) following the $^{7}$Li($^7$Li,$α$)$^{10}$Be reaction at 8 and 10 MeV. The new experiments significantly reduce systematic uncertainties in the DSAM technique. The J$^π$ = $2^+_1$ state at 3.37 MeV has $τ$ = 205${\pm}$(5)$_{stat}\pm$(7)$_{sys}$ fs corresponding to a $B$($E2\downarrow$) of 9.2(3) $e^2$fm$^4$ in broad agreement with many calculations. The J$^π$ = $2^+_2$ state at 5.96 MeV was found to have a $B$($E2\downarrow$) of 0.11(2) $e^2$fm$^4$ and provides a more discriminating test of nuclear models. New Green's Function Monte Carlo (GFMC) calculations for these states and transitions with a number of Hamiltonians are also reported and compared to experiment.

preprint2008arXiv

Quantum Monte Carlo calculations of magnetic moments and M1 transitions in A<=7 nuclei including meson-exchange currents

Green's function Monte Carlo calculations of magnetic moments and M1 transitions including two-body meson-exchange current (MEC) contributions are reported for A<=7 nuclei. The realistic Argonne v18 two-nucleon and Illinois-2 three-nucleon potentials are used to generate the nuclear wave functions. The two-body meson-exchange operators are constructed to satisfy the continuity equation with the Argonne v18 potential. The MEC contributions increase the A=3,7 isovector magnetic moments by 16% and the A=6,7 M1 transition rates by 17--34%, bringing them into very good agreement with the experimental data.

preprint2004arXiv

Neutron Matter: A Superfluid Gas

We review recent progress in the theory of neutron matter with particular emphasis on its superfluid properties. Results of quantum Monte Carlo calculations of simple and realistic models of uniform superfluid neutron gas are discussed along with those of neutrons interacting in a potential well chosen to approximate neutron-rich oxygen isotopes. The properties of dilute superfluid Fermi gases that may be produced in atom traps, and their relations with neutron matter, are illustrated. The density dependence of the effective interaction between neutrons, used to describe neutron-rich systems with the mean field approximations, is also discussed.

preprint2001arXiv

Benchmark Test Calculation of a Four-Nucleon Bound State

In the past, several efficient methods have been developed to solve the Schroedinger equation for four-nucleon bound states accurately. These are the Faddeev-Yakubovsky, the coupled-rearrangement-channel Gaussian-basis variational, the stochastic variational, the hyperspherical variational, the Green's function Monte Carlo, the no-core shell model and the effective interaction hyperspherical harmonic methods. In this article we compare the energy eigenvalue results and some wave function properties using the realistic AV8' NN interaction. The results of all schemes agree very well showing the high accuracy of our present ability to calculate the four-nucleon bound state.