Engineering topics

Pieper, Steven C.

Publications and source records attributed to Pieper, Steven C..

10 records

DOE OSTI2022

Nuclear Physics Exascale Requirements Review: An Office of Science Review sponsored jointly by Advanced Scientific Computing Research and Nuclear Physics, June 15 - 17, 2016, Gaithersburg, Maryland

Imagine being able to predict — with unprecedented accuracy and precision — the structure of the proton and neutron, and the forces between them, directly from the dynamics of quarks and gluons, and then using this information in calculations of the structure and reactions of atomic nuclei and of the properties of dense neutron stars (NSs). Also imagine discovering new and exotic states of matter, and new laws of nature, by being able to collect more experimental data than we dream possible today, analyzing it in real time to feed back into an experiment, and curating the data with full tracking capabilities and with fully distributed data mining capabilities. Making this vision a reality would improve basic scientific understanding, enabling us to precisely calculate, for example, the spectrum of gravity waves emitted during NS coalescence, and would have important societal applications in nuclear energy research, stockpile stewardship, and other areas. This review presents the components and characteristics of the exascale computing ecosystems necessary to realize this vision.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Impact of O 16 ( γ , α ) C 12 measurements on the C 12 ( α , γ ) O 16 astrophysical reaction rate

The C 12 ( α , γ ) O 16 reaction, an essentialt component of stellar helium burning, plays a key role in nuclear astrophysics. It has direct impact on the evolution and final state of massive stars, while also influencing the elemental abundances resulting from nucleosynthesis in such stars. Providing a reliable estimate for the energy dependence of this reaction at stellar helium burning temperatures has been a major goal for the field. In this work, we study the role of potential new measurements of the inverse reaction, O 16 ( γ , α ) C 12 , in reducing the overall uncertainty. A multilevel R -matrix analysis is used to make extrapolations of the astrophysical S factor for this reaction to the stellar energy of 300 keV. The statistical precision of the S-factor extrapolation is determined by performing multiple fits to existing E1 and E2 ground-state capture data, including the impact of possible future measurements of the O 16 ( γ , α ) C 12 reaction. Specifically, we consider a proposed Jefferson Laboratory (JLab) experiment that will make use of a high-intensity low-energy bremsstrahlung beam that impinges on an oxygen-rich single-fluid bubble chamber in order to measure the total cross section for the inverse reaction. The importance of low-energy data as well as high-precision data is investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Local chiral interactions and magnetic structure of few-nucleon systems

The magnetic form factors of 2 H, 3 H, and 3 He, deuteron photodisintegration cross sections at low energies, and deuteron threshold electrodisintegration cross sections at backward angles in a wide range of momentum transfers are calculated with the chiral two-nucleon (and three-nucleon) interactions including Δ intermediate states that have recently been constructed in configuration space. The A=3 wave functions are obtained from hyperspherical-harmonics solutions of the Schrödinger equation. The electromagnetic current includes one- and two-body terms, the latter induced by one- and two-pion exchange (OPE and TPE, respectively) mechanisms and contact interactions. The contributions associated with Δ intermediate states are only retained at the OPE level and are neglected in TPE loop (tree-level) corrections to two-body (three-body) current operators. Expressions for these currents are derived and regularized in configuration space for consistency with the interactions. The low-energy constants that enter the contact currents are determined by reproducing the magnetic moments of these few-nucleon systems. The predicted form factors and deuteron electrodisintegration cross section are in excellent agreement with experiment for momentum transfers up to 2–3fm –1 . However, the experimental values for the deuteron photodisintegration cross section are consistently underestimated by theory, unless use is made of the Siegert form of the electric dipole transition operator. In conclusion, a complete analysis of the results is provided, including the clarification of the origin of the aforementioned discrepancy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Local chiral interactions, the tritium Gamow-Teller matrix element, and the three-nucleon contact term

