Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear forces”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

SU(3)-guided realistic nucleon-nucleon interactions for large-scale calculations

We examine nucleon-nucleon realistic interactions, based on their SU(3) decomposition to SU(3)-symmetric components. We find that many of these interaction components are negligible, which, in turn, allows us to identify a subset of physically relevant components that are sufficient to describe the structure of low-lying states in 12C and related observables, such as excitation energies, electric quadrupole transitions, and rms radii. Here, we find that paring down the interaction by half or more yields results that practically coincide with the corresponding ab initio calculations with the full interaction. In addition, we show that while various realistic interactions differ in their SU(3) decomposition, their renormalized effective counterparts exhibit a striking similarity and composition that can be linked to dominant nuclear features such as deformation, pairing, clustering, and spin-orbit effect.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Proton decays in 16 Ne and 18 Mg and isospin-symmetry breaking in carbon isotopes and isotones

We report that proton-rich nuclei possess unique properties in the nuclear chart. Due to the presence of both continuum coupling and Coulomb interaction, phenomena such as halos, Thomas-Ehrman shift, and proton emissions can occur. Relevant experimental data are difficult to obtain, so that theoretical calculations are needed to understand nuclei at drip lines and to guide experimentalists. In particular, the 16 Ne and 18 Mg isotopes are supposed to be one-proton and/or two-proton emitting nuclei, but associated experimental data are either incomplete or even unavailable. Consequently, we performed Gamow shell model calculations of carbon isotones bearing A = 15 - 18 . Isospin-symmetry breaking occurring in carbon isotones and isotopes is also discussed. It is hereby shown that the mixed effects of continuum coupling and Coulomb interaction at drip lines generate complex patterns in isospin multiplets. Added to that, it is possible to determine the one-proton and two-proton widths of 16 Ne and 18 Mg . Obtained decay patterns are in agreement with those obtained in previous experimental and theoretical works. Moreover, to our knowledge, this is the first theoretical calculation of binding energy and partial decay widths of 18 Mg in a configuration interaction picture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nuclear spin features relevant to ab initio nucleon-nucleus elastic scattering

Effective interactions for elastic nucleon-nucleus scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations, as well as a consistent treatment of the relevant operators at each order. Previous work using these interactions has shown good agreement with available data. Here, we study the physical relevance of one of these operators, which involves the spin of the struck nucleon, and examine the interpretation of this quantity in a nuclear structure context. Using the framework of the spectator expansion and the underlying framework of the no-core shell model, we calculate and examine spin-projected, one-body momentum distributions required for effective nucleon-nucleus interactions in $J=0$ nuclear states. The calculated spin-projected, one-body momentum distributions for $^4$He, $^6$He, and $^8$He display characteristic behavior based on the occupation of protons and neutrons in single particle levels, with more nucleons of one type yielding momentum distributions with larger values. Additionally, we find this quantity is strongly correlated to the magnetic moment of the $2^+$ excited state in the ground state rotational band for each nucleus considered. In conclusion, we find that spin-projected, one-body momentum distributions can probe the spin content of a $J=0$ wave function. This feature may allow future ab initio nucleon-nucleus scattering studies to inform spin properties of the underlying nucleon-nucleon interactions. The observed correlation to the magnetic moment of excited states illustrates a previously unknown connection between reaction observables such as the analyzing power and structure observables like the magnetic moment.

6 ≤ A ≤ 19↗
DOE OSTI2021

Auxiliary function approach for determining symmetry energy at suprasaturation densities

Nuclear symmetry energy E sym (ρ) at density ρ is normally expanded or simply parameterized as a function of χ = (ρ – ρ 0 )/3ρ 0 in the form of Esym(ρ) ≈ S + Lχ + 2 –1 K sym χ 2 + 6 –1 J sym χ 3 + … using its magnitude S, slope L, curvature K sym and skewness J sym at the saturation density ρ0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. Furthermore, such kind of expansions/parameterizations do not converge at supra-saturation densities where χ is not small enough, hindering an accurate determination of high-density E sym (ρ) even if its characteristic parameters at ρ0 are all well determined by experiments/observations. By expanding the E sym (ρ) in terms of a properly chosen auxiliary function Π sym (χ,Θ sym ) with a parameter Θsym fixed accurately by an experimental E sym (ρ r ) value at a reference density ρr, we show that the shortcomings of the χ-expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym (ρ). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ-order contributions and converges to the same E sym (ρ) much faster than the conventional χ-expansion at densities ≲3ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Nucleon-pair coupling scheme in Elliott's SU(3) model

