Engineering topics

Piarulli, M.

Publications and source records attributed to Piarulli, M..

15 records

DOE OSTI2024

Longitudinal form Factors of A ≤ 10 Nuclei in a Chiral Effective Field Theory Approach

Here, in this work, we present the elastic electron scattering longitudinal form factors of A ≤ 10 nuclei computed in a variational Monte Carlo approach. We employ the Norfolk family of local chiral interactions and a consistent electromagnetic charge operator. Our calculations are compared both to data and past theoretical evaluations. This work represents, to our knowledge, the first accurate calculation of longitudinal form factors using manybody methods based on interacting nucleon degrees of freedom in the 7 ≤ A ≤ 10 mass range. Finally, we identify 9 Be and 10 B as candidate targets for renewed experimental interest, as they exhibit the potential to provide more stringent constraints on the theoretical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Magnetic structure of A ≤ 10 nuclei using the Norfolk nuclear models with quantum Monte Carlo methods

Here we present quantum Monte Carlo calculations of magnetic moments, form factors, and densities of A ≤ 10 nuclei within a chiral effective field theory approach. We use the Norfolk two- and three-body chiral potentials and their consistent electromagnetic one- and two-nucleon current operators. We find that two-body contributions to the magnetic moment can be large (up to ≈ 33% in A = 9 systems). We study the model dependence of these observables and place particular emphasis on investigating their sensitivity to using different cutoffs to regulate the many-nucleon operators. Calculations of elastic magnetic form factors for A ≤ 10 nuclei show excellent agreement with the data out to momentum transfers q ≈ 3 fm -1 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Quantum Monte Carlo Calculations of Magnetic Form Factors in Light Nuclei

Here, we present Quantum Monte Carlo calculations of magnetic form factors in A = 6-10 nuclei, based on Norfolk two- and three-nucleon interactions, and associated one- and two-body electromagnetic currents. Agreement with the available experimental data for 6 Li, 7 Li, 9 Be and 10 B up to values of momentum transfer q ~ 3 fm -1 is achieved when two-nucleon currents are accounted for. We present a set of predictions for the magnetic form factors of 7 Be, 8 Li, 9 Li, and 9 C. In these systems, two body currents account for ~ 40-60% of the total magnetic strength. Measurements in any of these radioactive systems would provide valuable insights on the nuclear magnetic structure emerging from the underlying many-nucleon dynamics. A particularly interesting case is that of 7 Be, as it would enable investigations of the magnetic structure of mirror nuclei

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Quenching of the quantum p -strength in light exotic nuclei

One-nucleon spectroscopic overlaps, their strength or spectroscopic factors (SFs), and nucleon removal cross sections for light nuclei in the mass range A≤12 were evaluated using the fully correlated Quantum Monte Carlo (QMC) wave functions (WFs). Harder (AV18+UX) and softer (NV2+3) bare interactions were used providing consistent results. The SFs were also taken from simple Shell Model (SM) calculations. This structure information was incorporated in the standard three-body Faddeev/Alt-Grassberger-Sandhas reaction formalism to evaluate the (p,pN) cross sections. The results further our understanding of the quenching of the quantum p-shell strength obtained from structure and reactions. We have found the ratios of the total QMC sums of SFs to the SM ones to be uniform and ∼ 3/4. The corresponding ratios of the sums Below Particle Threshold (BPT) of SFs, as well as of total cross sections, deviate strikingly from the uniform trend in some special cases. We find these ratios to be close to unit for 9 Li → 8 Li + n and 9 C → 8 B + p. In contrast, they are strongly reduced for 11 C → 10 C + n and 11 B → 10 Be + p mirror transitions, resulting from the quenching of the strength for low-lying 2 + final state transitions. The QMC theoretical cross section BPT for 11 C(p,pn) is about two times smaller than the experimental data.

(p,pN) reactions↗

DOE OSTI2023

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗

DOE OSTI2023

Densities and momentum distributions in A ≤ 12 nuclei from chiral effective field theory interactions

Current and future electron and neutrino scattering experiments will be greatly aided by a better understanding of the role played by short-range correlations in nuclei. Two-body physics, including nucleon-nucleon correlations and two-body electroweak currents, is required to explain the body of experimental data for both static and dynamical nuclear properties. In this work, we focus on examining nucleon-nucleon correlations from a chiral effective field theory perspective and provide a comprehensive set of new variational Monte Carlo calculations of one- and two-body densities and momentum distributions based on the Norfolk many-body nuclear Hamiltonians for A ≤ 12 systems. Here, online access to detailed tables and figures is available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Muon capture on deuteron using local chiral potentials

The muon capture reaction μ – + d → n + n + ν μ in the doublet hyperfine state is studied using nuclear potentials and consistent currents derived in the chiral effective field theory, which are local and expressed in coordinate space (the so-called Norfolk models). Only the largest contribution due to the 1 S 0 nn scattering state is considered. Particular attention is given to the estimate of theoretical uncertainty, for which four sources have been identified: 1) the model dependence, 2) the chiral-order convergence for the weak nuclear current, 3) the uncertainty in the single-nucleon axial form factor, and 4) the numerical technique adopted to solve the bound and scattering A = 2 systems. This last source of uncertainty has turned out to be essentially negligible. For the 1 S 0 doublet muon capture rate Γ D ( 1 S 0 ), we obtain Γ D ( 1 S 0 ) = 255.8(0.6) (4.4)(2.9) s –1 , where the three errors come from the first three sources of uncertainty. The value for Γ D ( 1 S 0 ) obtained within this local chiral framework is compared with previous calculations and found in very good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Benchmark calculations of infinite neutron matter with realistic two- and three-nucleon potentials