The Gamow-Teller (GT) matrix element contributing to tritium β decay is calculated with trinucleon wave functions obtained from hyperspherical-harmonics solutions of the Schrödinger equation with the chiral two- and three-nucleon interactions including Δ intermediate states that have recently been constructed in configuration space. Predictions up to next-to-next-to-next-to-leading order (N3LO) in the chiral expansion of the axial current (with Δ's) overestimate the empirical value by 1–4%. By exploiting the relation between the low-energy constant (LEC) in the contact three-nucleon interaction and two-body axial current, we provide new determinations of the LECs c D and c E that characterize this interaction by fitting the trinucleon binding energy and tritium GT matrix element. Some of the implications that the resulting models of three-nucleon interactions have on the spectra of light nuclei and the equation of state of neutron matter are briefly discussed. We also provide a partial analysis, which ignores Δ's, of the contributions due to loop corrections in the axial current at next-to-next-to-next-to-next-to-leading order (N4LO). Lastly, explicit expressions for the axial current up to N4LO have been derived in configuration space, which other researchers in the field may find useful.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Quantum Monte Carlo calculation of neutral-current ν - C 12 inclusive quasielastic scattering

Quasielastic neutrino scattering is an important aspect of the experimental program to study fundamental neutrino properties including neutrino masses, mixing angles, the mass hierarchy and CP-violating phase. Proper interpretation of the experiments requires reliable theoretical calculations of neutrino-nucleus scattering. In this paper we present calculations of response functions and cross sections by neutral-current scattering of neutrinos off $^{12}$C. These calculations are based on realistic treatments of nuclear interactions and currents, the latter including the axial, vector, and vector-axial interference terms crucial for determining the difference between neutrino and anti-neutrino scattering and the CP-violating phase. Here in this paper, we find that the strength and energy-dependence of two-nucleon processes induced by correlation effects and interaction currents are crucial in providing the most accurate description of neutrino-nucleus scattering in the quasielastic regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Quantum Monte Carlo calculations of weak transitions in A = 6 – 10 nuclei

{\it Ab initio} calculations of the Gamow-Teller (GT) matrix elements in the $\beta$ decays of $^6$He and $^{10}$C and electron captures in $^7$Be are carried out using both variational and Green's function Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange phenomenology or chiral effective field theory. The agreement with experimental data is excellent for the electron captures in $^7$Be, while theory overestimates the $^6$He and $^{10}$C data by $\sim 2\%$ and $\sim 10\%$, respectively. We show that for these systems correlations in the nuclear wave functions are crucial to explain the data, while many-body currents increase by $\sim 2$--$3\%$ the one-body GT contributions. These findings suggest that the longstanding $g_A$-problem, {\it i.e.}, the systematic overprediction ($\sim 20 \%$ in $A\le 18$ nuclei) of GT matrix elements in shell-model calculations, may be resolved, at least partially, by correlation effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Light-Nuclei Spectra from Chiral Dynamics

In recent years local chiral interactions have been derived and implemented in quantum Monte Carlo methods in order to test to what extent the chiral effective field theory framework impacts our knowledge of few- and many-body systems. Here in this Letter, we present Green’s function Monte Carlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4–12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Quantum Monte Carlo Calculations in Nuclear Theory

To a remarkable extent, atomic nuclei can be described as collections of point-like particles whose dynamics are dictated by a nonrelativistic Hamiltonian involving two- and three-nucleon potentials.Ab initio approaches are aimed at solving the many-body Schr$ \overline{o} $dinger equation associated with the nuclear Hamiltonian, which is a highly non-trivial problem because of the nonperturbative nature and the strong spin-isospin dependence of the nuclear forces. Within this framework it is possible to disentangle the theoretical uncertainty coming from modeling the nuclear potential and currents from that due to the approximations that are usually inherent in many-body techniques.Our many-body method of choice is quantum Monte Carlo (QMC), in particular Green’s function Monte Carlo (GFMC), which allows solving the nuclear Schrodinger equation with the required 1%accuracy level for both the ground- and the low-lying excited states of A ≤ 12 (A is the number of nucleons) nuclei. Since 2000 we have been achieving excellent agreement with experiment for nuclei with increasing number of nucleons by using the Argonne V18 two-nucleon interaction and the Illinois three-nucleon interactions (AV18+IL7). During this period many other groups have been developing chiral effective field theory (χEFT) potentials; these are generally nonlocal and notsuitable for QMC methods. We and collaborators have recently developed a local χEFT potential including the Δ-isobar degrees of freedom which can be used with QMC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Variational Calculation of the Ground State of Closed-Shell Nuclei Up to $A$ = 40

Variational calculations of ground-state properties of 4He, 16O and 40Ca are carried out employing realistic phenomenological two- and three-nucleon potentials. The trial wave function includes twoand three-body correlations acting on a product of single-particle determinants. Expectation values are evaluated with a cluster expansion for the spin-isospin dependent correlations considering up to five-body cluster terms. The optimal wave function is obtained by minimizing the energy expectation value over a set of up to 20 parameters by means of a nonlinear optimization library. We present results for the binding energy, charge radius, point density, single-nucleon momentum distribution, charge form factor, and Coulomb sum rule. We find that the employed three-nucleon interaction becomes repulsive for A ≥ 16. In 16O the inclusion of such a force provides a better description of the properties of the nucleus. In 40Ca instead, the repulsive behavior of the three-body interaction fails to reproduce experimental data for the charge radius and the charge form factor. We find that the high-momentum region of the momentum distributions, determined by the short-range terms of nuclear correlations, exhibit a universal behavior independent of the particular nucleus. The comparison of the Coulomb sum rules for 4He, 16O, and 40Ca reported in this work will help elucidate in-medium modifications of the nucleon form factors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Evolution of a minimal parallel programming model

Here, we take a historical approach to our presentation of self-scheduled task parallelism, a programming model with its origins in early irregular and nondeterministic computations encountered in automated theorem proving and logic programming. We show how an extremely simple task model has evolved into a system, asynchronous dynamic load balancing (ADLB), and a scalable implementation capable of supporting sophisticated applications on today’s (and tomorrow’s) largest supercomputers; and we illustrate the use of ADLB with a Green’s function Monte Carlo application, a modern, mature nuclear physics code in production use. Our lesson is that by surrendering a certain amount of generality and thus applicability, a minimal programming model (in terms of its basic concepts and the size of its application programmer interface) can achieve extreme scalability without introducing complexity.

97 MATHEMATICS AND COMPUTING↗