SEARCH · Search NASA
Results for “Nuclear theory”
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.
Target mass corrections in lepton–nucleus DIS: Theory and applications to nuclear PDFs
Motivated by the wide range of kinematics covered by current and planned deep-inelastic scattering (DIS) facilities, we revisit the formalism, practical implementation, and numerical impact of target mass corrections (TMCs) for DIS on unpolarized nuclear targets. An important aspect is that we only use nuclear and later partonic degrees of freedom, carefully avoiding a picture of the nucleus in terms of nucleons. After establishing that formulae used for individual nucleon targets (p,n), derived in the Operator Product Expansion (OPE) formalism, are indeed applicable to nuclear targets, here we rewrite expressions for nuclear TMCs in terms of re-scaled (or averaged) kinematic variables. As a consequence, we find a representation for nuclear TMCs that is approximately independent of the nuclear target. We go on to construct a single-parameter fit for all nuclear targets that is in good numerical agreement with full computations of TMCs. We discuss in detail qualitative and quantitative differences between nuclear TMCs built in the OPE and the parton model formalisms, as well as give numerical predictions for current and future facilities.
Unexpected Rise in Nuclear Collectivity from Short-Range Physics
We discover a surprising relation between the collective motion of nucleons within atomic nuclei, traditionally understood to be driven by long-range correlations, and short-range nucleon-nucleon interactions. Specifically, we find that quadrupole collectivity in low-lying states of 6 Li and 12 C, calculated with state-of-the-art ab initio techniques, is significantly influenced by two opposing 𝑆-wave contact couplings that subtly alter the surface oscillations of one largely deformed nuclear shape, without changing that shape’s overall contribution within the nucleus. The results offer new insights into the nature of emergent nuclear collectivity and its link to the underlying nucleon-nucleon interaction at short distances.
Computational schemes for the Magnus expansion of the in-medium similarity renormalization group
The in-medium similarity renormalization group (IMSRG) is a popular many-body method used for computations of nuclei. It solves the many-body Schrödinger equation through a continuous unitary transformation of the many-body Hamiltonian. The IMSRG transformation is typically truncated at the normal-ordered two-body level, the IMSRG(2), but recently several approaches have been developed to capture the effects of normal-ordered three-body operators, the IMSRG(3). In particular, a factorized approximation to the IMSRG(3) proposes to capture the leading effects of three-body operators at the same computational cost as the IMSRG(2) approximation. This approach often employs an approximate scheme for solving the IMSRG equations, the so-called hunter-gatherer scheme. In this work, I study the uncertainty associated with this scheme. I find that the hunter-gatherer scheme differs by up to 7MeV for ground-state energies and 0.5MeV for excitation energies from standard IMSRG(2) approaches. These differences are in some cases comparable to the expected size of IMSRG(3) corrections.
Multinucleon structure and dynamics via quantum computing
We propose a framework for computing the structure and dynamics for second-quantized many-nucleon Hamiltonians on quantum computers. We develop an oracle-based Hamiltonian input model that computes the many-nucleon states and nonzero Hamiltonian matrix elements of the many-nucleon system. With our Fock-state based input model, we show how to implement the sparse matrix simulation algorithms to calculate the dynamics of the second-quantized many-nucleon Hamiltonian. Based on the dynamics simulation methods, we also present the methodology for structure calculations of the many-nucleon system. In this work, we provide an explicit circuit design of our input model of the second-quantized Hamiltonian within a direct encoding scheme that maps the occupation of each available single-particle state in the many-nucleon state to the state of specific qubit in a quantum register. Here, we analyze our method and provide the asymptotic cost in computing resources for structure and dynamics calculations of many-nucleon systems. For pedagogical purposes, we demonstrate our input model with two model problems in restricted model spaces.
Quasielastic lepton scattering and back-to-back nucleons in the short-time approximation
Understanding quasielastic electron- and neutrino-scattering from nuclei has taken on new urgency with current and planned neutrino oscillation experiments, and with electron scattering experiments measuring specific final states, such as those involving nucleon pairs in “back-to-back” configurations. Accurate many-body methods are available for calculating the response of light (A ≤ 12) nuclei to electromagnetic and weak probes, but they are computationally intensive and only applicable to the inclusive response. In the present work we introduce a novel approach, based on realistic models of nuclear interactions and currents, to evaluate the short-time (high-energy) inclusive and exclusive response of nuclei. The approach accounts reliably for crucial two-nucleon dynamics, including correlations and currents, and provides information on back-to-back nucleons observed in electron and neutrino scattering experiments. We demonstrate that in the quasielastic regime and at moderate momentum transfers both initial- and final-state correlations, and two-nucleon currents are important for a quantitatively successful description of the inclusive response and final state nucleons. Finally, the approach can be extended to include relativistic—kinematical and dynamical—effects, at least approximately in the two-nucleon sector, and to describe the response in the resonance-excitation region.
Isotopically resolved neutron total cross sections at intermediate energies
The neutron total cross sections σ tot of 16,18 O, 58 64 Ni, 103 Rh, and 112,124 Sn have been measured at the Los Alamos Neutron Science Center from low to intermediate energies (3 ≤ E lab ≤ 450 MeV) by leveraging wave-form-digitizer technology. The σ tot relative differences between isotopes are presented, revealing additional information about the isovector components needed for an accurate optical-model description away from stability. Digitizer-enabled σ tot -measurement techniques are discussed and a series of uncertainty-quantified dispersive optical model (DOM) analyses using these new data is presented, validating the use of the DOM for modeling light systems ( 16,18 O) and systems with open neutron shells ( 58,64 Ni and 112,124 Sn). Here, the valence-nucleon spectroscopic factors extracted for each isotope reaffirm the usefulness of high-energy proton reaction cross sections for characterizing depletion from the mean-field expectation.
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.
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.
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.
Role of quadrupole deformation and continuum effects in the “island of inversion” nuclei 28,29,31 F
The peculiar properties of nuclei in the so-called “island of inversion” around Z = 10 and N = 20 are the focus of current nuclear physics research. Recent studies showed that 28 F has a negative-parity ground state and thus lies within the southern shore of the island of inversion, and 29 F presents a halo structure in its ground state, but it is unclear which effects, such as deformation, shell evolution due to tensor forces, or couplings to the continuum, lead to this situation. Here, we investigate the role of quadrupole deformation and continuum effects on the single-particle structure of 28,29,31 F from a relativistic mean-field approach and show how both phenomena can lead to a negative-parity ground state in 28 F and halo structures in 29,31 F. We solve the Dirac equation in the complex-momentum (Berggren) representation for a potential with quadrupole deformation at the first order obtained from relativistic mean-field calculations using the NL3 interaction and calculate the continuum level densities using the Green’s function method. We extract single-particle energies and widths from the continuum level densities to construct the Nilsson diagrams of 28,29,31 F in the continuum and analyze the evolution of both the widths and occupation probabilities of relevant Nilsson orbitals in 28 F and find that some amount of prolate deformation must be present. In addition, we calculate the density distributions for bound Nilsson orbitals near the Fermi surface in 29,31 F and reveal that, for a quadrupole deformation 0.3 ≤ β 2 ≤ 0.45 (prolate), characteristic halo tails appear at large distances. Using the relativistic mean-field approach in the complex-momentum representation with the Green’s function method, we demonstrate that in neutron-rich fluorine isotopes, while in the spherical case the pf shells are already inverted and close to the neutron emission threshold, a small amount of quadrupole deformation can dramatically reduce the gap between positive- and negative-parity states and increase the role of continuum states, ultimately leading to the negative parity in the ground state of 28 F and the halo structures in 29,31 F.
Normal-ordering approximations and translational (non)invariance
Not Available
Using spin alignment of inelastically excited nuclei in fast beams to assign spins: The spectroscopy of 13 O as a test case
Excited states in O 13 were investigated using inelastic scattering of an E / A = 69.5 MeV O 13 beam off of a Be 9 target. The excited states were identified in the invariant-mass spectra of the decay products. Both single-proton and sequential two-proton decays of the excited states were examined. Furthermore, for a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave Born-approximation calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced O 13 level scheme is compared to ab initio no-core shell model with continuum predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1 s 1 / 2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use.
Differential analysis of incompressibility in neutron-rich nuclei
Both the incompressibility K A of a finite nucleus of mass A and that (K ∞ ) of infinite nuclear matter are fundamentally important for many critical issues in nuclear physics and astrophysics. While some consensus has been reached about K ∞ , accurate theoretical predictions and experimental extractions of K τ characterizing the isospin dependence of K A have been very difficult. We propose a differential approach to extract K τ and K ∞ independently from the K A data of any two nuclei in a given isotope chain. Applying this method to the K A data from isoscalar giant monopole resonances (ISGMR) in even-even Pb, Sn, Cd, and Ca isotopes taken by Garg et al. at the Research Center for Nuclear Physics (RCNP), Osaka University, Japan, we find that the 106 Cd– 116 Cd and 112 Sn– 124 Sn pairs having the largest differences in isospin asymmetries in their respective isotope chains measured so far provide consistently the most accurate up-to-date K τ value of K τ = –616 ± 59 MeV and K τ =–623 ± 86 MeV, respectively, largely independent of the remaining uncertainties of the surface and Coulomb terms in expanding K A , while the K ∞ values extracted from different isotopes chains are all well within the current uncertainty range of the community consensus for K ∞ . Moreover, the size and origin of the “soft Sn puzzle” is studied with respect to the “stiff Pb phenomenon.” Furthermore, it is found that the latter is favored due to a much larger (by ≈ 380 MeV) K τ for Pb isotopes than for Sn isotopes, while K ∞ from analyzing the K A data of Sn isotopes is only about 5 MeV less than that from analyzing the Pb data.
Quantum many-body calculations using body-centered cubic lattices
It is often computationally advantageous to model space as a discrete set of points forming a lattice grid. This technique is particularly useful for computationally difficult problems such as quantum many-body systems. For reasons of simplicity and familiarity, nearly all quantum many-body calculations have been performed on simple cubic lattices. Since the removal of lattice artifacts is often an important concern, it would be useful to perform calculations using more than one lattice geometry. In this paper we show how to perform quantum many-body calculations using auxiliary-field Monte Carlo simulations on a three-dimensional body-centered cubic (BCC) lattice. As a benchmark test we compute the ground state energy of 33 spin-up and 33 spin-down neutrons in the unitary limit, which is an idealized limit where the interaction range is zero and scattering length is infinite. As a fraction of the free Fermi gas energy E FG , we find that the ground state energy is E 0 /E FG =0.369(2),0.371(2), using two different definitions of the finite-system energy ratio. This is in excellent agreement with recent results obtained on a cubic lattice [He et al., Phys. Rev. A 101, 063615 (2020)]. We find that the computational effort and performance on a BCC lattice is approximately the same as that for a cubic lattice with the same number of lattice points. We discuss how the lattice simulations with different geometries can be used to constrain the size of lattice artifacts in simulations of continuum quantum many-body systems.
Spectroscopic factors in dripline nuclei
Single-nucleon knockout reaction studies of the proton-dripline nuclei 9 C and 13 O suggest an appreciable suppression of spectroscopic factors. In this work, we calculate the one-neutron and one-proton spectroscopic factors for the mirror pair 9 C – 9 Li and 13 O using two variants of the continuum shell model: the complex-energy Gamow shell model and the real-energy shell model embedded in the continuum. Furthermore, our results indicate that the continuum effects strongly suppress the spectroscopic factors of well-bound orbits in the dripline systems but have less impact on the spectroscopic factors of weakly bound states.
Electron scattering on A = 3 nuclei from quantum Monte Carlo based approaches
We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
Glassy quantum nuclear pasta in neutron star crusts
We conduct a comprehensive survey of the shape parameter space of the nuclear pasta phases in neutron star crusts by conducting three-dimensional Hartree-Fock + BCS calculations. Spaghetti, waffles, lasagna, bicontinuous phases and cylindrical holes occupy local minima in the resulting constant-pressure Gibbs energy surfaces, implying multiple geometries coexist at a given depth. Notably, the bicontinuous phase, in which both the neutron gas and nuclear matter extend continuously in all dimensions appears over a large depth range. Our results support the idea that nuclear pasta is a glassy system. At a characteristic temperature, of order 10 8 – 10 9 K, different phases may become frozen into domains whose sizes we estimate to be 1–50 times the lattice spacing and over which the local density and electron fraction can vary. Above this temperature, very little long-range order exists and matter is an amorphous solid. Electron scattering off domain boundaries may contribute to the disorder resistivity of the pasta phases. Annealing of the domains may occur during cooling; repopulating of local minima during crustal heating might lead to temperature-dependent transport properties in the deep crust layers. We identify four regions distinguished by whether pasta is the true ground state, and whether the pasta structure allows delocalization of protons. The whole pasta region can occupy up to 70% of the crust by mass and 25% by thickness, and the layer in which protons are delocalized could occupy 45% of the crust mass and 15% of its thickness.