Elliott's SU(3) model is at the basis of the shell-model description of rotational motion in atomic nuclei. Here we demonstrate that SU(3) symmetry can be realized in a truncated shell-model space if constructed in terms of a sufficient number of collective S, D, G,...pairs (i.e., with angular momentum zero, two, four,...) and if the structure of the pairs is optimally determined either by a conjugate-gradient minimization method or from a Hartree-Fock intrinsic state. We illustrate the procedure for six protons and six neutrons in the pf (sdg) shell and exactly reproduce the level energies and electric quadrupole properties of the ground-state rotational band with SDG (SDGI) pairs. The SD-pair approximation without significant renormalization, on the other hand, cannot describe the full SU(3) collectivity. A mapping from Elliott's fermionic SU(3) model to systems with s, d, g,... bosons provides insight into the existence of a decoupled collective subspace in terms of S, D, G,... pairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Extracting the number of short-range correlated nucleon pairs from inclusive electron scattering data

The extraction of the relative abundances of short-range correlated (SRC) nucleon pairs from inclusive electron scattering is studied using the generalized contact formalism (GCF) with several nuclear interaction models. GCF calculations can reproduce the observed scaling of the cross-section ratios for nuclei relative to deuterium at high xB and large Q 2 , a 2 =(σ A /A)/(σ d /2). In the nonrelativistic instant-form formulation, the calculation is very sensitive to the model parameters and only reproduces the data using parameters that are inconsistent with ab initio many-body calculations. Using a light-cone GCF formulation significantly decreases this sensitivity and improves the agreement with ab initio calculations. The ratio of similar mass isotopes, such as 40 Ca and 48 Ca, should be sensitive to the nuclear asymmetry dependence of SRCs, but is found to also be sensitive to low-energy nuclear structure. Thus the empirical association of SRC pair abundances with the measured a2 values is only accurate to about 20%. Finally, improving this will require cross-section calculations that reproduce the data while properly accounting for both nuclear structure and relativistic effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Evidence against the Efimov effect in 12 C from spectroscopy and astrophysics

Background: The Efimov effect is a universal phenomenon in physics whereby three-body systems are stabilized via the interaction of an unbound two-body subsystems. A hypothetical state in 12 C at 7.458-MeV excitation energy, comprising a loose structure of three α particles in mutual two-body resonance, has been suggested in the literature to correspond to an Efimov state in nuclear physics. The existence of such a state has not been demonstrated experimentally. Purpose: Using a combination of γ spectroscopy, charged-particle spectroscopy, and astrophysical rate calculations allowing for strict limits on the existence of such a state to been established here. Method: Using the combined data sets from two recent experiments, one with the TexAT (Texas Active Target) TPC (Time Projection Chamber) to measure α decay and the other with Gammasphere to measure γ decay of states in 12 C populated by 12 N and 12 B β decay, respectively, we achieve high sensitivity to states in close proximity to the α threshold in 12 C. Results: No evidence of a state at 7.458 MeV is seen in either data set. Using a likelihood method, the 95% confidence limit γ -decay branching ratio is determined as a function of the β-decay feeding strength relative to the Hoyle state. In parallel, calculations of the 3α reaction rate show the inclusion of the Efimov corresponds to a large increase in the reaction rate around 5 × 10 7 K. Conclusion: From decay spectroscopy—at the 95% confidence limit, the Efimov state cannot exist at 7.458 MeV with any γ-decay branching ratio unless the β strength is less than 0.7% of the Hoyle state. This limit is evaluated for a range of different excitation energies and the results are not favorable for existence of the hypothetical Efimov state in 12 C . Furthermore, the 3α reaction rate with the inclusion of a state between 7.43 and 7.53 MeV exceeds the rate required for stars to undergo the red giant phase.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Reduction of the neutron imaginary potential off the stability line and its possible impact on neutron capture rates

The effect of the isovector imaginary optical potential has been studied experimentally by using the particle-evaporation technique for the 11 B + 48 Ca reaction with a 21.8 MeV 11 B beam. Spectra of neutron, proton, and α particles emitted from the neutron-rich compound nucleus 59 Mn have been measured and analyzed with traditional optical model potentials with their original parametrizations as well as with adjusted isovector imaginary components. It is shown that the isovector component of the imaginary potential is indispensable in the reproduction of proton and α-particle yields from this reaction and even needs to be enhanced compared with the suggestions of the original model parametrizations. Furthermore, this can lead to important consequences for astrophysical reaction-rate calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Singular value decomposition and similarity renormalization group evolution of nuclear interactions

One of the main challenges for ab initio nuclear many-body theory is the growth of computational and storage costs as calculations are extended to heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform Singular Value Decompositions of nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization scheme and resolution scales. We develop and implement the Similarity Renormalization Group (SRG) evolution of the factorized interaction, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves, and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Lastly, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the In-Medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Ground-state properties of light 4 n self-conjugate nuclei in ab initio no-core Monte Carlo shell model calculations with nonlocal N N interactions