In this study, we present the equation of state of infinite neutron matter as obtained from highly realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne v 18 and the Argonne v' 6 nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ( ρ 0 =0.16 fm ⁻3 ) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and nd doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium β decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Partial muon capture rates in A = 3 and A = 6 nuclei with chiral effective field theory

Searches for neutrinoless-double beta decay rates are crucial in addressing questions within fundamental symmetries and neutrino physics. The rates of these decays depend not only on unknown parameters associated with neutrinos, but also on nuclear properties. In order to reliably extract information about the neutrino, one needs an accurate treatment of the complex many-body dynamics of the nucleus. Neutrinoless-double beta decays take place at momentum transfers on the order of 100 MeV/$c$ and require both nuclear electroweak vector and axial current matrix elements. Muon capture, a process in the same momentum transfer regime, has readily available experimental data to validate these currents. In this work, we present results of {\it ab initio} calculations of partial muon capture rates for $^3$He and $^6$Li nuclei using variational and Green's Function Monte Carlo computational methods. We estimate the impact of the three-nucleon interactions, the cutoffs used to regularize two-nucleon ($2N$) interactions, and the energy range of $2N$ scattering data used to fit these interactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Reply to “Comment on ‘Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter' ”

In their comment to our paper [1], Bertsch and Stroberg [2] provide three criticisms. Two of these concern the interpretation of our dispersive optical model (DOM) results and their relation to the liquid drop model (LDM) parameters. The third criticism focuses on the potential systematic uncertainties on our DOM results associated with missing three-body contributions. Before addressing these critiques, we want to state that the key message of our paper remains whether or not the DOM results agree with the LDM predictions in the nuclear interior. Here, the key point is that the standard determination of the saturation energy from the LDM is not ideal since the total binding energy has a minimal contribution from the core of the nucleus as pointed out in Figs. 1–3 of our paper.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Two- and three-nucleon contact interactions and ground-state energies of light- and medium-mass nuclei

Classes of two-nucleon (2N) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet np scattering length and deuteron binding energy at LO, and np and pp scattering data in the laboratory-energy ranges 0–15 MeV at NLO and 0–25 MeV at N3LO. These interactions are regularized by including two Gaussian cutoffs, one for T = 0 and the other for T = 1 channels. The cutoffs are taken to vary in the ranges R 0 = (1.5 –2.3) fm and R 1 = (1.5 –3.0) fm. The 780 (1100) data points up to 15 (25) MeV energy, primarily differential cross sections, are fitted by the NLO (N3LO) models with a χ 2 / datum about 1.7 or less (well below 1.5), when harder cutoff values are adopted. As a first application, we report results for the binding energies of nuclei with mass numbers A = 3 –6 and 16 obtained with selected LO and NLO 2N models both by themselves as well as in combination with a LO three-nucleon (3N) contact interaction. The latter is characterized by a single low-energy constant that is fixed to reproduce the experimental 3 H binding energy. The inclusion of the 3N interaction largely removes the sensitivity to cutoff variations in the few-nucleon systems and leads to predictions for the 3 He and 4 He binding energies that cluster around 7.8 and 30 MeV, respectively. However, in 16 O this cutoff sensitivity remains rather strong. Finally, predictions at LO only are also reported for medium-mass nuclei with A = 40 , 48, and 90.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Many-body factorization and position–momentum equivalence of nuclear short-range correlations

While mean-field approximations, such as the nuclear shell model, provide a good description of many bulk nuclear properties, they fail to capture the important effects of nucleon–nucleon correlations such as the short-distance and high-momentum components of the nuclear many-body wave function1. Here, we study these components using the effective pair-based generalized contact formalism and ab initio quantum Monte Carlo calculations of nuclei from deuteron to 40 Ca. We observe a universal factorization of the many-body nuclear wave function at short distance into a strongly interacting pair and a weakly interacting residual system. The residual system distribution is consistent with that of an uncorrelated system, showing that short-distance correlation effects are predominantly embedded in two-body correlations. Spin- and isospin-dependent ‘nuclear contact terms’ are extracted in both coordinate and momentum space for different realistic nuclear potentials. The contact coefficient ratio between two different nuclei shows very little dependence on the nuclear interaction model. These findings thus allow extending the application of mean-field approximations to short-range correlated pair formation by showing that the relative abundance of short-range pairs in the nucleus is a long-range (that is, mean field) quantity that is insensitive to the short-distance nature of the nuclear force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Chiral effective field theory calculations of weak transitions in light nuclei

In this work, we report quantum Monte Carlo calculations of weak transitions in A≤10 nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. Furthermore, we find that the contribution from two-body currents is at the 2–3% level, with the exception of matrix elements entering the rates of 8Li, 8B, and 8He β decay. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a 20–30% correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of A=3 to A=10 systems considered here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions

Here, we report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner–Bethe–Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne v 18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