We report J π = 0 + ground-state energies and point-proton radii of 4 He, 8 Be, 12 C, 16 O, and 20 Ne nuclei calculated by the ab initio no-core Monte Carlo shell model with the JISP16 and Daejeon16 nonlocal NN interactions. Ground-state energies are obtained in the basis spaces up to seven oscillator shells (N shell =7) with several oscillator energies (ℏω) around the optimal oscillator energy for the convergence of ground-state energies. These energy eigenvalues are extrapolated to obtain estimates of converged ground-state energies in each basis space using energy variances of computed energy eigenvalues. We further extrapolate these energy-variance-extrapolated energies obtained in the finite basis spaces to infinite basis-space results with an empirical exponential form. This form features a dependence on the basis-space size but is independent of the value of ℏω used for the harmonic-oscillator basis functions. Point-proton radii for these states of atomic nuclei are also calculated following techniques employed for the energies. From these results, it is found that the Daejeon16 NN interaction provides good agreement with experimental data up to approximately 16 O, while the JISP16 NN interaction provides good agreement with experimental data up to approximately 12 C. Beyond these nuclei, the interactions produce overbinding accompanied by radii that are too small. These findings suggest and encourage further revisions of nonlocal NN interactions towards the investigation of nuclear structure in heavier-mass regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Application of an efficient generator-coordinate subspace-selection algorithm to neutrinoless double- β decay

The generator coordinate method begins with the variational construction of a set of nonorthogonal mean-field states that span a subspace of the full many-body Hilbert space. These states are then often projected onto states with good quantum numbers to restore symmetries, leading to a set with members that can be similar to one another, and it is sometimes possible to reduce this set without greatly affecting results. Here, we propose a greedy algorithm that we call the energy-transition-orthogonality procedure (ENTROP) to select subsets of important states. As applied here, the approach selects on the basis of diagonal energy, orthogonality, and contribution to the matrix element that governs neutrinoless double-β decay. We present both shell-model and preliminary ab initio calculations of this matrix element for the decay of 76Ge, with quadrupole deformation parameters and the isoscalar pairing strength as generator coordinates. ENTROP converges quickly, reducing significantly the number of basis states needed for an accurate calculation.

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

Complete set of bound negative-parity states in the neutron-rich nucleus 18 N

High-resolution γ -ray spectroscopy of N 18 is performed with the Advanced GAmma Tracking Array, following deep-inelastic processes induced by an O 18 beam on a Ta 181 target. Six states are newly identified, which together with the three known excitations exhaust all negative-parity excited states expected in N 18 below the neutron threshold. Spin and parities are proposed for all located states on the basis of decay branchings and comparison with large-scale shell-model calculations performed in the p-sd space, with the YSOX interaction. Of particular interest is the location of the 0 1 - and 1 2 - excitations, which provide strong constrains for cross-shell p-sd matrix elements based on realistic interactions and help to simultaneously reproduce the ground and first-excited states in N 16 and N 18 , for the first time. Overall, understanding the N 18 structure may also have significant impact on neutron-capture cross-section calculations in r-process modeling including light neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Isospin mixing and the cubic isobaric multiplet mass equation in the lowest T=2, A=32 quintet

The isobaric multiplet mass equation (IMME) is known to break down in the first T=2,A=32 isospin quintet. In this work we combine high-resolution experimental data with state-of-the-art shell-model calculations to investigate isospin mixing as a possible cause for this violation. The experimental data are used to validate isospin-mixing matrix elements calculated with newly developed shell-model Hamiltonians. Our analysis shows that isospin mixing with nonanalog T=1 states contributes to the IMME breakdown, making the requirement of an anomalous cubic term inevitable for the multiplet.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Slope parameter of the symmetry energy and the structure of three-particle interactions in nuclear matter

In the first part of this paper, we present a study of the symmetry energy (a s ) and its slope parameter (L) for nuclear matter in the framework of the Fermi liquid theory of Landau and Migdal. We derive an exact relation between a s and L, which involves the nucleon effective masses and three-particle Landau-Migdal parameters. We present simple estimates which suggest that there are two main mechanisms to explain the empirical values of L: The proton-neutron effective mass difference in isospin asymmetric matter and the ℓ = 0 moment of the isovector in-medium three-particle scattering amplitude. In the second part of this paper, we discuss the general structure of three-particle interactions in nuclear matter in the framework of the Fermi liquid theory. The connections to the Bethe-Brueckner-Goldstone theory and other approaches are also discussed. Finally, we show explicitly how the first few terms in the Faddeev series, together with medium induced three-particle interactions, emerge naturally in the Fermi liquid theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